GEE在R中的负二项分布模型数据有更新吗？

目标：使用负二项分布运行GEE以分析R中的数据集

有人问过这个问题，但答案已经有5年了，我想知道是否有新的发展。我总结了我的数据，然后概述了我是如何试图找到解决方案的

我的模型（使用泊松分布而非负二项分布成功运行，并使用包geeglm）：

glmmPQL（Diract~DirPar*RecPar+offset（LnScan）+Year，random=~1 |中出错：找不到函数“corCompSymm”

---

如果有人有建议或解决方案，我愿意接受任何建议。非常感谢！

多年来，这个问题被问了很多次，我在遇到同样的问题并搜索了几个小时的答案后发现了它。两个解决方案：

geeM

包允许使用负二项族

通过R

networkite

包导入python

statsmodels

模块，并调用R中允许负二项式族的gee拟合函数

m1 <- geeglm(Diract ~ Dir*Rec + Year + offset(LnScan), 
             family = poisson("log"), data = Direct, id = ID, corstr = "exchangeable")

m2 <- glmmPQL(Diract ~ Dir*Rec + offset(LnScan) + Year, 
              random = ~ 1 | ID, 
              family = negative.binomial(1), 
              data = Direct, 
              correlation = corCompSymm(form=~1|ID))