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

目的
ワイブル分布のパラメータを最尤推定する。

使用法
weibull.par(x, epsilon=1e-7)

引数
x データベクトル
epsilon 収束判定値（省略時には 1e-7 が仮定される）
max.loop 収束計算の上限回数（省略時には 500 回が仮定される）

このプログラムで計算される尺度パラメータは，
分布関数が F(x) = 1 - exp(- (x/b)^a) と表現されるときの b である.

cf.
http://blog.goo.ne.jp/hideunuma/e/b9f92b3183c7eef5767d94867fc90b1b

References
ワイブル分布のパラメータの最尤推定
http://aoki2.si.gunma-u.ac.jp/R/weibull-par.html