Skip to contents

Returns the data from all random effects terms.

Usage

get_random(x)

Arguments

x

A fitted mixed model.

Value

The data from all random effects terms, as data frame. Or NULL

if model has no random effects.

Examples

data(sleepstudy)
# prepare some data...
sleepstudy$mygrp <- sample(1:5, size = 180, replace = TRUE)
sleepstudy$mysubgrp <- NA
for (i in 1:5) {
  filter_group <- sleepstudy$mygrp == i
  sleepstudy$mysubgrp[filter_group] <-
    sample(1:30, size = sum(filter_group), replace = TRUE)
}

m <- lmer(
  Reaction ~ Days + (1 | mygrp / mysubgrp) + (1 | Subject),
  data = sleepstudy
)
#> boundary (singular) fit: see help('isSingular')

head(get_random(m))
#>   mysubgrp mygrp Subject
#> 1       25     4     308
#> 2        6     3     308
#> 3       26     3     308
#> 4       24     5     308
#> 5       19     3     308
#> 6        3     1     308