summary.p2canon {dae} | R Documentation |
Produces a summary of the efficiency criteria computed from the
canonical efficiency factors for the joint decomposition of two
sets of projectors (Brien and Bailey, 2009) obtained using
projs.2canon
. It takes the form of a decomposition or skeleton
ANOVA table.
## S3 method for class 'p2canon'
summary(object, which.criteria = c("aefficiency", "eefficiency", "order"), ...)
object |
A |
which.criteria |
A character |
... |
further arguments affecting the summary produced. |
An object of classes summary.p2canon
and data.frame
, whose
rows correspond to the pairs of projectors, one from the
Q1
argument and the other from the Q2
argument from
projs.2canon
; only pairs with non-zero efficiency factors
are included. In addition, a line is included for each nonzero Residual
Q1
projector.
Chris Brien
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.
projs.2canon
, proj2.efficiency
,
efficiency.criteria
, proj2.combine
,
proj2.eigen
, pstructure
,
print.summary.p2canon
in package dae,
eigen
.
projector
for further information about this class.
## 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 <- designRandomize(allocated = trt,
recipient = PBIBD2.unit,
nested.recipients = PBIBD2.nest)
##obtain projectors using pstructure
unit.struct <- pstructure(~ Block/Unit, data = PBIBD2.lay)
trt.struct <- pstructure(~ trt, data = PBIBD2.lay)
##obtain combined decomposition and summarize
unit.trt.p2canon <- projs.2canon(unit.struct$Q, trt.struct$Q)
summary(unit.trt.p2canon)