Intervalles de crédibilité de densité a posteriori maximale
Vous savez que donner uniquement des estimations ponctuelles ne suffit pas. Il serait utile d’indiquer l’incertitude autour du taux d’efficacité du médicament, et vous avez tout ce qu’il faut pour le faire. Vous décidez d’ajouter ce qui suit à la note de service.
Les résultats expérimentaux indiquent qu’avec une probabilité de 90 %, le taux d’efficacité du nouveau médicament se situe entre ___ et ___, et qu’avec une probabilité de 95 %, il se situe entre ___ et ___.
Vous devez calculer deux intervalles de crédibilité : l’un avec 90 % et l’autre avec 95 % de probabilité. Le tableau drug_efficacy_posterior_draws est toujours disponible dans votre espace de travail.
Cet exercice fait partie du cours
Analyse de données bayésienne en Python
Instructions
- Importez le package
arvizsous le nomaz. - Calculez l’intervalle de crédibilité de densité a posteriori maximale à 90 % et assignez-le à
ci_90. - Calculez l’intervalle de crédibilité de densité a posteriori maximale à 95 % et assignez-le à
ci_95.
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de code.
# 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)}.")