使用lme4 glmer函數進行不平衡處理比較結果的可變長度錯誤

Question

我正在使用lme4包使用二進制響應運行比例數據的廣義線性混合模型。我有我的處理不公平的樣本量和正在以下錯誤，我的理解是，由於張女士我有不等的樣本大小：

Error in model.frame.default(data = POL3, drop.unused.levels = TRUE, formula = X2 ~ : variable lengths differ (found for '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可以處理不平衡的樣本，但不能得到它的工作。

Answer 1

無法肯定地說沒有一個重複的例子，但你可能需要確保該Trtmt變量包含在POL3（即不存在另一個Trtmt變量在全局工作區躺在附近）。

我可能會採用這一模式以這樣的方式

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