[R-sig-ME] BLUPs from MCMCglmm
Jarrod Hadfield
j.hadfield at ed.ac.uk
Thu Jan 19 12:34:43 CET 2012
Hi,
The random effects are actually stored in Sol (solutions) and will be
saved if you use the argument pr=TRUE. The marginal posterior modes
of the random effects should coincide with BLUPs if the variances are
fixed a priori to the value used when obtaining BLUPs and the fixed
effects are either fixed a priori or are given improper flat priors.
However, in practice (i.e. when the variances are not fixed etc.) the
correlation between BLUPs and marginal posterior modes is usually very
very high.
Be aware that if you have many (m) random effects and you store many
(n) iterations you end up with a lot (m*n) of numbers to store.
Cheers,
Jarrod
On 18 Jan 2012, at 14:51, Ben Bolker wrote:
> Eryn McFarlane <mcfarlas at ...> writes:
>
>>
>> Dear list,
>>
>> I was wondering if anyone knew of a way to estimate BLUPs
>> from an MCMCglmm model? I would just like to eyeball
>> the individuals with high and low BLUPs for my trait to
>> see if there are other relationships that I can see
>> (i.e. year effects, affect of territory). Does this make sense to
>> try to do
> from these models?
>
>
> I *think* you can just look at the $Liab component of the fit,
> which as stated is the posterior distribution of the latent
> variables --
> you need to set pl=TRUE.
>
> This should get you started (although HPDinterval() isn't
> behaving sensibly in this case -- not quite sure why not)
>
> data(PlodiaPO)
> model1<-MCMCglmm(PO~1, random=~FSfamily, data=PlodiaPO,
> verbose=FALSE, pl=TRUE)
> str(model1$Liab)
> mm <- data.frame(m=colMeans(model1$Liab),HPDinterval(model1$Liab))
> plot(mm[order(mm$m),"m")
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
--
The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.
More information about the R-sig-mixed-models
mailing list