作者phil5566 (5566)
看板Statistics
標題[問題] R內建的optim有無找不出最佳解?
時間Thu Nov 3 19:02:02 2016
如題
最近有個問題困擾蠻久
就是當我在建立模型來fit現有的data時
使用mle的方法來解模型參數時
由於函數太複雜
勢必要透電腦求算
洽好R內建的optim可以解決這個問題
但由於初始值給的不同,解出來的答案也不同
所以我就把現有的解出的解一個一個代入模型去畫圖
看看fit的程度好不好
無意中代入非現有的解值來畫圖
發現都比現有的解值畫圖的結果
fit的還要好.....
於是把該解當作optim的初始值代入
但跑出"In .... 產生了 NaNs"的訊息
想請教一下optim是否有可能解不出最佳解的情況?
又或是在解mle大家都如何找出最佳解的?
是一個一個試找出最合理的解?還是有其他比較好的方法?
指教了,謝謝
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.112.74.231
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Statistics/M.1478170924.A.B4E.html
1F：→ celestialgod: hello 標題 標題 提醒幾次看板規了？ 11/03 19:12
2F：→ celestialgod: MLE就乖乖算微分解吧，先確定是convex函數 11/03 19:13
3F：→ phil5566: 抱歉手滑按錯,函數很複雜要怎麼手算微分?謝謝 11/03 20:04
4F：→ LiamIssac: 有可能不是convex 先確認這個 再來算微分解 11/03 20:58
5F：推 JamesChen: 當然有可能... 11/04 06:44
6F：→ simon552614: 用PSO丟下去吧 11/08 00:19