# fitness example, p 81
dat <- c(7, 5, 9, 8, 3, 12, 10, 8, 7, 9,
5, 7, 7, 6, 4, 8, 12, 9, 7, 8,
3, 2, 7, 6, 8, 6, 5, 3, 4, 3,
4, 3, 3, 1, 2, 5, 7, 6, 8, 7)
gp <- as.factor(rep(1:4,each=10))
gp = gl(4,10,40, labels=c("6 hr", "3 hr", "Sports club", "Video games"))
contrast1 <- c(-3, -1, 1, 3) # linear
contrast2 <- c(1, -1, -1, 1) # quadratic
contrast3 <- c(1, 1, -1, -1)
L_1way_ANOVA(dat,gp,contrast3, contrast1)
Run the code above in your browser using DataLab