连续变量分箱Logit图

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

     在构建线形回归模型时，我们用观看x与y之间散点图的方式确定x以何种方式纳入回归模型，比如是取二次方形式，还是对数形式。

      而做逻辑回归时，也需要一种图形观看连续变量x与logit(y)的关系，可以SAS本身没有提供，需要我们自行绘制。以下就是做好的一个宏。

1. /****************************************************************************/
   
2. %Macro E_Logit(data,target,var,bins);
   
3. /*第一步:对变量进行RANK分组*/
   
4. proc rank data=&data groups=&bins out=out;
   
5.    var &var;
   
6.    ranks bin;
   
7. run;

复制代码

1. /*第二步:对每一组计算该变量的平均值;响应事件数和总事件数**/
   
2. 
   
3. /*数据集BINS 包含:          */
   
4. /* &target = 每个BIN里面响应事件数 */
   
5. /* _FREQ_ =每个BIN里面总事件数 */
   
6. /* &var =每个BIN里面&var平均值 */
   
7. 
   
8. proc means data=out noprint nway;
   
9.    class bin;
   
10.    var &target &var;
    
11.    output out=bins sum(&target)=&target mean(&var)=&var;
    
12. run;

复制代码

1. /*第三步:根据公式计算 empirical logit */
   
2. data bins;
   
3.    set bins;
   
4.    elogit=log((&target+(sqrt(_FREQ_ )/2))/
   
5.           ( _FREQ_ -&target+(sqrt(_FREQ_ )/2)));
   
6. run;

复制代码

1. /*第四步:画LOGIT与原变量平均值;LOGIT与BIN变量的线图*/
   
2. /*proc sgplot data = bins;*/
   
3. /*title "Empirical Logit against &var";*/
   
4. /*series y=elogit x=&var;*/
   
5. /*scatter y=elogit x=&var;*/
   
6. /*run;*/
   
7. proc sgplot data = bins;
   
8. title "Empirical Logit against Binned &var";
   
9. scatter y=elogit x=bin;
   
10. series y=elogit x=bin;
    
11. run;quit;

复制代码

1. /*第五步:用BIN变量替代原来的变量,并对BIN变量进行代码保存和改造*/
   
2. 
   
3. title;
   
4. 
   
5. proc means data = out noprint nway;
   
6.    class bin;
   
7.    var  &var;
   
8.    output out=endpts  max(&var)=max;
   
9. run;
   
10. 
    
11. filename rank "d:\rank.sas";
    
12. 
    
13. /*编写BIN代码*/
    
14. data _null_;
    
15.    file rank;
    
16.    set endpts end=last;
    
17.    if _n_ = 1 then put "select;";
    
18.    if not last then do;
    
19.      put "  when (&var <= " max ") B_&var =" bin ";";
    
20.      end;
    
21.    else if last then do;
    
22.      put "otherwise B_&var =" bin ";";
    
23.      put "end;";
    
24.    end;
    
25. run;
    
26. 
    
27. /* Use the code. */
    
28. data &data;
    
29.    set &data;
    
30.    %include rank /source2;
    
31. run;
    
32. 
    
33. /*proc means data = &data min max;*/
    
34. /*   class B_&var;*/
    
35. /*   var &var;*/
    
36. /*run;*/
    
37. proc delete data=out bins endpts;run;
    
38. %Mend E_Logit;
    
39. 
    
40. /****************************************************************************/
    
41. /*%Macro E_Logit(数据集,二分目标变量,连续自变量,分箱个数);*/
    
42. %E_Logit(gg,flag,score,5);
    
43. run;
    
44. %E_Logit(tel,churn,avPlan,4);
    
45. run;

复制代码

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：Logit图 logit 连续变量 Log Empirical groups target 平均值 散点图 模型 连续变量 分箱 Logit图 sas

好贴，谢谢分享。