各位大神,鄙人最近在研究广义倾向得分匹配,使用doseresponse 命令,
qui generate cut = 0.0103 if exrate0_w<=0.0103
qui replace cut = 0.0345 if exrate0_w>0.0103& exrate0_w<= 0.0345
qui replace cut = 0.0993 if exrate0_w> 0.0993
matrix define tp = (10\20\30\40\50\60\70\80\90\100)
doseresponse gq lnzj_w ldbl_w dshsize_w ddrs_w top10_w, outcome( v47) t( exrate0_w) gpscore(pscore) predict(hat_treat) sigma(hat_sd) cutpoints(cut) index(p50) nq_gps(5) t_transf(bcskew0) dose_response(dose_response) tpoints(tp) reg_type_t(quadratic) reg_type_gps(quadratic) interaction(1) bootstrap(yes) boot_reps(5) filename("output") analysis(yes) graph("graph_output") detail
前两步的回归结果是可以正常显示的,但是图显示不出来,提示
Bootstrapping of the standard errors
...variable treatment_level_plus not found
想问下各位大神,是不是tp设置的问题,因为exrate0_w较小,这边的tp(10\20\30\40\50\60\70\80\90\100)看案例应该是百分位的意思;还有按cut分类后,回归出现四个区间是怎么回事?
再或者是数据本身的问题?因为产出变量相对比较大,而处理变量特别小。
再或者进行Bootstrap本身有什么要求没有满足?
跪求在这方面有研究的大神帮帮我!