数据中有1153个变量，其中后面的1151个变量是我想要进行分析的数据。我想每一个变量都作一次逐步回归，即变量p1作为y时，其他剩余的p2-p1151都作为进行逐步回归的x，算完p1以后，再以p2作为y，p1和p3-p1151作为x的进行逐步回归。并且每个变量的观测值最多只有30个，而变量一共有一千多个，所以希望限制最终保留在方程中的最优的x的数量（例如10个）。
因为数据中并不是每个变量都有完整30个数据，有些是缺失的，所以我是想通过建立每个变量的逐步回归方程，然后数据缺失的地方用方程计算的预测值进行插值。
我想这应该需要写一个循环，而且我希望所有变量的参数可以都放在一张表上，所有变量的预测值放在一张表上。
我尝试写了一下，但是运行的时候说宏解释错误。
希望哪位大侠帮忙修改一下。我对宏和循环不太了解，所以希望可以帮忙写出完整的语句。附件是我说的数据
非常感谢！

%macro reg1(i,j) ;
%do i = 2 %to 1151 ;
j =1 %to (i-1) and j = (i+1) %to 1151;
proc reg data = p_stepwise outest = p_reg_log&j noprint;
model p&i = p1 - p&j/ edf selection = stepwise slstay = 0.10;
output out= p_stepwise_reg p = pre&i;
run;
%end;
%mend regl;
%regl;