-## compar.gee.R (2010-07-20)
+## compar.gee.R (2011-06-14)
## Comparative Analysis with GEEs
QIC = QIC,
coefficients = geemod$coefficients,
residuals = geemod$residuals,
+ fitted.values = MU,
family = geemod$family$family,
link = geemod$family$link,
scale = geemod$scale,
as.character(as.expression(fm)), "\n")
ans
}
+
+predict.compar.gee <-
+ function(object, type = c("link", "response"), ...)
+{
+ type <- match.arg(type)
+ pred <- object$fitted.values
+ if (type == "link") return(pred)
+ f <- match.fun(object$family)
+ f(link = object$link)$linkinv(pred)
+}