b-값 계산하기
b-값은 특정 지역의 지진 활동을 나타내는 대표적인 지표입니다. 지진 데이터를 다룰 때 이 값을 반복적으로 계산해야 하는 경우가 많습니다. 이처럼 자주 수행하는 작업은 함수로 만들어 두는 것이 좋습니다! 시그니처가 b_value(mags, mt, perc=[2.5, 97.5], n_reps=None)인 함수를 작성하세요. 이 함수는 주어진 진도 집합 mags에 대해 b-값과 (선택적으로, n_reps가 None이 아닌 경우) 신뢰 구간을 반환해야 합니다. 완전성 임계값은 mt로 주어지며, 키워드 인자 perc는 신뢰 구간의 하한과 상한에 해당하는 백분위수를 나타내고, n_reps는 신뢰 구간 계산에 사용할 부트스트랩 복제 횟수입니다.
이 연습은 강의의 일부입니다
통계적 사고 사례 연구
연습 안내
- 다음 동작을 수행하는
b_value(mags, mt, perc=[2.5, 97.5], n_reps=None)시그니처의 함수를 정의하세요.- 불리언 인덱싱을 사용하여
mags에서 완전성 임계값mt이상의 진도를 추출하고, 결과를 변수m에 저장하세요. - b-값의 최적 추정치를 계산하세요. b-값의 최적 추정치는 b = (m - mt)·ln(10)임을 기억하세요. 결과를 변수
b에 저장하세요. n_reps가 None이 아닌 경우, 다음을 수행하세요.m의 평균에 대한 부트스트랩 복제값을n_reps개 추출하고, 결과를 변수m_bs_reps에 저장하세요.m의 평균에 대한 부트스트랩 복제값을 b-값의 복제값으로 변환하고, 결과를b_bs_reps에 저장하세요.- b-값의 부트스트랩 복제값으로부터 신뢰 구간을 계산하고, 결과를
conf_int에 저장하세요.
b와conf_int를 반환하되,n_reps가None인 경우에는b만 반환하세요.
- 불리언 인덱싱을 사용하여
실습형 인터랙티브 연습
이 예제를 이 샘플 코드를 완성하여 풀어보세요.
____ ____:
"""Compute the b-value and optionally its confidence interval."""
# Extract magnitudes above completeness threshold: m
m = ____[____ >= ____]
# Compute b-value: b
b = ____
# Draw bootstrap replicates
if n_reps is None:
return b
else:
m_bs_reps = ____
# Compute b-value from replicates: b_bs_reps
b_bs_reps = (____ - ____) * ____
# Compute confidence interval: conf_int
conf_int = ____
return b, conf_int