[求助]麻烦高手看下我的极大似然方法估计分布问题出在哪

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我要用极大似然方法对收入的分布进行估计，假定服从gamma分布

    di=read.table("di.txt",header=TRUE)

    attach(di)   #di中有变量dio2和di03

    library(stats4)

    ll=function(lambda,alfa){n=500

                                       x=di02

    -n*alfa*log(lambda)+n*log(gamma(alfa))-(alfa-1)*sum(log(x))+lambda*sum(x)}

   est=mle(minuslog=ll,start=list(lambda=2,alfa=1))

  summary(est)

运行后提示大量出错信息：

Error in optim(start, f, method = method, hessian = TRUE, ...) :

  non-finite finite-difference value [1]

In addition: There were 50 or more warnings (use warnings() to see the first 50)

估计结果也不对，麻烦高手给看看问题出在哪，谢谢。

[此贴子已经被作者于2009-3-30 13:47:35编辑过]

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：极大似然 difference warnings function addition 高手 麻烦

你想过没有,收入是GAMMA分布,GAMMA分布的取值范围在什么地方