R语言glm.nb负二项回归设置offset参数报错问题咨询

报错的核心原因 #

你使用的glm.nb()函数来自R的MASS包，该函数的参数列表中没有独立的offset参数，offset项必须嵌入到模型公式内部，通过offset()函数声明，不能作为独立参数传入函数，这是你触发unused argument (offset = L4)报错的直接原因。
另外你原代码没有指定data = doctor，如果数据集变量没有被挂载到全局环境，也会额外触发变量找不到的报错。

修正后的代码示例 #

基础修正（匹配你当前写的模型结构） #

library(MASS)
nbinom <- glm.nb(V3 ~ V2 + V5 + V2:V5 + offset(L4), data = doctor)

匹配AIC筛选的最优模型结构 #

如果需要保留你之前通过AIC得到的三向交互项V2:L4:V5，仅将L4的主效应设为offset，调整为如下写法即可：

library(MASS)
nbinom <- glm.nb(V3 ~ V2 + V5 + V2:V5 + V2:L4:V5 + offset(L4), data = doctor)

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