R语言负二项回归stepAIC输出AIC与直接计算结果不一致咨询

• 核心原因是MASS::stepAIC处理负二项回归模型的参数固定机制
  运行逐步筛选时，stepAIC会默认沿用初始模型（你代码里的空模型out0）估计得到的负二项分布离散参数theta，计算所有候选模型的对数似然与AIC，不会在每一步增减变量时重新估计该变量组合下的最优theta。而你手动拟合全变量模型时，glm.nb会基于当前的变量集重新估计最优theta，两个场景下theta取值不同，最终计算得到的AIC自然会存在差异，你遇到的差值为2属于该机制下的典型情况。
• 额外注意代码调用顺序问题
  你给出的示例代码中，先调用glm.nb拟合空模型，再加载MASS包，而glm.nb是MASS包的专属函数，如果运行代码前没有提前加载过MASS，空模型拟合步骤会直接报错。建议统一将包加载代码放在所有对应包函数调用的最前端，避免环境混乱导致的参数估计偏差。

如果需要消除该差异，可以手动拟合所有候选变量组合的负二项回归模型，分别计算AIC后再比较，就可以得到和单独拟合一致的结果。

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