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 dấu và hạng Wilcoxon (Wilcoxon signed-rank test)

    Cho nnn cặp quan sát trước/sau (xi,yi)(x_i,y_i)(xi​,yi​), i=1,…,ni=1,\ldots,ni=1,…,n. Xét hiệu di=yi−xid_i = y_i - x_idi​=yi​−xi​. Loại bỏ các cặp có di=0d_i = 0di​=0; gọi mmm là số cặp còn lại (hiệu lực).

    Xếp hạng tăng dần các giá trị ∣di∣|d_i|∣di​∣ trong mmm cặp hiệu lực (dùng trung bình hạng khi trùng nhau như các bài trước). Tính: W+=∑i: di>0hang(∣di∣),W−=∑i: di<0hang(∣di∣),W=min⁡(W+,W−)W^+ = \sum_{i:\, d_i>0} \text{hang}(|d_i|), \qquad W^- = \sum_{i:\, d_i<0} \text{hang}(|d_i|), \qquad W=\min(W^+,W^-)W+=∑i:di​>0​hang(∣di​∣),W−=∑i:di​<0​hang(∣di​∣),W=min(W+,W−)

    Kỳ vọng và phương sai (hiệu chỉnh trùng): μW=m(m+1)4,σW2=m(m+1)(2m+1)24−∑g(tg3−tg)48\mu_W = \frac{m(m+1)}{4}, \qquad \sigma_W^2 = \frac{m(m+1)(2m+1)}{24} - \frac{\sum_g (t_g^3-t_g)}{48}μW​=4m(m+1)​,σW2​=24m(m+1)(2m+1)​−48∑g​(tg3​−tg​)​ với tổng chạy qua các nhóm trùng ggg của ∣di∣|d_i|∣di​∣ trong mmm giá trị hiệu lực.

    Thống kê chuẩn hoá (hiệu chỉnh liên tục 0.50.50.5), định nghĩa như bài Mann-Whitney (dựa trên dấu của W−μWW-\mu_WW−μW​), và p=2(1−Φ(∣z∣))p = 2(1-\Phi(|z|))p=2(1−Φ(∣z∣)).

    Quy ước biên: nếu m=0m=0m=0 (mọi di=0d_i=0di​=0), in khong xac dinh cho dòng z,pz,pz,p (không có dòng ppp riêng, gộp một dòng) và kết luận Khong bac bo H0.

    Kết luận: Bac bo H0 nếu p<αp<\alphap<α, ngược lại Khong bac bo H0.

    Ví dụ: n=4n=4n=4, x=(1,2,3,4)x=(1,2,3,4)x=(1,2,3,4), y=(2,3,4,5)y=(2,3,4,5)y=(2,3,4,5): mọi di=1>0d_i=1>0di​=1>0 nên W+=1+2+3+4=10W^+ = 1+2+3+4=10W+=1+2+3+4=10 (hạng của 4 giá trị ∣di∣|d_i|∣di​∣ bằng nhau, trung bình hạng =2.5=2.5=2.5 mỗi cái, tổng =10=10=10), W−=0W^-=0W−=0.

    • Định dạng đầu vào:
      • Dòng 1: số nguyên nnn (1≤n≤10001 \le n \le 10001≤n≤1000) và số thực α\alphaα (0<α<10<\alpha<10<α<1).
      • Dòng 2: nnn số thực x1,…,xnx_1,\ldots,x_nx1​,…,xn​ (trước).
      • Dòng 3: nnn số thực y1,…,yny_1,\ldots,y_ny1​,…,yn​ (sau).
    • Định dạng đầu ra:

      Nếu m=0m=0m=0 (không có cặp hiệu lực):

      • Dòng 1: 0 0
      • Dòng 2: khong xac dinh
      • Dòng 3: Khong bac bo H0

      Ngược lại (m>0m>0m>0):

      • Dòng 1: W+W^+W+ và W−W^-W−, mỗi giá trị làm tròn 1 chữ số thập phân, cách nhau bởi khoảng trắng.
      • Dòng 2: giá trị zzz làm tròn 4 chữ số thập phân.
      • Dòng 3: giá trị ppp làm tròn 6 chữ số thập phân.
      • Dòng 4: Bac bo H0 hoặc Khong bac bo H0.

    Ví dụ:

    Đầu vào:

    6 0.05
    10 12 14 16 18 20
    12 15 13 20 19 25
    

    Đầu ra:

    19.5 1.5
    -1.7870
    0.073937
    Khong bac bo H0
    

    Đầu vào:

    1 0.05
    5
    5
    

    Đầu ra:

    0 0
    khong xac dinh
    Khong bac bo H0
    

    Đang tải editor...