如何在mgcv中将平滑参数分别应用于主效应和交互作用？

Question

我使用下面的代码将平滑参数仅应用于交互项。

mgcv::gam(Y ~ s(X1, k=3, sp=-1) + s(X2, k=4, sp=-1) + 
            ti(X1, X2, k=4, sp=c(1,1)) + X3,
          method = "REML")

k 值对于我的模型来说很小，因为我只想获得非常平滑的变量样条曲线。

在此代码中，我尝试将 sp = 1 应用于交互项。

这种使用 sp 选项控制平滑度参数的方法是否正确？

而且，我可以在 ti() 交互项中单独设置 X1、X2 的 sp 参数值吗？（如下图）

ti(X1, X2, k=3, sp=c(0.5, 1))

Answer 1

是的。在第一个示例中，X1 和 X2 的边际平滑度将其平滑度参数固定为值 1。在第二个示例中，张量积中 X1 的边缘平滑有固定平滑参数 0.5，张量积中 X2 的边缘平滑有固定平滑参数 1交互流畅

How to apply smooth parameter to main effect and interaction separately in mgcv?