Stata的双向曲线与置信区间

Question

使用

clear 
score  group  test 
2   0   A 
3   0   B 
6   0   B 
8   0   A 
2   0   A 
2   0   A 
10   1   B 
7   1   B 
8   1   A 
5   1   A 
10   1   A 
11   1   B 
end 

我想散点图通过group每个test（同图）置信区间（实际数据有数以千计的意见）是指score。结果图将有两组两个点。 test==a（group==0 vs group==1）的一组点和test==b（group==0 vs group==1）的一组点。

我目前的做法很有效，但很费力。我使用egen来计算所有需要的统计数据：通过测试得到的每个组的平均值，观察次数，标准偏差...。我然后collapse的数据和情节。

还有另一种方式，不是吗？

我假设是塔塔将能够采取作为其输入scoregroup和test变量，然后计算和现在这个漂亮的标准曲线。在Google上花了很多时间后，我不得不问。

Answer 1

虽然有用户编写的程序，但我倾向于statsby作为这里的基本方法。讨论可访问in this paper。

本示例以您的数据示例（几乎可执行的代码）为例。一些选择取决于隐含的大置信区间。请注意，如果您的Stata版本不是最新的，则ci的语法将有所不同。 （只是省略means）

clear 
input score  group str1 test 
2   0   A 
3   0   B 
6   0   B 
8   0   A 
2   0   A 
2   0   A 
10   1   B 
7   1   B 
8   1   A 
5   1   A 
10   1   A 
11   1   B 
end 
save cj12 , replace 

* test A 
statsby mean=r(mean) ub=r(ub) lb=r(lb) N=r(N), by(group) clear : /// 
ci means score if test == "A" 
gen test = "A" 
save cj12results, replace 

* test B 
use cj12 
statsby mean=r(mean) ub=r(ub) lb=r(lb) N=r(N), by(group) clear : /// 
ci means score if test == "B" 
gen test = "B" 
append using cj12results 

* graph; show sample sizes too, but where to show them is empirical 
set scheme s1color 
gen where = -20 
scatter mean group, ms(O) mcolor(blue) || /// 
rcap ub lb group, lcolor(blue) /// 
by(test, note("95% confidence intervals") legend(off)) /// 
subtitle(, fcolor(ltblue*0.2)) /// 
ytitle(score) xla(0 1) xsc(r(-0.25 1.25)) yla(-10(10)10, ang(h)) || /// 
scatter where group, ms(none) mla(N) mlabpos(12) mlabsize(*1.5) 

我们不能比较的完整代码或您的图形，因为你既没有显示。