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 phù hợp phân phối Nhị thức với tham số ước lượng

    Một thí nghiệm gồm NNN lần lặp độc lập, mỗi lần gồm nnn phép thử con (ví dụ: tung nnn đồng xu một lượt), và ta đếm số lần "thành công" k∈{0,1,…,n}k \in \{0,1,\dots,n\}k∈{0,1,…,n} trong mỗi lần lặp. Cho tần số quan sát O0,O1,…,OnO_0, O_1,\dots,O_nO0​,O1​,…,On​ (với ∑k=0nOk=N\sum_{k=0}^n O_k = N∑k=0n​Ok​=N).

    Giả thuyết H0H_0H0​: số lần thành công tuân theo phân phối Nhị thức Binomial(n,p)\text{Binomial}(n, p)Binomial(n,p) với ppp chưa biết, được ước lượng bằng phương pháp hợp lý cực đại:

    p^=∑k=0nk Okn N\hat p = \frac{\sum_{k=0}^{n} k\,O_k}{n\,N}p^​=nN∑k=0n​kOk​​

    Tần số kỳ vọng của nhóm kkk: Ek=N(nk)p^k(1−p^)n−kE_k = N \binom{n}{k} \hat p^{k} (1-\hat p)^{n-k}Ek​=N(kn​)p^​k(1−p^​)n−k.

    Chỉ tính tổng trên các nhóm có Ek>0E_k > 0Ek​>0 (loại các nhóm suy biến khi p^∈{0,1}\hat p \in \{0,1\}p^​∈{0,1} khiến Ek=0E_k=0Ek​=0); gọi mmm là số nhóm có Ek>0E_k>0Ek​>0:

    χ2=∑k:Ek>0(Ok−Ek)2Ek,df=max⁡(m−2, 0)\chi^2 = \sum_{k: E_k>0} \frac{(O_k-E_k)^2}{E_k}, \qquad df = \max(m-2,\ 0)χ2=∑k:Ek​>0​Ek​(Ok​−Ek​)2​,df=max(m−2, 0)

    (trừ 111 do ràng buộc tổng tần số, trừ thêm 111 do ước lượng tham số ppp). So sánh χ2\chi^2χ2 với giá trị tới hạn χcrit2\chi^2_{crit}χcrit2​ cho trước: nếu χ2>χcrit2\chi^2 > \chi^2_{crit}χ2>χcrit2​ thì bác bỏ H0H_0H0​, ngược lại chưa đủ cơ sở bác bỏ.

    Ví dụ: n=4n=4n=4, O=5,20,40,25,10O = 5, 20, 40, 25, 10O=5,20,40,25,10 (N=100N=100N=100): p^=0⋅5+1⋅20+2⋅40+3⋅25+4⋅104⋅100=0.5375\hat p = \dfrac{0\cdot5+1\cdot20+2\cdot40+3\cdot25+4\cdot10}{4\cdot100} = 0.5375p^​=4⋅1000⋅5+1⋅20+2⋅40+3⋅25+4⋅10​=0.5375. Tính EkE_kEk​ theo Nhị thức rồi χ2≈1.1566\chi^2 \approx 1.1566χ2≈1.1566, df=3df=3df=3.

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

      Dòng 1: số nguyên nnn (1≤n≤301 \le n \le 301≤n≤30) và số thực χcrit2≥0\chi^2_{crit} \ge 0χcrit2​≥0, cách nhau bởi dấu cách. Dòng 2: n+1n+1n+1 số nguyên không âm O0,O1,…,OnO_0, O_1,\dots,O_nO0​,O1​,…,On​ (tổng N=∑Ok>0N=\sum O_k > 0N=∑Ok​>0), cách nhau bởi dấu cách.

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

      In ra 3 dòng: Dòng 1: giá trị χ2\chi^2χ2 làm tròn đến đúng 4 chữ số thập phân (%.4f). Dòng 2: số nguyên dfdfdf. Dòng 3: chuỗi Bac bo H0 nếu χ2>χcrit2\chi^2 > \chi^2_{crit}χ2>χcrit2​, ngược lại Khong bac bo H0.

    Ví dụ:

    Đầu vào:

    4 3.841
    100 0 0 0 0
    

    Đầu ra:

    0.0000
    0
    Khong bac bo H0
    

    Đầu vào:

    4 9.488
    5 20 40 25 10
    

    Đầu ra:

    1.1566
    3
    Khong bac bo H0
    

    Đang tải editor...