Intervale credibile Highest Posterior Density
Știi că raportarea unor estimări punctuale simple nu este suficientă. Ar fi util să oferi și o măsură a incertitudinii în estimarea ratei de eficacitate a medicamentului – și ai toate instrumentele necesare pentru asta. Decizi să adaugi următoarele în memo.
Rezultatele experimentale indică faptul că, cu o probabilitate de 90%, rata de eficacitate a noului medicament se situează între ___ și ___, iar cu o probabilitate de 95%, aceasta se situează între ___ și ___.
Va trebui să calculezi două intervale credibile: unul cu probabilitate de 90% și altul cu probabilitate de 95%. Vectorul drug_efficacy_posterior_draws este în continuare disponibil în spațiul tău de lucru.
Acest exercițiu face parte din cursul
Analiza bayesiană a datelor în Python
Instrucțiuni pentru exercițiu
- Importă pachetul
arvizcaaz. - Calculează intervalul credibil Highest Posterior Density de 90% și atribuie-l variabilei
ci_90. - Calculează intervalul credibil Highest Posterior Density de 95% și atribuie-l variabilei
ci_95.
Exercițiu interactiv practic
Încearcă acest exercițiu completând acest cod de exemplu.
# 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)}.")