Các lượt nộp
    Danh sách bài
    Trang chủ
    Báo lỗi

    solution

    Đề bài: [Xác suất - Thống kê] Kiểm định hậu định Fisher's LSD sau ANOVA

    Sau khi ANOVA một yếu tố cho kết luận bác bỏ H0H_0H0​ (có ít nhất một cặp nhóm khác biệt), người ta thường dùng kiểm định hậu định (post-hoc) Fisher's LSD (Least Significant Difference) để xác định cụ thể cặp nhóm nào khác biệt. Giả sử kkk nhóm có cùng cỡ mẫu nnn.

    Trước tiên tính tổng bình phương trong nhóm và trung bình bình phương trong nhóm như ANOVA một yếu tố:

    SSW=∑i=1k∑j=1n(xij−xˉi)2,dfW=k(n−1),MSW=SSWdfWSSW=\sum_{i=1}^k\sum_{j=1}^n (x_{ij}-\bar{x}_i)^2, \qquad df_W = k(n-1), \qquad MSW = \frac{SSW}{df_W}SSW=∑i=1k​∑j=1n​(xij​−xˉi​)2,dfW​=k(n−1),MSW=dfW​SSW​

    Khoảng cách khác biệt tối thiểu có ý nghĩa (LSD), với tαt_{\alpha}tα​ là giá trị ttt tới hạn cho trước:

    LSD=tα2 MSWnLSD = t_{\alpha}\sqrt{\frac{2\,MSW}{n}}LSD=tα​n2MSW​​

    Với mỗi cặp nhóm (i,j)(i, j)(i,j), i<ji < ji<j: nếu ∣xˉi−xˉj∣>LSD|\bar{x}_i - \bar{x}_j| > LSD∣xˉi​−xˉj​∣>LSD thì kết luận cặp đó khác biệt có ý nghĩa; ngược lại (kể cả bằng nhau) thì không khác biệt có ý nghĩa.

    Ví dụ: với 3 nhóm {1,2,3},{4,5,6},{7,8,9}\{1,2,3\},\{4,5,6\},\{7,8,9\}{1,2,3},{4,5,6},{7,8,9} và tα=2.776t_{\alpha}=2.776tα​=2.776, ta có MSW=1MSW=1MSW=1, LSD=2.7762/3≈2.2666LSD=2.776\sqrt{2/3}\approx 2.2666LSD=2.7762/3​≈2.2666. Tất cả 3 cặp đều có chênh lệch trung bình lớn hơn LSDLSDLSD nên đều khác biệt có ý nghĩa.

    • Định dạng đầu vào:

      Dòng đầu tiên chứa số nguyên kkk (2≤k≤202 \le k \le 202≤k≤20) — số nhóm.

      kkk dòng tiếp theo, mỗi dòng có định dạng n xi1 … xinn \ x_{i1} \ \dots \ x_{in}n xi1​ … xin​, trong đó nnn (2≤n≤10002 \le n \le 10002≤n≤1000) là cỡ mẫu — giống nhau ở tất cả các dòng — và xijx_{ij}xij​ là số thực.

      Dòng cuối cùng chứa một số thực tα>0t_{\alpha} > 0tα​>0.

    • Định dạng đầu ra:

      In ra 1+(k2)1 + \binom{k}{2}1+(2k​) dòng:

      • Dòng 1: giá trị LSDLSDLSD, làm tròn 4 chữ số thập phân.
      • (k2)\binom{k}{2}(2k​) dòng tiếp theo, mỗi dòng ứng với một cặp (i,j)(i, j)(i,j) với i<ji<ji<j theo thứ tự tăng dần của iii rồi jjj (tức (1,2),(1,3),…,(1,k),(2,3),…(1,2),(1,3),\dots,(1,k),(2,3),\dots(1,2),(1,3),…,(1,k),(2,3),…), có định dạng:

      i j diff verdict

      trong đó i, j là chỉ số nhóm (đánh số từ 1), diff là ∣xˉi−xˉj∣|\bar{x}_i-\bar{x}_j|∣xˉi​−xˉj​∣ làm tròn 4 chữ số thập phân, và verdict là chuỗi chính xác Khac biet hoặc Khong khac biet.

      Ví dụ với input mẫu ở đề bài, output là:

      2.2666
      1 2 3.0000 Khac biet
      1 3 6.0000 Khac biet
      2 3 3.0000 Khac biet
      

    Ví dụ:

    Đầu vào:

    3
    3 1 2 3
    3 4 5 6
    3 7 8 9
    2.776
    

    Đầu ra:

    2.2666
    1 2 3.0000 Khac biet
    1 3 6.0000 Khac biet
    2 3 3.0000 Khac biet
    

    Đầu vào:

    2
    4 10 12 11 13
    4 20 22 21 23
    2.447
    

    Đầu ra:

    2.2338
    1 2 10.0000 Khac biet
    

    Đang tải editor...