做安慰剂检验，最后总是   _b not found r(111)
use "D:\LHY\毕业论文\2-行业数据\全部\安慰剂检验\安慰剂检验\是否捐赠\是否捐赠原数据.dta",clear
xtset code
xtset Date
xtset Date_Donation
gen did=(Date>= Date_Donation)
global xlist " mreturn trading_volume C_Dead "
reghdfe CAR did $xlist ,absorb(Date )
mat b = J(500,1,0) //* 系数矩阵
mat se = J(500,1,0) //* 标准误矩阵
mat p = J(500,1,0) //* P值矩阵



forvalues i=1/500{
        use "D:\LHY\毕业论文\2-行业数据\全部\安慰剂检验\安慰剂检验\是否捐赠\是否捐赠1.dta", clear
        duplicates drop
        xtset Date //面板数据声明
        keep if Date ==201994   //保留一期数据
        sample 53, count   //随机抽取32个城市
        keep stockid   //得到所抽取样本的id编号
        save match_id.dta, replace   //另存id编号数据
        merge 1:m stockid using "D:\LHY\毕业论文\2-行业数据\全部\安慰剂检验\安慰剂检验\是否捐赠\是否捐赠原数据.dta" //与原数据匹配
        gen treat = (_merge == 3) //将所抽取样本赋值为1，其余为0，得到政策分组虚拟变量
        gen period = ( Date >= Date_Donation ) //生成政策时间虚拟变量
        gen did1 = treat*period
        reghdfe CAR did1 $xlist ,absorb( code Date ) vce(cluster Province_code )

        * 将回归结果赋值到对应矩阵的对应位置
        mat b[`i',1] = _b [did1] //* 系数矩阵
        mat se[`i',1] = _se[did1] //* 标准误矩阵

        * 计算P值并赋值于矩阵
        mat p[`i',1] = 2*ttail(e(df_r), abs(_b[did1]/_se[did1]))
}


循环之前总归是出现_b not found   r(111);
大神能帮我看看是为啥吗？救救孩子