##Kline (2004) Table 4.4
dat <- data.frame(y = c(9,12,13,15,16,8,12,11,10,14),
x = rep(factor(c("a","b")), each=5)
)
ind.t.test.second(m = tapply(dat$y, dat$x, mean),
sd = tapply(dat$y, dat$x, sd),
n = tapply(dat$y, dat$x, length)
)
Run the code above in your browser using DataLab