forvalues i=1/500{
        use "C:\Users\ASUS\Desktop\研究生模型所有数据\数据.dta", clear
        xtset id year  //面板数据声明
        keep if year==2012   //保留一期数据
        sample 29, count   //随机抽取32个城市
        keep id  //得到所抽取样本的id编号
        save match_id.dta, replace   //另存id编号数据
        merge 1:m id using "C:\Users\ASUS\Desktop\研究生模型所有数据\数据.dta" //与原数据匹配
        gen treat1 = (_merge == 3) //将所抽取样本赋值为1，其余为0，得到政策分组虚拟变量
        gen period = (year >= 2018) //生成政策时间虚拟变量
        gen did2 = treat1*period
        xtreg LnPL did2 $xlisti. i.id i.year

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

        * 计算P值并赋值于矩阵
        mat p[`i',1] = 2*ttail(e(df_r), abs(_b[did2]/_se[did2]))
}
之前做的时候可以，隔一段时间再去跑就显示. invalid name