Kredibilitetsintervall med högsta posteriora densitet
Du vet att det inte räcker att rapportera enskilda punktskattningar. Det vore bra att också ange ett mått på osäkerheten i skattningen av läkemedlets effektivitet – och du har alla verktyg som behövs. Du bestämmer dig för att lägga till följande i promemorian.
De experimentella resultaten visar att det med 90% sannolikhet ligger mellan ___ och ___, och med 95% sannolikhet mellan ___ och ___.
Du behöver beräkna två kredibilitetsintervall: ett med 90% sannolikhet och ett med 95% sannolikhet. Arrayen drug_efficacy_posterior_draws finns fortfarande tillgänglig i din arbetsyta.
Den här övningen är en del av kursen
Bayesiansk dataanalys i Python
Övningsinstruktioner
- Importera paketet
arvizsomaz. - Beräkna kredibilitetsintervallet med högsta posteriora densitet på 90% och tilldela det till
ci_90. - Beräkna kredibilitetsintervallet med högsta posteriora densitet på 95% och tilldela det till
ci_95.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# 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)}.")