A Latin Square is a one-factor design that extends blocking to two nuisance factors. For example, suppose we want to study the effects of five different formulations of rocket propellant. Suppose we also want to block (or average over) two nuisance factors:
1) Batch of raw material
2) Operator.
This is a 5 x 5 Latin Square and looks like this:
|
Raw material batch |
Operator |
||||
|
1 |
2 |
3 |
4 |
5 |
|
|
1 |
A |
B |
C |
D |
E |
|
2 |
B |
C |
D |
E |
A |
|
3 |
C |
D |
E |
A |
B |
|
4 |
D |
E |
A |
B |
C |
|
5 |
E |
A |
B |
C |
D |
This arrangement is a square and with each of the five treatments (a formulation of rocket propellant) represented by a Latin letter. Each letter occurs once for each operator and once for each batch of raw material. The Latin Square allows blocking of two nuisance factors. In this example we block on batch of raw material (the row) and operator (the column).
We will use this example to illustrate how to build and analyze a Latin square design using Design‑Expert® software.
* Douglas Montgomery, Design and Analysis of Experiments, 4th edition, John Wiley, Example 5-4 on page 196.
Procedure to create a 5x5 Latin Square design:
1. From the “File” menu choose “New Design…”
2. On the “Factorial” tab choose “General Factorial” and “3” Categorical Factors:
3. Name the first blocking factor (e.g. Row), enter the number of rows “5” and the block level names:
4. Name the second blocking factor (e.g. Column), enter the number of rows “5” and the block level names:
5. Name the treatment factor (e.g. Formulation), enter the number of rows “5” and the factor level names:
6.
7.
8. Right click on the column heading “Std” and choose “Sort by Standard Order”.
9. Right click on the column heading “Run” and choose “Re-order as currently displayed”.
10. Delete rows 26 – 125:
· Left click on the button to the left of row 26.
· While holding down the shift key, scroll down and left click the button to the left of the last row, i.e. row 125.
· Right click on a button to the left of any row in the selected area and choose “Delete Row[s]”:
“Are you sure you want to delete these rows??”
“If you delete or ignore these rows you will have only
one level of factor/component 3. Do you want to
delete/ignore these rows?”
11. Treatments (formulations) need to be assigned correctly. Type in the Latin letters that represent the five treatments (formulations) in a cyclic order. Starting at the top (row 1) type:
A, B, C, D, E B, C, D, E, A C, D, E, A, B D, E, A, B, C
and finally E, A, B, C, D
12. Right click on the column heading “Run” and choose “Randomize…”.
13. Save your design as “Example 5-4.dx6”.
14. Enter data.
· Right click on the column heading "Std" and choose "Sort by Standard Order"
· Type the following data into the response column. (You can copy / paste it if you are viewing this document online.)
| Std |
Response 1 |
| 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 |
24 17 18 26 22 20 24 38 31 30 19 30 26 26 20 24 27 27 23 29 24 36 21 22 31 |
15. Analyze the response:
· Click on the response node “Burn rate”.
· Click on the “Effects” button.
a. Right click on “A” and choose “Block”
b. Right click on “B” and choose “Block”
c. Right click on “C” and choose “Model”
16. Click on the “ANOVA” button:
Block term includes A, B
Analysis of variance table [Partial sum of squares]
Sum of Mean
F
Source Squares DF Square Value Prob
> F
Block 218.00 8 27.25
Model 330.00 4 82.50 7.73
0.0025
C 330.00 4 82.50 7.73 0.0025
Residual 128.00 12 10.67
Cor Total 676.00 24
17. Click on the “Diagnostics” button and check as usual:
18. Note the Box-Cox plot recommends an inverse transformation, this makes sense for a rate.
19. Reanalyze using the inverse transform:
Block term includes A, B
Analysis of variance table [Partial sum of squares]
Sum of
Mean F
Source Squares DF Square Value Prob
> F
Block 5.331E-004 8 6.664E-005
Model 9.072E-004 4 2.268E-004 11.73 0.0004
C 9.072E-004 4 2.268E-004 11.73 0.0004
Residual 2.321E-004 12 1.934E-005
Cor Total 1.672E-003 24
20. Click on the “Model Graphs” button:
Plot the treatment (Formulation) and average over the blocks (Rows and Columns):