データセットの形状と規模を推定しようとしています。
2つの異なる方法を使用しましたが、どちらの場合もエラーメッセージが表示されました。
surveyまず、パッケージを使って少しずつ試してみました。
survreg(Surv(all.ws)~1, dist="weibull")
エラーメッセージが表示されました：
この分布の無効な生存時間
次に、fitdistr()関数を使用してみました。
fitdistr(all.ws, densfun=dweibull, start=list(scale=1, shape=2))
エラーメッセージが表示されました：
optim（x = c（2.2、2.1,1.9 ....：
'vmmin'の初期値が有限ではありません
私がしていることの何が問題になっていますか？
Google 検索「fitdistr Weibull Error」は、この正確な質問が 1 年前に R-help メーリング リストで議論されたことを示しています: http://r.789695.n4.nabble.com/Problems-with-fitdistr-td1334772.html
そのリンクからのいくつかのポイント：
データにゼロがあると問題が発生します
パッケージ lmom の pelwei() 関数を使用する
ベータ版ディストリビューションで fitdistr() を使用すると、同様の問題が発生しました。その場合、データ内の 1 と 0 の両方でこのエラーが発生しました。
さらに、 limit 引数が fitdistr() 呼び出しで使用されると、別のエラーが生成されることがわかりました: たとえば (ここで、x は 1.0 または 0 を含むサンプルのベクトルです):
fitdistr(x, "beta", list(shape1 = 1, shape2 = 0.2),lower=0.001)
`Error in stats::optim: L-BFGS-B needs finite values of 'fn'`