[loveliness][loveliness]我这共有105个ID,30个实验组,743个日度数据,命令代码如下:
*========安慰剂检验=======
*生成备用矩阵
mat b = J(500,1,0)
mat se = J(500,1,0)
mat p = J(500,1,0)
*抽样过程
forvalues i = 1/500{
use D:\11月论文\安慰剂检验专用.dta, clear
xtset ID t
bsample 1, strata(ID) //根据**ID**分组,每组随机抽取一个年份
keep t
save matchdate.dta, replace
mkmat t, matrix(sampledate)
use D:\11月论文\安慰剂检验专用.dta, clear
xtset ID t
gen DID = 0
foreach j of numlist 1/743 {
replace DID = 1 if (ID == `j' & t >= sampledate[`j',1])
}
qui xtscc R DID duration turn lCap markettype t,fe
mat b[`i',1] = _b[DID]
mat se[`i',1] = _se[DID]
scalar df_r = e(N) - e(df_m) -1
mat p[`i',1] = 2*ttail(df_r,abs(_b[DID]/_se[DID]))
}
跑出来结果却是一堆的(0 real changes made)
请帮帮我这个可怜的少女吧!!!


雷达卡




京公网安备 11010802022788号







