最近在研究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论坛币答谢!