Computes the canonical efficiency factors for the joint
decomposition of two structures or sets of mutually orthogonally
projectors (Brien and Bailey, 2009), orthogonalizing projectors in the
list to those earlier in the
list of projectors with
which they are partially aliased. The results can be summarized in the
form of a skeleton ANOVA table.
Two loops, one nested within the other. are performed. The first cycles
over the components of
Q1 and the nested loop cycles over the
Q2. The joint decomposition of the two projectors
in each cycle, one from
Q1[[i]]) and the other
Q2[[j]]) is obtained using
In particular, the nonzero canonical efficiency factors for the joint
decomposition of the two projectors is obtained. The nonzero canonical
efficiency factors are the nonzero eigenvalues of
Q1[[i]] %*% Q2[[j]] %*% Q1[[i]] (James and Wilkinson, 1971).
An eigenvalue is regarded as zero if it is less than
daeTolerance, which is initially set to
.Machine$double.eps ^ 0.5 (about 1.5E-08). The function
set.daeTolerance can be used to change
However, a warning occurs if any pair of Q2 projectors (say
Q2[[k]]) do not have adjusted orthgonality
with respect to any Q1 projector (say
Q1[[i]]), because they are
partially aliased. That is, if
Q2[[j]] %*% Q1[[i]] %*%
Q2[[k]] is nonzero for any pair of different Q2 projectors and any
Q1 projector. When it is nonzero, the projector for the later term in
the list of projectors is orthogonalized to the projector that is
earlier in the list.
list of class
p2canon. It has a component for each
Q1. Each of the components for
Q1 is a
list; its components are one for each component of
Q2 and a component
Pres. Each of the
list of three components:
Qproj. These components are based on an eigenalysis of the
relationship between the projectors for the parent
Q2 components. Each
pairwise component is based on the
nonzero canonical efficiency factors for the joint decomposition of the
two parent projectors (see
component is based on the nonzero canonical efficiency factors for the
joint decomposition of the
Q1 component and the
component, the latter adjusted for all
Q2 projectors that have
occured previously in the
Qproj component is the adjusted projector for the parent
Q2 component. The
have the following components:
– for details see
Brien, C. J. and R. A. Bailey (2009). Decomposition tables for multitiered experiments. I. A chain of randomizations. The Annals of Statistics, 36, 4184 - 4213.
James, A. T. and Wilkinson, G. N. (1971) Factorization of the residual operator and canonical decomposition of nonorthogonal factors in the analysis of variance. Biometrika, 58, 279-294.
projector for further information about this class.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
## PBIBD(2) from p. 379 of Cochran and Cox (1957) Experimental Designs. ## 2nd edn Wiley, New York PBIBD2.unit <- list(Block = 6, Unit = 4) PBIBD2.nest <- list(Unit = "Block") trt <- factor(c(1,4,2,5, 2,5,3,6, 3,6,1,4, 4,1,5,2, 5,2,6,3, 6,3,4,1)) PBIBD2.lay <- fac.layout(unrandomized = PBIBD2.unit, nested.factors=PBIBD2.nest, randomized = trt) ##obtain projectors using projs.structure Q.unit <- projs.structure(~ Block/Unit, data = PBIBD2.lay) Q.trt <- projs.structure(~ trt, data = PBIBD2.lay) ##obtain combined decomposition and summarize unit.trt.p2canon <- projs.2canon(Q.unit, Q.trt) summary(unit.trt.p2canon)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.