泊松和负二项式回归何时拟合相同的系数?
我已经注意到,在R中,泊松和负二项式(NB)回归似乎总是适合相同的系数,以用于分类但非连续的预测变量。 例如,这是带有分类预测变量的回归: data(warpbreaks) library(MASS) rs1 = glm(breaks ~ tension, data=warpbreaks, family="poisson") rs2 = glm.nb(breaks ~ tension, data=warpbreaks) #compare coefficients cbind("Poisson"=coef(rs1), "NB"=coef(rs2)) 这是一个连续预测变量的示例,其中泊松和NB拟合不同的系数: data(cars) rs1 = glm(dist ~ speed, data=cars, family="poisson") rs2 = glm.nb(dist ~ speed, data=cars) #compare coefficients cbind("Poisson"=coef(rs1), "NB"=coef(rs2)) (当然,这些不是计数数据,模型也没有意义...) 然后,将预测变量重新编码为一个因子,然后两个模型再次拟合相同的系数: library(Hmisc) speedCat = cut2(cars$speed, g=5) #you can change …