Message-ID: <501FDDBA.3090803@gmail.com>
Date: 2012-08-06T15:07:38Z
From: Robert Long
Subject: Extracting the posterior distribution for a random, effect in MCMCglmm
In-Reply-To: <501FDD2A.3090305@gmail.com>
Thanks Ned, Thierry and Jarrod
My code is working nicely now.
Regards
Rob
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])
>