R语言glm.nb拟合负二项回归出现alternation limit reached报错求解

alternation limit reached警告本质是MASS包的glm.nb()函数在默认迭代次数上限内，没有完成负二项分布离散参数θ和回归系数的交替收敛估计，结合提供的数据集和代码，触发该问题的具体原因如下：

• 模型设定错误：L4是总暴露人年的对数，计数回归中暴露量应该作为偏移项offset()引入模型（理论上偏移项的系数固定为1，不需要模型估计），将L4作为普通协变量放入回归，导致模型需要额外估计一个本应固定的参数，参数估计空间异常，迭代无法在默认次数内收敛。
• 分类变量未正确编码：V2（1-5级年龄组）、V5（二分类吸烟状态）本质是分类变量，直接作为连续变量纳入模型，引入的设定偏误进一步加大了收敛难度。
• 样本量不足：数据集仅包含10条观测，而负二项回归需要额外估计离散参数θ，本身对样本量的要求高于普通泊松回归，小样本下迭代收敛的难度大幅提升。
• 默认迭代次数上限偏低：glm.nb()默认的最大迭代次数为25次，参数初始值与实际值偏差较大时，很容易在达到收敛阈值前耗尽迭代次数。

1. 优先修正模型设定，将L4作为偏移项放入，同时对分类变量做因子编码：

nb1 = glm.nb(V3 ~ factor(V2) + factor(V5) + offset(L4), data = doctor)

2. 若修正后仍出现警告，可手动调大迭代次数上限：

nb1 = glm.nb(V3 ~ factor(V2) + factor(V5) + offset(L4), data = doctor, control = glm.control(maxit = 1000))

3. 也可通过init.theta参数手动传入离散参数θ的初始估计值，降低迭代收敛的难度。

内容的提问来源于stack exchange，提问作者rrr