r - GEE中的负二项式

Question

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

1. 是否还有其他我不知道的 GEE R 包？

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

Answer 1

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

如果negative.binomial在 中不存在MASS，您可以定义自己的系列（诚然，这有点高级——我将从下载MASS包的源代码并查看文件开始R/neg.bin.R）。