다음 파크필드 대지진은 언제 발생할까요?
파크필드 지역에서 마지막으로 발생한 대지진은 현지 시각으로 2004년 9월 27일 저녁이었습니다. 이번 과제에서는 지수 모델과 가우시안 모델을 각각 적용하여 다음 파크필드 지진이 언제 발생할지 추정합니다. 두 모델 모두에서 최적 추정값은 평균 시간 간격으로 나타나며, 이전 연습 문제에서 계산한 값은 24.62년입니다. 즉, 다음 지진은 2029년경에 발생할 것으로 예상됩니다. 이전 연습 문제에서 계산한 mean_time_gap을 모수로 하는 지수 분포를 가정하여 다음 지진 발생 시점에 대한 95% 신뢰 구간을 계산하세요. 동일하게 mean_time_gap과 std_time_gap을 모수로 하는 정규 분포를 가정한 경우도 계산하세요.
이 연습은 강의의 일부입니다
통계적 사고 사례 연구
연습 안내
mean_time_gap을 평균으로 하는 지수 분포에서 샘플 100,000개를 추출하세요. 결과를exp_samples에 저장하세요.mean_time_gap을 평균으로 하고std_time_gap을 표준 편차로 하는 정규 분포에서 샘플 100,000개를 추출하세요. 결과를norm_samples에 저장하세요.- 오늘 현재까지 파크필드 지진이 발생하지 않았으므로,
today - last_quake보다 큰 샘플만 슬라이싱하세요. 여기서today는 오늘의 연도 소수 표현이고,last_quake = 2004.74는 마지막 파크필드 지진의 연도 소수 표현입니다. 슬라이싱한 배열로exp_samples와norm_samples변수를 각각 덮어쓰세요. np.percentile()을 사용하여 다음 파크필드 지진 발생 시점에 대한 95% 신뢰 구간을 계산하세요. 같은 함수 호출에서 50번째 백분위수를 포함하면 중앙값도 함께 계산할 수 있습니다.
실습형 인터랙티브 연습
이 예제를 이 샘플 코드를 완성하여 풀어보세요.
# Draw samples from the Exponential distribution: exp_samples
exp_samples = ____
# Draw samples from the Normal distribution: norm_samples
norm_samples = ____
# No earthquake as of today, so only keep samples that are long enough
exp_samples = ____[____ > ____ - ____]
norm_samples = ____[____ > ____ - ____]
# Compute the confidence intervals with medians
conf_int_exp = ____(____, [____, ____, ____]) + last_quake
conf_int_norm = ____(____, [____, ____, ____]) + last_quake
# Print the results
print('Exponential:', conf_int_exp)
print(' Normal:', conf_int_norm)