[R-sig-ME] Extracting the posterior distribution for a random, effect in MCMCglmm

W Robert Long longrob604 at gmail.com
Mon Aug 6 17:07:38 CEST 2012

Thanks Ned, Thierry and Jarrod

My code is working nicely now.


On 6/08/2012 4:05 PM, Ned Dochtermann wrote:
> The values you're referring to in the solution part of the output (Sol)
> are what for a linear model would be the BLUPs, hence some are negative.
> The random effects variances live in VCV. If I recall correctly
> posterior modes and HPD intervals for a model like yours can be
> extracted as:
> posterior.mode(m2a.7$VCV[,1])
> HPDinterval(m2a.7$VCV[,1])

More information about the R-sig-mixed-models mailing list