R code for simulating marginal predictions from random effects models