使用 lme4 glmer 函数进行不平衡处理比较会导致变长错误答案

【问题标题】：Using lme4 glmer function for unbalanced treatment comparison results in variable length error使用 lme4 glmer 函数进行不平衡处理比较会导致变长错误
【发布时间】：2015-05-14 01:47:10
【问题描述】：

我正在使用lme4 包来运行使用二元响应的比例数据的广义线性混合模型。我的治疗样本量不相等，并且收到以下错误，我理解这是因为我的样本量不相等：

model.frame.default 中的错误（数据 = POL3，drop.unused.levels = TRUE，公式 = X2 ~ : 可变长度不同（为 'Trtmt' 找到）

这是导致错误的代码：

#Excluding NA from the data set
POL3<-na.exclude(POL)
#Indicating the binary response
X2<-cbind(POL3$CHSd, POL3$TotSd-POL3$CHSd)
#Running the model
MMCHS4<-glmer(X2~Trtmt+(1|BSD)+(1|Hgt), family=binomial, data=POL3)

我读到lme4 可以处理不平衡的样本，但不能让它工作。

【问题讨论】：

标签： r lme4

【解决方案1】：

如果没有可重现的示例，就无法确定，但您可能需要确保 Trtmt 变量包含在 POL3 中（即，您的全局变量中没有另一个 Trtmt 变量工作区）。

我可能会以这种方式实现模型：

glmer(CHSd/TotSd~Trtmt+(1|BSD)+(1|Hgt),
      weights=TotSd,
      family=binomial,
      na.action=na.exclude,
      data=POL)

【讨论】：

谢谢！您的解决方案有效。现在我必须处理完全不同的非整数数据。非常感谢您的帮助。