Gam：为什么可以'；我是否在gamm中使用s（PTT，bs=“re”？）？

在

mgcv

软件包中，gamm函数可用于模拟随机效果

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

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

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

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

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

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

---

也在gamm内部？它似乎很好用

简单的回答是，

gamm

通过调用

nlme

包中的

lme

来模拟随机效应，因此语法与此相同。如果您使用

gam

语法调用它，则会发生一些不同的情况。

因此，gamm甚至不会像gam那样处理平滑术语？好吧，

summary（gam_model）

的输出看起来与

summary（gamm_model$gam）

的输出类似，但不完全相同。您应该使用一个或另一个适当的语法。