Displaying 1 result from an estimated 1 matches for "getranvar".
2006 Nov 24
2
low-variance warning in lmer
...ef = 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(y,treat,reef)
}
getranvar <- function(x) {
vc <- VarCorr(x)
c(diag(vc[[1]]),attr(vc,"sc")^2)
}
estfun <- function(reefeff,...) {
x <- simfun(reefeff=reefeff,...)
ow <- options(warn=2)
f1 <- try(lmer(y~treat+(1|reef),data=x))
w <- (class(f1)=="try-error" && length...