X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=R%2Fcompar.gee.R;h=e9075a1f724adf6a86ea9676bfa312ab19db385a;hb=0f717ab8186382152b581c9dd67bf85779112e7f;hp=f8a7b6ae43e217c51b79e8e79890f19021141358;hpb=d88302b4735b5b7c9132387090bb592d906fe1cb;p=ape.git diff --git a/R/compar.gee.R b/R/compar.gee.R index f8a7b6a..e9075a1 100644 --- a/R/compar.gee.R +++ b/R/compar.gee.R @@ -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)))