Plot the results of a simulation of the random effects

Plot the simulated random effects on a ggplot2 chart. Points that are distinguishable from zero (i.e. the confidence band based on level does not cross the red line) are highlighted. Currently, the plots are ordered according to the grouping factor.

Usage

plotREsim(
  data,
  level = 0.95,
  stat = "median",
  sd = TRUE,
  sigmaScale = NULL,
  oddsRatio = FALSE,
  labs = FALSE,
  facet = TRUE
)

Arguments

data: a data.frame generated by REsim with simulations of the random effects of a merMod
level: the width of the confidence interval
stat: a character value indicating the variable name in data of the midpoint of the estimated interval, e.g. "mean" or "median"
sd: a logical indicating whether or not to plot error bars around the estimates (default is TRUE). Calculates the width of the error bars based on level and the variable named "sd" in data
sigmaScale: a numeric value to divide the estimate and the standard deviation by in the case of doing an effect size calculation
oddsRatio: logical, should the parameters be converted to odds ratios before plotting
labs: logical, include the labels of the groups on the x-axis
facet: Accepts either logical (TRUE) or list to specify which random effects to plot. If TRUE, facets by both groupFctr and term. If list selects the panel specified by the named elements of the list

Value

a ggplot2 plot of the coefficient effects

Examples

# \donttest{
 fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)
 (p1 <- plotREsim(REsim(fm1)))

 #Plot just the random effects for the Days slope
 (p2 <- plotREsim(REsim(fm1), facet= list(groupFctr= "Subject", term= "Days")))

 # }