SAS基礎(chǔ)宏之8：GetMissNum

是 否 +2 論壇幣

k人 參與回答

經(jīng)管之家送您一份

應(yīng)屆畢業(yè)生專屬福利!

求職就業(yè)群

趙安豆老師微信：zhaoandou666

經(jīng)管之家聯(lián)合CDA

送您一個全額獎學(xué)金名額~ !

立即領(lǐng)取

感謝您參與論壇問題回答

經(jīng)管之家送您兩個論壇幣！

+2 論壇幣

這個宏用來統(tǒng)計數(shù)據(jù)集中指定變量的缺失值數(shù)量，支持數(shù)值變量和字符變量
最后都匯到一張結(jié)果表格中

之前我用的是把表格轉(zhuǎn)置后，用nmiss和cmiss求橫向的數(shù)組的缺失值
這樣代碼比較簡單，但是后來發(fā)現(xiàn)proc transpose實在是太慢了，而且耗資源
這下改了下，直接用sql對列進行計算了，不需要轉(zhuǎn)置了

1. %macro GetMissNum(SourceTable,TargetTable,InputVar);
   
2. /**********************************************************************/
   
3. /* 此宏的作用是統(tǒng)計原表中不同變量的缺失值數(shù)量。其中SourceTable是原始 */
   
4. /* 表格，SourceTable是結(jié)果表格；TargetTable是結(jié)果表格；InputVar是原始 */
   
5. /* 表格中的變量，可設(shè)多個變量，用空格分隔，也可如下設(shè)置：=_Numeric_表 */
   
6. /* 示統(tǒng)計全部數(shù)值型變量，=_Character_表示統(tǒng)計全部字符型變量，=_All_表 */
   
7. /* 示統(tǒng)計全部變量。 */
   
8. /* */
   
9. /* 最終結(jié)果表格中包含所有指定變量的名稱、類型和相應(yīng)的缺失值數(shù)量。 */
   
10. /* */
    
11. /* Created on 2013.5.8 */
    
12. /* Modified on 2013.5.8 */
    
13. /**********************************************************************/
    
14. 
    
15. /* 檢查TargetTable的存在性 */
    
16. %if &TargetTable EQ %STR() %then %do;
    
17. %put ERROR: The TargetTable should not be blank, please check it again.;
    
18. %goto exit;
    
19. %end;
    
20. 
    
21. /* 檢查InputVar的合法性 */
    
22. %if %UPCASE(&InputVar) NE _NUMERIC_ AND %UPCASE(&InputVar) NE _CHARACTER_ AND %UPCASE(&InputVar) NE _ALL_ %then %do;
    
23. %ChkVar(SourceTable=&SourceTable,InputVar=&InputVar,FlagVarExists=GMN_FlagVarExists);
    
24. 
    
25. %if %SYSFUNC(FIND(&GMN_FlagVarExists,0)) NE 0 %then %do;
    
26. %put ERROR: The InputVar should be _Numeric_, _Character_, _All_ or any variable name in SourceTable, case insensitive and without quotes.;
    
27. %put ERROR: The InputVar "%SCAN(&InputVar,%SYSFUNC(FIND(&GMN_FlagVarExists,0)))" does not exist in SourceTable, please check it again.;
    
28. %goto exit;
    
29. %end;
    
30. %end;
    
31. 
    
32. %if %UPCASE(&InputVar) EQ _NUMERIC_ %then %do;
    
33. proc contents data=&SourceTable position out=GMN_VarList(keep=name type varnum) noprint;
    
34. run;
    
35. 
    
36. /* 檢查數(shù)值型變量是否存在 */
    
37. %ChkValue(SourceTable=GMN_VarList,
    
38. InputVar=type,
    
39. Value=1,
    
40. FlagValueExists=GMN_FlagNumVarExists);
    
41. 
    
42. %if &GMN_FlagNumVarExists GT 0 %then %do;
    
43. proc sql noprint;
    
44. select NAME into :GMN_NumVarList separated by ' '
    
45. from GMN_VarList
    
46. where TYPE EQ 1;
    
47. quit;
    
48. 
    
49. /* 化簡原始表格 */
    
50. data &TargetTable(keep=&GMN_NumVarList drop=GMN_i);
    
51. set &SourceTable;
    
52. array VarList &GMN_NumVarList;
    
53. do GMN_i=1 to dim(VarList);
    
54. if VarList{GMN_i} NE . then VarList{GMN_i}=0;
    
55. else VarList{GMN_i}=1;
    
56. end;
    
57. run;
    
58. 
    
59. %GetStatsForTable(SourceTable=&TargetTable,
    
60. TargetTable=&TargetTable,
    
61. ByFactors=,
    
62. InputVar=&GMN_NumVarList,
    
63. InputVarType=,
    
64. OutputVarType=,
    
65. Weight=,
    
66. Statistic=SUM);
    
67. 
    
68. proc transpose data=&TargetTable out=&TargetTable;
    
69. var &GMN_NumVarList;
    
70. run;
    
71. 
    
72. data &TargetTable(rename=(_NAME_=VarName _LABEL_=VarLabel COL1=MissNum));
    
73. retain _NAME_ _LABEL_ VarType COL1;
    
74. set &TargetTable;
    
75. VarType='N';
    
76. run;
    
77. %end;
    
78. %else %do;
    
79. %put ERROR: There is no numeric variable existed in SourceTable, please check it again.;
    
80. %goto exit;
    
81. %end;
    
82. %end;
    
83. %else %if %UPCASE(&InputVar) EQ _CHARACTER_ %then %do;
    
84. proc contents data=&SourceTable position out=GMN_VarList(keep=name type varnum) noprint;
    
85. run;
    
86. 
    
87. /* 檢查字符型變量是否存在 */
    
88. %ChkValue(SourceTable=GMN_VarList,
    
89. InputVar=type,
    
90. Value=2,
    
91. FlagValueExists=GMN_FlagCharVarExists);
    
92. 
    
93. %if &GMN_FlagCharVarExists GT 0 %then %do;
    
94. proc sql noprint;
    
95. select NAME,STRIP(NAME)||'_Temp' into :GMN_CharVarList separated by ' ',:GMN_CharVarList_Temp separated by ' '
    
96. from GMN_VarList
    
97. where TYPE EQ 2;
    
98. quit;
    
99. 
    
100. /* 化簡原始表格 */
     
101. data &TargetTable(keep=&GMN_CharVarList_Temp drop=GMN_j);
     
102. set &SourceTable;
     
103. array VarList &GMN_CharVarList;
     
104. array VarList_Temp &GMN_CharVarList_Temp;
     
105. do GMN_j=1 to dim(VarList);
     
106. if VarList{GMN_j} NE "" then VarList_Temp{GMN_j}=0;
     
107. else VarList_Temp{GMN_j}=1;
     
108. end;
     
109. run;
     
110. 
     
111. %GetStatsForTable(SourceTable=&TargetTable,
     
112. TargetTable=&TargetTable,
     
113. ByFactors=,
     
114. InputVar=&GMN_CharVarList_Temp,
     
115. InputVarType=,
     
116. OutputVarType=,
     
117. Weight=,
     
118. Statistic=SUM);
     
119. 
     
120. proc transpose data=&TargetTable out=&TargetTable;
     
121. var &GMN_CharVarList_Temp;
     
122. run;
     
123. 
     
124. data &TargetTable(rename=(_NAME_=VarName _LABEL_=VarLabel COL1=MissNum));
     
125. retain _NAME_ _LABEL_ VarType COL1;
     
126. set &TargetTable;
     
127. VarType='C';
     
128. _NAME_=SUBSTR(_NAME_,1,LENGTH(_NAME_)-5);
     
129. run;
     
130. %end;
     
131. %else %do;
     
132. %put ERROR: There is no character variable existed in SourceTable, please check it again.;
     
133. %goto exit;
     
134. %end;
     
135. %end;
     
136. %else %do;
     
137. proc contents data=&SourceTable position out=GMN_VarList(keep=name type varnum) noprint;
     
138. run;
     
139. 
     
140. /* 篩選指定的變量 */
     
141. %if %UPCASE(&InputVar) NE _ALL_ %then %do;
     
142. %let InputVar_Comma=%PrxChange(InputString=&InputVar,PrxString=s/(\w+)/'$1'/); /* 給InputVar中的代碼加引號 */
     
143. %let InputVar_Comma=%SYSFUNC(TRANSLATE(&InputVar_Comma,%STR(,),%STR( ))); /* 替換InputVar中的空格為逗號 */
     
144. 
     
145. proc sql noprint;
     
146. create table GMN_VarList as
     
147. select * from GMN_VarList
     
148. where Name in (&InputVar_Comma)
     
149. order by Name;
     
150. quit;
     
151. %end;
     
152. 
     
153. /* 檢查數(shù)值型變量是否存在 */
     
154. %ChkValue(SourceTable=GMN_VarList,
     
155. InputVar=type,
     
156. Value=1,
     
157. FlagValueExists=GMN_FlagNumVarExists);
     
158. 
     
159. /* 檢查字符型變量是否存在 */
     
160. %ChkValue(SourceTable=GMN_VarList,
     
161. InputVar=type,
     
162. Value=2,
     
163. FlagValueExists=GMN_FlagCharVarExists);
     
164. 
     
165. %if &GMN_FlagNumVarExists GT 0 %then %do;
     
166. proc sql noprint;
     
167. select NAME into :GMN_NumVarList separated by ' '
     
168. from GMN_VarList
     
169. where TYPE EQ 1;
     
170. quit;
     
171. 
     
172. /* 化簡原始表格 */
     
173. data GMN_NumMiss(keep=&GMN_NumVarList drop=GMN_k);
     
174. set &SourceTable;
     
175. array VarList &GMN_NumVarList;
     
176. do GMN_k=1 to dim(VarList);
     
177. if VarList{GMN_k} NE . then VarList{GMN_k}=0;
     
178. else VarList{GMN_k}=1;
     
179. end;
     
180. run;
     
181. 
     
182. %GetStatsForTable(SourceTable=GMN_NumMiss,
     
183. TargetTable=GMN_NumMiss,
     
184. ByFactors=,
     
185. InputVar=&GMN_NumVarList,
     
186. InputVarType=,
     
187. OutputVarType=,
     
188. Weight=,
     
189. Statistic=SUM);
     
190. 
     
191. proc transpose data=GMN_NumMiss out=GMN_NumMiss;
     
192. var &GMN_NumVarList;
     
193. run;
     
194. 
     
195. data GMN_NumMiss(rename=(_NAME_=VarName _LABEL_=VarLabel COL1=MissNum));
     
196. retain _NAME_ _LABEL_ VarType COL1;
     
197. set GMN_NumMiss;
     
198. VarType='N';
     
199. run;
     
200. %end;
     
201. %if &GMN_FlagCharVarExists GT 0 %then %do;
     
202. proc sql noprint;
     
203. select NAME,STRIP(NAME)||'_Temp' into :GMN_CharVarList separated by ' ',:GMN_CharVarList_Temp separated by ' '
     
204. from GMN_VarList
     
205. where TYPE EQ 2;
     
206. quit;
     
207. 
     
208. /* 化簡原始表格 */
     
209. data GMN_CharMiss(keep=&GMN_CharVarList_Temp drop=GMN_l);
     
210. set &SourceTable;
     
211. array VarList &GMN_CharVarList;
     
212. array VarList_Temp &GMN_CharVarList_Temp;
     
213. do GMN_l=1 to dim(VarList);
     
214. if VarList{GMN_l} NE "" then VarList_Temp{GMN_l}=0;
     
215. else VarList_Temp{GMN_l}=1;
     
216. end;
     
217. run;
     
218. 
     
219. %GetStatsForTable(SourceTable=GMN_CharMiss,
     
220. TargetTable=GMN_CharMiss,
     
221. ByFactors=,
     
222. InputVar=&GMN_CharVarList_Temp,
     
223. InputVarType=,
     
224. OutputVarType=,
     
225. Weight=,
     
226. Statistic=SUM);
     
227. 
     
228. proc transpose data=GMN_CharMiss out=GMN_CharMiss;
     
229. var &GMN_CharVarList_Temp;
     
230. run;
     
231. 
     
232. data GMN_CharMiss(rename=(_NAME_=VarName _LABEL_=VarLabel COL1=MissNum));
     
233. retain _NAME_ _LABEL_ VarType COL1;
     
234. set GMN_CharMiss;
     
235. VarType='C';
     
236. _NAME_=SUBSTR(_NAME_,1,LENGTH(_NAME_)-5);
     
237. run;
     
238. %end;
     
239. 
     
240. data &TargetTable;
     
241. set
     
242. %if &GMN_FlagNumVarExists GT 0 %then %do;
     
243. GMN_NumMiss
     
244. %end;
     
245. %if &GMN_FlagCharVarExists GT 0 %then %do;
     
246. GMN_CharMiss
     
247. %end;
     
248. ;
     
249. run;
     
250. %end;
     
251. 
     
252. /* 刪除不必要的表格 */
     
253. proc datasets lib=work nolist;
     
254. delete GMN_:;
     
255. quit;
     
256. 
     
257. %exit:
     
258. %mend;
     
259. 
     
260. 
     
261. %macro Demo();
     
262. 
     
263. %let SourceTable=Cars;
     
264. %let TargetTable=MissNum;
     
265. %let InputVar=Cylinders; /* =_Numeric_表示統(tǒng)計全部數(shù)值型變量，=_Character_表示統(tǒng)計全部字符型變量，=_All_表示統(tǒng)計全部變量，大小寫不敏感 */
     
266. %GetMissNum(&SourceTable,&TargetTable,&InputVar);
     
267. 
     
268. %mend;

復(fù)制代碼

掃碼加我 拉你入群

請注明：姓名-公司-職位

以便審核進群資格，未注明則拒絕

分享0 收藏6 回帖

關(guān)鍵詞：sas基礎(chǔ) Miss Get Iss MIS 基礎(chǔ)

我擦，縮進正常了
老天有眼啊

這么長的宏，看著都眼花，想問一下，是怎么寫出來的？

考慮得很嚴謹啊。佩服。

Imasasor 發(fā)表于 2013-4-16 15:46
這么長的宏，看著都眼花，想問一下，是怎么寫出來的？

人肉寫的唄
只不過只要寫一次，以后就都方便了
這個宏用于檢查回歸中自變量的缺失數(shù)量
我發(fā)現(xiàn)proc reg里都沒有相應(yīng)的控制選項
所以只能人肉做

拜讀，收藏，學(xué)習(xí)，感謝。
人肉做的，了不起。

學(xué)習(xí)，學(xué)習(xí)，
這里還是有個疑問，如果數(shù)據(jù)行數(shù)太多，做transpose會不會受到影響？

00810112 發(fā)表于 2013-4-26 10:41
學(xué)習(xí)，學(xué)習(xí)，
這里還是有個疑問，如果數(shù)據(jù)行數(shù)太多，做transpose會不會受到影響？

會有影響
但是我想不到其他好辦法了

另外為了加快transpose時的速度
我在之前把原始表格化簡了下
即有值的換為'a'，沒值的換為''
這樣轉(zhuǎn)置的時候快多了

en ,

謝謝分享