# Kaps and Lamberson(2009)
data(data1)
data(data2)
data(data3)
data(data4)
# analysis in completely randomized design
r1<-ea1(data1)
names(r1)
r1
# analysis in randomized block design
r2<-ea1(data2)
# analysis in latin square design
r3<-ea1(data3)
# analysis in several latin squares design
r4<-ea1(data4)
r1[1]
r2[1]
r3[1]
r4[1]
# analysis in unbalanced randomized block design
response<-ifelse(data2$Gain>850, NA, data2$Gain)
ndata<-data.frame(data2[-3],response)
ndata
r5<-ea1(ndata)
r5['Adjusted means']
Run the code above in your browser using DataLab