蜜蜂精子计数的自助法假设检验
现在，您将检验如下假设：平均而言，接受新烟碱类杀虫剂处理的雄蜂与未处理雄蜂在每毫升精液中的活性精子数量相同。您将使用"均值之差"作为检验统计量。
供您参考，您在第 2 章编写的 draw_bs_reps() 函数的调用签名为 draw_bs_reps(data, func, size=1)，可在此处查看：链接。
本练习是课程的一部分
Python 统计思维（第 2 部分）
练习说明
- 计算
control的活性精子计数均值减去treated的活性精子计数均值。 - 计算所有活性精子计数的总体均值。为此，先连接
control和treated，再对连接后的数组取均值。 - 为
control和treated生成移位后的数据集，使两者具有相同的均值。此步骤已为您完成。 - 使用您的
draw_bs_reps()函数，分别为两个移位后的数组各生成 10,000 个均值的自助法重复值。 - 计算"均值之差"的自助法重复值。
- 计算并打印 p 值的代码已为您写好。点击提交即可查看结果！
交互式实操练习
通过完成这段示例代码来试试这个练习。
# Compute the difference in mean sperm count: diff_means
diff_means = ____
# Compute mean of pooled data: mean_count
mean_count = ____
# Generate shifted data sets
control_shifted = control - np.mean(control) + mean_count
treated_shifted = treated - np.mean(treated) + mean_count
# Generate bootstrap replicates
bs_reps_control = ____(____,
np.mean, size=10000)
bs_reps_treated = ____(____,
np.mean, size=10000)
# Get replicates of difference of means: bs_replicates
bs_replicates = ____
# Compute and print p-value: p
p = np.sum(bs_replicates >= np.mean(control) - np.mean(treated)) \
/ len(bs_replicates)
print('p-value =', p)