# NOT RUN {
#! Completely Randomized Design (CRD)
#! 1 factor - CRD
crd <- gexp(mu=1,
r=3,
fe=list(f1=c(1, 1, 5, 1, 1)),
fl=list(Treat=LETTERS[1:5]),
round=2)
summary(crd)
plot(crd) # Default
plot(crd, # Changing some arguments
main='',
sub='',
coltext='black',
colgrid='darkred',
ltygrid='solid',
lwdgrid=3)
#! Factorial Experiment (FE) - CRD
# 2 factors (f1,f2, level^factor): 3^1 * 2^1 = 6 experimental units
crd_fe <- gexp(mu=0,
r=2,
fe=list(f1=c(1, 1, 5),
f2=c(1, 1)),
fl=list(A=paste('a',
1:3,
sep=''),
B=paste('b',
1:2,
sep='')),
round=0,
inte = rep(1, 6),
type = 'FE')
summary(crd_fe)
plot(crd_fe)
#! Split-plot Experiment (SPE) - CRD
split_crd <- gexp(mu=30,
fe=list(f1=c(1, 1),
f2=c(2, 3)),
fl=list(P=paste('p',
1:2,
sep=''),
SP=paste('sp',
1:2,
sep='')),
inte=c(1, 15, 1, 1),
round=1,
type='SPE')
summary(split_crd)
plot(split_crd)
#! Randomized Complete Block Design (RCBD)
# 1 factor, 3 blocks
rcbd <- gexp(mu=0,
fe=list(f1=c(5, 1, 1)),
fl=list(TR=LETTERS[1:3]),
blke=c(1, 2, 3),
blkl=list(BLK=paste('B',
1:3,
sep='')),
round=1,
type='RCBD')
summary(rcbd)
plot(rcbd)
#! Factorial Experiment (FE) - RCBD
fe_rcbd <- gexp(mu=30,
fe=list(f1=c(1, 1, 1),
f2=c(2, 3)),
blke=c(1, 3),
inte=c(1, 15, 1, 1, 5, 1),
round=1,
type='FE')
summary(fe_rcbd)
plot(fe_rcbd)
fe_rcbd1 <- gexp(mu=30,
fe=list(f1=c(1, 1, 1),
f2=c(2, 3)),
blke=c(1, 3),
blkl=list(Blk=c('B1', 'B2')),
inte=c(1, 15, 1, 1, 5, 1),
round=1,
type='FE')
summary(fe_rcbd1)
plot(fe_rcbd1)
#! Split-plot Experiment (SPE) - RCBD
split_rcbd <- gexp(mu=30,
r = 2,
fe=list(f1=c(1, 1),
f2=c(2, 3),
f3=c(1, 1, 1)),
fl=list(A=paste('a',
1:2,
sep=''),
B=paste('b',
1:2,
sep=''),
C=paste('c',
1:3,
sep='')),
blke=c(1, 2, 3),
blkl=list(BLK=paste('B',
1:3,
sep='')),
inte=c(1, 15, 1, 1, 1, 3, 4, 2, 1, 1, 4, 1,
1, 2, 1, 1,
1, 1, 1, 1, 1, 1,
1, 1, 3, 3, 3, 3),
round=1,
type='SPE')
summary(split_rcbd)
plot(split_rcbd)
#! Latin Square Design (LSD)
lsd <- gexp(mu=30,
fe=list(f1=c(1, 1, 10)),
rowe=c(1, 1, 1),
cole=c(1, 1, 1),
rowl=list(Row=paste('r',
1:3,
sep='')),
coll=list(Col=paste('c',
1:3,
sep='')),
round=1,
type='LSD')
summary(lsd)
plot(lsd)
#! Factorial Experiment (FE) - LSD
fe_lsd <- gexp(mu=30,
fe=list(f1=c(1, 1),
f2=c(2, 3)),
rowe=c(1, 3, 2, 1),
cole=c(2, 2, 1, 1),
rowl=list(Row=paste('r',
1:4,
sep='')),
coll=list(Col=paste('c',
1:4,
sep='')),
inte=c(1, 15, 1, 1),
round=1,
type='FE')
summary(fe_lsd)
plot(fe_lsd)
#! Split-plot Experiment (SPE) - LSD
split_lsd <- gexp(mu=30,
fe=list(f1=c(1, 1, 2),
f2=c(2, 3, 1)),
fl=list(P=paste('p',
1:3,
sep=''),
SP=paste('sp',
1:3,
sep='')),
inte=c(1, 15, 1, 1, 1, 1, 1, 1, 1),
rowe = c(1, 1, 1),
cole = c(1, 1, 1),
rowl=list(Row=paste('r',
1:3,
sep='')),
coll=list(Col=paste('c',
1:3,
sep='')),
round=1,
type='SPE')
summary(split_lsd)
plot(split_lsd)
# }
# NOT RUN {
#! Using images in plannig
# CRD
crd <- gexp(mu=1,
r=2,
fe=list(f1=c(1, 1, 5, 1, 1)),
fl=list(Treat=LETTERS[1:5]),
round=2)
summary(crd)
plot(crd) # Default
# Dynamic
# Open picture 'crd.jpg' whem requested
crd_i <- update(crd,
r=3,
fe=list(f1=c(1, 1)),
fl=list(Treat=LETTERS[1:2]))
summary(crd_i)
plot(crd_i) # Default
# dynamic plot require 'jpeg' package!
plot(crd_i,
dynamic=TRUE,
xleftimg=0.6,
ybottomimg=.6,
xrightimg=1.4,
ytopimg=1.4)
# RCBD
rcbd <- gexp(mu=0,
fe=list(f1=c(5, 1, 1)),
fl=list(TR=LETTERS[1:3]),
blke=c(1, 2, 3),
blkl=list(BLK=paste('B',
1:3,
sep='')),
round=1,
type='RCBD')
summary(rcbd)
plot(rcbd)
rcbd_i <- update(rcbd,
r=1,
fe=list(f1=c(1, 1, 1, 1)),
blke=c(1, 1, 1))
plot(rcbd_i,
dynamic=TRUE)
# LSD
lsd <- gexp(mu=30,
fe=list(f1=c(1, 1, 10)),
rowe=c(1, 1, 1),
cole=c(1, 1, 1),
rowl=list(Row=paste('r',
1:3,
sep='')),
coll=list(Col=paste('c',
1:3,
sep='')),
round=1,
type='LSD')
summary(lsd)
#TODO: dynamic plot require 'png' package!
plot(lsd)
plot(lsd,
dynamic=TRUE)
# }
Run the code above in your browser using DataLab