我试图在具有一个固定效应的GAMM模型中指定随机截距和随机斜率项.
我已经使用mgcv库中的以下代码成功地使用随机截距拟合了一个模型,但现在无法确定gamm()函数中随机斜率的语法:
M1 = gamm(dur ~ s(dep, bs="ts", k = 4), random= list(fInd = ~1), data= df)
如果我在线性混合效果模型中同时使用随机截距和斜率,我会按以下方式编写它:
M2 = lme(dur ~ dep, random=~1 + dep|fInd, data=df)
在gamm()支持文档指出随机条款需要在给予list的形式为lme(),但我找不到任何解释的例子,其中包括斜率和截距项.任何建议/解决方案将不胜感激.
小智 6
包中的gamm4函数gamm4包含执行此操作的方法。您可以按照与lmer样式中相同的方式指定随机截距和斜率。在你的情况下：
M1 = gamm4(dur~s(dep,bs="ts",k=4), random = ~(1+dep|fInd), data=df)
这是 gamm4 文档：https ://cran.r-project.org/web/packages/gamm4/gamm4.pdf