在Stata twoway传奇中结合标记符号和错误栏

时间:2016-07-26 06:09:30

标签: stata

我经常使用-twoway scatter-and -twoway rcap-绘制使用95%误差条绘制平均值或系数的图表。下面的代码生成一个带有两个条目的图例:一个用于平均标记符号,另一个用于错误栏。但是我想让图例显示一个条目,显示标记符号和错误条的组合。下面是我通常如何制作图表的示例。

sysuse auto

gen b = .
gen se = .

mean mpg if foreign == 1
replace b = _b[mpg] in 1
replace se = _se[mpg] in 1

mean mpg if foreign == 0
replace b = _b[mpg] in 2
replace se = _se[mpg] in 2

gen lb = b - (1.96 * se)
gen ub = b + (1.96 * se)
gen index = _n in 1/2

twoway scatter b index || rcap lb ub index, legend(order(1 "Mean" 2 "95% Interval"))

-legend-中是否有选项允许我以我想要的方式覆盖两个图例条目?

1 个答案:

答案 0 :(得分:2)

我真的不知道如何做到你想要的。这似乎很难。

我也讨厌浪费传说中的房地产,因此另一种方法是标记手段而不是使用图例(并添加" 95%CI和#34;标题):

sysuse auto
reg mpg i.foreign
margins foreign, post
estimates store means

marginsplot, recast(scatter) xscale(reverse)
coefplot means

另一种方法是在没有任何回归/摘要的情况下使用ciplot

ciplot mpg, by(foreign) xscale(reverse)

coefplotciplot都是用户编写的。