1、關(guān)于常用函數(shù)簡介第一張，PPT共四十九頁，創(chuàng)作于2022年6月SAS函數(shù) SAS系統(tǒng)提供的SAS函數(shù)比一般高級語言的標準函數(shù)多得多，尤其是概率函數(shù)、分位數(shù)函數(shù)、樣本統(tǒng)計函數(shù)、隨機數(shù)函數(shù)。 如此豐富的SAS函數(shù),對用戶統(tǒng)計分析計算以及編寫SAS程序帶來極大的方便。第二張，PPT共四十九頁，創(chuàng)作于2022年6月SAS函數(shù)定義 SAS函數(shù)是一個子程序，對0個或幾個參數(shù)（自變量）返回一個結(jié)果值。SAS函數(shù)的形式：函數(shù)名（參數(shù)1，參數(shù)2，）例如： MEAN(x,y,z,a)MIN(SUM(x,y,z),1000)第三張，PPT共四十九頁，創(chuàng)作于2022年6月SAS函數(shù)中的參數(shù)1）參數(shù)類型：變量名；常數(shù)；

2、函數(shù)；表達式。例1 函數(shù)及其參數(shù)類型舉例。 第四張，PPT共四十九頁，創(chuàng)作于2022年6月2）參數(shù)個數(shù)： 有些函數(shù)不需要自變量，有些只需要一個自變量，有些函數(shù)有多個自變量。3）參數(shù)表示法：當函數(shù)有多個參數(shù)時，必須用逗號分隔開。4）參數(shù)名縮寫方法： 函數(shù)名(OF 變量名1 變量名n ) 第五張，PPT共四十九頁，創(chuàng)作于2022年6月參數(shù)的表示法舉例以逗號(，)分隔參數(shù),如 sum(x1, x2, x3, x4, x5, x6, x7, x8)去掉逗號(，)，最前面加上“of”,如 sum(of x1 x2 x3 x4 x5 x6 x7 x8)去掉逗號(，)，最前面加上“of”,變量簡寫：sum(

3、of x1 - x8)第六張，PPT共四十九頁，創(chuàng)作于2022年6月例2 正確與錯誤的表示方法例3 數(shù)組元素自變量的縮寫方法array y10 y1-y10; /*定義數(shù)組y10 */x=sum (of y*);z=sum (of y1-y10);例中，兩個SUM函數(shù)的表示法是等價的。 第七張，PPT共四十九頁，創(chuàng)作于2022年6月函數(shù)結(jié)果通常函數(shù)的結(jié)果由其自變量的屬性決定： 自變量是字符時結(jié)果變量為字符； 自變量是數(shù)值時結(jié)果變量為數(shù)值。PUT函數(shù)是一個例外，它不管自變量是什么類型，其結(jié)果總是字符值。第八張，PPT共四十九頁，創(chuàng)作于2022年6月用PUT語句在LOG窗口顯示函數(shù)值例4 顯示概率

4、值和分位數(shù)。data;Y=probnorm(1.96); /* 標準正態(tài)分布Z小于1.96的概率 */put 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第九張，PPT共四十九頁，創(chuàng)作于2022年6月利用HELP查找SAS系統(tǒng)提供的各種函數(shù)索引中鍵入：SAS FunctionHelp下拉菜單SAS 幫助與文件第十張，PPT共四十九頁，創(chuàng)作于2022

5、年6月SAS函數(shù)種類在SAS系統(tǒng)9.1版本,SAS函數(shù)分為二十五種類型,分別是算術(shù)函數(shù)、截取函數(shù)、數(shù)學(xué)函數(shù)、概率和密度函數(shù)、分位數(shù)函數(shù)、非中心函數(shù)、樣本統(tǒng)計函數(shù)、隨機數(shù)函數(shù)、財政金融函數(shù)、日期和時間函數(shù)、換算函數(shù)和特殊函數(shù)等。在SAS系統(tǒng)6.12版本，SAS函數(shù)分為十七種類型共有178個,SAS系統(tǒng)8.1版本，SAS函數(shù)分為二十二種類型共有319個函數(shù)。第十一張，PPT共四十九頁，創(chuàng)作于2022年6月第十二張，PPT共四十九頁，創(chuàng)作于2022年6月1.日期時間函數(shù)簡介第十三張，PPT共四十九頁，創(chuàng)作于2022年6月第十四張，PPT共四十九頁，創(chuàng)作于2022年6月例1.1：計算兩個日期之間相隔的

6、天數(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, 01jan2003d, act/act);days360=datdif(01jan2002d, 01jan2003d, 30/360);put actual= days360=; run;第十五張，PPT共四十九頁，創(chuàng)作于2

7、022年6月例1.2: 以日為單位計算當前日期的天數(shù)。data;x=date (); y=today();format x y yymmdd10.;put x= y=;run;第十六張，PPT共四十九頁，創(chuàng)作于2022年6月例1.3: 將時間值換算為以秒計數(shù)的計算方法。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 yymmdd10.;run;第十七張，PPT共四十九頁，創(chuàng)作于2022年6月練習(xí)1計算以下數(shù)據(jù)表中日期值所在的年季月以及處

8、于某月的第幾天、某周的第幾天。當日期值date=14JUL2007d時,試問這一天是星期幾? 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第十八張，PPT共四十九頁，創(chuàng)作于2022年6月解答：data a;Input date YYMMDD10.;Cards;1991-01-021991-01-03 1991-01-04 1991-01-05 1991-01-07 1991-01-08 1991-01-09 1991

9、-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;第十九張，PPT共四十九頁，創(chuàng)作于2022年6月解答：data _null_;weekday=weekday(12JUL2007d);put weekday=; run;第二十張

10、，PPT共四十九頁，創(chuàng)作于2022年6月2. 數(shù)學(xué)函數(shù)ABS (x) ：求x的絕對值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ù).第二十一張，PPT共四十九頁，創(chuàng)作于2022年6月例2：求-314的絕對值，31除以4的商的余數(shù)以及-314 、-87、0、1127這4個數(shù)的最大值data _null_;x=abs(-314);y=mod(31,4);z=max(

11、-314 ,-87,0,1127);put x= y= z=; run;第二十二張，PPT共四十九頁，創(chuàng)作于2022年6月練習(xí)2ex2.1 計算0、1.57、 3.14對應(yīng)的正弦函數(shù)值，并將結(jié)果顯示在輸出窗口。ex2.2 計算e的0次冪,并計算1的自然對數(shù)。第二十三張，PPT共四十九頁，創(chuàng)作于2022年6月解答：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

12、.2 c 6.2;Run;Data;X=exp(0);Y=log(1);Put x y;Run;第二十四張，PPT共四十九頁，創(chuàng)作于2022年6月3. 數(shù)組函數(shù)第二十五張，PPT共四十九頁，創(chuàng)作于2022年6月數(shù)組函數(shù)應(yīng)用舉例data;array mult5,10,2 mult1-mult100;a=DIM1(MULT);b=DIM2(MULT);c=DIM3(MULT);put a b c;run;例3：產(chǎn)生一個5102的三維數(shù)組，并提取各維度上元素的個數(shù)。第二十六張，PPT共四十九頁，創(chuàng)作于2022年6月練習(xí) 3產(chǎn)生一個行變量名為2:6，列變量名為4:13的二維數(shù)組，輸出各維度上的上界和下

13、界。第二十七張，PPT共四十九頁，創(chuàng)作于2022年6月解答：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;第二十八張，PPT共四十九頁，創(chuàng)作于2022年6月4. 樣本統(tǒng)計量函數(shù)第二十九張，PPT共四十九頁，創(chuàng)作于2022年6月簡單統(tǒng)計量函數(shù)應(yīng)用舉例例4：請求出2、4、6、3、1這五個數(shù)的最大值、最小值、算術(shù)平均值、標準差、變異系數(shù)。data test4;input x1-x5;cards;2 4 6 3 1;r

14、un;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;第三十張，PPT共四十九頁，創(chuàng)作于2022年6月用簡便的方法實現(xiàn)例1中目的data simple_statistics;max_

15、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;第三十一張，PPT共四十九頁，創(chuàng)作于2022年6月練習(xí)4某農(nóng)村地區(qū)1998年隨機抽查測量10名13歲女孩的身高，資料如下，試求全距、均值、標準差、標準誤、變異系數(shù)。 145.1 143.3 153.

16、2 152.3 150.6 148.1 144.0 137.1 158.2 146.7第三十二張，PPT共四十九頁，創(chuàng)作于2022年6月解答：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;第三十三張，PPT共四十九頁，創(chuàng)作于2022年6月5

17、. 截取函數(shù)簡介第三十四張，PPT共四十九頁，創(chuàng)作于2022年6月SAS截取函數(shù)應(yīng)用舉例例5：請把STD=1.92354和CV=0.601106用3個不同的截取函數(shù)(ceil、int、round)進行運算，并分析這些函數(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 c

18、lose;第三十五張，PPT共四十九頁，創(chuàng)作于2022年6月練習(xí)5求小于等于31.77的最小整數(shù)；將1984.314保留2位小數(shù)。第三十六張，PPT共四十九頁，創(chuàng)作于2022年6月解答：data;x=floor(31.07);y=round(1984.314,0.01);file print;put x y;run;第三十七張，PPT共四十九頁，創(chuàng)作于2022年6月6. SAS概率函數(shù)第三十八張，PPT共四十九頁，創(chuàng)作于2022年6月SAS概率函數(shù)應(yīng)用舉例例6：試求出服從標準正態(tài)分布的隨機變量z取值分別小于等于-2.576、-1.960 、-1.645 、-1.282 、0 、1.282 、1

19、.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),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

20、(probnorm(2.576),0.001);Format a1-a9 5.3;ods html;proc print data=test3;Var a1-a9;run;ods html close;第三十九張，PPT共四十九頁，創(chuàng)作于2022年6月練習(xí) 6設(shè)隨機變量X服從自由度為3的中心卡方分布，計算X7.94的概率。第四十張，PPT共四十九頁，創(chuàng)作于2022年6月解答data;p=1-probchi(7.94,3);put p=;run;第四十一張，PPT共四十九頁，創(chuàng)作于2022年6月7. SAS概率分布下分位數(shù)函數(shù)第四十二張，PPT共四十九頁，創(chuàng)作于2022年6月例7：試求t臨界值表

21、。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;Title; title t

22、臨界值表;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;第四十三張，PPT共四十九頁，創(chuàng)作于2022年6月練習(xí)7設(shè)隨機變量t(df),計算t分布的分位數(shù)tp(df),其中自由度和對應(yīng)概率值分別為：n=1, 5, 10, 20; p=0.10, 0.90, 0.95, 0.975;第