GEE中的负二项式

Question

对于实施GEE的R包，例如gee，geepack，似乎不包括负二项式族。我有两个问题：

1. GEE还有其他我不知道的R套餐吗？

2. 如果没有，是否有一个简单的步骤允许创建一个系列，即提供链接函数（log mu）和方差函数（mu + mu^2/theta），假设{{1指定（否则NB不是GLM），然后让theta或gee代码以与geepack类似的方式开展业务？

Answer 1

您应该能够使用negative.binomial包中定义的MASS系列来执行此操作（设置具有指定theta值的NB系列）。看起来geepack::geese（至少）会接受此表单中的系列规范。要估算θ，您可以尝试使用固定的theta将GEE拟合嵌入循环中，或者制作geefit_NB(theta)函数并优化theta。

如果negative.binomial中尚未存在MASS，您可以定义自己的家庭（这无疑是有点先进的 - 我首先要下载MASS的源代码打包并查看文件R/neg.bin.R）。