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()的公式中。