GAM:收缩平滑样条曲线是否也适用于弯曲?

机器算法验证 模型选择 aic 广义加法模型 样条 弯曲度
2022-03-20 15:08:29

我有一个基于三次样条 (bs = cr) 和 SELECT == T 或收缩三次样条 (bs = cs) 和 SELECT == F 的自动预测器选择的 gam 模型。
现在我想知道受弯曲度影响的预测器是否自动通过收缩方法从我的模型中删除,或者我是否应该在检查弯曲度后自己将它们从模型中删除。

我的模型如下所示:

genericModel <- mgcv::bam(formula = TT2_ScaledTransformed ~
            s(NDVI, k = k, bs = cr) + ... +
            s(BuildingHeight_10m_10std, k = k, bs = cr),
            data = trainSet,
            family = gaussian, method = "fREML", select = TRUE, 
            control = ctrl,
            cluster = cl, gamma = 1.4)

例如,在此示例中,BuildingHeight 和 NDVI(植被指数)部分弯曲,因为这两个参数部分相反。

我检查这样的concurvity:

print(concurvity(genericModel, full = TRUE))

我基于 AIC 评估我的模型。

1个回答

简而言之,不,使用select = TRUE不会自动删除弯曲项。您仍然应该检查结果模型中项的弯曲度,并决定是否删除高度弯曲的项,检查当您删除一个弯曲项时其他估计项如何变化。

也就是说,适合method = "REML"(或"ML""fREML"取决于上下文)并且select = TRUE如果您不想或不能删除弯曲条款,则可能是针对我们拥有的模型中的弯曲条款提出的问题的最佳保护。