НачатьНачать бесплатно

Достоверные интервалы наибольшей апостериорной плотности

Вы знаете, что приводить лишь точечные оценки недостаточно. Было бы хорошо указать меру неопределённости в оценке эффективности препарата — и у вас есть все необходимые инструменты для этого. Вы решаете добавить в служебную записку следующее.

Экспериментальные результаты показывают, что с вероятностью 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)}.")
Редактировать и запускать код