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.

# S3 method for MCMCglmm
ranef(object, use = c("all", "mean"), ...)

Arguments

object

An MCMCglmm model object to extract the effects from

use

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.

Value

A matrix of the fixed effects

See also

Examples

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, ])
}