生成置换重复值
正如视频中所讨论的,置换重复值是从一次置换样本计算得到的统计量的一个取值。类似于您在第 2 章编写的用于生成自助法重复值的 draw_bs_reps() 函数(https://campus.datacamp.com/courses/statistical-thinking-in-python-part-2/bootstrap-confidence-intervals?ex=6),我们同样希望有一个类似的函数 draw_perm_reps() 来生成置换重复值。本练习将编写这个实用函数。
该函数的调用签名为 draw_perm_reps(data_1, data_2, func, size=1)。需要注意,func 必须是一个以两个数组为参数的函数。在大多数情况下,func 将是您自己编写的函数。
本练习是课程的一部分
Python 统计思维(第 2 部分)
练习说明
- 定义具有以下签名的函数:
draw_perm_reps(data_1, data_2, func, size=1)。- 使用
np.empty()初始化用于存放置换重复值的数组。 - 编写
for循环以:- 使用您的
permutation_sample()函数计算一次置换样本; - 将所得样本传入
func()来计算单次重复值,并将结果存入重复值数组中。
- 使用您的
- 返回该重复值数组。
- 使用
交互式实操练习
通过完成这段示例代码来试试这个练习。
def draw_perm_reps(data_1, data_2, func, size=1):
"""Generate multiple permutation replicates."""
# Initialize array of replicates: perm_replicates
perm_replicates = ____
for i in ____:
# Generate permutation sample
perm_sample_1, perm_sample_2 = ____
# Compute the test statistic
perm_replicates[i] = ____
return perm_replicates