オートバイ衝突データ：非線形アプローチ
それでは、同じ mcycle データに対して mgcv パッケージの gam() 関数を使い、非線形モデルを当てはめます。
この演習はコースの一部です
Rで学ぶGeneralized Additive Models (GAMs) による非線形モデリング
演習の手順
mgcvパッケージを読み込みます。gam()関数を使い、mcycleデータでaccelがtimesに対して滑らかな非線形関係を持つモデルを当てはめます。- 提示された
plot()の呼び出しでモデル当てはめを可視化します。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Load mgcv
___
# Fit the model
gam_mod <- gam(___ ~ s(___), data = mcycle)
# Plot the results
plot(gam_mod, residuals = TRUE, pch = 1)