Displaying 1 result from an estimated 1 matches for "nreef".
Did you mean:
reef
2006 Nov 24
2
low-variance warning in lmer
...aven't found anything relevant 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 =...