从 centile 命令保存百分位数

如何解决从 centile 命令保存百分位数

我在 Stata 中使用命令 centile。我想将结果存储在本地矩阵或向量中。

我知道我可以使用 r(c_#) 获得单个百分位数的结果。但我正在计算 50 个百分位数，并希望将它们全部存储在一个本地向量/矩阵中，而不必将它们全部写出。

更新： 我找到了以下解决方案，但我仍然希望得到更优雅的答案。

centile var,centile(0(5)100)
forval i=1/21{
    local cut_`i'= r(c_`i')
}
display `cut_1' `cut_2'

解决方法

这里有一些小技巧。请注意，矩阵或向量不被视为“本地”，尽管其名称可以保存在本地宏中。

sysuse auto,clear
centile mpg,centile(0(5)100)
matrix results = J(21,1,.)

forval i = 1/21 {
    matrix results[`i',1] = r(c_`i')
}