1、常用SAS函數(shù)簡介,阮紅蓮 廣州醫(yī)科大學(xué) 統(tǒng)計學(xué)系,SAS函數(shù),SAS系統(tǒng)提供的SAS函數(shù)比一般高級語言的標準函數(shù)多得多，尤其是概率函數(shù)、分位數(shù)函數(shù)、樣本統(tǒng)計函數(shù)、隨機數(shù)函數(shù)。 如此豐富的SAS函數(shù),對用戶統(tǒng)計分析計算以及編寫SAS程序帶來極大的方便。,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),SAS函數(shù)中的參數(shù),1）參數(shù)類型： 變量名；常數(shù)；函數(shù)；表達式。,例1 函數(shù)及其參數(shù)類型舉例。,2）參數(shù)個數(shù)： 有些函數(shù)不需要自變量，有

2、些只需要一個自變量，有些函數(shù)有多個自變量。 3）參數(shù)表示法： 當(dāng)函數(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 正確與錯誤的表示方法 例3 數(shù)組元素自變量的縮寫方法 array y10 y1-y10; /*定義數(shù)組y10 */ x=sum (of

3、 y*); z=sum (of y1-y10); 例中，兩個SUM函數(shù)的表示法是等價的。,函數(shù)結(jié)果,通常函數(shù)的結(jié)果由其自變量的屬性決定： 自變量是字符時結(jié)果變量為字符； 自變量是數(shù)值時結(jié)果變量為數(shù)值。 PUT函數(shù)是一個例外，它不管自變量是什么類型，其結(jié)果總是字符值。,用PUT語句在LOG窗口顯示函數(shù)值,例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.

4、95分位數(shù)*/ put q1= q2=; run;,0.9750021049 q1=2.9199855804 q2=13.894376071,利用HELP查找SAS系統(tǒng)提供的各種函數(shù),索引中鍵入： SAS Function,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ù)。,1.日期時間函數(shù)簡介

5、,例1.1：計算兩個日期之間相隔的天數(shù),例1.2: 以日為單位計算當(dāng)前日期的天數(shù)。 data; x=date (); y=today(); format x y yymmdd10.; put x= y=; run;,例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;,練習(xí)1,計算以下數(shù)據(jù)表中日期值所在的年季月以及處于某月的第幾天、某周的第幾天。 當(dāng)日期值

6、date=14JUL2007d時,試問這一天是星期幾?,Date 1991-01-02 1991-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,解答：,解答：,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ù)ex SIN(x)、COS(

7、x)、TAN(x)：正弦、余弦、正切函數(shù) .,例2：,求-314的絕對值，31除以4的商的余數(shù)以及-314 、-87、0、1127這4個數(shù)的最大值,練習(xí)2,ex2.1 計算0、1.57、 3.14對應(yīng)的正弦函數(shù)值，并將結(jié)果顯示在輸出窗口。 ex2.2 計算e的0次冪,并計算1的自然對數(shù)。,解答：,3. 數(shù)組函數(shù),數(shù)組函數(shù)應(yīng)用舉例,例3：產(chǎn)生一個5102的三維數(shù)組，并提取各維度上元素的個數(shù)。,練習(xí) 3,產(chǎn)生一個行變量名為2:6，列變量名為4:13的二維數(shù)組，輸出各維度上的上界和下界。,解答：,4. 樣本統(tǒng)計量函數(shù),簡單統(tǒng)計量函數(shù)應(yīng)用舉例,例4：請求出2、4、6、3、1這五個數(shù)的最大值、最小值、算

8、術(shù)平均值、標準差、變異系數(shù)。,用簡便的方法實現(xiàn)例1中目的,data simple_statistics; 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年隨機抽查測量10名13歲女孩的身高

9、，資料如下，試求全距、均值、標準差、標準誤、變異系數(shù)。 145.1 143.3 153.2 152.3 150.6 148.1 144.0 137.1 158.2 146.7,解答：,5. 截取函數(shù)簡介,SAS截取函數(shù)應(yīng)用舉例,例5：請把STD=1.92354和CV=0.601106用3個不同的截取函數(shù)(ceil、int、round)進行運算，并分析這些函數(shù)的作用。,練習(xí)5,求小于等于31.77的最小整數(shù)； 將1984.314保留2位小數(shù)。,解答：,6. SAS概率函數(shù),SAS概率函數(shù)應(yīng)用舉例,例6：試求出服從標準正態(tài)分布的隨機變量z取值分別小于等于-2.576、-1.960 、-1.645

10、、-1.282 、0 、1.282 、1.645 、1.960 、2.576所對應(yīng)的概率。,練習(xí) 6,設(shè)隨機變量X服從自由度為3的中心卡方分布，計算X7.94的概率。,解答,7. SAS概率分布下分位數(shù)函數(shù),例7：試求t臨界值表。,練習(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;,解答：,8. SAS隨機數(shù)函數(shù),UNIFORM(seed)和RANUNI(seed)：產(chǎn)生均勻分布隨機數(shù)函數(shù)； NORMAL(seed)和RANNOR(seed)：產(chǎn)生正態(tài)分布隨機數(shù)函數(shù)； RANEXP(seed)：指數(shù)分布隨機數(shù)函數(shù)； RANGAM(seed,alpha)：伽馬分布隨機數(shù)函數(shù)； RANTRI(seed,h)：三角分布隨機數(shù)函數(shù)； RANCAU(seed) ：柯西分布隨機數(shù)函數(shù)； RANBIN(seed,n,p) ：二項分布隨機數(shù)