%let lag_n=60;
data final_data2(drop=index);
   array pp[0:%eval(&lag_n-1)] _temporary_;   /* 创建数组，下标0-59 */
   index = 0; /* 为下标变量赋初值 */
      set raw_data;   
      pp[index] = p; /* 为数组元素赋p值 */
      index = mod(index+1, &lag_n);   /* 下标变量+1，取模保证值域为0-59 */
     /* 计算部分 */
      if _n_>=&lag_n then do;
         mean_p = mean(of pp[*]);
         var_p = var(of pp[*]);
      end;
run;

看看上面的对不？？？

%let lag_n=60;
data final_data2(drop=index);
   array pp[0:%eval(&lag_n-1)] _temporary_;   /* 创建数组，下标0-59 */
retain  index  0; /* 为下标变量赋初值 */
      set raw_data;   
      pp[index] = p; /* 为数组元素赋p值 */
      index = mod(index+1, &lag_n);   /* 下标变量+1，取模保证值域为0-59 */
     /* 计算部分 */
      if _n_>=&lag_n then do;
         mean_p = mean(of pp[*]);
         var_p = var(of pp[*]);
      end;
run;

看看最后这个版本行不？？？

/*假设原数据集有变量date和p两字段*/
proc sort data= raw_data;
by date;
run;

%let lag_n=60;
data final_data2(drop=index);
   array pp[0:%eval(&lag_n-1)] _temporary_;   /* 创建数组，下标0-59 */
retain  index  0; /* 为下标变量赋初值 */
      set raw_data;   
      pp[index] = p; /* 为数组元素赋p值 */
      index = mod(index+1, &lag_n);   /* 下标变量+1，取模保证值域为0-59 */
     /* 计算部分 */
      if _n_>=&lag_n then do;
         mean_p = mean(of pp[*]);
         var_p = var(of pp[*]);
      end;
run;