这真的很奇怪，我想不出解决办法。大约几个月前，我正在运行这样的 GAM 交互代码。我已经为 GAM 安装了 mgcv 和 gamair 包。它运行良好，没有错误。
Model1 <- gam(response ~ Location + s(predictor1) + s(predctor2, by = Location, m=1) + s(predictor3), data=Data, method = "REML", bs='cr', family=Gamma(link=identity))
然后，突然，相同的代码停止工作。我做了一些测试，发现 by= 术语是问题所在。输出返回此错误消息：
'Error in smoothCon(split$smooth.spec[[i]], data, knots, absorb.cons, scale.penalty = scale.penalty, : Can't find by variable".
如果我删除 by=，代码运行正常。即使只有一个预测变量，如果我添加 by=，模型也不起作用。
我做了 RStudio 更新、软件包更新、病毒检查等等。我不知道为什么具有相同数据的相同代码突然停止工作。有什么建议吗？谢谢。