始める無料で始める

置換レプリケートの生成

動画で説明したように、置換レプリケートとは、置換サンプルから計算される統計量の単一の値のことです。第2章で作成した draw_bs_reps() 関数はブートストラップのレプリケートを生成するのに便利でしたが、同様に置換レプリケートを生成するための draw_perm_reps() 関数もあると便利です。この演習では、その関数を実装します。

この関数の呼び出しシグネチャは draw_perm_reps(data_1, data_2, func, size=1) です。重要なポイントとして、func は引数として配列を「2つ」受け取る関数である必要があります。多くの場合、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
コードを編集して実行