powerSim() R function from [simr]

Estimate power by simulation.

Perform a power analysis for a mixed model.


powerSim(
  fit,
  test = fixed(getDefaultXname(fit)),
  sim = fit,
  fitOpts = list(),
  testOpts = list(),
  simOpts = list(),
  seed,
  ...
)

fit: a fitted model object (see doFit).
test: specify the test to perform. By default, the first fixed effect in fit will be tested. (see: tests ).
sim: an object to simulate from. By default this is the same as fit (see doSim).
fitOpts: extra arguments for doFit.
testOpts: extra arguments for doTest.
simOpts: extra arguments for doSim.
seed: specify a random number generator seed, for reproducible results.
...: any additional arguments are passed on to simrOptions. Common options include:
- nsim:: the number of simulations to run (default is 1000).
- alpha:: the significance level for the statistical test (default is 0.05).
- progress:: use progress bars during calculations (default is TRUE).


fm1 <- lmer(y ~ x + (1|g), data=simdata)
powerSim(fm1, nsim=10)