g _b_random_gtp= _b[random_gtp]  //提取x的回归系数
        g _se_random_gtp= _se[random_gtp] //提取x的标准误
                gen tvalue=_b_random_gtp/_se_random_gtp
                gen pvalue = 2*ttail(e(df_r), abs(tvalue))