library(agricolae)
data(heterosis)
str(heterosis)
site1<-subset(heterosis,heterosis[,1]==1)
site2<-subset(heterosis,heterosis[,1]==2)
site3<-subset(heterosis,heterosis[,1]==3)
attach(site1)
model1<-lineXtester(Replication, Female, Male, v1)
DFe <- df.residual(model1)
CMe <- deviance(model1)/DFe
test1 <- HSD.test(v1, Factor,DFe,CMe)
test2 <- HSD.test(v1, Treatment,DFe,CMe)
detach(site1)
attach(site2)
model22<-lineXtester(Replication, Female, Male, v3)
detach(site2)
attach(site3)
model3<-lineXtester(Replication, Female, Male, v4)
detach(site3)
Run the code above in your browser using DataLab