搞了好几天，终于想明白了。
修正的琼斯模型是这样的：

TA=NI-CFO

TA/A_1=a0(1/A_1)+a1(Δsale-ΔAR)/A_1+a2*FA/A_1+e

NDA=a0(1/A_1)+a1(Δsale-ΔAR)/A_1+a2*FA/A_1

DA=TA/A_1-NDA

公式编辑得不好，但是应该可以看懂吧。A_1是上一期资产

本来想的是先分行业和年份回归，提取系数a0,a1,a2，再带入第二个公式回归。后来发现！！！残差项e就是DA值！所以不用那么麻烦了！

第一步计算TA，这个不用写了。

第二步，按年份和行业回归：（之前要先把年份和行业都编号）

sortyear indus

bys year indus: regata_wa_w asalear_wafa_w *我先把变量进行了winsorize处理*

*生成各回归的预测值yp与残差e：

predict ata_wp
predict e,r

最后这个e就是DA，这样就完了。

如果大家认为我做错了，希望指出来，大神们都来帮我看看，谢谢啦！