Hello,

I have a dataset derived from a multi-species experiment that I carried out.

I have 30 different plant species from several families and 3 geographical groups. I applied 6 different treatments to 5 replicates from each species.

For the analysis I was thinking of a GLMM with germination % as response and group and treatment as predictors. However, I'd like to take into account the phylogeny of the species as a random factor.

I managed to create a phylogenetic tree (and vcv matrix) using the R package V.PhyloMaker, but now I can't seem to find a way to implement it in a GLMM. I'm aware of the package pez, but it doesn't really work for me as i don't have community composition data. Another solution seems the MCMCglmm package, but I'm not familiar with bayesian statistics.ù

Any suggestions would be highly appreciated!!