r - r中的负二项式函数

Question

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

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