부트스트랩 신뢰 구간
통계적 사고 사례 연구
Justin Bois
Lecturer, Caltech
EDA는 첫 번째 단계
"탐색적 데이터 분석이 모든 것을 담을 수는 없지만, 그 어떤 것도 기초 단계인 첫 번째 단계를 대체할 수는 없다."
--John Tukey
활동 구간 길이 ECDF
1 데이터 제공: Avni Gandhi, Grigorios Oikonomou, David Prober, Caltech
최적 매개변수 값
- 최적 매개변수 값: 데이터를 가장 잘 설명하는 확률 분포의 매개변수 값
- 지수 분포의 최적 매개변수: 데이터 평균으로부터 계산
np.mean(nuclear_incident_times)
87.140350877192986
1 데이터 출처: Wheatley, Sovacool, and Sornette, Nuclear Events Database
부트스트랩 표본
데이터를 재샘플링한 배열
# Resample nuclear_incident_times with replacement
bs_sample = np.random.choice(
nuclear_incident_times,
replace=True,
size=len(inter_times)
)
부트스트랩 복제본
1 데이터 출처: Wheatley, Sovacool, and Sornette, Nuclear Events Database
부트스트랩 복제본
1 데이터 출처: Wheatley, Sovacool, and Sornette, Nuclear Events Database
부트스트랩 복제본
1 데이터 출처: Wheatley, Sovacool, and Sornette, Nuclear Events Database
부트스트랩 복제본
1 데이터 출처: Wheatley, Sovacool, and Sornette, Nuclear Events Database
부트스트랩 복제본
부트스트랩 복제본: 부트스트랩 표본에서 계산된 통계량
dcst.draw_bs_reps()
데이터셋에서 부트스트랩 복제본을 추출하는 함수
# Draw 10000 replicates of the mean from
# nuclear_incident_times
bs_reps = dcst.draw_bs_reps(
nuclear_incident_times, np.mean, size=10000
)
부트스트랩 신뢰 구간
1 데이터 출처: Wheatley, Sovacool, and Sornette, Nuclear Events Database
부트스트랩 신뢰 구간
측정을 반복했을 때, 관측값의 p%가 p% 신뢰 구간 내에 포함됩니다
부트스트랩 신뢰 구간
np.percentile(bs_reps, [2.5, 97.5])
array([ 73.31505848, 102.39181287])
