Example 1: RCBD Factorial (Interaction is not significant)
This article is a continuation of the RCBD Factorial.
The following is Experimental Data on the Effect of Soil Processing and Organic Fertilizers on the Aggregate Stability Index. Soil Processing consists of 3 levels and 4 levels of Organic Fertilizer. The experiment was arranged using the basic design of Completely Randomized Block Design (RCBD) . The following are the steps for the calculation of the Analysis of Variety followed by Post Hoc Test: Fisher's LSD .
A full discussion of the Example of Factorial RCBD (Interaction is not significat) can be read in the following document. Tutorial on data processing using SmartstatXL Add-In can be learned on the following link:
Example 1: Factorial RCBD (Not significant interaction)
Table 26. Experiment on the Effect of Tillage and Organic Fertilizers on the Aggregate Stability Index.
Tillage (A) | Organic Fertilizer | Block (K) | Grand Total ∑AB | ||
1 | 2 | 3 | |||
1 | 0 | 154 | 151 | 165 | 470 |
10 | 166 | 166 | 160 | 492 | |
20 | 177 | 178 | 176 | 531 | |
30 | 193 | 189 | 200 | 582 | |
2 | 0 | 143 | 147 | 139 | 429 |
10 | 149 | 156 | 171 | 476 | |
20 | 160 | 164 | 136 | 460 | |
30 | 190 | 166 | 169 | 525 | |
3 | 0 | 139 | 134 | 145 | 418 |
10 | 162 | 147 | 166 | 475 | |
20 | 181 | 161 | 149 | 491 | |
30 | 161 | 172 | 182 | 515 | |
Grand Total | ∑K | 1975 | 1931 | 1958 | 5864 |
Calculation Analysis of Variance:
Step 1: Calculate the Correction Factor
$$ CF=\frac{Y...^2}{abr}=\frac{(5864)^2}{3\times4\times3}=955180.44$$
Step 2: Calculate the Sum of The Total Squares
$$\begin{matrix}SSTOT=\sum_{i,j,k}{Y_{ijk}}^2-CF\\=(154)^2+(151)^2+...+(182)^2-955180.44\\=9821.56\\\end{matrix}$$
Step 3: Calculate the Sum of Squares of Groups
$$\begin{matrix}SSR=\frac{\sum_{k}{(r_k)^2}}{ab}-CF\\=\frac{(1975)^2+(1931)^2+(1958)^2}{3\times4}-955180.44\\=82.06\\\end{matrix}$$
The Sum of Squares of Treatment needs to be decomposed into the Sum of the Squares of its components.
Create a Table For Total Treatment
Organic Fertilizer (B) | |||||
Tillage (A) | 0 | 10 | 20 | 30 | ΣA = Yi.. |
1 | 470 | 492 | 531 | 582 | 2075 |
2 | 429 | 476 | 460 | 525 | 1890 |
3 | 418 | 475 | 491 | 515 | 1899 |
ΣB=Y. j. | 1317 | 1443 | 1482 | 1622 | 5864 |
Step 4: Calculate the Sum of Squares of Factor A
$$\begin{matrix}SS(A)=\frac{\sum_{i}{(a_i)^2}}{rb}-CF\\=\frac{(2075)^2+(1890)^2+(1899)^2}{3\times4}-955180.4444\\=1813.39\\\end{matrix}$$
Step 5: Calculate the Sum of Squares of Factor B
$$\begin{matrix}SS(B)=\frac{\sum_{j}{(b_j)^2}}{ra}-CF\\=\frac{(1317)^2+(1443)^2+(1482)^2+(1622)^2}{3\times3}-955180.44\\=5258.00\\\end{matrix}$$
Step 6: Calculate the Sum of Squares of AB Interactions
$$\begin{matrix}SS(AB)=\frac{\sum_{i,j}{(a_ib_j)^2}}{r}-CF-SS(A)-SS(B)\\=\frac{(470)^2+(492)^2+...+(491)^2+(515)^2}{3}-955180.44-1813.39-5258.00\\=463.50\\\end{matrix}$$
Note: SST = SS(A) + SS(B) + SS(AB)
Step 7: Calculate the Sum of Squares of Errors
$$\begin{matrix}SSE=\ SSTOT\ -\ SSB\ -\ SS(A)\ -\ SS(B)\ -SS(AB)\\=9821.56-82.06-1813.39-5258.00-463.50\\=2204.61\\\end{matrix}$$
Step 8: Create a Variance Analysis Table along with its F-tables.
Table 27. ANOVA Table of Two Factor Factorial Design in Completely Randomized Block Design
Source Variance | Degree of Freedom (df) | Sum Squares | Squares Mean | F-stat | F0. 05 | F0. 01 |
Block (R) | r-1 = 2 | 82.06 | 41.0277778 | 0.41 ns | 3.443 | 5.719 |
Treatment |
|
|
|
|
|
|
A | a-1 = 2 | 1813.39 | 906.6944444 | 9.05 ** | 3.443 | 5.719 |
B | b-1 = 3 | 5258.00 | 1752.666667 | 17.49 ** | 3.049 | 4.817 |
AB | (a-1) (b-1) = 6 | 463.50 | 77.25 | 0.77 ns | 2.549 | 3.758 |
Error | ab(r-1) = 22 | 2204.61 | 100.209596 | - |
|
|
Total | abr-1 = 35 | 9821.56 |
|
|
|
|
F(0.05,2,22) =3.443
F(0.0 1,2,22) = 5,719
F(0.05,3,22) = 3,049
F(0.0 1,3,22) = 4,817
F(0.05,6,22) = 2,549
F(0.0 1,6,22) = 3,758
Step 9: Make a Conclusion
Interaction Effect: not significant
Since F-stat (0.77) ≤ 2.549 then we fail to reject H0:μ1 = μ2 = .... at a confidence level of 95%. This means that at a confidence level of 95%, there is no difference in the effect of the interaction on the observed response.
Effect of Factor A: significant
Since F-stat (9.05) > 3,443 then we reject H0: μ1 = μ2 = .... at a confidence level of 95%. This means that at a confidence level of 95%, there is one or more of the average different treatments from others. Or in other words, a decision to reject Ho can be made, meaning that there is a difference in the effect of Factor A on the observed response.
Effect of Factor B: significant
Since F-stat (9.05) > 3,443 then we reject H0: μ1 = μ2 = .... at a confidence level of 95%. This means that at a confidence level of 95%, there is one or more of the average different treatments from others. Or in other words, a decision can be made to reject Ho, meaning that there is a difference in the effect of Factor B on the observed response.
Since Interaction is not significant (significant), then we proceed to the examination of its main effect. Both of its main effect are significant, so we need to further investigate which treatments are the same and which are different. Perform further testing to compare treatment averages, both treatment mean differences for Factor A and Factor B.
Post-Hoc
Based on the analysis of variance, the effect of the interaction between Factor A and Factor B is not significant, while both main effect are significant so that further testing is only carried out on the main effect of the two factors we are trying.
In this further test, the differences between the average pairs of treatments were carried out using the LSD test.
Main Effects of Tillage Factors (A)
Calculate LSD with the following Formula:
- $ LSD=t_{\alpha/2;db}\sqrt{\frac{2MSE}{rb}}$
- MSE = 100.21
- error-degree of freedom = 22
- Block (r) = 3; Level of Factor B (b) = 4
- t(α/2.22) = t(α/2.22) = 2,074 (See t-student table at a significant level, α = 0.05, and df = 22, or if using functions in MS Excel, write the formula "=tinv(0.05,22)"
- The above parameters are entered into the formula:
- $\begin{matrix}LSD=t_{\alpha/2;db}\sqrt{\frac{2MSE}{rb}}\\=t_{0.05/2;22}\sqrt{\frac{2(100.21)}{3\times4}}\\=2.074\times4.087\\=8.475\\\end{matrix}$
- Create a Table of treatment averages for Factor A (the main effect of Factor A), and then sort from small to large values (ascending order).
Tillage (O) | Average |
1 | 172.92 |
2 | 157.50 |
3 | 158.25 |
- Sorted by tillage:
Tillage (O) | Average |
2 | 157.50 |
3 | 158.25 |
1 | 172.92 |
- Compare the difference in the average treatment with the LSD value = 8,475. To make the work easier, create a matrix table of average differences as in the following example:
If Soil (O) | 2 | 3 | 1 |
| |
Average | 157.50 | 158.25 | 172.92 |
| |
2 | 157.50 | 0.00 | a | ||
3 | 158.25 | 0.75 | 0.00 | a | |
1 | 172.92 | 15.42 | 14.67 | 0.00 | b |
- After being given a letter notation, return the order based on the order of treatment (not the average order). The end result is as follows:
Tillage (O) | Average |
1 | 172.92 b |
2 | 157.50 a |
3 | 158.25 a |
Main Effect of Organic Fertilizer Factors (B)
Calculate LSD with the following Formula:
- $ LSD=t_{\alpha/2;db}\sqrt{\frac{2MSE}{ra}}$
- MSE = 100.21
- error-degree of freedom = 22
- Block (r) = 3; Factor Level A (a) = 3
- t(α/2.22) = t(α/2.22) = 2,074 (See t-student table at a significant level, α = 0.05, and df = 22, or if using functions in MS Excel, write the formula "=tinv(0.05,22)"
- The above parameters are entered into the formula:
- $\begin{matrix}LSD=t_{\alpha/2;db}\sqrt{\frac{2MSE}{ra}}\\=t_{0.05/2;22}\sqrt{\frac{2(100.21)}{3\times3}}\\=2.074\times4.719\\=9.787\\\end{matrix}$
- Create a Table of treatment averages for Factor B (the main effect of Factor B), then sort from small to large values (ascending order). It just so happens that in this example, the average value of the treatment has been sorted from small to large.
Organic Fertilizer (P) | Average |
0 | 146.3333 |
10 | 160.3333 |
20 | 164.6667 |
30 | 180.2222 |
- Compare the difference in the average treatment with the LSD value = 9.787. To make the work easier, create a matrix table of average differences as in the following example:
Fertilizer Organic (P) | 0 | 10 | 20 | 30 |
| |
Average | 146.33 | 160.33 | 164.67 | 180.22 |
| |
0 | 146.33 | 0.00 | a | |||
10 | 160.33 | 14.00 | 0.00 | b | ||
20 | 164.67 | 18.33 | 4.33 | 0.00 | b | |
30 | 180.22 | 33.89 | 19.89 | 15.56 | 0.00 | c |
- After being given a letter notation, return the order based on the order of treatment (not the average order). It just so happens that the order is as simple as according to the order of treatment. The final result is as follows:
Organic Fertilizer (P) | Average |
0 | 146.3333 a |
10 | 160.3333 b |
20 | 164.6667 b |
30 | 180.2222 c |
Calculation Using SmartstatXL Excel Add-In
Example 2: Factorial RCBD
The data is provided as follows:
A | B | Block | Yij. | |||
1 | 2 | 3 | 4 | |||
a0 | b0 | 12 | 15 | 14 | 13 | 54 |
a0 | b1 | 19 | 22 | 23 | 21 | 85 |
a1 | b0 | 29 | 27 | 33 | 30 | 119 |
a1 | b1 | 32 | 35 | 38 | 37 | 142 |
Y.. k | 92 | 99 | 108 | 101 | Y...= 400 |
Analysis of Variance Calculation:
$$ CF=\frac{Y...^2}{abr}=\frac{(400)^2}{2\times2\times4}=10000$$
$$\begin{matrix}SSTOT=\sum_{i,j,k}{Y_{ijk}}^2-CF\\=(12)^2+(15)^2+...+(37)^2-10000\\=1170\\\end{matrix}$$
$$\begin{matrix}SSR=\frac{\sum_{k}{(r_k)^2}}{ab}-CF\\=\frac{(92)^2+(99)^2+(108)^2+(101)^2}{2\times2}-10000\\=32.5\\\end{matrix}$$
Create a Table For Total Treatment
a0 | a1 | ΣB = Y.j. | |
b0 | 54 | 119 | 173 |
b1 | 85 | 142 | 227 |
ΣA=Yi.. | 139 | 261 | 400 |
$$\begin{matrix}SS(A)=\frac{\sum_{i}{(a_i)^2}}{rb}-CF\\=\frac{(139)^2+(261)^2}{4\times2}-10000\\=930.25\\\end{matrix}$$
$$\begin{matrix}SS(B)=\frac{\sum_{j}{(b_j)^2}}{ra}-CF\\=\frac{(173)^2+(227)^2}{4\times2}-10000\\=182.25\\\end{matrix}$$
$$\begin{matrix}SS(AB)=\frac{\sum_{i,j}{(a_ib_j)^2}}{r}-CF-SS(A)-SS(B)\\=\frac{(54)^2+(85)^2+(119)^2+(142)^2}{4}-10000-930.25-182.25\\=4\\\end{matrix}$$
Note: SST = SS(A) + SS(B) + SS(AB)
$$\begin{matrix}SSE=\ SSTOT\ -\ SSB\ -\ SS(A)\ -\ SS(B)\ -SS(AB)\\=1170-32.5-930.25-182.25-4\\=21\\\end{matrix}$$
Table 28. Anova Table of Two-Factorial Designs In RCBD
Source Variance | Degree of | Sum Squares | Squares Mean | F-stat | F0. 05 | F0. 01 |
Block (R) | r-1 = 3 | 32.5 | 10.833 | 4.64* | 3.86 | 6.99 |
Treatment |
|
|
|
|
|
|
A | a-1 = 1 | 930.25 | 930.25 | 398.679** | 5.11 | 10.56 |
B | b-1 = 1 | 182.25 | 182.25 | 78.107** | 5.11 | 10.56 |
AB | (a-1) (b-1) = 1 | 4 | 4 | 1.714 | 5.11 | 10.56 |
Error | ab(r-1) = 9 | 21 | 2.33 |
|
|
|
Total | abr-1 = 15 | 1170 |
|
|
|
|
Post-Hoc
Based on the analysis of variance, the effect of the interaction between Factor A and Factor B is not significant, while both main effects are significant so that further testing is only carried out on the main effect of the two factors we are trying.
In this follow-up test, the differences between the average pairs of treatments were carried out using the Duncan test.