Rを使って負の二項分布で回帰する Negative binomial regression

裾の長いカウントデータを扱いたい（外れ値が有る）ときに、負の二項分布を仮定したモデリングを考える。

MASS packageにはいっていたので、覚書として残しておく。

R: Fit a Negative Binomial Generalized Linear Model

以下で実行できる。

library(MASS)
quine.nb1 <- glm.nb(Days ~ Sex+Age, data = quine)
confint(quine.nb1)