Dear R-Sig-Me and R-Sig-Phylo users,
I have recently started using MCMCglmm with the aim of conducting
phylogenetic meta-analysis. I have managed to run a basic model using
MCMCglmm(), with weak uninformative priors specified, but I am
struggling to work out how to incorporate a phylogeny.
Hi all, Jarrod made a speedy response to my phylogenetic meta-analysis
query in MCMCglmm on the R sig-me list. Here was his response below.
I'll try both of his recommendations, and will post on how I get on.
Cheers,
Wayne
--
Dr. Wayne Dawson
Institute of Plant Sciences
University of Bern
Dear Wayne,
This is my fault. With phylogenies the ancestral nodes are treated as
missing data and so I set their measurement error to an arbitrary
value. The code for working out how many new measurement errors
there are was incorrect.
L98 of MCMCglmm.R should read
mev-c(mev, rep(1,