Displaying 1 result from an estimated 1 matches for "natlevel".
Did you mean:
atlevel
2003 Dec 10
0
nested analysis with lme - odd result?
...expecting near zero variance in those other levels except say, 5% of
the time.
Here's some example code with variance simulated in level 2 and spurious
variance appearing frequently in level 1:
library(nlme)
# 4 level nesting
simVar <- c(0,1,0,.1) # first is level one, last becomes error
nAtLevel <- c(5,5,5,5) # number of replicates at each level
F1V <- F2V <- F3V <- residV <- numeric(0)
for (rep in 1:100){
F1f <- F2f <- F3f <- value <- numeric(0)
mn <- numeric(length(nAtLevel))
# data generator
for (F1 in 1:nAtLevel[1]) {
mn[1] <- rnorm(...