b값 정량화
이전 연습 문제에서 생성한 플롯을 바탕으로, 완전성 임계값 mt = 3을 안전하게 사용할 수 있습니다. 이 임계값을 활용하여 1980년부터 2009년까지의 기간과 2010년부터 2017년 중반까지의 기간에 대한 b값을 계산해 보세요. 이전 챕터에서 작성한 b_value() 함수(규모 집합과 완전성 임계값으로 b값과 신뢰 구간을 계산하는 함수)가 네임스페이스에 준비되어 있으며, 지난 연습 문제의 numpy 배열 mags_pre, mags_post와 mt도 사용할 수 있습니다.
이 연습은 강의의 일부입니다
통계적 사고 사례 연구
연습 안내
- 부트스트랩 반복 횟수 10,000회를 사용하여 1980년부터 2009년까지의 지진에 대한 b값과 95% 신뢰 구간을 계산하세요.
- 부트스트랩 반복 횟수 10,000회를 사용하여 2010년부터 2017년 중반까지의 지진에 대한 b값과 95% 신뢰 구간을 계산하세요.
- '답변 제출'을 클릭하여 결과를 화면에 출력하세요.
실습형 인터랙티브 연습
이 예제를 이 샘플 코드를 완성하여 풀어보세요.
# Compute b-value and confidence interval for pre-2010
b_pre, conf_int_pre = ____(____, ____, ____=[____, ____], ____=____)
# Compute b-value and confidence interval for post-2010
b_post, conf_int_post = ____
# Report the results
print("""
1980 through 2009
b-value: {0:.2f}
95% conf int: [{1:.2f}, {2:.2f}]
2010 through mid-2017
b-value: {3:.2f}
95% conf int: [{4:.2f}, {5:.2f}]
""".format(b_pre, *conf_int_pre, b_post, *conf_int_post))