全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
13035 21
2011-01-28
我有一些关于用R计算MLE的问题:我知道的命令有optim,nlm,constrOptim,但是有很多困难。

1.有非线性的constraint怎么办,比如我的参数其中有6个是d11,d21,d22,d31,d32,d33,条件是:用这6个参数组成一个矩阵3*3的对称矩阵Q,要求Q是一个正定矩阵,如何在估计中加入这个非线性的constraint?
2.有什么比较科学的方法设定initial value么?
3.log likelihood function比较复杂,而且要进行simulation,那个命令是最节省时间的?

请高手不吝赐教,万分感谢!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2011-1-28 04:55:23
lz不要只关注那些一步计算出来的函数嘛。我们需要的东西大部分不是那么容易算出来的,需要自己处理,R只是一个高级的工具,但不是万能钥匙。

对于quadratic的可行域,LZ可以用矩阵分解自己推导一下,一般来说可以通过特征值分解或者Cholesky分解将可行域化为比较好的形式(不如高维球),然后进行相应处理(比如Schwaz Inequality)。或者寻找Nonlinear optimization的package(不过R在这方面的功能还很不完善)。

loglikelihood 可以自己编程用 gradient descent 或者 牛顿法做,十分方便,code也就5-10行。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2011-1-28 05:20:12
2# ltx5151

呵呵,谢谢,很有启发。
我一直懒得自己编,现在发现用现成的package总是被动的,一劳永逸的方法是提高能力能自己编出来。到处找package的时间还不如自己写code呢。
需要开始自己想办法了
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2011-1-28 09:42:05
Nonlinear optimization的package
有一些package可以考虑试一下
比如 Rdonlp2, Rsolnp和 alabama
当然可能不一定能够直接用
楼主要加的约束似乎不常见
初值是一个永远的痛呀
每次都要费好多力气
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2011-1-28 21:03:08
如果想自己动手编写MLE,不妨考虑一下ADMB: http://admb-project.org/
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2011-2-12 14:49:28
学习下,多谢了
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群