在Stata双向图例中结合标记符号和错误栏

Question

我经常使用 -twoway scatter- 和 -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- 中是否有允许我以我想要的方式覆盖两个图例条目的选项？

Answer 1

我真的不知道如何做你想要的。好像很难。

我也讨厌浪费图例空间，因此一种替代方法是标记手段而不是使用图例（并在标题中添加“With 95%CIs”）：

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)

coefplot 和 ciplot 都是用户编写的。