手指上的细菌。具有随机斜率但不具有MASS::glmmPQL截获的交叉随机效果的语法

Question

我有不正常的数据(戴和不戴手套接触表面后手指上的细菌)，所以使用MASS包中的glmmPQL。我有一个分类预测变量(Gloves)，一个重复测量变量(NumberContacts)，实验参与者戴着手套和不戴手套都被交叉了。我想使用参与者变量作为随机斜率的随机效果(但不是截取，因为他们一开始就没有细菌)。我想不出使用随机斜率而不是随机截距的随机效果的语法。你能教我怎么做吗？

到目前为止，我有：

require(MASS)
PQL <- glmmPQL(bacteria ~ Gloves+ NumberContacts, ~1|Participant,
               family =    gaussian(link = "log"),
               #weights = varIdent(form = ~1 | NumberContacts),
               #correlation = corAR1(NumberContacts),
               data = na.omit(Ksub), 
               verbose = F) 

每次接触后手指上的细菌

​

​

每次接触后手指上的细菌密度图

​

Answer 1

参见https://bbolker.github.io/mixedmodels-misc/glmmFAQ.html#model-specification，其中指出(0+x|group)或(-1+x|group)指定了“组内x的随机斜率:截距没有变化”。

以下示例中的机型规格是等效的：

library(MASS)
library(lme4)

fm1 <- lmer(Reaction ~ Days + (0 + Days | Subject), sleepstudy)
fm2 <- glmmPQL(Reaction ~ Days, random = ~ 0 + Days | Subject, 
  family = gaussian, data = sleepstudy)