SAS数组计算日增量变量问题：生成的b变量全部缺失求助

看起来你遇到的问题核心是数组下标和变量命名的错位，加上对b_0的处理缺失，导致所有b变量都变成了缺失值。我来帮你拆解问题并修正代码：

问题根源分析 #

1. 数组下标默认规则的坑：SAS中数组默认的下标是从1开始的，但你的变量是a_0到a_45（后缀从0到45）。当你写array a a_0-a_45;时，SAS会把a_0当作数组的第1个元素（a(1)=a_0），a_1当作第2个元素（a(2)=a_1）……直到a_45是数组的第46个元素（a(46)=a_45）。这就导致你后续的下标引用完全错位，自然计算不出正确值。
2. b_0未被赋值：你的需求是生成b_0到b_45，但原代码里只处理了b(1)到b(45)，b_0完全没被赋值，也会是缺失值。

修正后的代码 #

data test;
    set dataset;
    /* 明确指定数组下标范围为0到45，和变量名的后缀一一对应 */
    array a[0:45] a_0-a_45;
    array b[0:45] b_0-b_45;
    
    /* 处理b_0：根据业务需求定义，这里假设初始持有量就是b_0的取值 */
    b[0] = a[0]; 
    /* 如果b_0不需要初始值，可改为 b[0] = .; 设为缺失 */
    
    /* 循环计算从b_1到b_45的日增量 */
    do i = 1 to 45;
        b[i] = a[i] - a[i-1];
    end;
run;

关键说明 #

• 下标对齐：通过array a[0:45]强制让数组下标从0开始，这样a[0]对应a_0，a[1]对应a_1，以此类推，彻底避免下标错位问题。
• b_0的逻辑：因为第0天没有前一天，所以你需要根据业务场景定义b_0的值——比如初始持有量就是b_0，或者设为缺失值，根据你的实际需求调整即可。
• 循环范围：循环从i=1到45，确保每个b[i]都能正确计算为当日持有量减去前一日的持有量。

另外，建议你检查原数据集dataset中是否确实存在a_0到a_45这些变量，有没有拼写错误（比如把a_0写成a0），这也是导致缺失值的常见原因。

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