X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;ds=sidebyside;f=R%2Fcompar.gee.R;h=1df21079ca726daed70b8db034ea6e42eb8d42ce;hb=2e7fdb027f189e7c8eb771dadaf195c43bd41d5f;hp=8c9bf868fafd23580dc6eabdcd1bb8e62a81aaec;hpb=29792fb8cef8ddab4553474398369fb6c24c1d8b;p=ape.git

diff --git a/R/compar.gee.R b/R/compar.gee.R
index 8c9bf86..1df2107 100644
--- a/R/compar.gee.R
+++ b/R/compar.gee.R
@@ -1,4 +1,4 @@
-## compar.gee.R (2010-07-20)
+## compar.gee.R (2011-06-14)
 
 ##   Comparative Analysis with GEEs
 
@@ -80,6 +80,7 @@ compar.gee <-
                 QIC = QIC,
                 coefficients = geemod$coefficients,
                 residuals = geemod$residuals,
+                fitted.values = MU,
                 family = geemod$family$family,
                 link = geemod$family$link,
                 scale = geemod$scale,
@@ -151,3 +152,13 @@ you should be careful when interpreting the significance of the main effects.")
                                   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)
+}