Hi Diederik,
you can use MCMCglmm. The package allows for inclusion of phylogenetic
information, random effects and zero-inflated response variables. However, it
may take some time to get familiar with the package.
Best,
J
—
Jörg Albrecht, PhD
Postdoctoral researcher
Institute of Nature
l",
prior = prior1, data = dat)
# fixed effects should be around zero and one (+/- monte carlo error)
summary(m1$Sol)
# phylogenetic ICC should be 1/(2 + pi^2/3) = 0.189 (+/- monte carlo error)
summary(m1$VCV[, 1] / (rowSums(m1$VCV) + pi^2/3))
Hope this helps,
Jörg
—
Jörg Alb