做安慰剂检验,最后总是 _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);
大神能帮我看看是为啥吗?救救孩子


雷达卡


京公网安备 11010802022788号







