Dear R users,

Is there a simple way to use the effect function (library(effects)) with
a GEE estimated model?

 

library(gee)

library(MASS)

library(effects)

attach(epil)

 

b = gee(y ~ lbase*trt + lage + V4, family=poisson, id=subject,
corstr="exchangeable")

 

plot(effect("lbase*trt", b))

# Errore in effect("lbase*trt", b) : nessun metodo applicabile per
"effect"

 

 


	[[alternative HTML version deleted]]

