2006-02-18

最尤推定や最小自乗法等、ある目的関数の最適化に用いる汎用関数 RWikiの該当ページはこちら Rの Htmlヘルプ"nlm{stats}"も参考のこと ごくごく簡単に 今、に近似できるデータを持っているとする。aの値がわからないので、それを推定したいとする。データは datay,dataxというペアデータとしてもっているとする。 今、このdatay,dataxのデータセットについて、が、yの値の差であり、その自乗をすべてのデータペアについて足し合わせた値がもっとも小さくなるような a の値を求めたいものとする 最小にしたい値を計算する関数 f を次のように書く f<-function(x,y,a) { sum((y-a*x^2)^2) } aの推定にあたっては、この式と、推定を開始するためのaの初期値と、データセットをnlm関数に次のように渡す nlm(f,3,x=datax,y=data

[ma_ko]

nlm／optim みたいなやつ