Function designed to extract the random effects from an
MCMCglmm
model object. Can either extract all samples from the
random effects posteriors or return the posterior means.
An MCMCglmm
model object to extract the effects from
A character string indicating whether to extract all posterior samples or the mean of the posteriors. Defaults to "all".
Arguments passed on to the worker function.
A matrix of the fixed effects
if (FALSE) {
# a simple MCMCglmm model
data(PlodiaPO)
m <- MCMCglmm(PO ~ 1, random= ~ FSfamily, data=PlodiaPO, pr=TRUE, verbose=FALSE)
# only extract average fixed effects
head(ranef(m, use = "mean"))
# histogram of posterior samples of fixed effects
hist(ranef(m)[1, ])
# matches the mean
rowMeans(ranef(m)[1:6, ])
}