开始使用免费开始使用

假设检验:2012 年喙更深吗?

通过 ECDF 图和置信区间,您大概已经看到达芙妮主岛上 G. scandens 的喙更深了。但这会不会只是随机性的结果?换句话说,如果两年的均值其实相同,我们观察到当前均值差的概率是多少?

请注意!我们检验的假设并不是"喙深度来自同一分布"。若检验这一点,可以使用置换检验。这里的原假设是"均值相等"。要进行该假设检验,需先将两个数据集平移到具有相同的均值,然后用自助法抽样来计算均值之差。

本练习是课程的一部分

Python 统计思维(第 2 部分)

查看课程

练习说明

  • 将 1975 年和 2012 年的喙深度拼接为一个数组,并计算并存储其均值。
  • 平移 bd_1975bd_2012,使它们的均值等于您刚为合并数据计算的均值。
  • 分别对 1975 年和 2012 年的喙深度各生成 10,000 次均值的自助法重复。
  • 用 2012 年的重复减去 1975 年的重复,得到差值的自助法重复。
  • 计算并打印 p 值。上一练习中得到的观测均值差仍在命名空间中,变量名为 mean_diff

交互式实操练习

通过完成这段示例代码来试试这个练习。

# Compute mean of combined data set: combined_mean
combined_mean = ____(____((bd_1975, bd_2012)))

# Shift the samples
bd_1975_shifted = ____
bd_2012_shifted = ____

# Get bootstrap replicates of shifted data sets
bs_replicates_1975 = ____
bs_replicates_2012 = ____

# Compute replicates of difference of means: bs_diff_replicates
bs_diff_replicates = ____

# Compute the p-value
p = np.sum(____ >= ____) / len(____)

# Print p-value
print('p =', p)
编辑并运行代码