trt <- rep(c("A","V","C"),each=5)
rabbit <- seq(1:15)

y.mat <- matrix(c(-0.3,-0.2,1.2,3.1,-0.5,2.2,3.3,3.7,-1.1,2.4,2.2,2.7,
1.0,1.7,2.1,2.5,-0.3,0.8,0.6,0.9,
-1.1,-2.2,0.2,0.3,-1.4,-0.2,-0.5,-0.1,-0.1,-0.1,-0.5,-0.3,
-0.2,0.1,-0.2,0.4,-0.1,-0.2,0.7,-0.3,
-1.8,0.2,0.1,0.6,-0.5,0.0,1.0,0.5,-1.0,-0.3,-2.1,0.6,
0.4,0.4,-0.7,-0.3,-0.5,0.9,-0.4,-0.3),byrow=T,ncol=4)

y1 <- y.mat[,1]
y2 <- y.mat[,2]
y3 <- y.mat[,3]
y4 <- y.mat[,4]

multi_y <- data.frame(trt,rabbit,y1,y2,y3,y4)
attach(multi_y)

trt <- factor(trt)
rabbit <- factor(rabbit)

eartemp.mod1 <- manova(cbind(y1,y2,y3,y4) ~ trt)
summary(eartemp.mod1,test="Pillai")
summary(eartemp.mod1,test="Wilks")
summary(eartemp.mod1,test="Hotelling-Lawley")
summary(eartemp.mod1,test="Roy")

rm(multi_y)



trt <- rep(c("A","V","C"),each=20)
rabbit <- rep(seq(1:15),each=4)
time <- rep(seq(0,90,30),times=15)

y <- c(-0.3,-0.2,1.2,3.1,-0.5,2.2,3.3,3.7,-1.1,2.4,2.2,2.7,
1.0,1.7,2.1,2.5,-0.3,0.8,0.6,0.9,
-1.1,-2.2,0.2,0.3,-1.4,-0.2,-0.5,-0.1,-0.1,-0.1,-0.5,-0.3,
-0.2,0.1,-0.2,0.4,-0.1,-0.2,0.7,-0.3,
-1.8,0.2,0.1,0.6,-0.5,0.0,1.0,0.5,-1.0,-0.3,-2.1,0.6,
0.4,0.4,-0.7,-0.3,-0.5,0.9,-0.4,-0.3)

uni_y <- data.frame(trt,rabbit,time,y)
attach(uni_y)

trt <- factor(trt)
rabbit <- factor(rabbit)
time <- factor(time,ordered=T)

eartemp.mod2 <- aov(y ~ trt + trt/rabbit + time + trt:time)

summary(eartemp.mod2)


eartemp.mod3 <- aov(y ~ trt*time + Error(rabbit))

summary(eartemp.mod3)

summary(eartemp.mod3,split=list(time=list(linear=1, quadratic=2,
cubic=3)))