全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件
1248 3
2014-02-13
最近在研究RFM随机模型的时候碰到了个问题。用gamma-gamma模型对每次平均交易金额建模后,无法成功估计出三个参数。不是参数超大就是为负数。我建立对数拟然函数后分别用excel中的规划求解与eviews的最大对数拟然法来进行参数的估计都不行(PS:假设NBD分布的购买频次无论用规划求解还是最大对数拟然法都可以成功估计出参数,并且这两种方法的估参结果一致)。
excel规划求解对数拟然函数为:=LOG((EXP(GAMMALN($H$2+$I$2))/EXP(GAMMALN($H$2))*EXP(GAMMALN($I$2)))*(D2/($J$2+D2))^($H$2)*($J$2/($J$2+D2))^$I$2*(1/D2))。每个拟然函数的贡献值最后会被加法求和,求和单元格为$G$2。$H$2、$I$2、$J$2为待估参数的单元格地址,D2为每次平均购买金额。将$G$2作为目标然后,将$H$2、$I$2、$J$2作为变动参数进行规划求解,参数无法成功算出。国内的几篇论文都是用的这个方法并且都成功估计出参数。

eviews最大对数拟然法的命令为:
@logl logl1
logl1=log((@gamma(c(1)+c(2))/(@gamma(c(1))*@gamma(c(2))))*((z/(c(3)+z))^c(1))*((c(3)/(c(3)+z))^c(2))*(1/z))
@param c(1) 1 c(2) 2 c(3) 100
虽然试了各种参数的初始值,模型都无法收敛成功估计出参数。

不知道是我的数据问题(500个数据点、每次购买金额的频次图显示近似于正态分布)、还是数据处理时候的问题,请求做过这方面的高人给与指导与建议!对提出建设性建议的朋友送上300论坛币答谢!
二维码

扫码加我 拉你入群

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

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

全部回复
2014-2-15 12:57:20
我谈谈我用EXCEL试算(即求最大值)的发现:
1. 你那个EXCEL公式可能有错,你的“=LOG((EXP(GAMMALN($H$2+$I$2))/EXP(GAMMALN($H$2))*EXP(GAMMALN($I$2)...”我觉得第3个EXP前面应该是除号(即/)不是乘号(即*),应该是“=LOG((EXP(...)/EXP(...)/EXP(GAMMALN($I$2)...”。我这么想不是因为我对你模型的理解(我对你的模型和原理一概不知,我只顾用EXCEL求最值),而是如果不这么改,那么$i$2取正无穷小(即$I$2从右边无限接近于0),你那个函数就趋向正无穷大(即不会有最大值),这是GAMMA函数的性质造成的。
2. 建议将EXP(A)/EXP(B)/EXP(C)这样的形式改成EXP(A-B-C)这样的形式,免得EXP(A)的计算超过取值范围出现#NUM!的错误。
进行这两个修改后,我发现EXCEL的规划求解能算出最值(至少看上去像最大值),而且3个参数取值也看上去正常了。你也试试吧。
二维码

扫码加我 拉你入群

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

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

2014-2-15 13:07:14
TimeT 发表于 2014-2-15 12:57
我谈谈我用EXCEL试算(即求最大值)的发现:
1. 你那个EXCEL公式可能有错,你的“=LOG((EXP(GAMMALN($H$2+ ...
谢谢答复,我后来用极大拟然估计解决掉了。不过还是谢谢你的建议,请设一个300论坛币的附件我来购买。
二维码

扫码加我 拉你入群

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

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

2014-2-15 15:19:38
602dxz 发表于 2014-2-15 13:07
谢谢答复,我后来用极大拟然估计解决掉了。不过还是谢谢你的建议,请设一个300论坛币的附件我来购买。
谢谢你的诚意。既然你自己解决了问题,我就试着申请版主奖励了,不麻烦你出币了。
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

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