-## compar.gee.R (2010-07-20)
+## compar.gee.R (2011-06-14)
## Comparative Analysis with GEEs
}
id <- rep(1, dim(R)[1])
- geemod <<- geemod <- do.call("gee", list(formula, id, data = data, family = family, R = R,
+ geemod <- do.call("gee", list(formula, id, data = data, family = family, R = R,
corstr = "fixed", scale.fix = scale.fix,
scale.value = scale.value))
W <- geemod$naive.variance
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)
+}