文盲率/生育率数据的成对自助法
使用您刚刚编写的函数,执行成对自助法,绘制一个直方图来描述基于文盲率/生育率数据得到的斜率估计。同时报告斜率的 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()