加重最小二乗法を使用して、GEV (一般化極値) 分布のパラメーターを推定したいと考えています。私は R を使用しており、この目的に使用できると思われる nls という関数を見つけました。数式とオプションのデータセットを要求します。GEV 式と年次最大値シリーズがここにあるはずですが、どうすればよいかわかりません。誰も nls を使用していて、これを行う方法について何か考えがありますか?
#Vector of ranged annual maxima
x <- c(21,24,29,32,32,34,35,35,35,36,37,37,38,40,40,41,43,47,47,52)
w <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2)
data <- list(x=x,w=w)
nls(y ~ exp(-(1+((x-location)/scale))^(-1/shape)),data=data, weights=w,start=list(location=5,scale=2,shape=0.10))
y が見つからないというエラーが表示されます。y は、GEV パラメーターを最適化したときに得られるものであり、すべての x に対して y が可能な限り x に近づくようにします (これも重みによって異なります)。したがって、GEV パラメータを推定するまで y は不明です...