Věrohodnostní intervaly nejvyšší posteriorní hustoty
Víš, že samotné bodové odhady nestačí. Bylo by skvělé doplnit míru nejistoty odhadu účinnosti léku – a máš k tomu vše potřebné. Rozhodneš se do zprávy přidat následující:
Experimentální výsledky ukazují, že s 90% pravděpodobností leží míra účinnosti nového léku mezi ___ a ___, a s 95% pravděpodobností mezi ___ a ___.
Budeš muset vypočítat dva věrohodnostní intervaly: jeden s 90% a druhý s 95% pravděpodobností. Pole drug_efficacy_posterior_draws je stále dostupné v tvém pracovním prostoru.
Toto cvičení je součástí kurzu
Bayesian Data Analysis in Python
Pokyny k cvičení
- Importuj balíček
arvizjakoaz. - Vypočítej věrohodnostní interval nejvyšší posteriorní hustoty s 90% pravděpodobností a přiřaď ho do
ci_90. - Vypočítej věrohodnostní interval nejvyšší posteriorní hustoty s 95% pravděpodobností a přiřaď ho do
ci_95.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
# 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)}.")