gepuro.net
gepulog

データ分析エンジニアによる備忘録的ブログ

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

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

参考


似てる記事

似てない記事