全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
9522 11
2010-01-24
我用stats4包中的mle求负指数方程的2个参数,但是总出问题。
负指数方程f(x)=a*exp(-bx)
library(stats4)
x.fuzhi=c(112,80,49,45,42,35,34,29,23,17,14,12,13,5,3,3,1)
ll<-function(a,b){n<-17
x<-x.fuzhi
-n*log(a)+b*sum(x)}##
负对数似然方程
est<-mle(minuslog=ll, start=list(a=2,b=0.05))
summary(est)


结果总是Error in solve.default(oout$hessian) :
Lapack routine dgesv: system is exactly singular


请问,哪里出了问题呢?
二维码

扫码加我 拉你入群

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

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

全部回复
2010-1-24 20:01:59
hessian 矩阵式奇异的
二维码

扫码加我 拉你入群

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

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

2010-1-24 20:04:48
请问,怎么解决呢?我的似然方程不对吗,还是这个数据用极大似然估计不出来?
二维码

扫码加我 拉你入群

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

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

2010-1-25 10:50:03
可以用 modified Marquardt 方法,sata我不会用,但以前用matlab做过优化碰到过类似的情况,但愿对楼主有所帮助



sata网上查到此信息,最优化设置选项

moptimize_init_singularHmethod(M, singularHmethod) specifies the            action to be taken during optimization if the Hessian is found to            be singular and the technique requires the Hessian be of full            rank.  Allowed values are
               singularHmethod   description               --------------------------------------------------------               "m-marquardt"     modified Marquardt algorithm               "hybrid"          mixture of steepest descent and Newton               --------------------------------------------------------               The default is "m-marquardt".               "hybrid" is equivalent to ml's difficult option; see [R] ml.
二维码

扫码加我 拉你入群

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

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

2010-1-25 11:04:18
syslon
     谢谢你的回复,还不知道在R中如何处理这种情况,但我觉得数据应该没问题,只是找不到原因。
二维码

扫码加我 拉你入群

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

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

2010-1-25 11:20:39
遇到这种问题很正常,二阶导数出问题了,用修正的 Marquardt方法就可解决此问题
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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