Парный бутстрэп для данных о грамотности и рождаемости
Используйте написанную вами функцию, чтобы выполнить парный бутстрэп и построить гистограмму оценки наклона по данным о неграмотности и рождаемости. Также вычислите 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(), чтобы получить1000бутстрэп-реплик наклона и точки пересечения. Данные по оси X —illiteracy, по оси Y —fertility. - Вычислите и выведите на экран 95%-й бутстрэп-доверительный интервал для наклона.
- Постройте и отобразите гистограмму реплик наклона. Не забудьте подписать оси. Этот код уже готов — нажмите Отправить ответ, чтобы увидеть гистограмму!
Интерактивное практическое упражнение
Попробуйте выполнить это упражнение, дополнив этот пример кода.
# 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()