#LMM 03 - SPLIT-PLOT ANOVA
setwd("c:/DATA/Models/")
#Oats in Venables & Ripley 2002 p. 282
library(MASS)
oats

options(digits=6)
LM1=aov(Y~V*N+Error(B/V),data=oats)
summary(LM1)


attach(oats)
#GRAND MEAN:
GM=mean(Y)
GM

#BLOCKS:
B1=Y[B=="I"]
B2=Y[B=="II"]
B3=Y[B=="III"]
B4=Y[B=="IV"]
B5=Y[B=="V"]
B6=Y[B=="VI"]
BLOCKS=cbind(B1,B2,B3,B4,B5,B6)
BLOCKS

Bb=tapply(oats$Y,oats$B,mean) #blocking factor means
Bb

#PLOT FACTOR V:
V1=Y[V=="Victory"]
V2=Y[V=="Golden.rain"]
V3=Y[V=="Marvellous"]
PLOTS=cbind(V1,V2,V3)
PLOTS
Vb=tapply(Y,oats$V,mean) #PLOT factor level means
Vb

#SUBPLOT FACTOR N:
N1=Y[N=="0.0cwt"]
N2=Y[N=="0.2cwt"]
N3=Y[N=="0.4cwt"]
N4=Y[N=="0.6cwt"]
SUBPLOTS=cbind(N1,N2,N3,N4)
SUBPLOTS
tapply(Y,oats$N,mean)

#CELLS C:
C11=subset(oats,V=="Victory"&N=="0.0cwt")$Y
C12=subset(oats,V=="Victory"&N=="0.2cwt")$Y
C13=subset(oats,V=="Victory"&N=="0.4cwt")$Y
C14=subset(oats,V=="Victory"&N=="0.6cwt")$Y
C21=subset(oats,V=="Golden.rain"&N=="0.0cwt")$Y
C22=subset(oats,V=="Golden.rain"&N=="0.2cwt")$Y
C23=subset(oats,V=="Golden.rain"&N=="0.4cwt")$Y
C24=subset(oats,V=="Golden.rain"&N=="0.6cwt")$Y
C31=subset(oats,V=="Marvellous"&N=="0.0cwt")$Y
C32=subset(oats,V=="Marvellous"&N=="0.2cwt")$Y
C33=subset(oats,V=="Marvellous"&N=="0.4cwt")$Y
C34=subset(oats,V=="Marvellous"&N=="0.6cwt")$Y
cbind(C11,C12,C13,C14)
cbind(C21,C22,C23,C24)
cbind(C31,C32,C33,C34)

C1=cbind(mean(C11),mean(C12),mean(C13),mean(C14))
C2=cbind(mean(C21),mean(C22),mean(C23),mean(C24))
C3=cbind(mean(C31),mean(C32),mean(C33),mean(C34))
rbind(C1,C2,C3)

detach(oats)