如何为含截断长尾数据的直方图坐标轴添加加号(+)

右侧长尾截断场景（原需求） #

针对你提供的价格数据右侧截断到10000的情况，只需在绘制直方图时通过xlabel()自定义X轴刻度标签，将最后一个刻度的文本改为10000+即可：

完整代码 #

sysuse auto, clear
gen price_graph = price
replace price_graph = 10000 if price_graph > 10000

// 自定义X轴刻度与标签，最后一个刻度添加加号
hist price_graph, start(0) width(1000) ///
    xtick(0(1000)10000) /// 保留原刻度位置
    xlabel(0 "0" 1000 "1000" 2000 "2000" 3000 "3000" 4000 "4000" ///
           5000 "5000" 6000 "6000" 7000 "7000" 8000 "8000" ///
           9000 "9000" 10000 "10000+")

灵活适配不同截断值 #

如果需要更换截断值，可通过循环自动生成标签，避免手动逐个输入：

sysuse auto, clear
local cut_right = 10000 // 自定义右侧截断值
gen price_graph = price
replace price_graph = `cut_right' if price_graph > `cut_right'

// 自动生成X轴标签
local labels ""
forvalues i=0(1000)`cut_right' {
    if `i' == `cut_right' {
        local labels "`labels' `i' "`i'+'""
    }
    else {
        local labels "`labels' `i' "`i'""
    }
}

hist price_graph, start(0) width(1000) xtick(0(1000)`cut_right') xlabel(`labels')

左侧截断场景适配 #

如果是左侧长尾（小值端）截断，只需修改左侧第一个刻度的标签，用<截断值>标识该组包含所有小于截断值的数据：

完整代码 #

sysuse auto, clear
local cut_left = 1000 // 自定义左侧截断值
gen price_graph_left = price
replace price_graph_left = `cut_left' if price_graph_left < `cut_left'

// 自动生成X轴标签
local max_val = ceil(max(price_graph_left)/1000)*1000 // 取最大刻度为1000的倍数
local labels ""
forvalues i=`cut_left'(1000)`max_val' {
    if `i' == `cut_left' {
        local labels "`labels' `i' "<`i'""
    }
    else {
        local labels "`labels' `i' "`i'""
    }
}

hist price_graph_left, start(`cut_left') width(1000) xtick(`cut_left'(1000)`max_val') xlabel(`labels')

说明 #

核心逻辑是通过xlabel()覆盖默认刻度文本，针对截断边界的刻度添加标识，明确该组包含超出截断范围的所有数据，让直方图的信息展示更准确。

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