getVMat.twoPhase {infoDecompuTE} | R Documentation |
Construct the matrix for each variance components for the single-phase or two-phase experiment.
getVMat.twoPhase(Z.Phase1, Z.Phase2, design.df, var.comp = NA)
Z.Phase1 |
a list of block design matrix from |
Z.Phase2 |
a list of block design matrix from |
design.df |
a data frame containing the experimental design. Requires
every column be a |
var.comp |
a vector of characters containing the variance components
of interest this allows the user to specify the variance components to be
shown on the ANOVA table. This also allows the user to specify artificial
stratum to facilitate decomposition. Default is |
A list of matrices.
Kevin Chang
design2 <- local({
Run = as.factor(rep(1:4, each = 4))
Ani = as.factor(LETTERS[c(1,2,3,4,
5,6,7,8,
3,4,1,2,
7,8,5,6)])
Tag = as.factor(c(114,115,116,117)[rep(1:4, 4)])
Trt = as.factor(letters[c(1,2,1,2,
2,1,2,1,
1,2,1,2,
2,1,2,1)])
data.frame(Run, Ani, Tag, Trt, stringsAsFactors = TRUE )
})
blk.str1 = "Ani"
blk.str2 = "Run"
rT1 = terms(as.formula(paste("~", blk.str1, sep = "")), keep.order = TRUE)
#random terms phase 1
rT2 = terms(as.formula(paste("~", blk.str2, sep = "")), keep.order = TRUE)
#random terms phase 2
blkTerm1 = attr(rT1,"term.labels")
blkTerm2 = attr(rT2,"term.labels")
Z1 = makeBlkDesMat(design2, rev(blkTerm1))
Z2 = makeBlkDesMat(design2, rev(blkTerm2))
V = getVMat.twoPhase(Z1, Z2, design2, var.comp = NA)