[程式] R內建的optim有無找不出最佳解?

看板Statistics作者 (5566)時間9年前 (2016/11/03 19:02), 編輯推噓1(105)
留言6則, 5人參與, 最新討論串1/1
如題 最近有個問題困擾蠻久 就是當我在建立模型來fit現有的data時 使用mle的方法來解模型參數時 由於函數太複雜 勢必要透電腦求算 洽好R內建的optim可以解決這個問題 但由於初始值給的不同,解出來的答案也不同 所以我就把現有的解出的解一個一個代入模型去畫圖 看看fit的程度好不好 無意中代入非現有的解值來畫圖 發現都比現有的解值畫圖的結果 fit的還要好..... 於是把該解當作optim的初始值代入 但跑出"In .... 產生了 NaNs"的訊息 想請教一下optim是否有可能解不出最佳解的情況? 又或是在解mle大家都如何找出最佳解的? 是一個一個試找出最合理的解?還是有其他比較好的方法? 指教了,謝謝 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.112.74.231 ※ 文章網址: https://www.ptt.cc/bbs/Statistics/M.1478170924.A.B4E.html

11/03 19:12, , 1F
hello 標題 標題 提醒幾次看板規了?
11/03 19:12, 1F

11/03 19:13, , 2F
MLE就乖乖算微分解吧,先確定是convex函數
11/03 19:13, 2F

11/03 20:04, , 3F
抱歉手滑按錯,函數很複雜要怎麼手算微分?謝謝
11/03 20:04, 3F

11/03 20:58, , 4F
有可能不是convex 先確認這個 再來算微分解
11/03 20:58, 4F

11/04 06:44, , 5F
當然有可能...
11/04 06:44, 5F

11/08 00:19, , 6F
用PSO丟下去吧
11/08 00:19, 6F
文章代碼(AID): #1O6naijE (Statistics)