1. 学ぶ
  2. /
  3. コース
  4. /
  5. Rで学ぶGeneralized Additive Models (GAMs) による非線形モデリング

Connected

演習

マルチスケール相互作用をもつ汚染モデル

meuse データセットには、同じスケールの予測変数(x、y)と、異なるスケールの予測変数(elev、dist、om)が含まれています。前の演習では、位置と標高を説明変数としてカドミウム汚染を予測するモデルを当てはめました。

mod <- gam(cadmium ~ s(x, y) + s(elev), 
           data = meuse, method = "REML")

この演習では、テンソルスムース te() を使って、スケールが異なっていても複数の変数が相互作用できるモデルを作成します。

指示

100 XP
  • x、y、elev がそれぞれのスケールで変動しつつ、1つの te() 項の中で相互作用するモデルに書き換えてください。
  • その後、モデルを要約し、plot() で可視化しましょう。