r中的负二项式函数

Question

我正在比较四个不同地区的发病率，以确定它们是否相等。首先我拟合泊松模型，但我不确定如何在使用偏移时测试方差=平均假设。因此，我想要拟合负二项式函数，并将两者与似然比检验进行比较。但是，glm.bn中的MASS函数会引发错误。有什么想法是错的吗？

library(MASS)
b <- data.frame(
  s=c(1800,539,490,301),
  pop=c(2900000,1327000,880000,268000),
  reg=c("A","B","C","D")
)

glm.nb(s~reg+offset(log(pop)),data=b)

结果

Error in while ((it <- it + 1) < limit && abs(del) > eps) { : 
  missing value where TRUE/FALSE needed