如何在SAS PROC SQL中匹配最近前置事件日期并合并数据集？

现有数据集 #

我们有两个SAS数据集：

• Have_1：存储各ID的月度数据

ID   Month 
1     1
1     2
1     3
1     4
1     5
1     6
1     7
1     8
1     9
1     10

• Have_2：记录每个ID的事件发生月份

ID   Month_Event 
1     4
1     9

需求说明 #

需要将Have_2中的Month_Event匹配到Have_1，规则是：每个ID的Month对应大于等于该Month的最近Month_Event；如果没有符合条件的事件（比如Month超过所有事件月份），则Month_Event置为空。目标输出如下：

ID   Month   Month_Event
1     1      4
1     2      4
1     3      4
1     4      4
1     5      9
1     6      9
1     7      9
1     8      9
1     9      9
1     10     .

当前代码问题 #

你提供的PROC SQL左连接代码会产生重复行，因为同一Month会匹配多个满足a.Month <= b.Month_Event的Month_Event，无法直接得到目标结果：

proc sql;
create table want as
select a.*, b.Month_Event
from have_1 as a
left join have_2 as b
on a.ID = b.ID and a.Month <= b.Month_Event;
quit;

解决方案 #

方法1：PROC SQL子查询（简洁高效） #

通过子查询为每个ID和Month筛选出最小的、大于等于当前Month的Month_Event，也就是最近的后续事件：

proc sql;
create table want as
select 
    a.ID,
    a.Month,
    (select min(b.Month_Event) 
     from have_2 as b 
     where b.ID = a.ID and b.Month_Event >= a.Month) as Month_Event
from have_1 as a;
quit;

方法2：数据步合并（适合复杂逻辑扩展） #

先对两个数据集按ID和月份排序，再用数据步将事件收集到临时数组，逐个匹配每个Month对应的最近事件：

/* 先排序数据集 */
proc sort data=have_1; by ID Month; run;
proc sort data=have_2; by ID Month_Event; run;

/* 数据步匹配逻辑 */
data want;
    merge have_1(in=a) have_2(in=b);
    by ID;
    array events[100] _temporary_; /* 临时数组存储当前ID的所有事件月份，可根据实际调整长度 */
    retain event_count;

    /* 每个ID开始时初始化数组和计数 */
    if first.ID then do;
        call missing(of events[*]);
        event_count = 0;
    end;

    /* 将当前ID的事件月份存入数组 */
    if b then do;
        event_count + 1;
        events[event_count] = Month_Event;
    end;

    /* 为每个月度数据匹配最近事件 */
    if a then do;
        Month_Event = .; /* 默认置空 */
        /* 遍历事件数组，找到第一个大于等于当前Month的事件 */
        do i = 1 to event_count;
            if events[i] >= Month then do;
                Month_Event = events[i];
                leave; /* 找到后跳出循环，保证是最近的 */
            end;
        end;
        output; /* 输出当前行 */
    end;

    drop event_count i;
run;

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