7.4 - Latin Square Design (LSD)
7.4 - Latin Square Design (LSD)The fundamental idea of blocking can be extended to more dimensions. However, the full use of multiple blocking variables in a complete block design usually requires many experimental units. Latin Square design (LSD) can be useful when we want to achieve blocking simultaneously in two directions with a limited number of experimental units.
The limitation of the Latin Square experimental layout is that the design is only possible when
number of Row blocks = number of Column blocks = number of treatments.
The experimental design process begins with a Standard Latin Square. These have the treatment levels ordered across the first row and first column. For example, consider a single factor with three levels (A, B, C), blocked in two directions resulting in this standard \(3 \times 3\) square:
| A | B | C |
| B | C | A |
| C | A | B |
To randomize, first randomly permute the order of the rows and produce a new square.
| B | C | A |
| C | A | B |
| A | B | C |
Then randomly permute the order of the columns to yield the final square for the experimental layout.
| C | A | B |
| A | B | C |
| B | C | A |
Notice that each treatment occurs only once in each row and once in each column; the fundamental requirement of a LSD. This process assures that any row or column will have all treatment levels only once.
Using Technology
To obtain the design in SAS we can use:
proc plan;
factors Row=4 ordered Col=4 ordered / noprint;
treatments Treatment=4 cyclic;
output out=LatinSquare
Row cvals=('RowBlock 1' 'RowBlock 2' 'RowBlock 3' 'RowBlock 4') random
Col cvals=('ColBlock 1' 'ColBlock 2' 'ColBlock 3' 'ColBlock 4') random
Treatment nvals=(1 2 3 4) random;
run;
The ANOVA for the Latin Square is a direct extension of the RCBD with random blocking effects. The SAS random statement has to be modified accordingly to incorporate both blocking factors and with the assumption of no interaction between them (because of only one observation for each cell). For example, we could use the following SAS code to estimate the model:
proc mixed data=LatinSquare method=type3;
class Row Col Treatment;
model Response = Treatment;
random Row Col;
run;
To obtain a Latin Square Design for four treatments we can use the function rlatin from the magic package.
library(magic)
set.seed(1)
LSD4 = rlatin(4)
LSD4
[,1] [,2] [,3] [,4]
[1,] 4 2 3 1
[2,] 2 1 4 3
[3,] 1 3 2 4
[4,] 3 4 1 2