咨询可计算指定概率对应分位值的SAS函数

嘿，这个需求在SAS里完全能轻松搞定，我给你两种常用的实现方法，你可以根据自己的使用场景来选：

如果你想在DATA步里直接算出这个90%分位数（也就是你要的“90%数据低于它，10%高于它”的数值），可以用SAS的QUANTILE函数，它支持基于经验分布的分位数计算，完美适配你的样本数据场景。

举个例子，假设你的样本数据集叫sample，要分析的数值变量是value，可以这么写代码：

data _null_;
    set sample end=eof;
    /* 定义临时数组存储所有样本值，大小根据你的实际样本量调整 */
    array temp_vals[1000] _temporary_;
    temp_vals[_n_] = value;
    
    /* 当读取完所有数据后计算分位数 */
    if eof then do;
        /* 'EMPIRICAL'指定用经验分布，0.9对应90%分位数 */
        q90_value = quantile('EMPIRICAL', 0.9, of temp_vals[*]);
        put "90%数据低于的数值为：" q90_value;
    end;
run;

运行后，日志里会直接输出你要的那个数值。

如果你需要同时查看更多统计信息，或者想把结果保存到数据集里，PROC UNIVARIATE会更方便，它是SAS里做单变量统计分析的常用过程，分位数计算是它的基础功能之一。

代码示例：

/* 计算分位数并输出到新数据集 */
proc univariate data=sample;
    var value;
    /* 指定输出90%分位数，P_是结果变量的前缀，最终变量名会是P_90 */
    output out=quantile_results pctlpts=90 pctlpre=P_;
run;

/* 查看结果 */
proc print data=quantile_results;
run;

运行后，quantile_results数据集里的P_90变量就是你要的数值。

小提示：SAS的分位数计算有几种不同的定义（比如插值方式），如果需要调整，可以在PROC UNIVARIATE里加PCTLDEF=选项（比如pctldef=5是默认值，pctldef=2是另一种常用的定义），不过一般默认设置就满足大多数场景的需求啦。

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