In some experiment designs you must restrict the randomization. Otherwise it wouldn’t be practical to perform the runs. For example, temperature cannot be easily changed in many applications. The solution may be a “split-plot” design, which originated in the field of agriculture. Experimenters divided large areas of land, called “whole plots,” into smaller “subplots” that could be treated separately. For example, they planted several crop varieties—one each per whole plot, and then applied different fertilizers on each subplot.

The analysis of a split-plot design is tricky, even for statisticians. But
Design Expert^{®} version 13 software makes it easy.

To illustrate how Design-Expert does design and analysis of split plots, let’s follow an example from Montgomery’s Design and Analysis of Experiments. Feel free to skip the boxed bits—these are sidebars for those who want to know all the details and spend more time to explore things.

Note

**Details on this case study**: Refer to Section 14.4 “The
Split-Plot Design” in the 8th Edition of the book referenced above—published by
Wiley in 2012—for more details on the experiment’s structure and this particular
case.

A manufacturer wants to determine the effect of three methods of pulp preparation and four levels of cooking temperature on the tensile strength of paper. These twelve combinations can be conducted in one day.

To increase statistical power, the paper chemist decides to perform three replicates of this twelve-run general factorial design over three days, which produces a total of 36 runs. However, it is not feasible to perform the experiment in a random fashion. Instead, on each of the three days, this experimenter produces a batch of pulp and divides it into four samples to be cooked at four different temperatures. The process is repeated for pulp prepared by the two other methods.

A flowchart of the experiment is shown below.

In statistical terms, the split plot experiment can be structured as:

- Whole plots for the three batches of pulp (hard-to-change factor)
- Subplots for the four samples cooked at four different temperatures (easy to change factor).

You will set up this design as a blocked (by day) split-plot general factorial.

Let’s build the design. This is done via the usual procedure for general
factorials, with days accounted for as blocks. From the Design-Expert welcome
screen, click the **New Design Button** (or click the blank-sheet icon
on your toolbar, or choose File, New Design).

Then, from the default **Factorial** tab under **Split-Plot** click **Multilevel
Categoric**. Choose **2** as the number of factors. Click the **Vertical** radio
button. Then enter **Pulp Prep** as your Factor a name. It’s already set up as a
Hard-to-Change factor (HTC), so leave that be (note that HTC factors are lower
cased). Mouse down to the Levels column and type in **3**. Then enter **a1, a2**,
and **a3** as your **L(1), L(2)**, and **L(3)** treatment names.

Still on the same screen, name your second factor **Temp**, with **deg F** for
units, **4** for the number of levels, and **200, 225, 250**, and **275** as your
treatments. This factor is Easy to Change as already set by the program. For
**Type** select **Ordinal**, thus denoting its numeric nature. Your screen should
now match that shown below.

Note the warning due to lack of power on the HTC factor. To remedy this, for
**Replicates** enter **3** and check on (√) **Assign one block per replicate**
to produce 36 runs (36=3x3x4 levels from the three factors) in 9 groups.

Press **Next**. Now for **Block Names** enter **Day 1, Day 2** and
**Day 3**.

Press **Next** again and for the one response **Name** enter **Tensile** in
**Units** of **strength**. Then to assess the power of this design, type in
**Delta** of **6** and **Sigma** of **2**. Also, change the **variance ratio**
to **0.25**.

Next click **Edit Model…** and change the **Process Order** to **2FI** as shown
below and press **OK**.

Note

Design-Expert defaults to a variance ratio of 1 for group (whole plot) to residual variance. But in this case, the whole-plot factor, although hard to change, can be controlled with little variation relative to the sub-plot factor of temperature.

P.S. The changing of the order presumes that the experimenters must be interested in the interaction of the pulp prep with temperature, otherwise why bother doing it at all. Therefore they want to know if this model term (aB) can be estimated with a ‘decent’ (80% being generally acceptable) level of power.

Press **Next** and review the power for the design. Notice that, due to the
restriction in randomization in the split plot, the power for seeing effects of
the hard-to-change factor “a” (Pulp prep) is reduced, whereas the power for
seeing the effects of the easy-to-change factor “B” (Temperature) and its
interaction “aB” are increased, albeit only slightly.

Note

Due to the structure of this design, interactions of hard to change (whole plot) factors with ones that are easy to change (sub-plot) must be tested against the sub-plot error.

Notice the loss in power for the hard-to-change factor “a” (Pulp Prep) still
leaves it above the recommended minimum. **Finish** to complete the
design-building wizard. Press **OK** on the warning to reset factors
group-to-group. Design-Expert now displays your design in groups by level of the
hard-to-change factor “a” (Pulp prep), within which the easy-to-change factor “B”
(Temperature) is laid out in completely random order.

At this stage the design would normally be printed as a run-sheet for the actual experiment, with resulting responses of tensile strength recorded.

Load the design by clicking **Help, Tutorial Data -> Paper Pulp**.

The analysis of a split-plot design must be done differently than a completely randomized experiment. To get the correct statistical tests for each, Design-Expert creates two separate ANOVAs for:

- Whole-plot treatment (pulp prep), and
- Subplot treatment (temperature) and the interaction between whole plot and subplot individually.

Then it fits the full model in order to get meaningful diagnostics and model graphs.

Start by clicking the Analysis node labeled **R1:Tensile** which can be found in
the tree structure along the left of the main window. Then press ahead to the
**Model** tab displayed in the progressive toolbar. Design-Expert preselects
terms for the designed-for two-factor interaction (2FI) model—these are labeled
with a green . Enable the program to check if any terms can be safely
taken out of the model by clicking the **Auto Select…** button. Then, change the
**criterion** to **p-values** and the **Selection** to **Backward** as shown in
the screen shot below.

Click the **Start** button and the selection algorithm is detailed with any terms
that exceed the “Alpha out” criterion (p-value of 0.1 by default) being removed
from the model. They will be re-designated to Error () after you click
**Accept**. If after looking over the statistics on the ANOVA, you want to
remove a term, simply go back and double-click it (or right-click for a menu of
designation choices).

Note

**Heads-up on effects selection for multilevel categoric split plots**:
Due to the structure of these designs, the handy graphical half-normal plot of
effects is not produced like it would be for an experiment with factors only
at two levels. Instead, run one of the Auto Select… algorithms (we used
backward in this case)—details for which can by clicking Help button.

Continue to the **ANOVA (REML)** tab and review the p-values for significance.
In this case all the designed-for model terms survive the elimination process,
that is, they are all statistically significant.

Note

**A note about terminology**: As you probably know, ANOVA is the
acronym for analysis of variance—the fundamental approach to assessing model
significance for all DOE where runs can be completely randomized. However, in
a split plot the components of variance are estimated by a different approach
called REML, which stands for restricted maximum likelihood estimation.
Sometimes, depending on the structure of the design, REML agrees with
ordinary-least-squares ANOVA exactly, but often the results differ somewhat.

Note

**Details on statistics presented in ANOVA (REML) report**: Search
Help, Contents for much more information, including references, on everything
you see presented on this screen. Since these tutorials only provide the
highlights, we must move on.

Move on to the next phase of analysis—verifying the statistical assumptions—by
pressing the **Diagnostics** tab. First up is the normal plot of the residuals—
externally studentized to provide the most effective view of any obvious
discrepancies. There are no obvious irregularities in this case—the points fall
in line as expected from normal data.

Examine all the graphs available via the Diagnostics tool. You should see nothing abnormal.

Press ahead to **Model Graphs**. On the **Factors Tool** right-click the
**B: Temperature** bar and change it to the **X1 axis**.

You now see an interaction plot of the continuous factor temperature at three discrete levels of pulp preparation. From this graph we can see pulp preparation a2 (symbolized by the green triangles) will maximize tensile strength. Evidently it will be most robust to temperature changes, especially in the temperature range of 225 to 275 degrees F.

This is an ingenious application of designed experimentation using a split-plot structure to handle the hard-to-change pulp preparation. Design-Expert provides valuable insights on how power is reduced by restricting randomization. In this case the actual effect of pulp prep (factor “a”) far exceeded that of minimal importance so the experiment succeeded.