Displaying 1 result from an estimated 1 matches for "rfrac".
Did you mean:
frac
2006 Nov 24
2
low-variance warning in lmer
...("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.55),axes=FALSE)
axis(side=2)
box()
boxplot(rfrac...