Skip to content

Phân tích ANOVA cho thí nghiệm phân bón và năng suất lúa

Chúng ta sẽ từng bước thực hiện kiểm định ANOVA để xác định xem có sự khác biệt về năng suất lúa giữa ba loại phân bón A, B, C hay không.

1. Thiết lập giả thuyết kiểm định

Chúng ta sử dụng kiểm định phương sai một chiều (ANOVA một yếu tố) để so sánh giá trị trung bình giữa ba nhóm.

  • Giả thuyết gốc (H₀): Không có sự khác biệt về năng suất lúa giữa ba loại phân bón A, B và C.
  • Giả thuyết thay thế (H₁): Ít nhất một loại phân bón có ảnh hưởng đáng kể đến năng suất lúa.

2. Dữ liệu thí nghiệm

Chúng ta có năng suất lúa (kg/ha) được đo trên 15 ô ruộng, chia đều cho ba loại phân bón:

Ô thí nghiệmPhân bón APhân bón BPhân bón C
1500052004800
2510053004900
3495052504700
4505054004850
5490051504750
  • Số nhóm (k) = 3 (A, B, C)
  • Tổng số quan sát (n) = 15

3. Tính các tham số quan trọng

Bước 1: Tính giá trị trung bình của mỗi nhóm

    \[ \bar{X}_A = \frac{5000 + 5100 + 4950 + 5050 + 4900}{5} = 5000\]


    \[ \bar{X}_B = \frac{5200 + 5300 + 5250 + 5400 + 5150}{5} = 5260\]


    \[ \bar{X}_C = \frac{4800 + 4900 + 4700 + 4850 + 4750}{5} = 4800\]

Bước 2: Tính trung bình tổng thể (Overall Mean)

    \[ \bar{X}_{ overall}} = \frac{5000 + 5100 + 4950 + 5050 + 4900 + 5200 + 5300 + 5250 + 5400 + 5150 + 4800 + 4900 + 4700 + 4850 + 4750}{15}\]


    \[ = 5020\]

4. Tính tổng bình phương (Sum of Squares – SS)

(1) Tổng bình phương giữa các nhóm (SSB – Between Groups)

    \[ SSB = n \sum (\bar{X_i} - \bar\bar{{X}})^2\]


    \[ = 5 \left[ (5000 - 5020)^2 + (5260 - 5020)^2 + (4800 - 5020)^2 \right]\]


    \[ = 5 \left[ (-20)^2 + (240)^2 + (-220)^2 \right]\]


    \[ = 5 \times (400 + 57600 + 48400) = 5 \times 106400 = 532000\]

(2) Tổng bình phương trong nhóm (SSW – Within Groups)

    \[ SSW = \sum (X_{ij} - \bar{X}_i)^2\]

Tính riêng cho từng nhóm:

  • Nhóm A:

        \[ (5000 - 5000)^2 + (5100 - 5000)^2 + (4950 - 5000)^2 + (5050 - 5000)^2 + (4900 - 5000)^2\]


        \[ = 0 + 10000 + 2500 + 2500 + 10000 = 25000\]

  • Nhóm B:

        \[ (5200 - 5260)^2 + (5300 - 5260)^2 + (5250 - 5260)^2 + (5400 - 5260)^2 + (5150 - 5260)^2\]


        \[ = 3600 + 1600 + 100 + 19600 + 12100 = 57000\]

  • Nhóm C:

        \[ (4800 - 4800)^2 + (4900 - 4800)^2 + (4700 - 4800)^2 + (4850 - 4800)^2 + (4750 - 4800)^2\]


        \[ = 0 + 10000 + 10000 + 2500 + 2500 = 25000\]

Tổng SSW:

    \[ SSW = 25000 + 57000 + 25000 = 107000\]

5. Tính giá trị F

    \[ MSB = \frac{SSB}{k-1} = \frac{532000}{3-1} = \frac{532000}{2} = 266000\]

    \[ MSW = \frac{SSW}{n-k} = \frac{107000}{15-3} = \frac{107000}{12} = 8916.67\]

    \[ F = \frac{MSB}{MSW} = \frac{266000}{8916.67} = 36.69\]

6. So sánh với giá trị tới hạn

  • Với độ tự do giữa nhóm (df1) = 2độ tự do trong nhóm (df2) = 12, giá trị F tới hạn ở mức α = 0.05 khoảng 3.89 (tra bảng F-distribution).
  • F = 36.69 lớn hơn rất nhiều 3.89, ta bác bỏ giả thuyết H₀.

p-value rất nhỏ (≈ 0.00000771), nhỏ hơn 0.05, nghĩa là sự khác biệt giữa các nhóm có ý nghĩa thống kê.

7. Kết luận

  • Có sự khác biệt đáng kể về năng suất giữa các loại phân bón A, B và C.

Leave a Reply

error: Content is protected !!