【发布时间】：2020-11-11 08:08:36
【问题描述】：
我正在尝试使用 mgcv 包拟合模型，并且我需要包含分类术语之间的交互。我知道如何实现连续变量和分类变量之间的交互
s(x1, by=c1)
但我不知道如何更改这两个分类术语，我找不到任何东西。
【问题讨论】：
-
如果它只是两个分类变量并且不包括任何平滑项然后以通常的方式添加交互
c1*c2
【发布时间】：2020-11-11 08:08:36
【问题描述】：
您可以像在线性模型或 GLM 中一样将参数项添加到公式中。
gam(y ~ f1 + f2 + f1:f2 + s(x1), data = df, method = "REML")
或通过简写
gam(y ~ f1 * f2 + s(x1), data = df, method = "REML")
将是一个具有主效应加上两个分类变量f1 和f2 的交互作用加上连续变量x1 的平滑效应的模型。
您可以对lm() 或glm() 中的参数项执行的任何操作都可以用于gam()、gamm()、bam() 和gamm4::gamm4() 的公式中。
【讨论】：