最高後驗密度可信區間
你知道只回報單一點估計不夠。最好也提供此藥物效力率估計的不確定性度量，而你已具備所需的一切工具。你決定在備忘錄中加入以下內容。
實驗結果顯示，以 90% 的機率，新藥的效力率介於 ___ 與 ___ 之間；以 95% 的機率，介於 ___ 與 ___ 之間。
你需要計算兩個可信區間：一個是 90%，另一個是 95% 的機率。drug_efficacy_posterior_draws 陣列仍可在你的工作環境中使用。
本練習屬於課程
Bayesian Data Analysis in Python
練習說明
- 將
arviz套件以az的名稱匯入。 - 計算 90% 的最高後驗密度可信區間，並指定給
ci_90。 - 計算 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)}.")