Re: [問題] 使用R來解估計方程式
※ 引述《caraf.bbs@bbs.wretch.cc ()》之銘言：
: ※ 引述《qqwerty.bbs@ptt.cc (love ocean)》之銘言：
: > 最近剛剛接觸R，之前都是使用Gauss來寫程式
: > 之前的程式老師希望我能用R重新在寫一次
: > 但我碰到一個問題是，我寫出了一個估計方程式來解參數(a,b)
: > 例如
: > g1(a,b,x,y)=0
: > g2(a,b,x,y)=0
: > (x,y) is data
: > 之前在Gauss可以使用eqSolve來解這個連立方程式
: > 但在R裡卻不知道該用什麼指令來解
: > 不知道有沒有了解的人可以幫忙一下?...感謝了!
: 先令你的 objective function 為
: g1(a,b,x,y)^2 + g2(a,b,x,y)^2
: 再利用 optim() 來解未知參數 a, b
嗯～首先感謝大家熱心的幫忙，提供了相關的指令能參考
不過我有一些問題想請教一下
大部份推文所提到的指令都是對於某個f(a,b,x,y)找出使
這個函數最小的a,b的值
例如在R中的Example所提的f(x1,x2)=100(x2-x1^2)^2+(1-x1)^2
找出x1,x2使f(x1,x2)最小
用optim()可以求出為(1,1)
轉成微積分的問題可視為在解
f1(x1,x2)=x2-x1^2=0
f2(x1,x2)=200x1(x2-x1^2)+(1-x1)=0
若是以上面回覆所說的
我令g(x1,x2)=f1^2+f2^2去求x1,x2使g有最小值
是否跟原來要找f有最小值的問題不同？
除此之外，我也能夠令g*=|f1|+|f2|來解
這樣看來似乎與我原先的問題又不大相同了
我的問題比較像是解f1,f2的連立方程解
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 220.133.117.103
討論串 (同標題文章)