# Get Variance Components' Coefficients and Mean Squares for Single-Phase or Two-Phase Experiments Compute the variance components' coefficients and corresponding to random effects in the expected mean squares of ANOVA table in single-phase or two-phase experiments. These coefficients are then inserted to a matrix where the rows correspond to each source of variation and column correspond to DF and every variance component. The mean squares is calculated if the `response` argument is used. ```r getCoefVC.onePhase( Pb, design.df, v.mat, response, table.legend, decimal, digits ) ``` ## Arguments - `Pb`: a list of matrices generated by `infoDecompMat` function. - `design.df`: a data frame containing the experimental design. Requires every column be a `factor`. - `v.mat`: a list of matrix generated by `getVMat.onePhase` or `getVMat.twoPhase`. - `response`: a numeric vector contains the responses from the experiment. - `table.legend`: a logical allows users to generate a legend for the variance components of the ANOVA table for large designs. Default is `FALSE`, resulting in the use of original block factor names. - `decimal`: a logical allows users to display the coefficients as the decimals. Default is `FALSE`, resulting in the use of `fractions`. - `digits`: a integer indicating the number of decimal places. Default is 2, resulting in 2 decimal places. ## Returns A matrix containing the characters. ## Details The main purpose of this function is to combine the matrices presenting every source of variation of the ANOVA table and the variance matrix to compute the coefficients of the variance components. The complication arise in giving the row names of the matrix for the source of variation in the ANOVA table. ## Examples ```r design1 <- local({ Ani <- as.factor(LETTERS[c(1,2,3,4, 5,6,7,8)]) Trt <- as.factor(letters[c(1,1,1,1, 2,2,2,2)]) data.frame(Ani, Trt, stringsAsFactors = TRUE ) }) blk.str <- "Ani" rT <- terms(as.formula(paste("~", blk.str, sep = "")), keep.order = TRUE) blkTerm <- attr(rT,"term.labels") Z <- makeBlkDesMat(design1, blkTerm) trt.str <- "Trt" fT <- terms(as.formula(paste("~", trt.str, sep = "")), keep.order = TRUE) #fixed terms trtTerm <- attr(fT, "term.labels") effectsMatrix <- attr(fT, "factor") T <- makeContrMat(design1, trtTerm, effectsMatrix, contr.vec = NA) N = makeOverDesMat(design1, trtTerm) Replist = getTrtRep(design1, trtTerm) Rep <- Replist$Rep trt.Sca <- Replist$Sca effFactors = lapply(makeOrthProjectors(Z), function(z) getEffFactor(z, T, N, Rep, trt.Sca)) #Now construct variance matrices Pb <- effFactors[sort(1:length(effFactors), decreasing=TRUE)] v.mat <- getVMat.onePhase(Z.Phase1 = Z, design.df = design.df, var.comp = NA) getCoefVC.onePhase(Pb = Pb, design.df = design1, v.mat = v.mat, response = NA, table.legend = FALSE, decimal = FALSE, digit = 2) ``` ## Author(s) Kevin Chang

getCoefVC.onePhase function