# 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)
print(crd)
crd
#! Factorial Experiment (FE) - CRD
# 5 factors (f1..f5, level^factor): 3^1 * 2^1 * 3^1 * 2^1 * 4^1 * 2 = 360 experimental units
# 5 factors (f1..f5, level^factor): 3^2 * 2^2 * 5^1 = 180 * 2 = 360 experimental units
crd_fe <- gexp(mu=0,
r=2,
fe=list(f1=c(1, 1, 5),
f2=c(1, 1),
f3=c(2, 2, 1),
f4=c(1, 5),
f5=c(1, 2, 3, 4, 5)),
fl=list(A=paste('a',
1:3,
sep=''),
B=paste('b',
1:2,
sep=''),
C=paste('c',
1:3,
sep=''),
D=paste('d',
1:2,
sep=''),
E=paste('e',
1:5,
sep='')),
inte = rep(1,848),
type = 'FE',
round=0)
print(crd_fe)
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')
print(split_crd)
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')
print(rcbd)
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')
print(fe_rcbd)
fe_rcbd
#! Split-plot Experiment (SPE) - RCBD
split_rcbd <- gexp(mu=30,
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),
blkl=list(BLK=paste('B',
1:2,
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')
print(split_rcbd)
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')
print(lsd)
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')
print(fe_lsd)
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')
print(split_lsd)
split_lsd
# }
Run the code above in your browser using DataCamp Workspace