執行模型
你剛才在產量隨時間變化的圖中看到的平滑趨勢線,是用廣義加性模型(GAM)來決定曲線的位置。這類模型非常適合擬合非線性的曲線。為了能對未來產量進行預測,現在就明確地跑一次模型。這個 GAM 的語法如下:
gam(response ~ s(explanatory_var1) + explanatory_var2, data = dataset)
這裡的 s() 表示「對變數做平滑」,而「平滑」大致上可以理解為非線性。
已經載入 mgcv 和 dplyr;corn 與 wheat 資料集可用,且已擴充為包含公制單位與普查區域資訊。
本練習屬於課程
R 函式撰寫入門
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
# Run a generalized additive model of yield vs. smoothed year and census region
___(___ ~ s(___) + ___, data = ___)