モデルを実行する

これまで見てきた収量の経年プロットにある滑らかなトレンドラインは、一般化加法モデル（GAM）を使って線の位置を推定しています。この種のモデルは非線形な曲線の当てはめに最適です。将来の収量を予測できるように、ここでは明示的にモデルを実行してみましょう。GAM を実行するための構文は次のとおりです。

gam(response ~ s(explanatory_var1) + explanatory_var2, data = dataset)

ここで s() は「変数をスムーズ化する」ことを表し、スムーズ化はごく大まかに言うと非線形化を意味します。

mgcv と dplyr は読み込まれており、corn と wheat の各データセットはメートル法の単位と国勢調査の地域情報を含むように拡張されています。

corn データセットを使って、yield_kg_per_ha を目的変数、年をスムーズ化し国勢調査の地域を加えたモデルで GAM を実行します。