如何使用BRR重复权重在SAS中创建五分位数

针对带有BRR重复权重的数据集，要为连续变量Var1生成加权五分位数，可以分两步操作：先通过PROC SURVEYMEANS计算加权分位点阈值，再将阈值匹配到原数据生成分组。

1. 计算加权五分位数的分割阈值 #

使用PROC SURVEYMEANS指定BRR权重和需要的分位点，输出分割阈值：

/* 计算Var1的加权20%、40%、60%、80%分位点（五分位数的四个分割点） */
PROC SURVEYMEANS DATA=your_dataset NOPRINT;
    VAR Var1;
    WEIGHT base_weight; /* 替换为你的基础权重变量名 */
    REPWEIGHTS brr_w1-brr_w500; /* 替换为500个BRR重复权重的变量范围，比如brr_weight1到brr_weight500 */
    QUANTILE=P20 P40 P60 P80;
    OUTPUT OUT=quintile_thresholds 
           P20=cut_p20 
           P40=cut_p40 
           P60=cut_p60 
           P80=cut_p80;
RUN;

• NOPRINT：避免输出大量重复权重的冗余统计结果，只保留分位点输出。
• REPWEIGHTS：指定所有BRR重复权重变量，让程序考虑复杂抽样的加权结构。
• QUANTILE参数：指定五分位数的四个分割点（20%、40%、60%、80%分位点），将数据划分为五等份。
• OUTPUT：把计算得到的阈值存入quintile_thresholds数据集，每个阈值对应一个变量。

2. 为原数据分配五分位数分组 #

将第一步得到的阈值匹配到原数据，通过条件判断生成五分位数分组变量：

/* 给原数据的Var1添加加权五分位数分组 */
DATA dataset_with_quintile;
    MERGE your_dataset quintile_thresholds;
    BY _TYPE_; /* PROC SURVEYMEANS输出的_TYPE_=0对应总体统计，确保每条数据都匹配到阈值 */
    
    /* 处理缺失值：如果Var1缺失，分组设为缺失 */
    IF MISSING(Var1) THEN Var1_quintile = .;
    ELSE IF Var1 <= cut_p20 THEN Var1_quintile = 1;
    ELSE IF Var1 <= cut_p40 THEN Var1_quintile = 2;
    ELSE IF Var1 <= cut_p60 THEN Var1_quintile = 3;
    ELSE IF Var1 <= cut_p80 THEN Var1_quintile = 4;
    ELSE Var1_quintile = 5;
    
    LABEL Var1_quintile = "Var1的加权五分位数分组（1=最低组，5=最高组）";
RUN;

• MERGE：将原数据和阈值数据集合并，由于阈值数据集只有一条总体记录，BY _TYPE_可以保证所有原数据都能匹配到阈值。
• 条件判断：根据Var1的值落在哪个区间，分配1到5的分组标签，同时处理缺失值避免逻辑错误。

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