重金悬赏：帮忙给写一个SAS macro来处理重复的数据文件操作

恳请各位大侠，帮忙给写一个SAS macro来处理重复的数据文件操作。如果在写这个macro的时候果被卡到了什么地方，没有关系，即使只是提供部分方案，我也感激不尽，100个论坛币悬赏照送。此外，写这个macro完全没有必要了解任何的统计知识。这个macro完全是个数据文件操作问题(data step, proc sql)。

正文如下：
我有一个数据文件，my_data，包含一个y变量和上百个x变量。为了举例方便，这里用3个x变量来演示。下边这个程序生成我的数据文件，my_data。

DATA my_data;
INPUT y x1 x2 x3;
DATALINES;
10        1       1        1
20        1       2        2
30        2       3        3
40        2       4        4
50        3       4        4
60        3       5        4
70        4       5        5
80        4       6        5
90        5       7        6
;
RUN;

我把（x1 x2 x3）这个组合放入到下边这段程序中，这段程序会生成一个数据文件，model_parms。（你只需要知道这些就足够了）。

PROC REG DATA=my_data;
MODEL y=x1 x2 x3; *对于不同x变量组合，这是唯一需要改动的一行程序；
ODS OUTPUT ParameterEstimates=model_parms;
RUN;QUIT;

对于生成的这个数据文件，model_parms，稍作如下处理。
DATA model_parms;
SET model_parms;
KEEP Variable Estimate Probt;
IF Variable NE "Intercept";
RUN;

对于生成的这个数据文件，model_parms，我只看两个东西：

1是看x变量的符号对不对（所有x变量的估计参数都应该是负号）。通过观察Parameter Estimate，我发现2/3(或者66%)的参数符号不对（x1,x2符号不对，x3符号对了）。

2是看P-Value是不是小于0.05。通过观察Pr > |t|，我发现2/3(或者66%)的参数的P-Value小于0.05（x1,x2的P-Value小于0.05，x3的P-Value大于0.05）。

目前为止，我只看了一种x变量的组合情况，即（x1 x2 x3）这个组合。

下面，

我还需要看（x1 x2）这个组合（只需把上面的MODEL y=x1 x2 x3 这行程序变动成为MODEL y=x1 x2）。

我还需要看（x2 x3）这个组合（只需把上面的MODEL y=x1 x2 x3这行程序变动成为MODEL y=x2 x3）。

最后，我要把这些x变量的组合按照这个数据文件，model_parms，符号正确率和P-Value小于0.05率，来排一下序，选出一个最佳x变量的组合。对于这个最佳x变量的组合来说，其估计参数的符号正确率最高，而且P-Value小于0.05的率同样的高。

请大侠帮帮忙，把上面的这个例子写成一个macro。这个macro不需要遍历所有x变量的组合, 用户可以指定x变量的各种组合。在这个例子了，我要看看，下面这3个x变量的组合，哪个最好：（x1 x2 x3）这个组合，（x1 x2）这个组合，还是（x2 x3）这个组合？

如果我有什么地方解释的不清楚，尽情按照你的理解来写这个macro即可。说一千，道一万，我的根本问题是，给定一些x变量的组合，按照符号正确率和P-Value小于0.05率找到个最好的那个组合：（x1 x2 x3）这个组合，（x1 x2）这个组合，还是（x2 x3）这个组合？

你这里就一个遍历回归，再加一个结果的筛选
后者没有必要和前者放在一个宏里
我这里就大概说下前者

你不需要全部遍历，如果你把需要遍历自变量放在一个表里，如下所示：
x1 x2 x3 x4 ... x88
1   1   0  0  ...  0
即一行就是一种情况，=1代表取这个变量，否则=0，那么可以用下面的宏来解决你的问题

%macro IterReg(SourceTable,VarTable,YVar,TargetTable);

data IR_VarList(drop=i);
        set &VarTable;
        array xs{*} x:;
        length Var_Name $100;
        Var_Name='';
        do i=1 to dim(xs);
                if xs{i} EQ 1 then Var_Name=catx(' ',Var_Name,vname(xs{i}));
        end;
        N=_N_;
run;

proc sql noprint;
        select COUNT(*) into :IR_IterTime from &IR_VarList;
quit;

%do i=1 %to &IR_IterTime;
        proc sql noprint;
                select Var_Name into :IR_VarName
                        from IR_VarList
                        where N=&i;
        quit;

        proc reg data=&SourceTable outest=IR_Result noprint;
                model &YVar = &IR_VarName / EDF ADJRSQ;
        quit;

        %if &i EQ 1 %then %do;
               data &TargetTable;
                       set IR_Result;
               run;
        %end;
        %else %do;
               proc append base=&TargetTable data=IR_Result;
               quit;
         %end;
%end;

%mend;

每一行做y与x1 x2....的回归？

对，每一行都做y与x1 x2....的回归，
y=x1 x2 x3 x4 ...
y=x2 x3 x5 x8 ...y=x2 ...
y= x5 x8 ...
y=x2 x3 x8 ...
.
.
.
不要试图使用proc reg的某些功能，因为这里只是以proc reg 为例。按照我在例子里讲的，需要设计一个几乎基于数据文件操作(data step, proc sql)的一个macro来实现。没有必要使用任何回归方面的知识。
谢谢关注和回复。

为啥不用selection=forward选项呢？

不要试图使用proc reg的某些功能，因为这里只是以proc reg 为例。最终我使用其他回归程序，这里只是以proc reg 为例来说明我的问题。

mark