最高后验密度可信区间
您知道只报告点估计并不够。最好还能给出该药物疗效率估计的不确定性度量，而您已经具备所需的一切。您决定在备忘录中补充如下内容：
实验结果表明，以 90% 的概率，新药的疗效率位于 ___ 与 ___ 之间；以 95% 的概率，位于 ___ 与 ___ 之间。
您需要计算两个可信区间：一个是 90% 概率，另一个是 95% 概率。数组 drug_efficacy_posterior_draws 仍可在您的工作区中使用。
本练习是课程的一部分
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)}.")