我做的是日度面板数据的多时点DID回归，因为只有工作日有数据所以日期不连续，我就生成了时间虚拟变量t，参考了黄河泉老师的命令代码，将year改成t，但是一直出现(0 real changes made)，因为大多数文献都是用年度数据的，抠破了头皮我也不知道自己错在了哪里，请路过的老师同学帮忙看看，感激不尽！！！[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)
请帮帮我这个可怜的少女吧！！！