The yields (kg/plot) of three fertilizers, each applied to four plots, are recorded below. Test whether the three fertilizers differ significantly in mean yield.
| F1 | 22 | 26 | 24 | 28 |
|---|---|---|---|---|
| F2 | 30 | 34 | 32 | 36 |
| F3 | 20 | 22 | 24 | 22 |
To test \(H_0:\mu_1 = \mu_2 = \mu_3\) against \(H_1:\) not all equal, using one-way analysis of variance with equal group sizes.
Applying it:
Blank ANOVA table:
| Source | SS | df | MS | F |
|---|---|---|---|---|
| Treatments | ||||
| Error | ||||
| Total | — | — |
Totals: \(T_1 = 100,\; T_2 = 132,\; T_3 = 88\); \(G = 320,\; N = 12,\; CF = 320^2/12 = 8533.33\).
\(\sum x^2 = 8840\Rightarrow SS_T = 306.67\); \(SS_{Tr} = (100^2+132^2+88^2)/4 - CF = 8792 - 8533.33 = 258.67\); \(SS_E = 48\).
| Source | SS | df | MS | F |
|---|---|---|---|---|
| Treatments | 258.67 | 2 | 129.33 | 24.25 |
| Error | 48.00 | 9 | 5.33 | — |
| Total | 306.67 | 11 | — | — |
\(F = 24.25 > F_{0.05}(2,9) = 4.26\), so \(H_0\) is rejected: the three fertilizers differ significantly in mean yield.
Three varieties are grown on different numbers of plots: V1 (3 plots): 20, 24, 22; V2 (4 plots): 28, 30, 32, 34; V3 (2 plots): 18, 20. Test whether the variety means differ.
To perform one-way ANOVA with unequal group sizes and test \(H_0:\mu_1=\mu_2=\mu_3\).
Applying it:
Blank ANOVA table:
| Source | SS | df | MS | F |
|---|---|---|---|---|
| Treatments | ||||
| Error | ||||
| Total | — | — |
Totals: \(T_1 = 66\,(n_1=3),\; T_2 = 124\,(n_2=4),\; T_3 = 38\,(n_3=2)\); \(G = 228,\; N = 9,\; CF = 228^2/9 = 5776\).
\(\sum x^2 = 6048\Rightarrow SS_T = 272\); \(SS_{Tr} = (66^2/3 + 124^2/4 + 38^2/2) - 5776 = 6018 - 5776 = 242\); \(SS_E = 30\).
| Source | SS | df | MS | F |
|---|---|---|---|---|
| Treatments | 242 | 2 | 121 | 24.2 |
| Error | 30 | 6 | 5 | — |
| Total | 272 | 8 | — | — |
\(F = 24.2 > F_{0.05}(2,6) = 5.14\), so \(H_0\) is rejected: the varieties differ significantly.
Yields are classified by three treatments (rows) and four blocks (columns), one observation per cell. Test whether treatment means and block means differ.
| B1 | B2 | B3 | B4 | |
|---|---|---|---|---|
| T1 | 22 | 26 | 28 | 24 |
| T2 | 30 | 36 | 34 | 32 |
| T3 | 20 | 24 | 22 | 26 |
To perform two-way ANOVA (without replication) and test the equality of treatment (row) effects and of block (column) effects.
Applying it:
Blank ANOVA table:
| Source | SS | df | MS | F |
|---|---|---|---|---|
| Treatments (rows) | ||||
| Blocks (columns) | ||||
| Error | — | |||
| Total | — | — |
Row totals \(R = 100, 132, 92\); column totals \(C = 72, 86, 84, 82\); \(G = 324,\; N = 12,\; CF = 324^2/12 = 8748\).
\(\sum x^2 = 9032\Rightarrow SS_T = 284\); \(SS_R = 35888/4 - 8748 = 224\); \(SS_C = 26360/3 - 8748 = 38.67\); \(SS_E = 284 - 224 - 38.67 = 21.33\).
| Source | SS | df | MS | F |
|---|---|---|---|---|
| Treatments (rows) | 224.00 | 2 | 112.00 | 31.50 |
| Blocks (columns) | 38.67 | 3 | 12.89 | 3.63 |
| Error | 21.33 | 6 | 3.56 | — |
| Total | 284.00 | 11 | — | — |
Treatments: \(F_R = 31.50 > F_{0.05}(2,6) = 5.14\) ⇒ treatment effects differ significantly. Blocks: \(F_C = 3.63 < F_{0.05}(3,6) = 4.76\) ⇒ block effects are not significant.
Using the completely randomised design of Experiment 1 (3 fertilizers × 4 replicates, \(MS_E = 5.33,\; df_E = 9\)), determine the critical difference and identify which pairs of fertilizer means differ significantly.
To compute the critical difference (CD) for a CRD and compare all pairs of treatment means.
Applying it:
Blank working table:
| Pair | \(|\bar y_i - \bar y_j|\) | vs CD | Decision |
|---|---|---|---|
| F1 vs F2 | |||
| F1 vs F3 | |||
| F2 vs F3 |
\(t_{0.025,9} = 2.262\), so \(CD = 2.262\sqrt{2(5.33)/4} = 2.262\times 1.633 = 3.69\).
Means: \(\bar y_1 = 25,\; \bar y_2 = 33,\; \bar y_3 = 22\).
| Pair | \(|\bar y_i - \bar y_j|\) | vs CD (3.69) | Decision |
|---|---|---|---|
| F1 vs F2 | 8 | > | Differ |
| F1 vs F3 | 3 | < | No difference |
| F2 vs F3 | 11 | > | Differ |
\(CD = 3.69\). F2 differs from both F1 and F3, but F1 and F3 do not differ significantly. F2 is the best-yielding fertilizer.
A randomised block design has 3 treatments in 4 blocks with \(SS_{Tr} = 60,\; SS_{Bl} = 24,\; SS_E = 12,\; df_E = 6\). Test the treatment effect, obtain the CD, and find the relative efficiency of RBD over CRD.
To analyse an RBD, compute the critical difference, and evaluate the gain from blocking via the relative efficiency of RBD over CRD.
Applying it:
Blank ANOVA table:
| Source | SS | df | MS | F |
|---|---|---|---|---|
| Treatments | ||||
| Blocks | — | |||
| Error | — |
| Source | SS | df | MS | F |
|---|---|---|---|---|
| Treatments | 60 | 2 | 30 | 15 |
| Blocks | 24 | 3 | 8 | — |
| Error | 12 | 6 | 2 | — |
\(F_{Tr} = 30/2 = 15\). \(CD = t_{0.025,6}\sqrt{2(2)/4} = 2.447\times 1 = 2.447\).
\(MS_E^{CRD} = (24 + 12)/(3 + 6) = 36/9 = 4\), so \(RE = 4/2 = 2.0\).
\(F_{Tr} = 15 > F_{0.05}(2,6) = 5.14\) ⇒ treatments differ significantly; \(CD = 2.447\). RBD is twice as efficient as CRD, so blocking was worthwhile.
In a 3 treatments × 3 blocks RBD, the observation at T2–B3 is missing. Estimate it.
| B1 | B2 | B3 | |
|---|---|---|---|
| T1 | 10 | 12 | 14 |
| T2 | 14 | 16 | ? |
| T3 | 8 | 10 | 12 |
To estimate a single missing value in an RBD by minimising the error sum of squares.
Applying it:
Blank working table:
| Quantity | Value |
|---|---|
| \(T'\) (treatment total) | |
| \(B'\) (block total) | |
| \(G'\) (grand total) | |
| \(\hat y\) |
| Quantity | Value |
|---|---|
| \(T'_2 = 14+16\) | 30 |
| \(B'_3 = 14+12\) | 26 |
| \(G' \) (all available) | 96 |
| \(t = b\) | 3 |
\(\hat y = \dfrac{3(30) + 3(26) - 96}{(2)(2)} = \dfrac{90 + 78 - 96}{4} = \dfrac{72}{4} = 18\).
The missing value is estimated as \(\hat y = 18\). Insert it and analyse the RBD with the error df reduced by one.
A 5 × 5 Latin square design gives \(SS_R = 30,\; SS_C = 25,\; SS_{Tr} = 200,\; SS_E = 60,\; df_E = 12\). Test the treatment effect, find the CD, and compute the efficiency of LSD relative to CRD and to RBD.
To analyse a Latin square design and evaluate its relative efficiency against CRD and RBD.
Applying it:
Blank working table:
| Comparison | \(MS_E\) of alternative | RE |
|---|---|---|
| LSD vs CRD | ||
| LSD vs RBD (cols as blocks) | ||
| LSD vs RBD (rows as blocks) |
\(MS_E^{LSD} = 60/12 = 5\); \(F_{Tr} = (200/4)/5 = 10 > F_{0.05}(4,12) = 3.26\) ⇒ treatments differ.
\(CD = t_{0.025,12}\sqrt{2(5)/5} = 2.179\sqrt{2} = 3.08\).
| Comparison | \(MS_E\) of alternative | RE |
|---|---|---|
| LSD vs CRD | \((30+25+60)/20 = 5.75\) | 1.15 (15 %) |
| LSD vs RBD (cols as blocks) | \((30+60)/16 = 5.625\) | 1.125 (12.5 %) |
| LSD vs RBD (rows as blocks) | \((25+60)/16 = 5.31\) | 1.06 (6 %) |
Treatments differ significantly (\(F = 10\)); \(CD = 3.08\). The LSD is 15 % more efficient than CRD and 6–12.5 % more efficient than RBD, so controlling both row and column variation was beneficial.
In a 4 × 4 Latin square, one observation is missing. From the available data the row total \(R' = 65\), column total \(C' = 70\), treatment total \(T' = 60\) and grand total \(G' = 240\) (\(t = 4\)). Estimate the missing value.
To estimate a single missing value in a Latin square design.
Applying it:
Blank working table:
| Quantity | Value |
|---|---|
| \(R'\) | |
| \(C'\) | |
| \(T'\) | |
| \(G'\) | |
| \(\hat y\) |
| Quantity | Value |
|---|---|
| \(R'\) | 65 |
| \(C'\) | 70 |
| \(T'\) | 60 |
| \(G'\) | 240 |
\(\hat y = \dfrac{4(65 + 70 + 60) - 2(240)}{(3)(2)} = \dfrac{780 - 480}{6} = \dfrac{300}{6} = 50\).
The missing value is estimated as \(\hat y = 50\). Insert it and analyse the LSD with the error df reduced by one.