R言語でガンマ分布のパラメータを推定する方法を2通り示す。
データの生成
> x <- rgamma(1000, shape=2, rate=0.5)最尤法による推定
> fitdistr(x, "gamma")
shape rate
1.81070179 0.45850801
(0.07471446) (0.02177244)モーメント法による推定
> med.gam <- mean(x)
> var.gam <- var(x)
> (shape <- med.gam/var.gam)
[1] 0.45973
> (rate <- ((med.gam)^2)/var.gam)
[1] 1.815528参考