Displaying 1 result from an estimated 1 matches for "vc2m".
Did you mean:
vc2
2005 Aug 04
1
Where the error message comes from?
...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 <- vc2*(1+(m-1)*rc2)/(k*m)
# EFECT #
ae2 <- anova(lm(e~g,d2))
re2 <- (ae2[1,3]-ae2[2,3])/(ae2[1,3]+(m-1)*ae2[2,3])
if (re2 < 0) re2 <- 0
ve2 <- var(d2[,3])
me2 <- as.vector(by(d2[,3],as.numeric(d2[,1]),mean))
re2p <- 1- m*sum(me2*(1-me2))/(k*(m-1)*mean...