# NOT RUN {
data<-data.frame(Treatments=c("T1","T2","T3","T4","T5","T6","T7","T1","T2","T3","T4","T5","T6",
"T7","T1","T2","T3","T4","T5","T6","T7"),
yield=c(25,21,21,18,25,28,24,25,24,24,16,21,20,17,16,19,14,15,13,11,25),
height=c(130,120,125,135,139,140,145,136,129,135,150,152,140,148,130,135,145,160,145,130,160))
#CRD analysis with LSD test for yield only
crd(data[2],data$Treatments,1)
#CRD analysis with LSD test for both yield and height
crd(data[2:3],data$Treatments,1)
# }
Run the code above in your browser using DataLab