infoDecompuTE (version 0.5.1)

getVMat.twoPhase: Get Variance Matrix for Single-phase Experiment

Description

Construct the matrix for each variance components for the single-phase experiments.

Usage

getVMat.twoPhase(Z.Phase1, Z.Phase2, design.df, var.comp = NA)

Arguments

Z.Phase1

a list of block design matrix from makeBlkDesMatrix function from Phase 1 block structure.

Z.Phase2

a list of block design matrix from makeBlkDesMatrix function from Phase 2 block structure.

design.df

a data frame containing the design of the experiment. Every column must be factors.

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 NA, which uses every random factor as the variance components with the first phase's variance components in random.terms1 appear before the second phase's variance components in random.terms2.

Value

This function returns a list of matrices.

Examples

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)
})

    blk.str1 = "Ani"
     blk.str2 = "Ani"
   
	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 = makeBlkDesMatrix(design2, rev(blkTerm1))
		Z2 = makeBlkDesMatrix(design2, rev(blkTerm2))

   V = getVMat.twoPhase(Z1, Z2, design2, var.comp = NA)