CRD, RCBD, LS, Factorial, Split plot, Rectangular lattice designs และ Correletion

#CRD

d<-read.table('clipboard',header=TRUE) 

d

attach(d)

trt<-factor(trt)

md<-aov(yield~trt)

summary(md)

library(agricolae)

cv.model(md)

print(LSD.test(md,'trt',alpha=0.05))

detach(d)

detach(d)

rm(list=ls(all=TRUE))

#RCBD

d<-read.table('clipboard',header=TRUE) 

d

attach(d)

names (d)

blk<-factor(blk)

trt<-factor(trt)

md<-aov(yield~blk+trt)

summary(md)

library(agricolae)

cv.model(md)

print(LSD.test(md,'trt',alpha=0.05))

detach(d)

detach(d)

rm(list=ls(all=TRUE))

#LS

d<-read.table('clipboard',header=TRUE) 

d

attach(d)

row<-factor(row)

col<-factor(col)

trt<-factor(trt)

md<-aov(yield~row+col+trt)

summary(md)

library(agricolae)

print(LSD.test(md,'trt',alpha=0.05))

detach(d)

detach(d)

rm(list=ls(all=TRUE))

#Factorial

d<-read.table('clipboard',header=TRUE) 

d

attach(d)

names (d)

blk<-factor(blk)

var<-factor(var)

fer<-factor(fer)

md<-aov(yield~blk+fer*var)

summary(md)

library(agricolae)

print(LSD.test(yield,var:fer,10,0.4,alpha=0.05))

 #Split plot

<-read.table('clipboard',header=TRUE) 

d

attach(d)

names (d)

blk<-factor(blk)

var<-factor(var)

fer<-factor(fer)

md1<-aov(yield~blk+var*fer+Error(blk:var))

summary(md1)

library(agricolae)

print(LSD.test(yield,var,4,0.417,alpha=0.05))

print(LSD.test(yield,fer,6,0.3,alpha=0.05))

detach(d)

detach(d)

rm(list=ls(all=TRUE))


#Rectangular lattice designs

d<-read.table('clipboard',header=T)
d
attach(d)
names(d)
library(agricolae)
print(BIB.test(blk,trt,yield))
detach(d)
detach(d)
rm(list=ls(all=TRUE))

#Correletion

d<-read.table('clipboard',header=T)

d

nx<-ncol(d)-1

E.x<-matrix(NA,nrow=nx,ncol=nx)

cor(d)

detach(d)

detach(d)

rm(list=ls(all=TRUE))