시작하기무료로 시작하기

최고 사후밀도(HPD) 신뢰 구간

점 추정만 보고하는 것은 충분하지 않다는 것을 알고 계시죠. 신약의 효능 비율 추정치에 대한 불확실성도 함께 제시하면 좋습니다. 이를 위해 메모에 다음 내용을 추가하기로 합니다.

실험 결과에 따르면, 90%의 확률로 신약의 효능 비율은 ___와 ___ 사이에 있으며, 95%의 확률로는 ___와 ___ 사이에 있습니다.

두 개의 신뢰 구간(credible interval)을 계산해야 합니다. 하나는 90%, 다른 하나는 95% 확률입니다. drug_efficacy_posterior_draws 배열은 워크스페이스에 그대로 있습니다.

이 연습은 강의의 일부입니다

Python으로 배우는 Bayesian 데이터 분석

강의 보기

연습 안내

  • arviz 패키지를 az로 임포트하세요.
  • 90% 최고 사후밀도(HPD) 신뢰 구간을 계산해 ci_90에 할당하세요.
  • 95% 최고 사후밀도(HPD) 신뢰 구간을 계산해 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)}.")
코드 편집 및 실행