many changes!

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

diff --git a/R/compar.gee.R b/R/compar.gee.R
index caead0776b8dc6161a14034805071bdc89fbcc45..b2a9c1eb7eefa9acead17c37b1519b4507848b90 100644 (file)
--- a/R/compar.gee.R
+++ b/R/compar.gee.R
@@ -57,7 +57,7 @@ compar.gee <-
     ## <FIXME>
     ## maybe need to refine below in case of non-Brownian corStruct
     if (!missing(corStruct)) phy <- attr(corStruct, "tree")
-    dfP <- sum(phy$edge.length)*N / sum(diag(R))
+    dfP <- sum(phy$edge.length)*N / sum(diag(vcv(phy))) # need the variances
     ## </FIXME>
 
     ## compute QIC: