单变量分析的sas宏

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

如题所示，下面给出一个单变量分析的sas宏，对建立评分卡的探索分析工作非常有帮助，只需输入一个数据集，数据集的每一个变量的缺失情况，0值情况，排序后按groups个数百分比占比分析等等，简单实用，去掉工作中的重复劳动，其中sasin=是输入数据集，其他为空则分析全部变量，好用请点赞哦！

1. %MACRO CROSSTAB(SASIN=,VARLIST=,ADDVAR=,DROPLIST=,TARGET=,GROUPS=10,
   
2.                 PRTORDER=PVALUE,RISKTITL=1 -- HIGH RISK);
   
3. 
   
4.   /*------------------------------------------------------------------
   
5.     THIS MACRO DIVIDES EACH VARIABLE IN THE 'VARLIST' INTO
   
6.     DECILES, COMPUTE THE 'TARGET' RATES (EG. RESPONSE RATE, BAD RATE,
   
7.     ETC.) IN EACH DECILE, AND MEAN OF ALL 'ADDVAR' VARIABLES
   
8. 
   
9.     PRTORDER: TAKES VALUE OF 'PVALUE' OR 'VARIABLE', IT DETERMINES THE
   
10.               PRINT ORDER: BY P-VALUE OR ALPHABETIC ORDER OF VARIABLE NAMES
    
11.     DROPLIST: OVERRIDES VARLIST
    
12.   -------------------------------------------------------------------*/
    
13. 
    
14. OPTIONS NOCENTER PS=120 PAGENO=1;
    
15. 
    
16. * DM 'LOG'    CLEAR;
    
17. * DM 'OUTPUT' CLEAR;
    
18. 
    
19.   %IF &TARGET= %THEN %DO;
    
20.     %LET TARGET=DUMMYZZZ;
    
21.   %END;
    
22. 
    
23. 
    
24.   DATA SASINZZZ;
    
25.     SET &SASIN;
    
26.     %IF &TARGET=DUMMYZZZ %THEN %DO;
    
27.       &TARGET=1;
    
28.     %END;
    
29.     %IF &DROPLIST NE %THEN %DO;
    
30.       DROP &DROPLIST;
    
31.     %END;
    
32.     %ELSE %IF &VARLIST NE %THEN %DO;
    
33.       KEEP &VARLIST &TARGET &ADDVAR;
    
34.     %END;
    
35.     RUN;
    
36. 
    
37. 
    
38.   PROC CONTENTS DATA=SASINZZZ (DROP=&TARGET &ADDVAR)
    
39.     OUT=DTCONTS (KEEP=NAME LABEL TYPE) NOPRINT;
    
40.     RUN;
    
41. 
    
42.   DATA DTCONTS;
    
43.     SET DTCONTS;
    
44.     NAME=UPCASE(NAME);
    
45.     RUN;
    
46. 
    
47.   PROC SORT DATA=DTCONTS;
    
48.     BY NAME;
    
49.     RUN;
    
50. 
    
51.   DATA DTCONTS;
    
52.     SET DTCONTS;
    
53.     RANK='RK' || LEFT(_N_);
    
54.     RUN;
    
55. 
    
56. 
    
57.   DATA _NULL_;
    
58.     SET DTCONTS END=FINAL;
    
59.     I+1;
    
60.     CALL SYMPUT('VAR_A'||LEFT(I),TRIM(NAME));
    
61.     IF TYPE=1 THEN CALL SYMPUT('RK_A'||LEFT(I),TRIM(RANK));
    
62.     ELSE           CALL SYMPUT('RK_A'||LEFT(I),TRIM(NAME));
    
63.     CALL SYMPUT('LBL'|| LEFT(I),TRIM(LABEL));
    
64.     CALL SYMPUT('TYPE'||LEFT(I),TYPE);
    
65.     IF FINAL THEN CALL SYMPUT('N_TZ',LEFT(I));
    
66.     RUN;
    
67. 
    
68. 
    
69.   %DO I=1 %TO &N_TZ;
    
70.     %LET SASMEAN=SASINZZZ;
    
71.     %IF &&TYPE&I=1 %THEN %DO;
    
72.       %LET SASMEAN=SASINTMP;
    
73.       PROC RANK DATA=SASINZZZ (KEEP=&&VAR_A&I &TARGET &ADDVAR)
    
74.                 OUT=SASINTMP GROUPS=&GROUPS;
    
75.         VAR &&VAR_A&I;
    
76.         RANKS RANKTMP;
    
77.         RUN;
    
78. 
    
79. 
    
80.       DATA SASINTMP;
    
81.         SET SASINTMP;
    
82.         &&RK_A&I=PUT(RANKTMP,2.);
    
83.         DROP RANKTMP;
    
84.         RUN;
    
85.     %END;
    
86. 
    
87. 
    
88.     PROC MEANS DATA=&SASMEAN (KEEP=&&VAR_A&I &TARGET &&RK_A&I &ADDVAR)
    
89.                MISSING NOPRINT;
    
90.       CLASS &&RK_A&I;
    
91.       %IF &&TYPE&I = 1 %THEN %DO;
    
92.         VAR &&VAR_A&I &TARGET &ADDVAR;
    
93.         OUTPUT OUT=MEANOZZZ MEAN(&TARGET &ADDVAR)=
    
94.                             MEAN(&&VAR_A&I)=MEAN_
    
95.                             MIN(&&VAR_A&I)=MIN_
    
96.                             MAX(&&VAR_A&I)=MAX_;
    
97. 
    
98.       %END;
    
99.       %ELSE %DO;
    
100.         VAR &TARGET &ADDVAR;
     
101.         OUTPUT OUT=MEANOZZZ MEAN=;
     
102.       %END;
     
103.       RUN;
     
104. 
     
105. 
     
106.     *** CALCULATE P-VALUE FOR ROW AND COLUMN INDEPENDENCE;
     
107. 
     
108.     %LET DOCHISQ=Y;
     
109.     DATA _NULL_;
     
110.       SET MEANOZZZ END=FINAL;
     
111.       IF _TYPE_=0 AND (&TARGET=0 OR &TARGET=1) THEN
     
112.         CALL SYMPUT('DOCHISQ','N');
     
113.       IF FINAL AND _N_=2 THEN
     
114.         CALL SYMPUT('DOCHISQ','N');
     
115.       RUN;
     
116. 
     
117. 
     
118.     DATA CHISQA;
     
119.       SET MEANOZZZ;
     
120.       IF _TYPE_>0;
     
121.       GROUPZ='A';
     
122.       WGT=ROUND(_FREQ_*&TARGET,1);
     
123.       RUN;
     
124. 
     
125. 
     
126. 
     
127.     DATA CHISQB;
     
128.       SET MEANOZZZ;
     
129.       IF _TYPE_>0;
     
130.       GROUPZ='B';
     
131.       WGT=_FREQ_-ROUND(_FREQ_*&TARGET,1);
     
132.       RUN;
     
133. 
     
134. 
     
135.     DATA CHISQ;
     
136.       SET CHISQA CHISQB;
     
137.       KEEP &&RK_A&I GROUPZ WGT;
     
138.       RUN;
     
139. 
     
140.     %IF &DOCHISQ=N %THEN %DO;
     
141.       DATA FREQOUT;
     
142.         LENGTH VARNAME $15;
     
143.         _PCHI_=0;
     
144.         P_PCHI=1;
     
145.         VARNAME="&&VAR_A&I";
     
146.         RUN;
     
147. 
     
148.     %END;
     
149.     %ELSE %DO;
     
150.       PROC FREQ DATA=CHISQ NOPRINT;
     
151.         WEIGHT WGT;
     
152.         TABLES &&RK_A&I*GROUPZ/CHISQ MISSING;
     
153.         OUTPUT OUT=FREQOUT (KEEP=_PCHI_ P_PCHI) PCHI;
     
154.         RUN;
     
155. 
     
156.       DATA FREQOUT;
     
157.         LENGTH VARNAME $15;
     
158.         SET FREQOUT;
     
159.         VARNAME="&&VAR_A&I";
     
160.         RUN;
     
161.     %END;
     
162. 
     
163. 
     
164.     PROC APPEND BASE=PVALUE DATA=FREQOUT;
     
165.     RUN;
     
166. 
     
167.     **** CREATE PRINT OUT FILES;
     
168. 
     
169.     DATA MEANOZZZ;
     
170.       SET MEANOZZZ;
     
171.       RETAIN TOTAL AVGODDS;
     
172. 
     
173. 
     
174.       IF &TARGET=1 THEN ODDS=0;
     
175.       ELSE ODDS=&TARGET/(1 - &TARGET);
     
176.       IF _TYPE_=0 THEN DO;
     
177.         TOTAL=_FREQ_;
     
178.         AVGODDS=ODDS;
     
179.       END;
     
180.       IF _TYPE_=1 AND &&RK_A&I=' ' THEN &&RK_A&I='-1';
     
181. 
     
182. 
     
183.       FREQPCT=_FREQ_/TOTAL*100;
     
184.       &TARGET=&TARGET*100;
     
185.       IF AVGODDS=0 THEN STDODDS=0; ELSE STDODDS=ODDS/AVGODDS;
     
186.       RUN;
     
187. 
     
188.     PROC SORT DATA=MEANOZZZ OUT=MEANO&I;
     
189.       BY DESCENDING _TYPE_;
     
190.       RUN;
     
191. 
     
192. 
     
193.     %IF %EVAL(&I/50)*50 = &I %THEN %DO;
     
194.       OPTIONS NOTES;
     
195.       PROC DATASETS;
     
196.         DELETE SASINTMP MEANOZZZ CHISQA CHISQB CHISQ FREQOUT;
     
197.         QUIT;
     
198.     %END;
     
199.     %ELSE %DO;
     
200.       OPTIONS NONOTES;
     
201.     %END;
     
202.   %END;
     
203. 
     
204. 
     
205.   DATA PVALUE;
     
206.     SET PVALUE;
     
207.     ORD_ALPH=_N_;
     
208.     CALL SYMPUT('PCHI'|| LEFT(_N_),TRIM(PUT(_PCHI_,8.)));
     
209.     RUN;
     
210. 
     
211.   PROC SORT DATA=PVALUE; BY DESCENDING _PCHI_; RUN;
     
212. 
     
213.   DATA PVALUE;
     
214.     SET PVALUE;
     
215.     ORD_PCHI=_N_;
     
216.     CALL SYMPUT('OALPH'|| LEFT(_N_),TRIM(ORD_ALPH));
     
217.     RUN;
     
218. 
     
219. 
     
220.   %MACRO PRINTIT(I);
     
221.     PROC PRINT DATA=MEANO&I (DROP=_TYPE_) SPLIT='*';
     
222.       VAR &&RK_A&I
     
223.       %IF &&TYPE&I=1 %THEN %DO;
     
224.         MIN_ MAX_ MEAN_
     
225.       %END;
     
226.       _FREQ_ FREQPCT
     
227.       %IF &TARGET NE DUMMYZZZ %THEN %DO;
     
228.         &TARGET STDODDS
     
229. 
     
230.       %END;
     
231.       &ADDVAR;
     
232.       LABEL
     
233.         MIN_='MIN'
     
234.         MAX_='MAX'
     
235.         MEAN_='MEAN'
     
236.         &&RK_A&I='RANK/VALUE'
     
237.         _FREQ_='FREQ'
     
238.         FREQPCT='%FREQ'
     
239.         &TARGET='%1'
     
240.         STDODDS='STANDARDIZED*ODDS RATIO(%1)'
     
241.         ;
     
242. 
     
243.       FORMAT &&RK_A&I $RANKFT. FREQPCT &TARGET STDODDS
     
244.              MEAN_ &ADDVAR 8.4;
     
245.       RUN;
     
246.   %MEND PRINTIT;
     
247. 
     
248.   TITLE3 "TABLE FOR DEPENDENT VARIABLE: &TARGET, &RISKTITL";
     
249.   TITLE4 " ";
     
250. 
     
251. 
     
252.   %DO J=1 %TO &N_TZ;
     
253.     %IF %UPCASE(%SUBSTR(&PRTORDER,1,1))=P %THEN %LET I=&&OALPH&J;
     
254.     %ELSE %LET I=&J;
     
255.     %IF &J=1 %THEN %DO;
     
256.       TITLE5
     
257.     %END;
     
258.     %ELSE %DO;
     
259.       TITLE1 " ";
     
260.       TITLE2
     
261.     %END;
     
262.     "&&VAR_A&I --- &&LBL&I, _PCHI_ --- &&PCHI&I";
     
263.     %PRINTIT(&I);
     
264.   %END;
     
265. 
     
266. 
     
267.   %IF %UPCASE(%SUBSTR(&PRTORDER,1,1))=P %THEN %DO;
     
268.     PROC SORT DATA=PVALUE; BY ORD_ALPH; RUN;
     
269. 
     
270.     TITLE1 " ";
     
271.     TITLE2 "ALPHABETIC ORDER OF VARIABLES";
     
272. 
     
273.     DATA _NULL_;
     
274.       SET PVALUE;
     
275.       FILE PRINT LS=110;
     
276.       PUT
     
277.       VARNAME $15. +2
     
278.       ORD_PCHI 4.  +1
     
279.       @@;
     
280.       RUN;
     
281. 
     
282. 
     
283.     PROC SORT DATA=PVALUE; BY ORD_PCHI; RUN;
     
284. 
     
285.     TITLE1 " ";
     
286.     TITLE2 "P-VLUE ORDER OF VARIABLES";
     
287. 
     
288.     DATA _NULL_;
     
289.       SET PVALUE;
     
290.       FILE PRINT LS=110;
     
291.       PUT
     
292.       VARNAME $15. +2
     
293.       ORD_PCHI 4.  +1
     
294.       @@;
     
295.       RUN;
     
296.   %END;
     
297. 
     
298. 
     
299.   OPTIONS NOTES;
     
300. 
     
301.   PROC DATASETS;
     
302.     DELETE DTCONTS SASINZZZ MEANO1 - MEANO&N_TZ PVALUE;
     
303.     QUIT;
     
304.   TITLE;
     
305. %MEND CROSSTAB;

复制代码

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：SAS宏 单变量 independence Standardized standardize

sas code贴上来之后发现格式有点乱了，不过复制到txt文档里阅读无障碍，使用也无障碍！