July 22, 2014

ガンマ分布のパラメータ推定をする

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

参考

© gepuro 2013