作者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