Парний bootstrap для даних про неписьменність/народжуваність
Скориставшись щойно написаною функцією, виконайте парний bootstrap і побудуйте гістограму, що описує оцінки коефіцієнта нахилу за даними про неписьменність/народжуваність. Також наведіть 95%-й довірчий інтервал для нахилу. Дані доступні у масивах NumPy illiteracy та fertility.
Нагадаємо, draw_bs_pairs_linreg() має сигнатуру draw_bs_pairs_linreg(x, y, size=1) і повертає два значення: bs_slope_reps та bs_intercept_reps.
Ця вправа є частиною курсу
Статистичне мислення в Python (Частина 2)
Інструкції до вправи
- Використайте вашу функцію
draw_bs_pairs_linreg()для отримання1000bootstrap-реплік коефіцієнта нахилу та вільного члена. Дані по осі x — цеilliteracy, а по осі y —fertility. - Обчисліть і виведіть 95% bootstrap-довірчий інтервал для нахилу.
- Побудуйте та покажіть гістограму реплік нахилу. Обов'язково підпишіть осі. Це вже зроблено за вас, тож натисніть «Надіслати», щоб побачити свою гістограму!
Інтерактивна практична вправа
Спробуйте виконати цю вправу, доповнивши цей зразок коду.
# Generate replicates of slope and intercept using pairs bootstrap
bs_slope_reps, bs_intercept_reps = ____
# Compute and print 95% CI for slope
print(np.percentile(____, ____))
# Plot the histogram
_ = plt.hist(bs_slope_reps, bins=50, normed=True)
_ = plt.xlabel('slope')
_ = plt.ylabel('PDF')
plt.show()