Stata中centile命令百分位计算结果批量存储方法问询

下面提供两种比逐一声明局部宏更简洁的实现方式，适配不同的后续使用需求： #

方案1：存储为单个局部宏向量 #

适合需要快速调用百分位值做截断、分组等场景，不需要额外依赖矩阵操作：

* 计算50个百分位（0(2)100刚好生成51个分位点，对应50个百分位区间，可根据需求调整步长）
centile var, centile(0(2)100)
local pct_vector ""
* 直接调用返回值r(n_cent)获取分位点总数，无需手动写死循环次数
forval i = 1/`=r(n_cent)' {
    local pct_vector `pct_vector' `=r(c_`i')'
}

* 调用示例：
* 输出所有分位点
display "所有百分位结果：`pct_vector'"
* 提取第k个分位点（比如第10个）
local p10 `:word 10 of `pct_vector''
display "第10个百分位值：`p10'"

方案2：存储为矩阵 #

适合后续需要做数值计算、批量导出等场景，调用更规范：

centile var, centile(0(2)100)
* 初始化和分位点数量等长的单列矩阵
mat pct_matrix = J(`=r(n_cent)', 1, .)
forval i = 1/`=r(n_cent)' {
    mat pct_matrix[`i', 1] = r(c_`i')
}

* 调用示例：
* 查看完整矩阵内容
mat list pct_matrix
* 提取第k个分位点（比如第20个）
local p20 = el(pct_matrix, 20, 1)
display "第20个百分位值：`p20'"

两种方案的共同优势 #

• 无需手动指定循环上限，修改centile的分位点参数时不需要同步调整循环代码，适配性更强
• 统一存储为单个对象，不会生成大量零散的局部宏，运行环境更整洁
• 拓展性更强，宏向量可以直接传入分组命令的cut参数，矩阵可以直接用svmat命令转为变量保存到数据集里

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