отрицательная биномиальная функция в r

Я сравниваю заболеваемость в четырех разных регионах, чтобы увидеть, равны они или нет. Сначала я подогнал модель Пуассона, но я не уверен, как я могу проверить предположение о дисперсии = среднему, когда я использую смещение. Следовательно, я хотел подогнать отрицательную биномиальную функцию и сравнить их с помощью теста отношения правдоподобия. Однако функция 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