various updates for ape 2.7-2 (fixed)

[ape.git] / R / compar.gee.R

diff --git a/R/compar.gee.R b/R/compar.gee.R
index 8c9bf868fafd23580dc6eabdcd1bb8e62a81aaec..1df21079ca726daed70b8db034ea6e42eb8d42ce 100644 (file)
--- 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)
+}