r语言 - 如何将平滑参数分别应用于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。在第二个例子中，你有一个固定的光滑参数0.5用于X1的边缘光滑，一个光滑参数1用于张量积相互作用光滑中X2的边缘光滑。