用 GAM 建立黃豆生長模型
在這個練習中,你要將黃豆植株的平均葉片重量建模為時間(播種後)的函數。 你會看到,黃豆並不是以固定速率成長,而是會出現一段「生長衝刺」,之後逐漸趨緩。因此,葉片重量不適合用線性模型來描述。
回想一下,你可以在公式中用 s()(docs)指定要以非線性方式建模的變數:
y ~ s(x)
另外要記得,套件 mgcv 的 gam()(docs)呼叫介面為:
gam(formula, family, data)
一般迴歸可使用 family = gaussian(預設值)。
已預先載入黃豆的訓練資料 soybean_train。它有兩個欄位:目標 weight 與變數 Time。作為比較,使用公式 weight ~ Time 擬合的線性模型 model.lin 也已載入。
本練習屬於課程
R 中的監督式學習:回歸
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
# soybean_train is available
summary(soybean_train)
# Plot weight vs Time (Time on x axis)
ggplot(soybean_train, aes(x = ___, y = ___)) +
geom_point()