Generating predictive simulations from a multilevel model with random intercepts
问题 I am struggling to understand how, in R, to generate predictive simulations for new data using a multilevel linear regression model with a single set of random intercepts. Following the example on pp. 146-147 of this text, I can execute this task for a simple linear model with no random effects. What I can't wrap my head around is how to extend the set-up to accommodate random intercepts for a factor added to that model. I'll use iris and some fake data to show where I'm getting stuck. I'll