Gam：为什么我不能在 gamm 中使用 s(PTT, bs = "re")？

Question

在mgcv包中，gamm函数可以用来模拟随机效应。

在通常的 gam 函数中，我们会使用

gam(y ~ s(Group, bs = "re")

为此，但在 gamm 中，我们使用

gamm(y ~ 1, random = list(Group = ~1))

.....我的问题是：为什么我们不能只使用

gamm(y ~ s(Group, bs = "re"))

也在 gamm 里面？它似乎工作正常？

Answer 1

简短的回答是 gamm 通过调用 nlme 包中的 lme 来模拟随机效应，因此语法与此相同。如果您使用 gam 语法调用它，则会发生一些不同的事情。