Displaying 1 result from an estimated 1 matches for "npertreat".
2006 Nov 24
2
low-variance warning in lmer
...imulation
code/demonstration:
---------------
library(lme4)
library(lattice)
simfun <- function(reefeff,ntreat=2,nreef=12,
nreefpertreat=3,
t.eff=10,
totvar=25,seed=NA) {
if (!is.na(seed)) set.seed(seed)
ntot = nreef*nreefpertreat
npertreat=ntot/ntreat
reef = gl(nreef,nreefpertreat)
treat = gl(ntreat,npertreat)
r.sd = sqrt(totvar*reefeff)
e.sd = sqrt(totvar*(1-reefeff))
y.det = ifelse(treat==1,0,t.eff)
r.vals = rnorm(nreef,sd=r.sd)
e.vals = rnorm(ntot,sd=e.sd)
y <- y.det+r.vals[as.numeric(reef)]+e.vals
data.frame(...