ช่วงความน่าเชื่อถือแบบ Highest Posterior Density
การรายงานเพียงค่าประมาณจุดเดียวนั้นไม่เพียงพอ การระบุความไม่แน่นอนของค่าประมาณประสิทธิภาพของยาจะทำให้ผลการวิเคราะห์มีความน่าเชื่อถือมากขึ้น คุณจึงตัดสินใจเพิ่มข้อความต่อไปนี้ลงในบันทึก
ผลการทดลองระบุว่า ด้วยความน่าจะเป็น 90% อัตราประสิทธิภาพของยาใหม่อยู่ระหว่าง ___ ถึง ___ และด้วยความน่าจะเป็น 95% อยู่ระหว่าง ___ ถึง ___
ต้องคำนวณช่วงความน่าเชื่อถือสองช่วง ได้แก่ ช่วงที่มีความน่าจะเป็น 90% และ 95% โดยที่อาร์เรย์ drug_efficacy_posterior_draws ยังคงพร้อมใช้งานใน workspace ของคุณ
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
การวิเคราะห์ข้อมูลแบบเบย์เซียนด้วย Python
คำแนะนำการฝึกหัด
- Import แพ็กเกจ
arvizในชื่อaz - คำนวณช่วงความน่าเชื่อถือแบบ Highest Posterior Density ที่ความน่าจะเป็น 90% แล้วกำหนดให้กับตัวแปร
ci_90 - คำนวณช่วงความน่าเชื่อถือแบบ Highest Posterior Density ที่ความน่าจะเป็น 95% แล้วกำหนดให้กับตัวแปร
ci_95
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Import arviz as az
____
# Calculate HPD credible interval of 90%
ci_90 = ____.____(____, ____=____)
# Calculate HPD credible interval of 95%
ci_95 = ____
# Print the memo
print(f"The experimental results indicate that with a 90% probability \nthe new drug's efficacy rate is between {np.round(ci_90[0], 2)} and {np.round(ci_90[1], 2)}, \nand with a 95% probability it is between {np.round(ci_95[0], 2)} and {np.round(ci_95[1], 2)}.")