0

4つの異なる地域での発生率を比較して、それらが等しいかどうかを確認しています. 最初にポアソン モデルを当てはめましたが、オフセットを使用した場合に分散 = 平均の仮定をどのようにテストできるかわかりません。したがって、負の二項関数を当てはめ、2 つを尤度比検定で比較したいと思いました。ただし、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
4

0 に答える 0