差分
このページの2つのバージョン間の差分を表示します。
r:mass:fitdistr [2018/12/11 17:50] – created watalu | r:mass:fitdistr [2018/12/11 19:38] (現在) – watalu | ||
---|---|---|---|
行 52: | 行 52: | ||
< | < | ||
fitdistr(c(1, | fitdistr(c(1, | ||
+ | </ | ||
+ | |||
+ | < | ||
+ | library(MASS) | ||
+ | </ | ||
+ | |||
+ | [[r: | ||
+ | < | ||
+ | fitdistr(X$time, | ||
+ | </ | ||
+ | その結果、次の数値を得る。 | ||
+ | < | ||
+ | | ||
+ | 1.6822280 | ||
+ | | ||
+ | </ | ||
+ | ワイブル分布には形状母数shapeと尺度母数scaleがある。それぞれの点推定値は1.6822280と6.8219220であり、標準誤差は括弧の中に表示される。 | ||
+ | |||
+ | これらは、次のようにしても表示される。 | ||
+ | < | ||
+ | X.fit.weibull <- fitdistr(X$time, | ||
+ | names(X.fit.weibull) | ||
+ | X.fit.weibull$estimate | ||
+ | X.fit.weibull$sd | ||
+ | X.fit.weibull$loglik | ||
+ | </ | ||
+ | 最後の数値は、最尤推定法で最大化した対数尤度関数の最大値である。この値を参考に、寿命分布を選ぶこともある。 | ||
+ | |||
+ | 対数正規分布も推定しておく。 | ||
+ | < | ||
+ | X.fit.lognormal <- fitdistr(X$time, | ||
+ | names(X.fit.lognormal) | ||
+ | X.fit.lognormal$estimate | ||
+ | X.fit.lognormal$sd | ||
+ | X.fit.lognormal$loglik | ||
+ | </ | ||
+ | |||
+ | 同様に、ガンマ分布。 | ||
+ | < | ||
+ | X.fit.gamma <- fitdistr(X$time, | ||
+ | names(X.fit.gamma) | ||
+ | X.fit.gamma$estimate | ||
+ | X.fit.gamma$sd | ||
+ | X.fit.gamma$loglik | ||
</ | </ | ||