作者Wush978 (柴人)
看板Statistics
標題Re: [問題] 用R作MLE時 有辦法加入限制條件嗎?
時間Sat Jul 11 16:35:36 2009
※ 引述《KirinGuess (Kirin)》之銘言:
: 假設我想運用MLE估計3個參數,p1、p2與p3,
: 如果是使用R來作這個MLE,
: 怎麼才能對參數估計值設定(p1<p2 & p2<p3)的限制條件?
: 我目前只會模仿R官網提供的教學手冊,
: 用nlm估計沒有限制條件的p1、p2、p3。
: 另外,我有試過optim並選用L-BFGS-B的方法,
: 可是它的上下界設定好像只能使用實際的值,
: 估計的參數不能被當成限制條件。
: 剛才試了一下,
: nlminb好像也有類似的情況。
: (也可能是我不會用)
: 麻煩大家解答了,謝謝。
nlminb 可以在找minimal 時可以設定lower, upper
限制參數的範圍
以你的需求(p1 < p2 & p2 < p3)
我建議你先做一下變數變換
把p1, p2-p1, p3-p2 當成變數
這樣變成後兩者的值要>0
接下來再套用nlminb應該就行了
--
文章po了, 沒有人關心是頂難受的事情.
但是刻意去關心又是一件不自在的事情.
所以就保持緘默吧... 但要是真的喜歡,
你就摸摸右耳; 不然, 就摸左耳好了. 不知道? 那就摸鼻子吧
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.114.34.64
1F:推 KirinGuess:非常謝謝,我去試試看。 07/11 18:41