Stata中带有均值和CI的折线图

时间:2015-06-10 03:49:32

标签: confidence-interval stata

我有一个数据集,其中包含不同时间点(1个月,2个月,6个月,12个月和18个月)的效果估计值及其标准误差。现在我想绘制每个时期的均值和手段周围的相应CI。

我的示例如下:

effect  horizon   se
0.03    1         0.2
0.02    6         0.01
0.01    6         0.3
0.00    1         0.4
0.04    18        0.2
0.02    2         0.05
0.01    2         0.02
...     ...

每个地平线的效果均值导致5个数据点,我想在折线图中与置信区间一起绘制。我试过这个:

egen means = mean(effect), by(horizon)
line means horizon

但是如何添加对称置信带?这样我得到的东西看起来像这样: enter image description here

1 个答案:

答案 0 :(得分:1)

不完全确定这在统计上是有道理的,但是我可以这样做:

gen variance = se^2
collapse (mean) effect (sum) SV = variance (count) noobs = effect, by(horizon)
gen se_mean = sqrt(SV*(1/noobs)^2)
gen LB = effect - 1.96*se_mean
gen UB = effect + 1.96*se_mean

twoway (rline LB UB horizon, lpattern(dash dash)) (line effect horizon, lpattern(solid)), yline(0, lcolor(gray))   

哪个收益率:

enter image description here

要获得平均效果的SE T ̅,我使用公式

V T ̅)= 1 /( n 2 )Σ n V T i

(假设效果的协方差都为零)。然后我拿平方根得到 T SE的SE。