Extract random effects from an MCMCglmm object

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.

Usage

# S3 method for class '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

fixef.MCMCglmm

Examples

if (FALSE) { # \dontrun{
  # 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, ])
} # }