【请教】SAS回归模型保存

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

问题是这样的：我在A数据集中得到了回归模型，
怎么样把这个回归模型运用到其他数据集B上呢？
两个数据集的 列名设置为一样的，
但是，由于自变量比较多，不能手工输入，
所以请教一下大家~非常谢谢~

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：回归模型 数据集 自变量 怎么样 请教 模型 SAS 保存

这是一个较为复杂的宏
等待高手解答

谢谢！1！！！！

你用得到什么呢？
如果是多重线性回归，通过A数据集的回归模型来来获得B数据集的预测值，其实思路可以这样
先用APPEND，当然前提两者的列名设置为一样的
在程序运行之后回自动计算预测值的

回LS的~想知道这个模型在B数据集上的表现如何~是好是坏~

简单线性回归+SAS9.1.3

data height;
    input height;
        datalines;
63
61
59
62
;



%macro reg(modelds=sashelp.class
      ,model=%str(weight=height)
      ,testds=height
          ,outds=result
       ,Yvar=weight
          );
   proc reg data=&modelds;
     model &model;
         ods output ParameterEstimates=para_est;
run;
quit;

data _null_;
     set para_est;
         call symputx(variable,estimate,l);
run;

proc sql noprint;
     select distinct variable into : vars separated by " "
           from para_est
       where upcase(variable)^="INTERCEPT";
quit;

data &outds;
    set &testds;
        &yvar=&intercept
        %do i=1 %to &sqlobs;
            %let var&i=%scan(&vars,&i);
                        +&&var&i*&&&&&&var&i
                %end;
                ;
run;
             
%mend;

options mprint;
%reg()

手上正好有代码，剥出来看看
&Var_Lst是所有变量的列表，用空格分开，不包括ID变量以及目标变量，例如%let Var_Lst=x1 x2 x3;
或者proc contents来获得
&Target 是目标变量 例如 %let Target=y;
Coefficient数据集包括回归模型结果的系数,并且变量名与B中相对应

1. %MACRO T;
   
2. %let Var_Cnt=%eval(%sysfunc(count(&Var_Lst,%str( )))+1);
   
3. %put &Var_Cnt;
   
4. proc sql;
   
5.         create table SCORE as
   
6.                 select
   
7.                         %do i=1 %to &Var_Cnt %by 1;%if &i>=2 %then %do;,%end;a.%scan(&Var_Lst,&i)%end;
   
8.                         ,a.&Target
   
9.                         sum(%do i=1 %to &Var_Cnt %by 1;%if &i>=2 %then %do;,%end;a.%scan(&Var_Lst,&i)*b.%scan(&Var_Lst,&i)%end;,b.Intercept) as SCORE
   
10.                 from B a
    
11.                 inner join  Coefficient b
    
12.                 on        1=1
    
13.                 order by SCORE descending
    
14.         ;
    
15. quit;
    
16. %MEND T;

复制代码

另外EM中的SCORE节点也可以,proc score也可以......
我习惯这样做了.....

恩恩~非常谢谢LS的两位~

进来学习了

crackman 发表于 2010-5-20 11:31
你用得到什么呢？
如果是多重线性回归，通过A数据集的回归模型来来获得B数据集的预测值，其实思路可以这样
先用APPEND，当然前提两者的列名设置为一样的
在程序运行之后回自动计算预测值的

Agree.
Moreover, I think, it also depends on the criteria by which the good model is defined.
Assume the ‘Residual’, a difference between predicted value by model from set A and the real Y in set B, is of your interest.

1. 
   
2. DATA AB;
   
3. SET  A B(IN = b RENAME = (y = yb)); *Let's say response as Y;
   
4. IF b THEN y = .A; *to differ the possible missing (.) in set A;
   
5. RUN;
   
6. ods output OutputStatistics=os(where = (depVar = .A));
   
7.    proc reg data=AB;   
   
8.       model y = x1 x2 / r cli clm; *assume fit x1 x2 as independent vars;
   
9.    run;
   
10. quit;
    
11. data C;
    
12. retain type 'Predicted ';
    
13. set  B; set os;
    
14. Residual = PredictedValue-y;
    
15. run;

复制代码