用模型診斷修正問題
你可以使用 gam.check(),根據其結果更新模型以加以改進。
本練習屬於課程
R 中的廣義加法模型(GAM)之非線性建模
練習說明
- 對
mod執行模型診斷。 - 根據診斷結果,重新配適模型為
mod2,並針對未通過的平滑項調整基底函式數量(k)。 - 對
mod2執行模型診斷,確認你已修正問題。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
library(mgcv)
# Fit the model
mod <- gam(y ~ s(x0, k = 3) + s(x1, k = 3) + s(x2, k = 3) + s(x3, k = 3),
data = dat, method = "REML")
# Check the diagnostics
___
# Refit to fix issues
mod2 <- gam(y ~ s(x0, k = ___) + s(x1, k = ___) + s(x2, k = ___) + s(x3, k = ___),
data = dat, method = "REML")
# Check the new model
___