SAS編程技術(shù)數(shù)據(jù)步循環(huán)與轉(zhuǎn)移控制課件

1、第8章 數(shù)據(jù)步循環(huán)與轉(zhuǎn)移控制清華大學經(jīng)管學院 朱世武ZhushwResdat樣本數(shù)據(jù)：SAS論壇： DO語句 DO語句必須由一個END語句來結(jié)束。 DO和END語句之間的這些語句稱為一個DO組。 DO組可以嵌套任意次。DO語句有5類： 簡單DO語句； 循環(huán)DO語句； DO OVER語句； DO WHILE語句； DO UNTIL語句。簡單DO語句 語句格式DO； 一些SAS語句END;應用舉例 例8.1 IF/THEN語句用中DO.data a;set ResDat.class;if age14 then do;h_cm=30.5*height/12;put name= sex= age= h

2、_cm=;end;run;329 data a;330 set fdata.class;331 if age14 then do;332 h_cm=30.5*height/12;333 put name= sex= age= h_cm=;334 end;335 run;Name=Janet Sex=F Age=15 h_cm=158.85416667Name=Mary Sex=F Age=15 h_cm=169.02083333Name=Philip Sex=M Age=16 h_cm=183Name=Ronald Sex=M Age=15 h_cm=170.29166667Name=Will

3、iam Sex=M Age=15 h_cm=169.02083333NOTE: 從數(shù)據(jù)集 FDATA.CLASS 讀取了 19 個觀測。NOTE: 數(shù)據(jù)集 WORK.A 有 19 個觀測和 6 個變量。NOTE: “DATA 語句”所用時間（總處理時間）: 實際時間 0.01 秒 CPU 時間 0.01 秒上例與下面的if語句效果一樣，但是IF語句效率較低。If age14 then H_cm=30.5*height/12;If age14 then Put name= sex= age= h_cm=; 循環(huán)DO語句 DO組中的語句需要被重復執(zhí)行時要用到循環(huán)DO語句。循環(huán)DO語句根據(jù)下標變量重

4、復執(zhí)行DO和END語句之間的語句 語句格式DO index_variable=specification-1;一些SAS語句END;其中：index-variable選項定義下標變量，用于控制DO組的執(zhí)行方式和重復次數(shù)。下標變量包含在創(chuàng)建的數(shù)據(jù)集中，可以用DROP語句刪除。 SPECIFICATION選項格式 Start ; Specification選項說明： 例8.2 缺省的步長為1。do I=1 to 1000;do I=1 to y+3;do I=1 to exit;例8.3 起始值例句。do i=5;do i=1 to n ;do i=n to 1 by 1;do i=k+1 to

5、n-1;do i=1 to k-1, k+1 to n;do i=2, 3, 5, 7, 11, 13, 17;do i=0.1 to 0.9 by 0.1, 1 to 10 by 1, 20 to 100 by 10;do i=saturday, sunday;do i=01jan99d, 25feb99d;do i=01jan99d to 01jan2000d by 1;例句子中，起始值必須全部是數(shù)值或全部是字符常數(shù)，也可以是變量。字符常數(shù)必須用引號括起來。 例8.4 使用表達式選項。do i=1 to 10 while(xy);do i=10 to 0 by 1 while(month=

6、jan);應用舉例 例8.5 用GO TO語句跳出循環(huán)體。data a;input x y; if xy then goto skip; /*skip是循環(huán)體外的語句標號*/y=log(y-x);yy=y-20;skip: if yy then goto skip; /*skip是循環(huán)體外的語句標號*/339 y=log(y-x);340 yy=y-20;341 skip: if y=5時停止。data a;n=0;do until(n=5);put n=;n+1;end;run; END語句 END語句是組成DO組或SELECT組的最后一個語句。 語句格式END; 應用舉例 DO;若干個SA

7、S語句END；SELECT (expression); When (expression) SAS語句； Otherwise SAS語句； END; SELECT語句 SELECT語句選擇執(zhí)行SAS語句。 語句格式SELECT; WHEN (when-expression-1) statement; WHEN(when-expression-1) statement; END; 選項說明： 應用舉例 例8.13 有SELECT表達式的SELECT語句。data a;set ResDat.st_list;obs=_n_;x=uniform(0);select (obs);when (1) x=x

8、*10;when (2,4,6);when (3,5,7,13,15,17) x=x*100;otherwise x=1;end;run;例中，根據(jù)OBS來計算新變量X的值。當OBS為1時，執(zhí)行x=x*10。當OBS為2，4，6時，X值不變。當OBS為3，5，7，13，15，17時，執(zhí)行X=X*100。OBS為其它值時，執(zhí)行x=1。 例8.14 沒有SELECT表達式的SELECT語句。data _null_;do mon=jan,feb,mar,apr,may,jun,jul,aug,sep;select; when (mon in (jun,jul,aug) put summer mon=

9、;when (mon in (mar,apr,may) put spring mon=;otherwise put fall or winter mon=;end;end;run; 361 data _null_;362 do mon=jan,feb,mar,apr,may,jun,jul,aug,sep;363 select;364 when (mon in (jun,jul,aug) put summer mon=;365 when (mon in (mar,apr,may) put spring mon=;366 otherwise put fall or winter mon=;367

10、 end;368 end;369 run;fall or wintermon=janfall or wintermon=febspring mon=marspring mon=aprspring mon=maysummer mon=junsummer mon=julsummer mon=augfall or wintermon=sepNOTE: “DATA 語句”所用時間（總處理時間）: 實際時間 0.03 秒 CPU 時間 0.01 秒例8.15 SELECT組的正確及錯誤用法舉例。select(x);when (x=2) put two; /*錯誤用法*/end;select(x);whe

11、n(2) put two; /*正確用法*/end;select;when (x=2) put two; /*正確用法*/end;run;例中，后兩個程序的效果相同。 IF語句IF語句的兩種類型：條件IF語句，含有一個THEN子句或ELSE子句；子集IF語句，沒有子句。IF-THEN與IF-THEN/ELSE語句表達式為真時執(zhí)行THEN后面的語句, 表達式為假執(zhí)行ELSE后面的語句。語句格式：IF expression THEN expression;例8.16 應用舉例。data;if x then delete;if status=OK and type=3 then count+1;if

12、 age ne agecheck then delete;if x=0 thenif y ne 0 then put x zero, y nozero;else put x zero, y zero;else put x nozero;if ans=9 thendo; ans=. ;put invalid answer for id=;end;else do; /*當IF條件不成立時，執(zhí)行下面另一DO組 */ans=ans1;valid+1;end;子集IF語句 語句格式：IF expression;其中：expression是任意有效的表達式。如果表達式真的，SAS語句對正被創(chuàng)建的觀測繼續(xù)執(zhí)

13、行DATA步的語句。如果表達式是假的，SAS立即返回到DATA步的開始對其他觀測執(zhí)行。 例8.17 應用舉例。data a;set ResDat.Idx000001;if _n_100;run;例中，選擇前99條觀測值。GO TO語句與語句標號 GO TO語句使SAS跳到本程序步帶有標號的語句繼續(xù)往下執(zhí)行。GO TO語句和指定的目標語句必須在同一DATA步。 GO TO語句語句格式：GO TO label; GOTO label;Label選項規(guī)定語句標號來指示GO TO的目標。 例8.18 應用舉例。data a;input x ;if 1=x=5 then goto ok;put x; c

14、ount+1; ok: sumx+x;cards;1 2 7 2 12 24 22;例中，COUNT僅對大于5的數(shù)計數(shù)，SUMX+X對每個觀測都執(zhí)行。 xcountsumx110122033711042112512224624348722470data a;input x ;if 1=x=5 then goto ok;put x; count+1;return;ok: sumx+x;cards;1 2 7 2 12 24 22;run;例中，COUNT僅對大于5的數(shù)計數(shù)，SUMX+X僅對X在1和5之間的那些觀測執(zhí)行。 xcountsumx1101220337134215512256243572

15、245例8.19 上面兩例中用DO-END語句替代GOTO語句。data a;input x ;if x5 then do;put x; count+1;end;sum+x;cards;1 2 7 2 12 24 22;data a;input x ;if x5 then do;put x; count+1;end;else sumx+x;cards;1 2 7 2 12 24 22;語句標號 語句標號給出GO TO，LINK, FILE語句中的選項HEADER=，或在INFILE語句中選項EOF=的位置。語句格式：LABEL: statement; LINK語句 LINK語句告訴SAS系統(tǒng)立

16、即轉(zhuǎn)到由LINK語句指示的語句標號，并從那里繼續(xù)執(zhí)行直到一個RETURN語句。RETURN語句讓SAS立即返回到LINK語句后面的那個語句并從那里繼續(xù)執(zhí)行。LINK語句和語句標號所在的目標語句必須在同一個DATA步。 語句格式LINK label; 例8.20 執(zhí)行一組語句。 data bydro;input type $ wd station $;label type=station type wd=depth to water;elev=.;if type=aluy then link calcu;year=1985;return;calcu: if station=site_1 then

17、 elev=6650-wd; if station=site_2 then elev=5500-wd;return; /*返回到link語句下面的year=1985;語句*/cards;aluv 523 site_1uppa 234 site_2aluy 666 site_2;typewdstationelevyear1aluv523site_1.19852uppa234site_2.19853aluy666site_248341985RETURN語句RETURN語句告訴SAS系統(tǒng)在DATA步當前位置上停止執(zhí)行語句，返回到一個預定位置上繼續(xù)執(zhí)行。語句格式RETURN;RETURN語句的作用：在

18、LINK語句后面，RETURN語句讓SAS系統(tǒng)立即返回到LINK后面的語句，并繼續(xù)往下執(zhí)行；在FILE語句的選項HEADER=中，RETURN語句讓SAS系統(tǒng)立即返回到開始新頁之前執(zhí)行的最后一個語句后面的語句，并繼續(xù)執(zhí)行；在DATA步的其它地方，RETURN語句讓SAS系統(tǒng)返回到DATA步開頭。例8.23 讓SAS系統(tǒng)返回到DATA步開頭。data survey;input x y z;if x=y then return;x=y+z; a=x*2;cards;1 2 3 3 3 45 6 78 8 9;例中，當X=Y時，RETURN語句被執(zhí)行。SAS系統(tǒng)添加這個觀測到數(shù)據(jù)集SURVEY，并返

19、回到DATA步的開頭。IF語句后面的二個賦值語句沒有被執(zhí)行。 xyza1523252334 .313671694889 .CONTINUE語句 CONTINUE語句停止當前的這次循環(huán)過程，繼續(xù)進行下一次循環(huán)。語句格式：CONTINUE; 例8.25 CONTINUE語句停止當前的循環(huán)過程，繼續(xù)進行下一次循環(huán)。data a;do i=1 to 5;input name $ idno status $;if status=PT then continue; /*返回到下次循環(huán)的開頭*/input benefits $10.;output;end;cards;/*接上頁 */Jones 9011 P

20、TThomas 876 PTRichards 1002 FTEye/DentalKelly 85111 PTSmith 433 FTHMOKelly 85111 PTSmith 433 FTHMOJones 9011 PTThomas 876 PTRichards 1002 FTEye/Dental;Run; inameidnostatusbenefits13Richards1002FTEye/Dental25Smith433FTHMO32Smith433FTHMO45Richards1002FTEye/DentalLEAVE 語句 LEAVE 語句停止當前DO組循環(huán)或SELECT組的處理過程

21、，并繼續(xù)執(zhí)行DO組或SELECT組后面的語句。語句格式：LEAVE; 例8.26 給定條件下停止DO循環(huán)的處理過程data week;input name $ idno start status $ dept $;bonus=0;do year=start to 1991;if bonus ge 500 then leave;bonus+50;end;cards;/*接上頁 */Jones 9011 1990 PT PUBThomas 876 1976 PT HRBarnes 7899 1991 ft techHarrell 1250 1975 ft hrRichards 1002 1990

22、ft devKelly 85 1981 pt pubStone 091 1990 pt mait;Run;例中，首先由IF語句檢查BONUS的值，當BONUS的值超過允許的最大值500時，用LEAVE語句來停止DO循環(huán)的處理過程nameidnostartstatusdeptbonusyear1Jones90111990PTPUB10019922Thomas8761976PTHR50019863Barnes78991991fttech5019924Harrell12501975fthr50019855Richards10021990ftdev10019926Kelly851981ptpub500

23、19917Stone911990ptmait1001992第13章 變量輸入輸出格式清華大學經(jīng)管學院 朱世武ZhushwResdat樣本數(shù)據(jù)：SAS論壇： 本章內(nèi)容包括： 輸入格式； 輸出格式； SAS日期時間存貯方式； 日期時間輸入格式； 日期時間輸出格式； 缺失值處理。輸入格式 輸入格式是SAS系統(tǒng)用來讀入數(shù)據(jù)值的一個指令。輸入格式形式 INFORMAT .選項說明：通用規(guī)則： 所有輸入格式必須包含一個點(.)作為名字的一部分； 對于省略W和D值的輸入格式，使用系統(tǒng)缺省值。輸入格式使用方法 INPUT語句； INPUT函數(shù)； DATA步中用INFORMAT或ATTRIB語句； PROC步中

24、用INFORMAT或ATTRIB語句。INPUT語句 例13.1 列格式輸入。input 6 date ddmmyy10. +1 stocd $ 8. oppr 8.2;例中，從第6列開始以DDMMYY10.格式讀入變量DATE, 第7列開始以長度為8的字符格式讀入變量STOCD,然后以數(shù)值格式8.2讀入變量OPPR。INPUT函數(shù) 例13.2 字符轉(zhuǎn)換數(shù)值。data;x=98.6;y=input(x, 4.);put x= y=;run;例中，INPUT函數(shù)和w.d輸入格式將X的字符值轉(zhuǎn)換為一個數(shù)值Y.INPUT函數(shù)對于數(shù)據(jù)的轉(zhuǎn)換是很有用的，當然，有更簡單的方法將字符值轉(zhuǎn)換為數(shù)值。data;

25、x=98.6;y=x-0;put x= y=;run; 例13.3 數(shù)值轉(zhuǎn)換為字符時會產(chǎn)生不正確結(jié)果。data;x=2557898;y=input(x,$8.);put y;run;例中，產(chǎn)生的結(jié)果為255, 不正確。例13.4 PUT函數(shù)將數(shù)值轉(zhuǎn)換為字符。data;x=2557898;y=put(x, $8.);put y;run;例中，用PUT函數(shù)得到正確結(jié)果，即字符型2557898.data;x=2557898;y=put(x, $3.);put y;run;例中，Y的結(jié)果為3E6.在DATA步中用INFORMAT語句或ATTRIB語句 例13.5 DATA步中用INFORMAT語句。d

26、ata a;input x;informat x comma12.4;cards;12,345.123422,345.1234;run;例中，最好就是只用INFORMAT語句，不要用ATTRIB，這樣可以少記些東西。 永久性與臨時性聯(lián)系 DATA步規(guī)定的輸入格式是永久聯(lián)系的； PROC步規(guī)定的輸入格式是臨時聯(lián)系的。 輸入格式類型 數(shù)值輸入格式； 字符輸入格式； 日期時間輸入格式； 豎式二進制數(shù)輸入格式； 使用FORMAT過程自定義的輸入格式。 w.d輸入格式 讀入標準數(shù)值數(shù)據(jù)。例13.6 INPUT語句對變量使用列指針控制和w.d的輸入格式。INPUT X 8.2;INPUT X 1-10 .

27、2;W.D輸入格式把一個點(.)作為缺失值。用W.D輸入格式時，數(shù)值域尾部的空白不表示0，用BZ.輸入格式可以把尾部空白讀為0.數(shù)值變量輸入格式COMMAw.d 輸入格式 COMMAw.d輸入格式讀入包含字符的數(shù)值。數(shù)值數(shù)據(jù)內(nèi)可能包含的字符有： 逗號; 空格; $; %; 破折號; 圓括號等。例13.8 讀入帶千分號的數(shù)據(jù)。data a;input x comma12.2;cards;12,522.2315;run;例13.10 數(shù)值變量輸入格式應用舉例。data ;informat x 8.2;x=12345.1234;put x=; run;data a;x=12345.1234;info

28、rmat x 12.4;format x 8.2;put x=;run;data a;input x;informat x comma12.4;cards;12,345.123422,345.1234;run;結(jié)果顯示x=12345.1234.結(jié)果顯示x=12345.12.$w.輸入格式 讀入標準字符數(shù)據(jù)。在讀入字符值之前，$w.輸入格式清除字符值開頭的空格。$w.輸入格式把一個點(.)轉(zhuǎn)換為空格，因為它把一個點看作缺失值。例13.11 讀入數(shù)據(jù)時，清除數(shù)據(jù)開頭空格并將（.）轉(zhuǎn)換為缺失值。data;input name $5.;cards;xyz.uvw;options nocenter;pr

29、oc print noobs;run;輸出結(jié)果：xyzuvw字符變量輸入格式 $CHARw輸入格式 讀含有空格的字符數(shù)據(jù)。$CHARw.輸入格式除了不清除字符值開頭的空作格外，它等同于標準的$w.輸入格式。例13.12 保留開頭和結(jié)尾的空格。data;input name $char10.;cards;xyz .uvw;options nocenter;proc print noobs;run;INPUT語句采取自由格式輸入時，INFORMAT或ATTRIB語句中不能使用$CHAR.輸入格式，因為SAS把空格看作數(shù)據(jù)行中數(shù)值間的分隔符。輸出結(jié)果：xyz.uvw例13.13 不能放放在一起使用的

30、語句。informat x $char12.;input x y z;例13.14直接賦值為字符型變量。 data a; x1=st92;x2=st92;informat x $8. ;x3=12,345.1234; run;輸出格式 輸出格式是SAS系統(tǒng)用來輸出數(shù)據(jù)值的一個指令。 輸出格式形式 FORMAT .選項說明： 通用規(guī)則： 所有輸出格式必須包含一個點(.)作為名字的一部分； 省略W和D的值時，使用系統(tǒng)的缺省值； 無論怎樣規(guī)定輸出格式中的小數(shù)位，輸出格式都不會影響存貯的數(shù)據(jù)值； 規(guī)定的輸出格式寬度太窄小時，對字符格式截去右邊的字符，對數(shù)值格式轉(zhuǎn)換為BESTw.的格式； 使用一個不協(xié)調(diào)

31、的輸出格式時，SAS系統(tǒng)首先試著使用其它類型的類似格式。如果行不通，將輸出一個錯誤信息在SAS日志。輸出格式使用方法 PUT語句； PUT函數(shù)； DATA步中用FORMAT或ATTRIB語句； PROC步中用FORMAT或ATTRIB語句。PUT語句 例13.15 PUT語句中使用輸出格式。data;x=1145.32;put x dollar10.2;run;結(jié)果為$1，145.32。輸出格式帶千分號和美元號。PUT函數(shù) PUT函數(shù)對于將數(shù)值轉(zhuǎn)換為字符，或者改變變量的字符輸出格式是很有用的。例13.16 PUT函數(shù)中使用輸入格式。data;cc=16;cchex=put(cc,hex3.);

32、put cc hex3.;run;結(jié)果為010.例中，轉(zhuǎn)換數(shù)值變量cc的值為三個字符的十六進制表達式。data;cc=16;cc=put(cc,3.);run;例中，將數(shù)值CC轉(zhuǎn)換為字符值。DATA步中用FORMAT或ATTRIB語句 例13.17 DATA步中用FORMAT或ATTRIB語句。format sales1-sales12 comma10.2;例中，F(xiàn)ORMAT語句對變量sales1至sales12規(guī)定輸出格式COMMAw.d.attrib sales1-sales12 format=comma10.2;例中，ATTRIB語句對變量sales1至sales12規(guī)定輸出格式COMM

33、Aw.d.輸出格式類型 數(shù)值輸出格式; 字符輸出格式; 日期時間輸出格式; 使用FORMAT過程創(chuàng)建的自定義輸出格式。w.d輸出格式 輸出標準數(shù)值數(shù)據(jù)。對超出輸出格式寬度的數(shù)值進行四舍五入。數(shù)值過大時，使用BESTw.輸出。例13.18 w.d輸出格式應用。data a; x=23.45; put x 6.3;run;結(jié)果是23.450。例13.19 使用列輸出法等價于使用指針控制和w.d的格式。put x 1-8 .2;put 1 x 8.2;例中，兩個語句等價。數(shù)值變量輸出格式BESTw.格式 沒有對變量規(guī)定輸出格式時，使用BEST.的格式。SAS在規(guī)定的寬度列數(shù)中選擇能夠給出該值最多信息

34、的表示法作為輸出格式。SAS仍存貯原來的完整數(shù)值。例13.20 用不同的BESTw.形式輸出時，顯示結(jié)果不一樣，但在SAS存貯的都是原來的完整數(shù)值。data;x=12570000;put 10 x best6.;run;輸出值為1.26E6。正確地輸出這個值需要8列，使用E表示法把這個值壓縮為6列.data;x=12570000;put 10 x best3.;run; 輸出值為1E6.data;x=12570000;put 10 x best2.;run; 輸出值為*.Ew.輸出格式 例13.22 用科學記數(shù)法表示輸出數(shù)值。data;x=1257;put 10 x e10.;run;輸出值為

35、1.257E+03，共占9列。輸入輸出控制流程 輸入-存儲-輸出控制流程 $w.輸出格式 $w.是最常用的字符數(shù)據(jù)輸出格，w值給出輸出字符值的列數(shù)。例13.24 $w.和列輸出格式產(chǎn)生同樣效果。data ;name=ABC;put 10 name $3.;put name $10-12;run;例中，都輸出值為ABC.字符變量輸出格式 日期時間存貯方式 SAS日期值存貯方式SAS系統(tǒng)存貯日期值為1960年1月1日到這個日期之間的天數(shù)。日期存貯方式舉例例13.25 表中天數(shù)的計算程序。data;date=mdy(1,1,1960);put date 8.;date=mdy(3,3,1962);p

36、ut date 8.;date=mdy(8,4,1985);put date 8.;date=mdy(10,1,2002);put date 8.;date=mdy(7,4,1776);put date 8.;run;SAS日期時間值存貯方式 SAS存貯時間值從0時開始以秒記數(shù)。日期時間值存貯方式舉例例13.26 表中秒數(shù)的計算程序。data;hms=hms (21,50,51); put hms=;x=21*60*60+50*60+51;put x=; date=mdy(10,1,2002);put date=;dhms=dhms (date,21,50,51); put dhms=;y=date*24*3600+x;put y=;run;日期時間輸入格式 輸入格式描述例子寬度范圍缺省范圍datew.Ddmmyy1jan200301jan031-jan-20037-327datetimew.Ddmmyy hh:mm:ss.ss01jan03:8:56:10.201jan2003/8:56:1013-4018yymmddw.03 01 0103/01/01030101200301016-326應用舉例 例13 27 使用SAS日期時間的輸入格式。data a;input date mmd