부트스트랩 신뢰 구간

Python으로 하는 통계적 사고 (2부)

Justin Bois

Lecturer at the California Institute of Technology

부트스트랩 복제 함수

def bootstrap_replicate_1d(data, func):
    """Generate bootstrap replicate of 1D data."""
    bs_sample = np.random.choice(data, len(data))
    return func(bs_sample)

bootstrap_replicate_1d(michelson_speed_of_light, np.mean)

299859.20000000001

bootstrap_replicate_1d(michelson_speed_of_light, np.mean)

299855.70000000001

bootstrap_replicate_1d(michelson_speed_of_light, np.mean)

299850.29999999999
Python으로 하는 통계적 사고 (2부)

다수의 부트스트랩 복제본

bs_replicates = np.empty(10000)

for i in range(10000):
    bs_replicates[i] = bootstrap_replicate_1d(
                  michelson_speed_of_light, np.mean)
Python으로 하는 통계적 사고 (2부)

부트스트랩 복제본의 히스토그램 그리기

_ = plt.hist(bs_replicates, bins=30, normed=True)
_ = plt.xlabel('mean speed of light (km/s)')
_ = plt.ylabel('PDF')
plt.show()
Python으로 하는 통계적 사고 (2부)

평균의 부트스트랩 추정

ch2-2.011.png

Python으로 하는 통계적 사고 (2부)

통계량의 신뢰구간

  • 측정을 무한히 반복한다면, 관측값의 p%가 p% 신뢰구간에 포함됩니다.
Python으로 하는 통계적 사고 (2부)

부트스트랩 신뢰 구간

conf_int = np.percentile(bs_replicates, [2.5, 97.5])

array([ 299837.,  299868.])

ch2-2.016.png

Python으로 하는 통계적 사고 (2부)

연습해 봅시다!

Python으로 하는 통계적 사고 (2부)

Preparing Video For Download...