##Kline (2004) Table 6.3
dat <- data.frame(y = c(9,12,13,15,16,
8,12,11,10,14,
10,11,13,11,15),
x = rep(factor(c("a","b","c")), each=5)
)
##contrast 1: a - c, contrast 2: 1/2(a + c) - b
my.cont <- matrix(c(1,0,-1,1/2,-1,1/2), ncol=3, nrow=2, byrow=TRUE)
ind.oneway.second(m = tapply(dat$y, dat$x, mean),
sd = tapply(dat$y, dat$x, sd),
n= tapply(dat$y, dat$x, length))
ind.oneway.second(m = tapply(dat$y, dat$x, mean),
sd = tapply(dat$y, dat$x, sd),
n= tapply(dat$y, dat$x, length),
contr = my.cont)
Run the code above in your browser using DataLab