始める無料で始める

複数のインピュテーションモデルをまとめて比較する

異なるインピュテーション手法を評価するには、単一のデータフレームにまとめる必要があります。ここでは、データセット oceanbuoys を使って欠損値への3つの対応方法を比較します。

  • 1つ目の方法は、完成ケースのみを使用するもので、ocean_cc として読み込まれています。
  • 2つ目の方法は、風の情報を用いた線形モデルで値をインピュートするもので、ocean_imp_lm_wind として読み込まれています。

3つ目のインピュート済みデータセット ocean_imp_lm_all を作成します。線形モデルを用い、wind_ew、wind_ns、year、latitude、longitude の変数を使って、sea_temp_c、air_temp_c、humidity をインピュートします。

その後、すべてのデータセット（ocean_cc、ocean_imp_lm_wind、ocean_imp_lm_all）を結合し、bound_models という名前にします。

この演習はコースの一部です

Rでの欠損データの扱い方

コースを見る

演習の手順

  • 線形モデルを使って ocean_imp_lm_all という名前のインピュート済みデータセットを作成し、wind_ew、wind_ns、year、latitude、longitude の変数を用いて sea_temp_c、air_temp_c、humidity をインピュートします。
  • すべてのデータセットを同じオブジェクトに結合し、bound_models と名付けます。

実践的なインタラクティブ演習

このサンプルコードを完成させて、この演習に挑戦してみましょう。

# Create an imputed dataset using a linear models
ocean_imp_lm_all <- bind_shadow(oceanbuoys) %>%
  add_label_shadow() %>%
  impute_lm(sea_temp_c ~ wind_ew + wind_ns + ___ + ___ + ___) %>%
  impute_lm(air_temp_c ~ wind_ew + wind_ns + ___ + ___ + ___) %>%
  impute_lm(humidity ~ wind_ew + wind_ns + ___ + ___ + ___)

# Bind the datasets
bound_models <- bind_rows(cc = ___,
                          imp_lm_wind = ___,
                          imp_lm_all = ___,
                          .id = "imp_model")
# Look at the models
bound_models
コードを編集して実行