Displaying 1 result from an estimated 1 matches for "ve2m".
Did you mean:
v2m
2005 Aug 04
1
Where the error message comes from?
.../(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(me2)*(1-mean(me2)))
ve2m <- ve2*(1+(m-1)*re2)/(k*m)
## Combining and Computing ICER ##
d <- data.frame(s = c(rep(1,k),rep(2,k)),
i = c(1:k,1:k), n =rep(m,2*k),
mc=c(mc1,mc2),me=c(me1,me2))
mmc1 <- mean(d[1:k,4])
mmc2 <- mean(d[(k+1):(2*k),4])
mme1 <- mean(d[1:k,5])
mm...