オートバイ衝突データ：線形アプローチ
この最初の演習では、データセットに線形モデルを当てはめ、結果を可視化して、データ内の関係をどれくらい捉えられているかを確認します。データフレーム mcycle に格納されたこのデータセットは、オートバイ衝突時にクラッシュテスト用ダミーの頭部加速度を測定したものです。accel 列には加速度 (g)、times 列には時間 (ミリ秒) が入っています。
この演習はコースの一部です
Rで学ぶGeneralized Additive Models (GAMs) による非線形モデリング
演習の手順
head()とplot()を使って、データフレームmcycleの中身を確認します。lm()を使って、mcycleデータに対し、accel変数をtimesの「線形」関数としてモデル化するモデルを当てはめます。- 用意された
termplot()の呼び出しで、モデルの当てはまりを可視化します。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Examine the mcycle data frame
head(___)
plot(___)
# Fit a linear model
lm_mod <- lm(___, data = mcycle)
# Visualize the model
termplot(lm_mod, partial.resid = TRUE, se = TRUE)