是否有必要在游戏中包含一个因子（拟合为平滑因子交互）作为参数项？

Question

我有兴趣研究数据集中的非线性时间趋势，因此我想使用 R 包 mgcv 来适应以下 GAM：

model1 <- gam(Variable ~ s(Date, by = Site.Factor), data = data)

其中Variable 是感兴趣的连续变量，Site.Factor 是具有两个水平的因子，Date 是连续变量。

我已经读过，知道由于在平滑函数中包含了 by 因子，因此两个因子水平的均值差异没有被考虑在内。因此，我应该将Site.Factor 包含为参数项，如下所示：

model2 <- gam(Variable ~ Site.Factor + s(Date, by = Site.Factor), data = data)

但是，虽然我可能预计 Site.Factor 对平滑的影响很大，但我并不认为该因子的每个级别的均值都很重要。我是否仍需要像 model1 那样在模型中单独包含该因子，还是 model2 可以？

Answer 1

除非您知道从中提取数据的总体具有完全相同相同的平均值，否则您应该将术语 Site.Factor 包含为固定效应术语，样本中的差异是否显着。