在 R 中使用 GAM 为最佳拟合模型选择变量？

Question

我也是 R 的新手，我正在尝试使用 mgcv 将鱼类丰度与 11 个环境因素之间的关系建模为 11 个解释变量。

y ~ s(x1) + s(x2) + ... + s(x11)

我尝试了太多模型（来自 1 个变量；2 个变量.... 11 个变量模型），然后使用 AIC 选择最佳拟合模型。但是，我无法同时运行所有模型。

你能帮我解决这个问题吗？如果您可以共享代码，那就太好了。谢谢，

NVN

Answer 1

这不是进行推理的好方法，但如果您只对预测感兴趣，可能还可以。

相反，考虑将select = TRUE 添加到模型中（它是gam() 等的一个参数。这个参数为基的零空间（完美平滑的部分）增加了额外的惩罚，它允许您进行模型选择因为每个平滑中的线性项也向零收缩（有点像套索），因此也可以从模型中选择出来。