วันจันทร์ที่ 11 สิงหาคม พ.ศ. 2557

ฟังก์ชันสุ่ม(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)