# ranef.mumm: Extract Random Effects In mumm: Multiplicative Mixed Models using the Template Model Builder

## Description

A function to extract the estimated random effects from a model object of class `mumm`.

## Usage

 ```1 2``` ```## S3 method for class 'mumm' ranef(object, ...) ```

## Arguments

 `object` an object of class "mumm" `...` Currently not used

## Value

A named list with the estimated random effects, where each element in the list is a numeric vector consisting of the estimated random effect coefficients for a random factor in the model.

## Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28``` ```set.seed(100) sigma_e <- 1.5 sigma_a <- 0.8 sigma_b <- 0.5 sigma_d <- 0.7 nu <- c(8.2, 6.2, 2.3, 10.4, 7.5, 1.9) nA <- 15 nP <- 6 nR <- 5 a <- rnorm(nA, mean = 0, sd = sigma_a) b <- rnorm(nA, mean = 0, sd = sigma_b) d <- rnorm(nA*nP, mean = 0, sd = sigma_d) e <- rnorm(nA*nP*nR, mean = 0, sd = sigma_e) Assessor <- factor(rep(seq(1,nA),each = (nP*nR))) Product <- factor(rep(rep(seq(1,nP),each = nR), nA)) AssessorProduct <- (Assessor:Product) y <- nu[Product] + a[Assessor] + b[Assessor]*(nu[Product]-mean(nu)) + d[AssessorProduct] + e sim_data <- data.frame(y, Assessor, Product) fit <- mumm(y ~ 1 + Product + (1|Assessor) + (1|Assessor:Product) + mp(Assessor,Product) ,data = sim_data) ranef(fit) ```

mumm documentation built on May 2, 2019, 9:44 a.m.