通过模型诊断修复问题
您可以使用 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
___