GAM（mgcv）中分类术语之间的相互作用

Question

我正在尝试使用 mgcv package 来拟合 model 并且我需要包括分类术语之间的交互。 我知道如何实现连续变量和分类变量之间的交互

s(x1, by=c1)

但我不知道如何为两个分类术语更改它，我找不到任何东西。

Answer 1

您可以将参数项添加到公式中，就像在线性 model 或 GLM 中一样。

gam(y ~ f1 + f2 + f1:f2 + s(x1), data = df, method = "REML")

或通过速记

gam(y ~ f1 * f2 + s(x1), data = df, method = "REML")

将是一个 model ，主效应加上两个分类变量f1和f2的交互作用加上连续变量x1的平滑效应。

您可以对lm()或glm()中的参数项执行的任何操作都可以用于gam() 、 gamm() 、 bam()和gamm4::gamm4()的公式中。