# NOT RUN {
#'## Microarray example from Jarrett & Ruggiero (2008) - see Brien (2019)
jr.lay <- fac.gen(list(Set = 7, Dye = 2, Array = 3))
jr.lay <- within(jr.lay,
{
Block <- factor(rep(1:7, each=6))
Plant <- factor(rep(c(1,2,3,2,3,1), times=7))
Sample <- factor(c(rep(c(2,1,2,2,1,1, 1,2,1,1,2,2), times=3),
2,1,2,2,1,1))
Treat <- factor(c(1,2,4,2,4,1, 2,3,5,3,5,2, 3,4,6,4,6,3,
4,5,7,5,7,4, 5,6,1,6,1,5, 6,7,2,7,2,6,
7,1,3,1,3,7),
labels=c("A","B","C","D","E","F","G"))
})
jr.anat <- designTwophaseAnatomies(formulae = list(array = ~ (Set:Array)*Dye,
plot = ~ Block/Plant/Sample,
trt = ~ Treat),
which.designs = c("first","cross"),
data = jr.lay)
## Three-phase sensory example from Brien and Payne (1999)
# }
# NOT RUN {
data(Sensory3Phase.dat)
Sensory.canon <- designTwophaseAnatomies(formulae = list(
eval= ~ ((Occasions/Intervals/Sittings)*Judges)/Positions,
field= ~ (Rows*(Squares/Columns))/Halfplots,
treats= ~ Trellis*Method),
data = Sensory3Phase.dat)
# }
Run the code above in your browser using DataCamp Workspace