optim関数で数値的に逆関数の計算したい - My Life as a Mock Quant

optim関数で数値的に逆関数の計算したい - My Life as a Mock Quant

これがオーソドックスな正しいやり方なのかわからないけど、とりあえずやってみた。 optim関数で数値的...概要を表示 これがオーソドックスな正しいやり方なのかわからないけど、とりあえずやってみた。 optim関数で数値的に逆関数の計算したい(一次元版) とりあえず一次元の関数をと書いた時の、が与えられた下でのの推定を考えたい。 これはもちろん解析解があるとハッピーなわけだが、ない場合には数値的にやるしかないわけです。ここでは適当にって関数を題材にやってみんとす。 この関数をRで書くと以下のような感じか。 #適当な関数f f <- function(x) { A <- 1:2 A[1]*sin(x) + A[2]*sin(x)^2 } 描画すると、こんな感じ。 a <- seq(-5, 5, 0.1) plot(a, sapply(a, f)) とりあえずこの関数に対する逆関数を計算してみたい。 そこで となるようなを計算するコードをoptim関数*1を使って以下のように書いた。要するに「最小二乗法でそれっ