ฟังก์ชันสุ่ม(Random)ด้วยR

#CRD

library(agricolae)
trt <- c ('A' , 'B' , 'C' , 'D')
trt
rep <- c(3, 3, 3, 3)
rep
library(agricolae)
plan <- design.crd(trt, rep, number = 1, seed =1123)
plan
detach(d)
rm(list=ls(all=TRUE))


#RCBD

trt <- c ('1' , '2' , '3' , '4' , '5')
blk <- 4
library(agricolae)
plan2 <- design.rcbd(trt, blk, number = 1, seed =111)
plan2
detach(d)
detach(d)
rm(list=ls(all=TRUE))



#LS

trt <- c ('Sw1' , 'Sw2' , 'Sw3' , 'Sw5')
library(agricolae)
plan3 <- design.lsd(trt, number = 1, seed =111)
plan3
detach(d)
detach(d)
rm(list=ls(all=TRUE))



#Factorial

var<-c('A','B','c')
fer<-c('0','10')
library(agricolae)
plant5<-design.ab(var, fer, 2, number=1, seed=600)
plant5
detach(d)
detach(d)
rm(list=ls(all=TRUE))



#Split plot

var<-c('A','B','c')
fer<-c('0','10')
library(agricolae)
plant7<-design.split(var, fer, r=2, number=1, seed=600)
plant7
detach(d)
detach(d)
rm(list=ls(all=TRUE))


#triple lattice

library(agricolae)
trt<-c(1:132)
trt
outdesign<-design.lattice(trt,r=3,serie=3)
outdesign

# triple lattice

library(agricolae)
trt<-LETTERS[1:9]
design.lattice(trt,r=3,serie=2)

# simple lattice
trt<-1:100
plan<-design.lattice(trt,r=2,serie=3)

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))