全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
13912 19
2016-02-09
悬赏 40 个论坛币 未解决
         最近在学习使用r语言进行极大似然估计的操作,知道可以使用optim、nlminb和MaxLix包来对未知参数进行估计,但发现这三个函数都只能给出参数的估计值和似然函数最优值,并不能给出待估参数的近似标准差,造成后续的参数检验无法进行。        请问有没有办法能直接输出参数标准差呢?或者怎样编写语言计算Fisher信息矩阵在参数估计值处的逆矩阵?
        在此请教大家,希望能有高手赐教!先行谢过!


二维码

扫码加我 拉你入群

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

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

全部回复
2016-2-9 19:56:04
没人回复,自己先来顶一下!!
另外刚才发现maxLik可以输出hess值,这个矩阵是我想要的信息阵吗?
二维码

扫码加我 拉你入群

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

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

2016-2-9 20:28:10
可以考虑使用bootstrap的方法,这种方法结合最大似然法就可以给出待估参数的标准差(其实是该参数的整个分布).但是,我不知什么软件来实现。如果是我就自己编程罗。
二维码

扫码加我 拉你入群

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

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

2016-2-9 21:01:48
参数标准差就是信息阵的相应对角元(例如,如果是第一个参数那么就是对角线第一个数)的开方
二维码

扫码加我 拉你入群

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

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

2016-2-9 22:43:36
oliyiyi 发表于 2016-2-9 21:01
参数标准差就是信息阵的相应对角元(例如,如果是第一个参数那么就是对角线第一个数)的开方
是的,这个我明白,但我不知道调用hess值得到的是否为信息阵呢,谢谢您!
二维码

扫码加我 拉你入群

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

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

2016-2-10 06:52:36
hopelesscat 发表于 2016-2-9 22:43
是的,这个我明白,但我不知道调用hess值得到的是否为信息阵呢,谢谢您!
求MLe时一般用N-L算法,迭代矩阵就是Hess阵,也就是信息阵
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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