[R-sig-ME] Simulating linear mixed model data for factorial design with 3 levels

Maarten Jung Maarten.Jung at mailbox.tu-dresden.de
Fri Mar 30 11:16:38 CEST 2018

It would be really helpful if anyone could give a hint.
I also added the exact model definition to make clear what I am trying to


On Sat, Mar 24, 2018, 11:54 Maarten Jung <Maarten.Jung at mailbox.tu-dresden.de>

> Dear list,
> I struggle to simulate linear mixed model data for a within-subject
> within-item factorial design (with crossed participant and item effects)
> where the categorical independent variable has three levels.
> The default parameters correspond to the estimates I got from an analysis
> of previous data.
> As I am using the mvrnorm() function with empirical = TRUE I expected
> that the estimated parameters from lmer and my specified parameters for
> the simulation would match closely.
> However, they differ more than I expected.
> I would be grateful if somebody more experienced in mixed model
> simulations could have a look at my R code on Cross Validate:
> https://stats.stackexchange.com/q/335643/136579
> Best Regards,
> Maarten

	[[alternative HTML version deleted]]

More information about the R-sig-mixed-models mailing list