Hi all,

I'm successfully using ape::binaryPGLMM to fit logistic regression models while taking phylogeny into account. I can obtain the estimates, standard errors and p-values, as they are directly returned in the model object.


However, can I calculate confidence intervals for the regression coefficients?

One way to get some numbers is to define the vcov and coef functions for binaryPGLMM objects and use the standard confint function:

m <- binaryPGLMM(Y ~ X1, phy=phy, data=sim.dat) # from ?binaryPGLMM
vcov.binaryPGLMM <- function(object, ...) { object$B.cov }
coef.binaryPGLMM <- function(object, ...) { object$B[, 1] }
confint(m)

This gives us CI's, but given the fact that confint.glm uses profiling it seems unlikely to me that these numbers are correct.

If I can't figure it out, I think I will change the phylopath package to use phylolm::phyloglm instead.

Thanks for your help in advance,

Wouter

_______________________________________________
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/

Reply via email to