求助安慰剂检验是否通过

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

求助大家，在安慰剂检验中随机改变政策发生时间。主要做法是将样本数据按照省份分组，然后在每个省份内的时间变量中随机抽取一个年份作为其政策时间。同样将上述过程重复500次后，得到这个图。它不是以0为轴线进行分布，想问下是否通过检验，谢谢大家了。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：安慰剂 样本数据 求助大家 时间变量 样本数

你这明显有问题，检查下代码

917968079 发表于 2023-4-1 12:16
你这明显有问题，检查下代码

你好，可以帮我看下代码么，救救我吧，感谢
mat b = J(500,1,0)
mat se = J(500,1,0)
mat p = J(500,1,0)

forvalues i = 1/500{
use C:\Users\Yu\Desktop\2.20数据.dta, clear
xtset id year
bsample 1, strata(id)
keep year
save matchyear.dta, replace
mkmat year, matrix(sampleyear)
use C:\Users\Yu\Desktop\2.20数据.dta, clear
xtset id year
gen DID = 0
foreach j of numlist 1/18 {
  replace DID = 1 if (id == `j' & year >= sampleyear[`j',1])
}
qui xtreg y DID x1-x5, 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]))
}

*绘图
svmat b, names(coef)
svmat se, names(se)
svmat p, names(pvalue)

drop if pvalue1 == .
label var pvalue1 p值
label var coef1 估计系数

# delimit ;
twoway (scatter pvalue1 coef1, xline(0.2142 , lwidth(0.2) lp(shortdash))
       xlabel(-0.4(0.1)1.0, grid) xtitle(估计系数) ytitle(p值)
           msymbol(smcircle_hollow)  legend(off))
           (kdensity coef1, title(安慰剂检验))
           ;
# delimit cr

我今天也在搞安慰剂检验，要是通过不了的话是不是之前做的都白费了