Latin Square Design
Latin Square Design
Latin Square Design
between a CRD and a RCB design is essentially due to blocking, the value
of the relative efficiency is indicative of the gain in precision due to
blocking.
For our example, the R.E. value is computed as:
Because the error d.f. is only 15, the adjustment factor is computed as:
k= [(3)(5) + 1][6(3) + 31 =0.982
[(3)(5) + 3][6(3) + 11
= 1.60
The results indicate that the use of the RCB design instead of a CRD design
increased experimental precision by 60%.
The process of randomization and layout for a LS design is shown below for
an experiment with five treatments A, B, C, D, and E.
A B C D E
B A E C D
C D A E B
D E B A C
E C D B A
846 2 4
475 3 2
902 4 5
452 5 1
Use the rank to represent the existing row number of the selected plan
and the sequence to represent the r,,w number of the new plan. For our
example, the third row of the selected plan (rank = 3) becomes the first
row (sequence = 1) of the new plan; the fourth row of the selected plan
becomes the second row of the new plan; and so on. The new plan, after
the row randomization is:
C D A E B
D E B A C
B A E C D
E C D B A
A B C D E
The rank will now be used to represent the column number of the plan
obtained in step 2 (i.e., with rearranged rows) and the sequence will be used
to represent the column number of the final plan.
For our example, the fourth column of the plan obtained in step 2
becomes the first column of the final plan, the first column of the plan of
step 2 becomes the second column of the final plan, and so on. The final
Latin Square Design 33
There are four sources of variation in a LS de-ign, two more than that for the
CRD and one more than that for the RCB design. The sources of variation are
row, column, treatment, and experimental error.
To illustrate the computation procedure for the analysis of variance of a LS
design, we use data on grain yield of three promising maize hybrids (A, B,and
D) and of a check (C) from an advanced yield trial with a 4 X 4 latin square
design (Table 2.7).
The step-by-step procedures in the construction of the analysis of variance are:
o STEP 1. Arrange the raw data according to their row and column designa
tions, with the corresponding treatment clearly specified for each observa
tion, as shown in Table 2.7.
o STEP
2. Compute row totals (R), column totals (C), and the grand total
(G) as shown in Table 2.7. Compute treatment totals (T) and treatment
Total d.f. = 12 - 1 = 16 - 1 = 15
Error d.f. = Totad d.f. - Row d.f - Column d.f .- Treatment d.f.
= 15 - 3-33-3 = 6
0 smrP 5. Compute the correction factor and the various sums of squares as:
C.F.G
"'(21.365)2 . 28.528952
16
Latin Square Design 35
Total SS _ZX 2-
C.F.
-1.413923
Row SS t - C.F.
Column SS = X - C.F.
I
= 0.827342
Treatment SS = T - C.F.
I
= 0.129585
1sup 6. Compute the mean square for each source of variation by dividing
the sum of squares by its corresponding degree of freedom:
Row MS = Row SS
t-1
0.030154_
0.03015 = 0.010051
3
36 Ningle.Factor Experiments
Column SS
Column MS C- t SS
0.827342 0.275781
___ --
3
Treatment SS
Treatment MS = T t- 1
0.426842
= 3 = 0.142281
Error MS = Error SS
(t- 1)(t- 2)
0.129585
= (3)(2)
= 0.021598
(3)(2)
3 STEP 7. Compute the F value for testing the treatment effect as:
Treatment MS
Error MS
0.142281
0.021598
3 STEP 8. Compare the computed F value with the tabular F value, from
Appendix E, with f, = treatment d.f. = t - 1 and f2 = error d.f. =
(t - 1)(t - 2) and make conclusions following the guidelines in step 9 of
Section 2.1.2.1.
For our example, the tabular F values, from Appendix E, with f, = 3 and
f2 = 6 degrees of freedom, are 4.76 at the 5% level of significance and 9.78
at the 1%level. Because the computed F value is higher than the tabular F
value at the 5% level of significance but lower than the tabular F value at the
1% level, the treatment difference is significant at the 5% level of signifi
cance.
/Error MS
Grand mean
tested, it does not identify the specific pairs or groups of varieties that
differed. For example, the F test is not able to answer the question of
whether every one of the three hybrids gave significantly higher yield than
that of the check variety or whether there is any significant difference among
the three hybrids. To answer these questions, the procedures for mean
comparisons discussed in Chapter 5 should be used.
As in the RCB design, where the efficiency of one-way blocking indicates the
gain in precision relative to the CRD (see Section 2.2.4), the efficiencies of both
row- and column-blockings in a LS design indicate the gain in precision
relative to either the CRD or the RCB design. The procedures are:
C SiuP 1. Test the level of significance of the differences among row- and
column-blocks:
A. Compute the F values for testing the row difference and column
difference as:
Row MS
F(row) = Error MS
0.010051
- <1
0.021598
F(column) = Column MS
Error MS
= 0.275781
0.021598
B. Compare each of the computed F values that is larger than I with
the tabular F values (from Appendix E) with f, = t - 1 and f2 =
It - lXt - 2) degrees of freedom. For our example, the computed
38 Single-Factor Experiments
F(row) value is smaller than 1 and, hence, is not significant. For the
computed F(column) value, the corresponding tabular F values with
f= 3 and 12 = 6 degrees of freedom are 4.76 at the 5% level of
significance and 9.78 at the 1%level. Because the computed F(column)
value is greater than both tabular F values, the difference among
column-blocks is significant at the 1%level. These results indicate the
success of column-blocking, but not that of row-blocking, in reducing
experimental error.
where E, is the row mean square, E, is the column mean square, and E, is
the error mean square in the LS analysis of variance; and t is the number
of treatments.
For our example, the R.E. is computed as:
RE(CRD) = 0.010051 + 0.275781 + (4 - 1)(0.021598)
(4 + 1)(0.021598)
- 3.25
where E,, Eo, Er, and t are as defined in the preceding formula.
Lattice Design 39
When the error d.f. in the LS analysis of variance is less than 20, the
R. E. value should be multiplied by the adjustment factor k defined as:
3]
[(t- 1)(t- 2) + 1][(t- 1)2 +
[(t- 1)(t- 2) + 31[(- 1)2 + 1]
For our example, the values of the relative efficiency of the LS design
compared to a RCB design with rows as blocks and with columns as
blocks are computed as:
R.E.(RCB, row) = 0.010051 + (4 - 1)(0.021598)
4(0.021598)
= 0.87
= 3.94
Because the error d.f. of the LS design is only 6, the adjustment factor k
is computed as: