如何在SAS数据集中计算并添加±1周、±2周的时间滞后变量？

针对你的SAS数据集，要添加**前1周（lag-1）、前2周（lag-2）、后1周（lag+1）、后2周（lag+2）**的变量，推荐以下两种高效方法，可根据数据集规模选择：

方法一：DATA步数组法（大数据集首选，效率最高） #

这种方法先将全量数据读入内存临时数组，再循环输出时直接通过索引访问前后周的数据，仅需读取数据集两次，IO开销极低，适合大数据量场景。

/* 先获取数据集的总观测数，用于数组定义 */
proc sql noprint;
select count(*) into :nobs from dataset;
quit;

data dataset_with_lags_leads;
    /* 定义临时数组存储所有周数、Var1、Var2的值 */
    array wk[&nobs] _temporary_;
    array v1[&nobs] _temporary_;
    array v2[&nobs] _temporary_;

    /* 第一步：读取所有数据到临时数组 */
    do i=1 to &nobs;
        set dataset;
        wk[i] = week_number;
        v1[i] = Var1;
        v2[i] = Var2;
    end;

    /* 第二步：循环输出每个观测，计算滞后/超前变量 */
    do i=1 to &nobs;
        week_number = wk[i];
        Var1 = v1[i];
        Var2 = v2[i];

        /* 计算前1周、前2周变量（滞后） */
        lag1_Var1 = ifn(i>1, v1[i-1], .);
        lag1_Var2 = ifn(i>1, v2[i-1], .);
        lag2_Var1 = ifn(i>2, v1[i-2], .);
        lag2_Var2 = ifn(i>2, v2[i-2], .);

        /* 计算后1周、后2周变量（超前） */
        lead1_Var1 = ifn(i<&nobs, v1[i+1], .);
        lead1_Var2 = ifn(i<&nobs, v2[i+1], .);
        lead2_Var1 = ifn(i<&nobs-1, v1[i+2], .);
        lead2_Var2 = ifn(i<&nobs-1, v2[i+2], .);

        output;
    end;

    drop i;
run;

方法二：PROC SQL自连接法（中小数据集首选，代码简洁） #

如果数据集规模不大，用SQL自连接逻辑更直观，通过week_number的加减匹配前后周数据，代码易读易维护。

proc sql;
create table dataset_with_lags_leads as
    select 
        a.week_number,
        a.Var1,
        a.Var2,
        /* 前1周：匹配当前周数-1的记录 */
        b.Var1 as lag1_Var1,
        b.Var2 as lag1_Var2,
        /* 前2周：匹配当前周数-2的记录 */
        c.Var1 as lag2_Var1,
        c.Var2 as lag2_Var2,
        /* 后1周：匹配当前周数+1的记录 */
        d.Var1 as lead1_Var1,
        d.Var2 as lead1_Var2,
        /* 后2周：匹配当前周数+2的记录 */
        e.Var1 as lead2_Var1,
        e.Var2 as lead2_Var2
    from dataset a
    left join dataset b on a.week_number = b.week_number + 1
    left join dataset c on a.week_number = c.week_number + 2
    left join dataset d on a.week_number = d.week_number - 1
    left join dataset e on a.week_number = e.week_number - 2
    order by a.week_number;
quit;

方法对比 #

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