Sidebar Menu

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
 (B)

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 H01 = μ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:

  1.  $ 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}$
  1. 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

  1. 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

  1. 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:

  1.  $ 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}$
  1. 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

  1. 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

  1. 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

Graphical user interface, application Description automatically generated

Table Description automatically generated

Table Description automatically generated with medium confidence

  

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
Freedom

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.