產生置換重抽樣值
如同影片所述,置換重抽樣值(permutation replicate)是由一個置換樣本計算出的統計量單一數值。就像你在第 2 章所寫的 draw_bs_reps() 函式(可幫助你產生自助法重抽樣值)一樣,擁有一個類似的函式 draw_perm_reps() 來產生置換重抽樣值也很實用。你會在本練習撰寫這個實用的函式。
此函式的呼叫介面為 draw_perm_reps(data_1, data_2, func, size=1)。重點是,func 必須是一個接受「兩個」陣列作為引數的函式。在多數情況下,func 會是你自己撰寫的函式。
本練習屬於課程
Statistical Thinking in 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