powered by
REsim
Simulate random effects from merMod REsim simulates random effects from merMod object posterior distributions
REsim(merMod, n.sims = 200, oddsRatio = FALSE, seed = NULL)
a data frame with the following columns
groupFctr
Name of the grouping factor
groupID
Level of the grouping factor
term
Name of random term (intercept/coefficient)
mean
Mean of the simulations
median
Median of the simulations
sd
Standard deviation of the simulations, NA if oddsRatio=TRUE
NA
oddsRatio=TRUE
a merMod object from the lme4 package
number of simulations to use
logical, should parameters be converted to odds ratios?
numeric, optional argument to set seed for simulations
Use the Gelman sim technique to build empirical Bayes estimates. Uses the sim function in the arm package
require(lme4) m2 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy) re2 <- REsim(m2, 25) head(re2)
Run the code above in your browser using DataLab