Displaying 1 result from an estimated 1 matches for "lev4".
Did you mean:
dev4
2003 Dec 17
2
variance estimates in lme biased?
...y(nlme)
F1V <- F2V <- F3V <- residV <- numeric(0)
for (rep in 1:500){
F1f <- F2f <- F3f <- value <- numeric(0)
# data generator
for (F1 in 1:5) {
for (F2 in 1:5) {
lev2 <- rnorm(1,sd=1.0)
for (F3 in 1:5) {
for (F4 in 1:5) {
lev4 <- rnorm(1,sd=0.1)
value <- c(value, lev2+lev4)
F1f <- c(F1f,F1);F2f <- c(F2f,F2);F3f <- c(F3f,F3)
}
}
}
}
L1 <- as.factor(F1f); L2 <- as.factor(F2f); L3 <- as.factor(F3f)
y.lme <- lme(value ~ 1,random = ~ 1 | L1/L2/L3)...