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ê] Cập nhật Bayes tuần tự qua nhiều lần xét nghiệm

    Một bệnh nhân có xác suất tiên nghiệm mắc bệnh là p0p_0p0​. Bệnh nhân được thực hiện kkk lần xét nghiệm liên tiếp bằng cùng một loại xét nghiệm có độ nhạy sesese và độ đặc hiệu spspsp.

    Sau mỗi lần xét nghiệm, xác suất hậu nghiệm được tính bằng định lý Bayes, và giá trị này được dùng làm xác suất tiên nghiệm cho lần xét nghiệm kế tiếp (cập nhật Bayes tuần tự), giả sử các lần xét nghiệm độc lập có điều kiện khi biết tình trạng bệnh thực sự:

    pj=pj−1⋅L(rj)pj−1⋅L(rj)+(1−pj−1)⋅Lˉ(rj)p_j = \dfrac{p_{j-1}\cdot L(r_j)}{p_{j-1}\cdot L(r_j) + (1-p_{j-1})\cdot \bar L(r_j)}pj​=pj−1​⋅L(rj​)+(1−pj−1​)⋅Lˉ(rj​)pj−1​⋅L(rj​)​

    trong đó L(duong)=seL(\text{duong})=seL(duong)=se, L(am)=1−seL(\text{am})=1-seL(am)=1−se, Lˉ(duong)=1−sp\bar L(\text{duong})=1-spLˉ(duong)=1−sp, Lˉ(am)=sp\bar L(\text{am})=spLˉ(am)=sp.

    Nếu tại bước nào đó mẫu số bằng 0 (tổ hợp dữ liệu không thể xảy ra), xác suất hậu nghiệm từ bước đó trở đi được coi là không xác định.

    Ví dụ: p0=0.01p_0=0.01p0​=0.01, se=0.9se=0.9se=0.9, sp=0.95sp=0.95sp=0.95, k=2k=2k=2, kết quả duong duong: bước 1 cho p1≈0.153846p_1\approx0.153846p1​≈0.153846, bước 2 dùng p1p_1p1​ làm tiên nghiệm cho p2≈0.765957p_2\approx0.765957p2​≈0.765957.

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

      Dòng 1: số thực p0p_0p0​.

      Dòng 2: hai số thực se spse\ spse sp.

      Dòng 3: số nguyên kkk (1≤k≤10001 \le k \le 10001≤k≤1000).

      Dòng 4: kkk từ, mỗi từ là duong hoặc am - kết quả của từng lần xét nghiệm theo thứ tự thực hiện.

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

      In ra kkk dòng; dòng thứ jjj là xác suất hậu nghiệm pjp_jpj​ sau khi cập nhật với kết quả xét nghiệm thứ jjj, làm tròn 6 chữ số thập phân.

      Nếu tại bước jjj mẫu số bằng 0, in -1.000000 cho dòng jjj và tất cả các dòng còn lại.

    Ví dụ:

    Đầu vào:

    0.5
    0.8 0.8
    1
    am

    Đầu ra:

    0.200000
    

    Đầu vào:

    0.01
    0.9 0.95
    2
    duong duong

    Đầu ra:

    0.153846
    0.765957
    

    Đang tải editor...