如何在Stata中绘制单值变量直方图并叠加背景分布？

核心解决思路 #

当变量仅取单一值时，Stata默认分箱逻辑会导致直方图显示异常，需手动指定分箱参数，再通过twoway命令实现两个直方图的叠加。

完整实现代码 #

假设参考分布为0到1区间的正态分布（可替换为你的实际数据），代码如下：

* 生成示例数据
gen example_var = 0.85
gen ref_dist = rnormal(0.5, 0.2)  // 生成参考分布的模拟数据

* 叠加两个直方图
twoway (hist ref_dist, percent xscale(range(0 1)) graphregion(color(white)) bin(20) lcolor(gray) fcolor(gray%30)) ///
       (hist example_var, percent bin(1) width(0.02) lcolor(black) fcolor(black%60)), ///
       xtitle("Example Data") legend(order(1 "参考分布" 2 "单一值数据"))

关键参数说明 #

• 针对单一值变量example_var：bin(1)强制生成1个分箱，width(0.02)设置窄分箱宽度，避免出现占满整个x轴的怪异柱子。
• 针对参考分布ref_dist：bin(20)调整分箱数量适配分布形态，fcolor(gray%30)设置半透明填充，防止遮挡前景的单一值柱子。
• twoway的括号结构分别定义两个直方图图层，实现背景与前景的叠加效果。

替代方案（基于已有统计结果） #

若参考分布是已有的统计结果而非原始数据，可先绘制参考分布直方图，再用addplot添加单一值的柱子：

* 绘制参考分布直方图并叠加单一值柱子
hist ref_dist, percent xscale(range(0 1)) graphregion(color(white)) bin(20) lcolor(gray) fcolor(gray%30) ///
       addplot(bar 100 0.85, barwidth(0.02) lcolor(black) fcolor(black%60)) ///
       xtitle("Example Data") legend(order(1 "参考分布" 2 "单一值数据"))

其中bar 100 0.85表示在x=0.85位置绘制高度为100%的柱子（因单一值占比为100%）。

内容的提问来源于stack exchange，提问作者garzi