]> git.donarmstrong.com Git - ape.git/blobdiff - R/compar.gee.R
fixes for ape 2.3-1
[ape.git] / R / compar.gee.R
index f8a7b6ae43e217c51b79e8e79890f19021141358..e9075a1f724adf6a86ea9676bfa312ab19db385a 100644 (file)
@@ -1,4 +1,4 @@
-## compar.gee.R (2008-01-14)
+## compar.gee.R (2008-02-21)
 
 ##   Comparative Analysis with GEEs
 
@@ -30,7 +30,9 @@ do not match: the former were ignored in the analysis.")
                                   corstr = "fixed", scale.fix = scale.fix,
                                   scale.value = scale.value))
     W <- geemod$naive.variance
-    if (family == "binomial")
+    fname <-
+        if (is.function(family)) deparse(substitute(family)) else family
+    if (fname == "binomial")
       W <- summary(glm(formula, family = quasibinomial, data = data))$cov.scaled
     N <- geemod$nobs
     dfP <- sum(phy$edge.length)*N / sum(diag(vcv.phylo(phy)))