Displaying 1 result from an estimated 1 matches for "reefvar".
Did you mean:
refvar
2006 Nov 24
2
low-variance warning in lmer
...s(f1)=="try-error" && length(grep("effectively zero",f1))>0)
options(ow)
f2 <- lmer(y~treat+(1|reef),data=x)
c(getranvar(f2),as.numeric(w))
}
rvec <- rep(c(0.01,0.05,0.1,0.15,0.2,0.3,0.5),each=100)
X <- t(sapply(rvec,estfun))
colnames(X) <- c("reefvar","resvar","warn")
rfrac <- X[,"reefvar"]/(X[,"reefvar"]+X[,"resvar"])
fracwarn <- tapply(X[,"warn"],rvec,mean)
est.mean <- tapply(rfrac,rvec,mean)
op <- par(mfrow=c(1,2))
plot(rvec,rfrac,type="n",xlim=c(-0.02,0...