Displaying 2 results from an estimated 2 matches for "residv".
Did you mean:
reside
2003 Dec 10
0
nested analysis with lme - odd result?
...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(1,sd=simVar[1]) # set mean for level 1
for (F2 in 1:nAtLevel[2]) {
mn[2] <- rno...
2003 Dec 17
2
variance estimates in lme biased?
...the
level 1 estimate is close to level 2's. Nothing like that happens in
level 3.
I use R 1.8.1 & Win2000. Since I last posted, I have also seen these
symptoms in SAS when I use PROC VARCOMP, although I haven't explored it
much there.
library(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 &...