在SAS中实现分位数回归：获取响应变量80分位数及预测区间

嘿，我完全懂你的需求——常规线性回归聚焦于响应变量的均值及其置信区间，但你想要的是特定分位数（比如你提到的80分位数）对应的响应值，以及这个分位数的预测区间（也就是给定预测变量值时，有指定比例的响应值低于该数值）。在SAS里，PROC QUANTREG就是专门解决这类问题的工具，下面给你详细的实现方法：

1. 拟合80分位数回归模型 #

首先用PROC QUANTREG针对80分位数构建回归模型，假设你的数据集名为my_data，响应变量是Response，预测变量是Predictor，代码如下：

proc quantreg data=my_data;
    model Response = Predictor / quantile=0.8;
run;

这个代码会输出80分位数回归的参数估计，逻辑和常规回归类似，但所有结果都是针对响应变量的80分位数，而非均值。

2. 获取指定预测点的分位数预测值与区间 #

如果要得到Predictor=3时的80分位数响应值，以及对应的预测区间，需要结合预测点数据集和OUTPUT语句来实现：

步骤1：创建包含目标预测点的数据集 #

先生成一个只包含你要预测的Predictor值的小数据集：

data pred_points;
    input Predictor;
    datalines;
3
;
run;

步骤2：执行预测并输出结果 #

使用PROC QUANTREG的PREDICT选项调用预测点数据集，同时用OUTPUT语句保存预测值和区间：

proc quantreg data=my_data;
    model Response = Predictor / quantile=0.8;
    predict data=pred_points out=pred_summary;
    /* 保存分位数预测值和其置信区间 */
    output out=interval_details p=pred_q80 lcl=lcl_q80 ucl=ucl_q80;
    /* 如果需要针对新观测的预测区间，用下面的语句 */
    /* output out=prediction_range p=pred_q80 lclp=lcl_pred_q80 uclp=ucl_pred_q80; */
run;

• pred_q80：对应Predictor=3时，80分位数的响应值（也就是你要的“80%的响应值会低于该数值”的目标值）
• lcl_q80/ucl_q80：该分位数估计值的95%置信区间（可通过ALPHA=参数调整置信水平，比如ALPHA=0.2就是80%置信区间）
• 如果你需要的是针对单个新观测的预测区间（类似你例子里的10到16），可以用注释里的lclp=和uclp=参数，输出的就是新观测80分位数的预测区间。

3. 查看结果 #

运行完代码后，查看pred_summary或interval_details数据集，就能直接获取你需要的所有结果：Predictor=3对应的80分位数响应值，以及该值的区间范围。

补充小技巧 #

• 若需要同时拟合多个分位数（比如50分位、80分位），可以在QUANTILE=参数里指定多个值，比如quantile=0.5 0.8
• PROC QUANTREG支持更多复杂模型（如带交互项、分类变量），只需在model语句里按常规回归的写法添加即可。

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