r - 如何使用 R 在 GLM 中指定负二项式误差分布？

Question

我正在使用 R 中的 glm() 函数构建模型。假设我知道我的数据具有符合负二项分布的误差分布。

当我在 R 手册中搜索各种系列时，family=binomial 作为一个选项提供，但 negative binomial 不是。

在 R 手册(系列)的同一部分中，NegBinomial 在“另请参阅”部分中链接，但它是在二项式系数的上下文中呈现的(我什至没有确定这是指什么)。

因此，总而言之，我希望找到类似于 glm(y~x, family=negbinomial, data=d,na.omit) 的语法。

Answer 1

对于未知的过度离散参数，负二项式不是负指数族的一部分，因此不能作为标准 GLM(或 glm())进行拟合。 MASS 包中有一个glm.nb() 函数可以帮助您...

library(MASS)
glm.nb(y~x, ...)

如果您碰巧有一个已知/固定的过度分散参数(例如，如果您想要拟合具有 theta=1 的几何分布模型)，您可以使用 negative.binomial 来自 MASS 的家族:

glm(y~x,family=negative.binomial(theta=1), ...)

如果 MASS::glm.nb 位于 ?glm 的“另见”部分，可能不会有什么坏处 ...