常用SAS函數(shù)簡介

1、常用SAS函數(shù)簡介阮紅蓮廣州醫(yī)科大學(xué) 統(tǒng)計(jì)學(xué)系SAS函數(shù) SAS系統(tǒng)提供的SAS函數(shù)比一般高級語言的標(biāo)準(zhǔn)函數(shù)多得多，尤其是概率函數(shù)、分位數(shù)函數(shù)、樣本統(tǒng)計(jì)函數(shù)、隨機(jī)數(shù)函數(shù)。 如此豐富的SAS函數(shù),對用戶統(tǒng)計(jì)分析計(jì)算以及編寫SAS程序帶來極大的方便。SAS函數(shù)定義 SAS函數(shù)是一個(gè)子程序，對0個(gè)或幾個(gè)參數(shù)（自變量）返回一個(gè)結(jié)果值。SAS函數(shù)的形式：函數(shù)名（函數(shù)名（參數(shù)1，參數(shù)2，）例如： MEAN(x,y,z,a)MIN(SUM(x,y,z),1000)SAS函數(shù)中的參數(shù)1）參數(shù)類型：變量名；常數(shù)；函數(shù)；表達(dá)式。例1 函數(shù)及其參數(shù)類型舉例。 2）參數(shù)個(gè)數(shù)： 有些函數(shù)不需要自變量，有些只需要一個(gè)自變

2、量，有些函數(shù)有多個(gè)自變量。3）參數(shù)表示法：當(dāng)函數(shù)有多個(gè)參數(shù)時(shí)，必須用逗號分隔開。4）參數(shù)名縮寫方法： 函數(shù)名(OF 變量名1 變量名n ) 參數(shù)的表示法舉例以逗號(，)分隔參數(shù),如 sum(x1, x2, x3, x4, x5, x6, x7, x8)去掉逗號(，)，最前面加上“of”,如 sum(of x1 x2 x3 x4 x5 x6 x7 x8)去掉逗號(，)，最前面加上“of”,變量簡寫：sum(of x1 - x8)例2 正確與錯(cuò)誤的表示方法例3 數(shù)組元素自變量的縮寫方法array y10 y1-y10; /*定義數(shù)組y10 */x=sum (of y*);z=sum (of y1-

3、y10);例中，兩個(gè)SUM函數(shù)的表示法是等價(jià)的。 正確的表示方法錯(cuò)誤的表示方法sum (x, y, z)sum (x y z)sum (of x1-x10)sum (x1-x10)sum (of x1-x10 d y1-y100)sum (x1-x10 d y1-y100)函數(shù)結(jié)果通常函數(shù)的結(jié)果由其自變量的屬性決定： 自變量是字符時(shí)結(jié)果變量為字符； 自變量是數(shù)值時(shí)結(jié)果變量為數(shù)值。PUT函數(shù)是一個(gè)例外，它不管自變量是什么類型，其結(jié)果總是字符值。用PUT語句在LOG窗口顯示函數(shù)值例4 顯示概率值和分位數(shù)。data;Y=probnorm(1.96); /* 標(biāo)準(zhǔn)正態(tài)分布Z小于1.96的概率 */pu

4、t Y;q1=tinv(.95, 2); /*自由度為2的t分布的0.95分位數(shù)*/q2=tinv(.95, 2, 3); /*自由為2，非中心參數(shù)為3的t分布的0.95分位數(shù)*/put q1= q2=;run;0.9750021049q1=2.9199855804 q2=13.894376071利用HELP查找SAS系統(tǒng)提供的各種函數(shù)索引中鍵入：SAS FunctionHelp下拉菜單SAS 幫助與文件SAS函數(shù)種類在SAS系統(tǒng)9.1版本,SAS函數(shù)分為二十五種類型,分別是算術(shù)函數(shù)、截取函數(shù)、數(shù)學(xué)函數(shù)、概率和密度函數(shù)、分位數(shù)函數(shù)、非中心函數(shù)、樣本統(tǒng)計(jì)函數(shù)、隨機(jī)數(shù)函數(shù)、財(cái)政金融函數(shù)、日期和時(shí)間

5、函數(shù)、換算函數(shù)和特殊函數(shù)等。在SAS系統(tǒng)6.12版本，SAS函數(shù)分為十七種類型共有178個(gè),SAS系統(tǒng)8.1版本，SAS函數(shù)分為二十二種類型共有319個(gè)函數(shù)。1.日期時(shí)間函數(shù)簡介例1.1：計(jì)算兩個(gè)日期之間相隔的天數(shù)data _null_;sdate=01jan2012d;edate=01jan2013d;actual=datdif(sdate, edate, act/act); days360=datdif(sdate, edate, 30/360); put actual= days360=; run;data _null_;actual=datdif(01jan2002d, 01jan20

6、03d, act/act);days360=datdif(01jan2002d, 01jan2003d, 30/360);put actual= days360=; run;例1.2: 以日為單位計(jì)算當(dāng)前日期的天數(shù)。data;x=date (); y=today();format x y yymmdd10.;put x= y=;run;例1.3: 將時(shí)間值換算為以秒計(jì)數(shù)的計(jì)算方法。data;hms=hms (21,50,51); put hms=;x=21*60*60+50*60+51; put x=;mdy=mdy (08,18,2001); put mdy=;format mdy yymm

7、dd10.;run;練習(xí)1計(jì)算以下數(shù)據(jù)表中日期值所在的年季月以及處于某月的第幾天、某周的第幾天。當(dāng)日期值date=14JUL2007d時(shí),試問這一天是星期幾? Date1991-01-021991-01-03 1991-01-04 1991-01-05 1991-01-07 1991-01-08 1991-01-09 1991-01-10 1991-01-11 1991-01-12解答：data a;Input date YYMMDD10.;Cards;1991-01-021991-01-03 1991-01-04 1991-01-05 1991-01-07 1991-01-08 1991-0

8、1-09 1991-01-10 1991-01-11 1991-01-12;Proc print;Format date YYMMDD10.;run;Data b;Set a;year=year (date);qtr=qtr (date);month=month (date);day=day (date);Weekday=weekday(date);proc print noobs; format date YYMMDD10.;Run;解答：data _null_;weekday=weekday(12JUL2007d);put weekday=; run;2. 數(shù)學(xué)函數(shù)ABS (x) ：求x的

9、絕對值Mod (x,y) ：求x/y的余數(shù)MAX (x1,x2,x3,xn）：求最大值 MIN (x1,x2,x3,xn）：求最小值SQRT(x)：求x的算術(shù)平方根 LOG(x)：求x 的自然對數(shù)EXP(x) ：求指數(shù)函數(shù)exSIN(x)、COS(x)、TAN(x)：正弦、余弦、正切函數(shù).例2：求-314的絕對值，31除以4的商的余數(shù)以及-314 、-87、0、1127這4個(gè)數(shù)的最大值data _null_;x=abs(-314);y=mod(31,4);z=max(-314 ,-87,0,1127);put x= y= z=; run;練習(xí)2ex2.1 計(jì)算0、1.57、 3.14對應(yīng)的正弦

10、函數(shù)值，并將結(jié)果顯示在輸出窗口。ex2.2 計(jì)算e的0次冪,并計(jì)算1的自然對數(shù)。解答：Data;a=sin(0);b=sin(1.57);c=sin(3.14);File print;Put a= 6.2 b= 6.2 c= 6.2;Run;Data;a=sin(0);b=sin(1.57);c=sin(3.14);proc print;format a 6.2 b 6.2 c 6.2;Run;Data;X=exp(0);Y=log(1);Put x y;Run;3. 數(shù)組函數(shù)數(shù)組函數(shù)應(yīng)用舉例data;array mult5,10,2 mult1-mult100;a=DIM1(MULT);b=

11、DIM2(MULT);c=DIM3(MULT);put a b c;run;例3：產(chǎn)生一個(gè)5102的三維數(shù)組，并提取各維度上元素的個(gè)數(shù)。練習(xí) 3產(chǎn)生一個(gè)行變量名為2:6，列變量名為4:13的二維數(shù)組，輸出各維度上的上界和下界。解答：data;array mult2:6,4:13 mult1-mult50;a=HBOUND(MULT);b=HBOUND2(MULT);c=lBOUND(MULT);d=lBOUND(MULT,2);put a b c d;run;4. 樣本統(tǒng)計(jì)量函數(shù)簡單統(tǒng)計(jì)量函數(shù)應(yīng)用舉例例4：請求出2、4、6、3、1這五個(gè)數(shù)的最大值、最小值、算術(shù)平均值、標(biāo)準(zhǔn)差、變異系數(shù)。data

12、 test4;input x1-x5;cards;2 4 6 3 1;run;data simple_statistics;set test1;max_x=max(of x1-x5);min_x=min(of x1-x5);mean_x=mean(of x1-x5);std_x=std(of x1-x5);cv_x=cv(of x1-x5);run;ods html;proc print data=simple_statistics;var max_x min_x mean_x std_x cv_x;run;ods html close;用簡便的方法實(shí)現(xiàn)例1中目的data simple_sta

13、tistics;max_x=max(2,4,6,3,1);min_x=min(2,4,6,3,1);mean_x=mean(2,4,6,3,1);std_x=std(2,4,6,3,1);cv_x=cv(2,4,6,3,1);run;ods html;proc print data=simple_statistics;var max_x min_x mean_x std_x cv_x;run;ods html close;練習(xí)4某農(nóng)村地區(qū)1998年隨機(jī)抽查測量10名13歲女孩的身高，資料如下，試求全距、均值、標(biāo)準(zhǔn)差、標(biāo)準(zhǔn)誤、變異系數(shù)。 145.1 143.3 153.2 152.3 150.6

14、 148.1 144.0 137.1 158.2 146.7解答：data aa;input x1-x10;cards;145.1 143.3 153.2 152.3 150.6 148.1 144.0 137.1 158.2 146.7;run;data bb;set aa;a=range(of x1-x10);b=mean(of x1-x10);c=std(of x1-x10);d=cv(of x1-x10);file print;put a b c d;run;5. 截取函數(shù)簡介SAS截取函數(shù)應(yīng)用舉例例5：請把STD=1.92354和CV=0.601106用3個(gè)不同的截取函數(shù)(ceil、

15、int、round)進(jìn)行運(yùn)算，并分析這些函數(shù)的作用。data test5;STD=1.92354;CV=0.601106;a1=ceil(STD); b1=ceil(CV);a2=int(STD);b2=int(CV);a3=round(STD, 0.001);b3=round(CV, 0.001);ods html;proc print data=test2;Var a1-a3 b1-b3;run;ods html close;練習(xí)5求小于等于31.77的最小整數(shù)；將1984.314保留2位小數(shù)。解答：data;x=floor(31.07);y=round(1984.314,0.01);fi

16、le print;put x y;run;6. SAS概率函數(shù)SAS概率函數(shù)應(yīng)用舉例例6：試求出服從標(biāo)準(zhǔn)正態(tài)分布的隨機(jī)變量z取值分別小于等于-2.576、-1.960 、-1.645 、-1.282 、0 、1.282 、1.645 、1.960 、2.576所對應(yīng)的概率。data test3;a1=round(probnorm(-2.576),0.001);a2=round(probnorm(-1.960),0.001);a3=round(probnorm(-1.645),0.001);a4=round(probnorm(-1.282),0.001);a5=round(probnorm(0)

17、,0.001);a6=round(probnorm(1.282),0.001);a7=round(probnorm(1.645),0.001);a8=round(probnorm(1.960),0.001);a9=round(probnorm(2.576),0.001);Format a1-a9 5.3;ods html;proc print data=test3;Var a1-a9;run;ods html close;練習(xí) 6設(shè)隨機(jī)變量X服從自由度為3的中心卡方分布，計(jì)算X7.94的概率。解答data;p=1-probchi(7.94,3);put p=;run;7. SAS概率分布下分位

18、數(shù)函數(shù)例7：試求t臨界值表。options LS=95 PS=55;data abc;array T(50,9);array df(50) df1-df50;g=1;k=1;do d=1 to 40,50,60,70,80,90,100,200, 500,1000,100000000;df(k)=d;k+1;w=1;do alpha=0.50,0.20,0.10,0.05,0.02, 0.01,0.005,0.002,0.001;p=1-alpha/2;b=tinv(p,d);b=round(b,0.0001);T(g,w)=b; w=w+1;Output; end; g=g+1; end;T

19、itle; title t臨界值表;ods html style=journal;file print;put #1 3df alpha: 140.50 230.20 310.10390.05 480.02 570.01 660.005 760.002 860.001;put #2 _;do L=1 to 50;c=L+2; df(50)=9999;put #c 1 df(L) 4.0 13 T(L,1) 6.4 21T(L,2) 6.4 29 T(L,3) 6.437 T(L,4) 7.4 46T(L,5) 7.4 55 T(L,6) 7.4 64 T(L,7) 8.4 74 T(L,8) 8.4 84 T(L,9) 8.4;end;run;ods html close;footnote;footnote note: df=9999 stands for df=100000000;練習(xí)7設(shè)隨機(jī)變量t(df),計(jì)算t分布的分位數(shù)tp(df),其中自由度和對應(yīng)概率值分別為：n=1, 5, 10, 20; p=0.10, 0.90, 0.95, 0.975;解答：data _null_;t1=tinv(0.10,1);t2=tinv(0.90,5);t3=tinv(0.95,10);t4=tinv(0.975,20);put t1=