手指上的细菌。在MASS :: glmmPQL中具有随机斜率但没有截距的交叉随机效应的语法

Question

我有非正常数据（戴着和不戴着手套接触表面后手指上的细菌），因此请使用MASS软件包中的glmmPQL。我有一个类别预测变量（Gloves），一个重复的测量变量（NumberContacts），参加实验的参与者戴着手套，没有手套，因此被交叉了。我想将Participant变量用作具有随机斜率的随机效果（但不要拦截，因为它们的细菌总数为0）。我无法弄清楚具有随机斜率但没有随机截距的随机效果的语法。你能告诉我怎么做吗？

到目前为止，我有：

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)