Displaying 1 result from an estimated 1 matches for "ve1m".
Did you mean:
ve1
2005 Aug 04
1
Where the error message comes from?
.../(k*m)
# EFFECT #
ae1 <- anova(lm(e~g,d1))
re1 <- (ae1[1,3]-ae1[2,3])/(ae1[1,3]+(m-1)*ae1[2,3])
if (re1 < 0) re1 <- 0
ve1 <- var(d1[,3])
me1 <- as.vector(by(d1[,3],as.numeric(d1[,1]),mean))
re1p <- 1- m*sum(me1*(1-me1))/(k*(m-1)*mean(me1)*(1-mean(me1)))
ve1m <- ve1*(1+(m-1)*re1)/(k*m)
## CONTROL ##
r2 <- cor(d2[,2:3])[1,2]
# COST #
ac2 <- anova(lm(c~g,d2))
rc2 <- (ac2[1,3]-ac2[2,3])/(ac2[1,3]+(m-1)*ac2[2,3])
if (rc2 < 0) rc2 <- 0
vc2 <- var(d2[,2])
mc2 <- as.vector(by(d2[,2],as.numeric(d2[,1]),mean))
vc2m...