*===============================================
*                                        会计稳健性计算
*===============================================
import excel 数据.xlsx, sheet("DA计算") firstrow clear


gen TA_A=应计利润总额/前一年期末资产总计
gen _1_A=1/前一年期末资产总计
gen REV_AR_A=((主营营业收入-前一年的主营营业收入)-(应收账款净额-前一年应收账款净额))/前一年期末资产总计
gen PPE_A=固定资产原值/前一年期末资产总计
gen ROA=前一年总资产收益率

//分行业分年度回归

egen ind=group(行业)
sum ind
local Ns=r(max)

gen DA=.
forvalues y=2014/2016{
   forvalues d=1/`Ns'{
      qui reg TA_A _1_A REV_AR_A PPE_A ROA if (year==`y' & ind==`d'), nocons
          qui predict temp if e(sample), resid
          qui replace DA=temp if e(sample)
          drop temp
   }
}