가설 검정: b값에 차이가 있을까요?
이전 연습 문제에서 설계한 가설 검정을 수행하세요. 변수 mags_pre, mags_post, 그리고 mt = 3은 이미 네임스페이스에 로드되어 있습니다.
이 연습은 강의의 일부입니다
통계적 사고 사례 연구
연습 안내
- 2010년 이전 지진 중 완전성 임계값 이상의 규모를 가진 지진만 추출하여
mags_pre를 덮어쓰세요.mags_post도 동일하게 처리하세요. - 2010년 이후 지진의 평균 규모에서 2010년 이전 지진의 평균 규모를 빼서 관측된 평균 규모 차이를 계산하세요.
dcst.draw_perm_reps()를 사용하여 순열 복제값 10,000개를 생성하세요.func인수로는dcst.diff_of_means를 사용하세요.- "최소한 동일하게 극단적인" 기준을 검정 통계량이 관측값보다 작은 경우로 정의하여 p-값을 계산하고 출력하세요.
실습형 인터랙티브 연습
이 예제를 이 샘플 코드를 완성하여 풀어보세요.
# Only magnitudes above completeness threshold
mags_pre = ____[____ >= ____]
mags_post = ____[____ >= ____]
# Observed difference in mean magnitudes: diff_obs
diff_obs = ____ - ____
# Generate permutation replicates: perm_reps
perm_reps = ____
# Compute and print p-value
p_val = ____(____ < ____) / 10000
print('p =', p_val)