※ 引述《[email protected] ()》之銘言： : ※ 引述《[email protected] (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