[問題] optim函數代不同初始值解出的解都不同?

作者: phil5566 (5566)   2016-10-22 22:01:36
[問題類型]:
程式諮詢(我想用R 做某件事情,但是我不知道要怎麼用R 寫出來)
[軟體熟悉度]:
請把以下不需要的部份刪除
新手(沒寫過程式,R 是我的第一次)
[問題敘述]:
我想解 likelihood 函數裡的變數,
打好了自定的function後,要解fr的X1,X2,X3
用了optim(C(a,b,c),fr)的指令
但解出來的值,根據給的初始值a,b,c不同而不同
甚至還會出現Error in optim 不可使用初始參數來評估函數
讓我感到困擾,究竟是哪一個解才是正確的
如果用牛頓疊代法解多變數,要怎麼用R來寫?
請高手指點迷津,謝謝
程式碼可貼於以下網站:
http://ideone.com/0MCowj
[環境敘述]:
Win7 Rx64 3.3.1
作者: carl090105 (Jing)   2016-10-22 22:20:00
就求最佳值問題需要注意的是所求的解是local還是global極值,這可能跟初始值、迭代次數等有關,可以參考https://en.m.wikibooks.org/wiki/R_Programming/Optimization 有一些相關的求極值套件
作者: Edster (Edster)   2016-10-23 11:09:00
做敏感度測試
作者: phil5566 (5566)   2016-10-23 12:20:00
c大的意思是從調optim裡的初始值、迭代次數嗎?如果希望X1解出來介於0~1之間有辦法辦到嗎?E大的意思我不太懂?要怎麼把敏感度測試用在解極值?

Links booklink

Contact Us: admin [ a t ] ucptt.com