Displaying 1 result from an estimated 1 matches for "nreefpertreat".
2006 Nov 24
2
low-variance warning in lmer
...evant in Pinheiro and Bates, but
I may not have looked in the right place ...)
thanks
Ben Bolker
self-contained but unnecessarily complicated simulation
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 = i...