ワイブル累積密度分布からの累積確率があり、基礎となる関数パラメーターを見つけたいと考えています。
データ：
At <- data.frame(cumu = c(0.004722245,0.004740720,0.008768300,0.010900330,0.194195851,0.215752640,0.226863587,0.237534826,0.240753194,0.324707908,0.354541558,0.421462196,0.421879734,0.423690298,0.454939143,0.597567194,0.647963670,0.688549850,0.729575737,0.742892615,0.758127212,0.760495725,0.762376494,0.763717789,0.764586120,0.985689160,0.993160505,0.995931775,1.000000000),yr = c(1,2,4,5,6,7,8,9,10,11,12,13,14,15,16,17,22,32,33,34,35,36,37,38,40,41,42,43,46))
nls 関数を使用して MLE によるパラメーター フィッティングを試みましたが、エラーがスローされます。
test <- nls(At$cumu ~ c * exp(-(At$yr/b)^a), start = list(a = .01, b = .2, c = 1))
エラー：
numericDeriv(form[[3L]], names(ind), env) のエラー: モデルの評価時に値が欠落しているか、無限大が生成されました
問題は開始値が間違っていることだと思います。私の質問は、私のデータに適した開始値はどれですか?
ありがとう、
サイモン