R -- ワイブル分布のパラメータの最尤推定

ワイブル分布のパラメータの最尤推定 Last modified: Dec 07, 2004 目的 ワイブル分布のパラメータを最尤推定する。 使用法 weibull.par(x, epsilon=1e-7) 引数 x データベクトル epsilon 収束判定値（省略時には 1e-7 が仮定される） max.loop 収束計算の上限回数（省略時には 500 回が仮定される） ソース 大きさ n の標本 x1, x2, ..., xn から，α および m の最尤推定量を求めるには， 1/α = n / Σ xim m = n / (Σ(xim * ln(xi))-Σln(xi)) を，1/α と m に関して反復法により求める。 インストールは，以下の 1 行をコピーし，R コンソールにペーストする source("http://aoki2.si.gunma-u.ac.jp/R/src