首先用stata跑MLE其实不是很方便。建议你用OLS估计系数,反正在iid正态分布的假设下二者区别就是自由度调整。那么假设你数据里面eps叫x, 然后前一年的eps叫x_l1y, 公司id叫cusip, 年度时间变量叫year,那么code这么写:
reg x x_l1y if cusip == "XXX" & year > YYY-10 & year <=YYY
就会做你上面的的估值。然后在stata里面每做一个这样的regression后台都会存储相关宏观变量,具体有哪些你可以打h reg看一下。你的那个predictability应该存在e(rmse) 里面, 当然这个还是有可能有自由度调整带来的一点点的区别。你跑完我上面写的一行代码写 di e(rmse)就可以看到predictability了。