已编知情交易概率ekop模型的stata编程，用极大似然方法求解，总是报错

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

请大神赐教ekop知情交易概率pin的eb es a u c5个参数的估计，用stata怎么写命令我自己写了一个命令，但是运行到ml check总是说出错
程序如下：
capture program drop log10
program log10
version 1.0
args lnf eb es a u c
quietly replace `lnf' = -`eb'+($ML_y1+$ML_y2)*log(`eb'/(`u'+`eb'))+$ML_y1*log(`u'+`eb')-`es'+($ML_y1+$ML_y2)*log(`es'/(`u'+`es'))+$ML_y2*log(`u'+`es')+log(`a'*(1-`c')*exp(-`u')*((`eb'/(`u'+`eb'))^(-($ML_y1+$ML_y2)))*((`es'/(`u'+`es'))^($ML_y2-($ML_y1+$ML_y2)))+`a'*`c'*exp(-`u')*((`eb'/(`u'+`eb'))^($ML_y1-($ML_y1+$ML_y2)))*((`es'/(`u'+`es'))^(-($ML_y1+$ML_y2)))+(1-`a')*((`eb'/(`u'+`eb'))^($ML_y1-($ML_y1+$ML_y2)))*((`es'/(`u'+`es'))^($ML_y2-($ML_y1+$ML_y2))))
end
ml model lf log10 (equation1:buy2010=) (equation2:sell2010=)
ml check
ml maximize


其中eb es a u c是计算PIN需要估计的5个参数，buy2010,sell2010都是已知的变量

我看别人的mle回归都只有$ml_y1，而我有两个被解释变量，是不是因为这一点出错啊，是不是我理解错了其实不是两个被解释变量...可是我已经找不到其他被解释变量了

大神快教我救命啊，重谢!!!!
数据在这里，真心请求救命，可以给论坛币~~
2010m0delete.xls (940 KB)

2017-2-12 23:19:26 上传

数据

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：stata编程 Stata 极大似然 tata EKOP stata 极大似然估计 MLE EKOP 编程

发现程序显示有点奇怪，再粘贴一次

capture program drop log10
program log10
version 1.0
args lnf eb es a u c
quietly replace `lnf' = -`eb'+($ML_y1+$ML_y2)*log(`eb'/(`u'+`eb'))+$ML_y1*log(`u'+`eb')-`es'+($ML_y1+$ML_y2)*log(`es'/(`u'+`es'))+$ML_y2*log(`u'+`es')+log(`a'*(1-`c')*exp(-`u')*((`eb'/(`u'+`eb'))^(-($ML_y1+$ML_y2)))*((`es'/(`u'+`es'))^($ML_y2-($ML_y1+$ML_y2)))+`a'*`c'*exp(-`u')*((`eb'/(`u'+`eb'))^($ML_y1-($ML_y1+$ML_y2)))*((`es'/(`u'+`es'))^(-($ML_y1+$ML_y2)))+(1-`a')*((`eb'/(`u'+`eb'))^($ML_y1-($ML_y1+$ML_y2)))*((`es'/(`u'+`es'))^($ML_y2-($ML_y1+$ML_y2))))
end
ml model lf log10 (equation1:buy2010=) (equation2:sell2010=)
ml search
ml check
ml maximize

还是有奇怪，这样子

1. capture program drop log10
   
2. program log10
   
3. version 1.0
   
4. args lnf eb es a u c
   
5. quietly replace `lnf' = -`eb'+($ML_y1+$ML_y2)*log(`eb'/(`u'+`eb'))+$ML_y1*log(`u'+`eb')-`es'+($ML_y1+$ML_y2)*log(`es'/(`u'+`es'))+$ML_y2*log(`u'+`es')+log(`a'*(1-`c')*exp(-`u')*((`eb'/(`u'+`eb'))^(-($ML_y1+$ML_y2)))*((`es'/(`u'+`es'))^($ML_y2-($ML_y1+$ML_y2)))+`a'*`c'*exp(-`u')*((`eb'/(`u'+`eb'))^($ML_y1-($ML_y1+$ML_y2)))*((`es'/(`u'+`es'))^(-($ML_y1+$ML_y2)))+(1-`a')*((`eb'/(`u'+`eb'))^($ML_y1-($ML_y1+$ML_y2)))*((`es'/(`u'+`es'))^($ML_y2-($ML_y1+$ML_y2))))
   
6. end
   
7. ml model lf log10 (equation1:buy2010=) (equation2:sell2010=)
   
8. ml search
   
9. ml check
   
10. ml maximize

复制代码

之前的程序都错了~算了~用截图吧

夏小雨1993 发表于 2017-2-12 23:21
请大神赐教ekop知情交易概率pin的eb es a u c5个参数的估计，用stata怎么写命令我自己写了一个命令，但是运 ...

你问题解决了吗？能把代码发给我看一下吗？

武天皓月 发表于 2017-5-12 12:03
你问题解决了吗？能把代码发给我看一下吗？

请问大佬们的问题是否解决！求联系

你好，请问问题解决了吗？