如何使用 GAM 应用两个变量之间的交互以进行多重非线性回归？

Question

我有一个 Y 和 X1 的数据，X2 有不同的维度，比如 X1 = xxxx 伏特，X2 = xx 小时。

为了使用这些数据建立回归模型，我使用了下面的代码。

MODEL <- gam(Y ~ s(X1) + s(X2), data = DATA, method = "REML")

它似乎运行良好，但我想在我的代码中应用 X1 和 X2 之间的交互。

我可以使用如下代码进行多重非线性回归：

MODEL <- gam(Y ~ s(X1) + s(X2) + ti(X1, X2), data = DATA, method = "REML")

或者我应该使用不同的方程式来完成这项工作吗？

Answer 1

首选方法是：

MODEL <- gam(Y ~ te(X1, X2), data = DATA, method = "REML")

因为这需要选择更少的平滑度参数，所以它是您展示的模型的一个稍微简单的版本。

但是，如果您想查看交互是否显着，或者想从交互中分离出主要影响，那么是的，您展示的模型就是要走的路

MODEL <- gam(Y ~ s(X1) + s(X2) + ti(X1, X2), data = DATA, method = "REML")