R：使用glmer命令

Question

我正在分析在不同学校进行的调查问卷中的分类数据，看看哪些因素可能会影响学生的反应。因此，我使用R＆lt; lme4包中的glmer命令构建混合模型。对于每个调查问题回答，我有六个预测变量，我希望将学校作为随机效应包括在内，因为拦截和斜率都因学校而异。我在网上和线下都经过长时间的努力搜索，并找到了有关正确编码方式的相互矛盾的帐户，并且作为R新手，我不确定哪个是正确的！以下是我提出的内容（Like是响应变量）：

LikeM1 <-glmer(Like ~ Treatment + Regularity + Learn + Age + Gender + Organisation_Membership_Summary + (1 + Like|School), data = MagpieData, na.action = "na.omit",family = binomial(logit))

我是否正确指定学校作为随机效果，以便截距和坡度因学校而异？我或许应该提到作为分类数据，我所有的变量都是R

中的因素

Answer 1

如果您希望斜率和截距都随组而变化，则一般形式为：y〜x +（1 + x |组）。括号中的1表示截距应按组变化，x表示预测变量x的系数应按组变化。您的模型中有很多预测变量。我将一次只使用一个预测变量，以使解释变得容易些。

Answer 2

我认为您想这样做：

  

LikeM1 <-glmer（例如〜治疗+规律+学习+年龄+性别+ Organisation_Membership_Summary +（1 |学校）+（0 +治疗+规律+学习+年龄+性别+ Organisation_Membership_Summary |学校），数据= MagpieData， na.action =“ na.omit”，family =二项式（logit））

括号中公式的第一部分是随机截距，第二部分是随机斜率。 This link提供了很好的解释。