求助，多期DID安慰剂检验如图，能帮我看下是什么问题么，感谢

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

代码如下：
search permute  //先安装此命令

permute did beta = _b[did] se = _se[did] df = e(df_r), seed(10000) ///
reps(500) saving("安慰剂检验.dta"): ///
    reghdfe y did $c ,a(id year)

use "安慰剂检验.dta" , clear
gen t_value = beta / se
gen p_value = 2 * ttail(df, abs(beta/se))
replace beta = beta

#delimit ;
                twoway (kdensity beta, yaxis(1) lc(black) ) (scatter p_value beta, yaxis(2) m(Oh) mc(black))  ,
                 xline(0, lc(black*0.5) lp(solid))
                 yline(0.05, lc(black*0.5) lp(dash) axis(2))
                 xtitle("Estimator", size(*0.8)) xlabel(, format(%4.1f) labsize(small))
                 ytitle("Density", size(*0.8)) ylabel(, nogrid format(%4.1f) labsize(small))
                 ytitle("P Value", size(*0.8) axis(2)) ///
                 ylabel(, nogrid format(%4.1f) labsize(small) axis(2))
                 legend(ring(2) order(1 "Estimator" 2 "P Value") pos(6) rows(1) )
                 graphregion(color(white)) scheme(s1mono) ;
#delimit cr

graph export "安检验.png", width(1000) replace

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：DID 安慰剂 estimator kdensity Density

看着没什么问题，不过需要再加一条垂直线标明你的真实估计值

哥哥海哥哥 发表于 2023-1-2 14:29
看着没什么问题，不过需要再加一条垂直线标明你的真实估计值

您好，谢谢您的回复，可以教我下加垂直线的代码么，刚入手，还不熟悉，万分感谢

用xline()这个选项就可以了，你的源代码里也有

哥哥海哥哥 发表于 2023-1-3 14:05
用xline()这个选项就可以了，你的源代码里也有

你好，我的回归系数是0.2142，我用xline(0.2142),出现not allowed.可以帮我写下这句代码么，我自己尝试了都没有出现结果，真的万分感谢。

把代码贴出来看一下吧，这样说我也不知道问题出在哪里

哥哥海哥哥 发表于 2023-1-4 13:44
把代码贴出来看一下吧，这样说我也不知道问题出在哪里

您好，我的代码如下，麻烦看下是什么问题，感谢万分。
search permute  //先安装此命令

permute did beta = _b[did] se = _se[did] df = e(df_r), seed(10000) ///
reps(500) saving("安1检验.dta"): ///
    reghdfe y did $c ,a(id year)

use "安1检验.dta" , clear
gen t_value = beta / se
gen p_value = 2 * ttail(df, abs(beta/se))
replace beta = beta

#delimit ;
                twoway (kdensity beta, yaxis(1) lc(black) ) (scatter p_value beta, yaxis(2) m(Oh) mc(black))  ,
                xline（0.2142）
                xline(0, lc(black*0.5) lp(solid))
                 yline(0.05, lc(black*0.5) lp(dash) axis(2))
                 xtitle("Estimator", size(*0.8)) xlabel(, format(%4.1f) labsize(small))
                 ytitle("Density", size(*0.8)) ylabel(, nogrid format(%4.1f) labsize(small))
                 ytitle("P Value", size(*0.8) axis(2)) ///
                 ylabel(, nogrid format(%4.1f) labsize(small) axis(2))
                 legend(ring(2) order(1 "Estimator" 2 "P Value") pos(6) rows(1) )
                 graphregion(color(white)) scheme(s1mono) ;
#delimit cr

graph export "安1检验.png", width(1000) replace

xline（0.2142）不要用中文括号

哥哥海哥哥 发表于 2023-1-5 09:32
xline（0.2142）不要用中文括号

您好，我调整了英文括号后，怎么还是没有出来垂直线，结果如下，打扰您看下是哪里的问题，很感谢帮我解答。


search permute  //先安装此命令

permute did beta = _b[did] se = _se[did] df = e(df_r), seed(10000) ///
reps(500) saving("安4检验.dta"): ///
    reghdfe y did $c ,a(id year)

use "安4检验.dta" , clear
gen t_value = beta / se
gen p_value = 2 * ttail(df, abs(beta/se))
replace beta = beta
#delimit ;
                twoway (kdensity beta, yaxis(1) lc(black) ) (scatter p_value beta, yaxis(2) m(Oh) mc(black))  ,
                xline(0.2142)
                xline(0, lc(black*0.5) lp(solid))
                 yline(0.05, lc(black*0.5) lp(dash) axis(2))
                 xtitle("Estimator", size(*0.8)) xlabel(, format(%4.1f) labsize(small))
                 ytitle("Density", size(*0.8)) ylabel(, nogrid format(%4.1f) labsize(small))
                 ytitle("P Value", size(*0.8) axis(2)) ///
                 ylabel(, nogrid format(%4.1f) labsize(small) axis(2))
                 legend(ring(2) order(1 "Estimator" 2 "P Value") pos(6) rows(1) )
                 graphregion(color(white)) scheme(s1mono) ;
#delimit cr

graph export "安1检验.png", width(1000) replace

您好，我调整了英文括号后，结果还是没有出现垂直线，麻烦您再帮我看下，很感谢您一直帮我解答。（代码与安慰剂图如下，感谢 ）

search permute  //先安装此命令

permute did beta = _b[did] se = _se[did] df = e(df_r), seed(10000) ///
reps(500) saving("安4检验.dta"): ///
    reghdfe y did $c ,a(id year)

use "安4检验.dta" , clear
gen t_value = beta / se
gen p_value = 2 * ttail(df, abs(beta/se))
replace beta = beta

#delimit ;
                twoway (kdensity beta, yaxis(1) lc(black) ) (scatter p_value beta, yaxis(2) m(Oh) mc(black))  ,
                xline(0.2142)
                xline(0, lc(black*0.5) lp(solid))
                 yline(0.05, lc(black*0.5) lp(dash) axis(2))
                 xtitle("Estimator", size(*0.8)) xlabel(, format(%4.1f) labsize(small))
                 ytitle("Density", size(*0.8)) ylabel(, nogrid format(%4.1f) labsize(small))
                 ytitle("P Value", size(*0.8) axis(2)) ///
                 ylabel(, nogrid format(%4.1f) labsize(small) axis(2))
                 legend(ring(2) order(1 "Estimator" 2 "P Value") pos(6) rows(1) )
                 graphregion(color(white)) scheme(s1mono) ;
#delimit cr

graph export "安4检验.png", width(1000) replace