Displaying 1 result from an estimated 1 matches for "vc1m".
Did you mean:
vc1
2005 Aug 04
1
Where the error message comes from?
...=0.5,z=1.281552) # for p=0.1
## TREATMENT ##
r1 <- cor(d1[,2:3])[1,2]
# COST #
ac1 <- anova(lm(c~g,d1))
rc1 <- (ac1[1,3]-ac1[2,3])/(ac1[1,3]+(m-1)*ac1[2,3])
if (rc1 < 0) rc1 <- 0
vc1 <- var(d1[,2])
mc1 <- as.vector(by(d1[,2],as.numeric(d1[,1]),mean))
vc1m <- vc1*(1+(m-1)*rc1)/(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...