Question

我正在R中使用lme4软件包来拟合我想在调解中使用的负二项式模型（比Poisson更好）。

我收到一个错误：“不支持用于结果模型的glmer系列”

为了简单起见，我只是将类切换为lmer。例如，来自

> class(model.Y)
[1] "glmerMod"
attr(,"package")
[1] "lme4"

到

> class(model.Ynew)
[1] "lmerMod"
attr(,"package")
[1] "lme4"

，并且有效。显然，这可能会产生问题，但是我检查了残差，拟合值和预测值，并且没有变化。结果也与泊松模型结果吻合。

> summary(residuals(model.Y)-residuals(model.Ynew))
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
      0       0       0       0       0       0 
> summary(fitted(model.Y)-fitted(model.Ynew))
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
      0       0       0       0       0       0 
> summary(predict(model.Y)-predict(model.Ynew))
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
      0       0       0       0       0       0

有人有没有洞悉这种影响或不是临时性的可能方法的结果？

Answer 1

如果您通过负二项式回归分析进一步探讨调解，并且可能从中学到的知识，我会很感兴趣

带有调解的负二项式混合模型

1 个答案: