Special considerations designing fractional factorials.
library(ggplot2)
theme424 <- theme_minimal() +
theme(
panel.grid.minor = element_blank(),
panel.background = element_rect(fill = "#f7f7f7"),
panel.border = element_rect(fill = NA, color = "#0c0c0c", size = 0.6),
axis.text = element_text(size = 14),
axis.title = element_text(size = 16),
legend.position = "bottom"
)
theme_set(theme424)
include_graphics("https://uwmadison.box.com/shared/static/b94hfeyugeun1xywlm2ccls9am5l38xg.jpg")
Projection makes it possible to turn a fractional factorial screening experiment into a full factorial experiment. Effectively, we delete the columns of our design matrix involving the factors outside the ones we want to project down to. Depending on what design we project, a \(2^{K - p}\) design will become either (1) a full factorial design or (2) a \(2^{K - r}\) fractional factorial design, where \(r < p\).
Beware that conclusions will be somewhat tentative, because effects in the full factorial could be aliased with high-order effects from factors that have been projected out.
Warm-up: (Projecting a \(2 ^ {4 - 1}\) Design). Consider the \(2 ^ {4 - 1}\) design with defining relation \(I = ABCD\).
A | B | C | D |
---|---|---|---|
- | - | - | - |
+ | - | - | + |
- | + | - | + |
- | - | + | + |
+ | + | - | - |
+ | - | + | - |
- | + | + | - |
+ | + | + | + |
What happens when we project out A?
B | C | D |
---|---|---|
- | - | - |
- | - | + |
+ | - | + |
- | + | + |
+ | - | - |
- | + | - |
+ | + | - |
+ | + | + |
This design includes all corners of the cube B-C-D, so it’s a full \(2^{3}\) design.
A | B | C | D | E = ABC | F = BCD |
---|---|---|---|---|---|
- | - | - | - | - | - |
+ | - | - | - | + | - |
- | + | - | - | + | + |
- | - | + | - | + | + |
- | - | - | + | - | + |
+ | + | - | - | - | + |
+ | - | + | - | - | + |
+ | - | - | + | + | + |
- | + | + | - | - | - |
- | + | - | + | + | - |
- | - | + | + | + | - |
+ | + | + | - | + | - |
+ | + | - | + | - | - |
+ | - | + | + | - | - |
- | + | + | + | - | + |
+ | + | + | + | + | + |
This study examined 6 factors over 16 runs. Its complete defining relations are \(I = ABCE = BCDF = ADEF\). What happens when we project out \(A\)?
B | C | D | E | F | |
---|---|---|---|---|---|
- | - | - | - | - | |
- | - | - | + | - | |
+ | - | - | + | + | |
- | + | - | + | + | |
- | - | + | - | + | |
+ | - | - | - | + | |
- | + | - | - | + | |
- | - | + | + | + | |
+ | + | - | - | - | |
+ | - | + | + | - | |
- | + | + | + | - | |
+ | + | - | + | - | |
+ | - | + | - | - | |
- | + | + | - | - | |
+ | + | + | - | + | |
+ | + | + | + | + |
C | D | E | F | ||
---|---|---|---|---|---|
- | - | - | - | ||
- | - | + | - | ||
- | - | + | + | ||
+ | - | + | + | ||
- | + | - | + | ||
- | - | - | + | ||
+ | - | - | + | ||
- | + | + | + | ||
+ | - | - | - | ||
- | + | + | - | ||
+ | + | + | - | ||
+ | - | + | - | ||
- | + | - | - | ||
+ | + | - | - | ||
+ | + | - | + | ||
+ | + | + | + |
B | C | D | F | ||
---|---|---|---|---|---|
- | - | - | - | ||
- | - | - | - | ||
+ | - | - | + | ||
- | + | - | + | ||
- | - | + | + | ||
+ | - | - | + | ||
- | + | - | + | ||
- | - | + | + | ||
+ | + | - | - | ||
+ | - | + | - | ||
- | + | + | - | ||
+ | + | - | - | ||
+ | - | + | - | ||
- | + | + | - | ||
+ | + | + | + | ||
+ | + | + | + |
Exercise: Try projecting onto \(A, B, C\). Is it a fractional or full factorial? How many replicates each? Could you have determined this, just by looking at the defining relations?
A design is said to have projectivity \(v\) if collapsing to any subset of \(v\) factors will result in a full factorial design. The example above has projectivity 3.
Let’s say we’ve subsetted to a set of runs corresponding to a fractional factorial design. But, we still have to break the runs into batches, we have too many runs to be able to reasonably run them in one block. How can we appropriately block a fractional factorial design?
The typical solution is to choose a high-order effect and alias it with the batches. This divides the runs into two blocks, those where the effect is + and those where it is -.
We should ensure that the alias group doesn’t contain any low-order effects of interest (e.g., we don’t want the blocks to be confounded with any main effects). We effectively sacrifice the alias group that is confounded with the block.
Example: For our \(2 ^{6 - 2}\) design above, we can confound \(ACD\) with the block, to avoid aliasing any main effects with the block effect. This let us divide our runs into two blocks of size 8.
A | B | C | D | E | F | Block = ACD |
---|---|---|---|---|---|---|
- | - | - | - | - | - | - (B1) |
+ | - | - | - | + | - | + (B2) |
- | + | - | - | + | + | - (B1) |
- | - | + | - | + | + | + (B2) |
- | - | - | + | - | + | + (B2) |
+ | + | - | - | - | + | + (B2) |
+ | - | + | - | - | + | - (B1) |
+ | - | - | + | + | + | - (B1) |
- | + | + | - | - | - | + (B2) |
- | + | - | + | + | - | + (B2) |
- | - | + | + | + | - | - (B1) |
+ | + | + | - | + | - | - (B1) |
+ | + | - | + | - | - | - (B1) |
+ | - | + | + | - | - | + (B2) |
- | + | + | + | - | + | - (B1) |
+ | + | + | + | + | + | + (B2) |