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 Kruskal-Wallis H

    Kiểm định Kruskal-Wallis H là kiểm định phi tham số mở rộng của Mann-Whitney cho k≥2k \ge 2k≥2 mẫu độc lập, dùng để kiểm tra xem các mẫu có cùng phân phối hay không.

    Cho kkk nhóm, nhóm thứ jjj có njn_jnj​ quan sát. Gộp toàn bộ N=∑jnjN=\sum_j n_jN=∑j​nj​ quan sát và xếp hạng (rank) tăng dần; giá trị trùng nhận hạng trung bình (average rank). Gọi RjR_jRj​ là tổng hạng của nhóm jjj.

    Thống kê H chưa hiệu chỉnh trùng: H0=12N(N+1)∑j=1kRj2nj−3(N+1)H_0 = \frac{12}{N(N+1)} \sum_{j=1}^{k} \frac{R_j^2}{n_j} - 3(N+1)H0​=N(N+1)12​∑j=1k​nj​Rj2​​−3(N+1)

    Hệ số hiệu chỉnh do giá trị trùng (tie correction): gọi các nhóm giá trị trùng nhau (kích thước ≥2\ge 2≥2) có kích thước lần lượt t1,t2,…t_1, t_2, \dotst1​,t2​,…, thì: C=1−∑i(ti3−ti)N3−NC = 1 - \frac{\sum_i (t_i^3 - t_i)}{N^3 - N}C=1−N3−N∑i​(ti3​−ti​)​

    nếu không có giá trị trùng nào thì C=1C=1C=1. Thống kê đã hiệu chỉnh: H=H0/CH = H_0 / CH=H0​/C; quy ước nếu C=0C=0C=0 thì H=H0H = H_0H=H0​.

    Yêu cầu: Tính HHH (đã hiệu chỉnh trùng) và bậc tự do df=k−1df = k-1df=k−1.

    Ví dụ: 3 nhóm A=(1,2,3)A=(1,2,3)A=(1,2,3), B=(4,5,6)B=(4,5,6)B=(4,5,6), C=(7,8,9)C=(7,8,9)C=(7,8,9) (không có giá trị trùng). N=9N=9N=9. Hạng: A→(1,2,3)A\to(1,2,3)A→(1,2,3), B→(4,5,6)B\to(4,5,6)B→(4,5,6), C→(7,8,9)C\to(7,8,9)C→(7,8,9). RA=6,RB=15,RC=24R_A=6, R_B=15, R_C=24RA​=6,RB​=15,RC​=24. H0=1290(363+2253+5763)−30=1290⋅279−30=37.2−30=7.2H_0 = \frac{12}{90}\left(\frac{36}{3}+\frac{225}{3}+\frac{576}{3}\right) - 30 = \frac{12}{90}\cdot 279 - 30 = 37.2-30=7.2H0​=9012​(336​+3225​+3576​)−30=9012​⋅279−30=37.2−30=7.2. Do không trùng nên C=1C=1C=1, H=7.2H=7.2H=7.2, df=2df=2df=2.

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

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

      kkk dòng tiếp theo, dòng thứ jjj bắt đầu bằng số nguyên njn_jnj​ (nj≥1n_j \ge 1nj​≥1) rồi tiếp theo là njn_jnj​ số thực — các quan sát của nhóm jjj (trị tuyệt đối không vượt quá 10610^6106). Tổng số quan sát N=∑nj≤5000N=\sum n_j \le 5000N=∑nj​≤5000.

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

      In trên một dòng số thực HHH (làm tròn đến 4 chữ số thập phân) và số nguyên dfdfdf, cách nhau bởi khoảng trắng. Ví dụ: 7.2000 2.

    Ví dụ:

    Đầu vào:

    2
    1 5
    1 5
    

    Đầu ra:

    0.0000 1
    

    Đầu vào:

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

    Đầu ra:

    7.2000 2
    

    Đang tải editor...