ブートストラップの実行
前の演習で calc_gender_coef() をうまく作成できましたね。この関数は、ブートストラップサンプルを生成して補完を行い、映画の主人公が女性であることが興行収入に与える影響を表す線形回帰係数を出力します。
この演習では、boot パッケージを使って、そうした係数のブートストラップ分布を求めます。この分布の広がりが、補完に伴う不確実性を反映します。また、ブートストラップ分布が、1回だけ補完して回帰を行った場合とどう異なるかも確認しましょう。さっそく始めましょう!
この演習はコースの一部です
R による欠損データの補完処理
演習の手順
bootパッケージを読み込みます。biopicsに対してブートストラップを実行します。statisticにcalc_gender_coefを渡し、Rを 50 に設定して、結果をboot_resultsに代入します。- ブートストラップの結果を
print()とplot()で確認します。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Load the boot library
___
# Run bootstrapping on biopics data
boot_results <- ___(___, statistic = ___, R = ___)
# Print and plot bootstrapping results
___
___