djhocking

Predict vs simulate in lme4

For this investigation we are going to use the sleepdata data set from the lme4 package. Here is the head of the data frame:

djhocking

# make fake dataset
df <- data.frame(x = runif(100, 0, 1), y = rnorm(100, 10, 3), z = rpois(100, 10))

# subset dataframe
df_sub <- df[which(df$x >= 0.75), ]

# subset using dplyr
library(dplyr)

df_sub2 <- df %>%