Stata绘制分类带可信区间的折线图
论文中折线图可以表示两个变量间的变化关系，带误差或可信区间的折线图表示其中一个变量的变化范围，既往我们已经使用R语言绘制了分类带误差和可信区间的折线图，今天我们使用Stata来绘制分类带误差和可信区间的折线图。
继续使用我们的汽车销售数据（公众号回复：汽车销售，可以获得该数据）来演示，先导入数据，我是直接黏贴
我们来看下数据，car就是汽车售价，age是年龄，gender是性别，inccat是收入，这里分成4个等级，ed是教育程度。
假设想知道不同教育水平中不同性别人群购车的平均价，我们先得对数据进行整理
collapse (mean) car= car age=age,by( ed gender )
Collapse函数一次可以分组两个变量，但是不能出两个结果，如果你还想知道标准误的结果，只能再做一次，然后合并，这一点就比不上R了。不过也就是多做一次，也不麻烦。
collapse (semean ) se= car ,by( ed gender )
合并数据，其实很简单直接粘贴在Excel合并就好了。
把黏贴好的数据重新导入stata
因为我们绘制的是带可信区间的折线图，所以要先生成可信区间
g up= car+1.96* se
g lp= car-1.96* se
生成可信区间就可以绘图了，我们先画女性f的，这里要用lcolor控制颜色，不然后期颜色不一样
twoway(rcap up lp ed if gender=="f",lcolor(blue))
然后画连接线，这里我们使用mlcolor(blue)lcolor(blue)控制颜色，控制的是先得颜色和便面环的颜色，填充的颜色还没有改，所以有点怪怪的
twoway(rcap up lp ed if gender=="f",lcolor(blue))(connected car ed if gender=="f",mlcolor(blue)lcolor(blue))
使用mfcolor控制填充颜色
twoway(rcap up lp ed if gender=="f",lcolor(blue))(connected car ed if gender=="f",mlcolor(blue)lcolor(blue)mfcolor(blue))
这样一条线就画好了，我们按照刚才的方法画另一条，非常简单把
twoway(rcap up lp ed if gender=="f",lcolor(blue))(connected car ed if gender=="f",mlcolor(blue)lcolor(blue)mfcolor(blue)) (rcap up lp ed if gender=="m",lcolor(red))(connected car ed if gender=="m",mlcolor(red)lcolor(red)mfcolor(red))
这个图例不好看，我们关掉它，然后添加上标题
twoway(rcap up lp ed if gender=="f",lcolor(blue))(connected car ed if gender=="f",mlcolor(blue)lcolor(blue)mfcolor(blue)) (rcap up lp ed if gender=="m",lcolor(red))(connected car ed if gender=="m",mlcolor(red)lcolor(red)mfcolor(red)),legend(off) xtitle("教育程度") ytitle("汽车售价")
最后在修整一下就弄好了
如果觉得有用就分享给有需要的人吧。