1、第一次實(shí)驗(yàn),1.了解SAS系統(tǒng)功能 2.熟悉SAS的頁(yè)面操作 3.熟悉SAS/BASE一些常用的語(yǔ)句 要求： (1)掌握數(shù)據(jù)的輸入輸出（data input put output cards等語(yǔ)句的配合使用）; (2)掌握常用的基礎(chǔ)統(tǒng)計(jì)計(jì)算（二項(xiàng)、正態(tài)等分布的概率函數(shù)值和分位數(shù)的計(jì)算）； (3)熟悉數(shù)據(jù)的圖形顯示功能,第一講 SAS系統(tǒng)概述,一、什么是 SAS系統(tǒng),二、 SAS/BASE簡(jiǎn)介,一、什么是 SAS系統(tǒng),SAS系統(tǒng)是一個(gè)模塊化的集成軟件系統(tǒng),提供數(shù)據(jù)訪問(wèn)、數(shù)據(jù)管理、數(shù)據(jù)分析、數(shù)據(jù)呈現(xiàn)。包括：,1. 基本部分（也是核心部分）： BASE SAS；,2. 統(tǒng)計(jì)分析： SAS/STAT；

2、,3. 高級(jí)繪圖： SAS/GRAPH；,4. 矩陣運(yùn)算： SAS/IML；,5. 運(yùn)籌學(xué)和線性規(guī)劃： SAS/OR；,6. 經(jīng)濟(jì)預(yù)測(cè)和時(shí)間序列分析： SAS/ETS；,7.面向任務(wù)的菜單驅(qū)動(dòng)界面：SAS/ASSIST,二、SAS系統(tǒng)的特點(diǎn),2. SAS語(yǔ)言是編程能力強(qiáng)且簡(jiǎn)潔易學(xué)的非過(guò)程語(yǔ)言,3. SAS系統(tǒng)把數(shù)據(jù)處理與統(tǒng)計(jì)分析融為一體,4. 適用性強(qiáng)、應(yīng)用面廣,二、SAS系統(tǒng)的特點(diǎn),1. 使用靈活方便、功能齊全,三、 SAS/BASE簡(jiǎn)介,主要功能是：,可以單獨(dú)使用，也可以同其它軟件產(chǎn)品一起組成一個(gè)用戶化的SAS系統(tǒng)。,3.報(bào)表生成和圖形顯示功能,1.數(shù)據(jù)管理,2.基礎(chǔ)統(tǒng)計(jì)計(jì)算功能,1.數(shù)據(jù)

3、管理功能：,（1）信息存貯和檢索：,可用任何格式讀入數(shù)據(jù)值，然后組成SAS數(shù)據(jù)集。具有很強(qiáng)的與外部文件交換信息的功能，可以用文件操作管理方法把不同數(shù)據(jù)集的數(shù)據(jù)組合在一起，供SAS過(guò)程分析處理，生成報(bào)表.,(2) 數(shù)據(jù)修改與程序設(shè)計(jì)：,(3) 文件操作：,用于對(duì)數(shù)據(jù)進(jìn)行編輯、整理、連接、合并及更新的功能，還能同時(shí)處理多個(gè)輸入文件，或者對(duì)一次輸入的數(shù)據(jù)能夠生成幾種報(bào)表。,3.報(bào)表生成和圖形顯示功能,2.基礎(chǔ)統(tǒng)計(jì)計(jì)算功能,.計(jì)算簡(jiǎn)單的描述性統(tǒng)計(jì)量,.SAS函數(shù)可用于計(jì)算概率分布函數(shù)、分位數(shù)、產(chǎn)生隨機(jī)數(shù)等,一、SAS顯示管理系統(tǒng)(Display Manager System，DMS),第二講 SAS使

4、用概述,（）編輯窗口（PROGRAM EDITOR）： 編輯程序和數(shù)據(jù)文件； （）日志窗口（LOG）： 記錄運(yùn)行情況，顯示信息； （）輸出窗口（OUTPUT）：輸出運(yùn)行的結(jié)果；（）圖形窗口（GRAPH）：輸出圖形。,點(diǎn)擊 Globals 菜單中的 Program editor、Log、Output、Graph 命令可以進(jìn)入編輯、日志、輸出及圖形窗口。 按功能鍵F5、F6、F7也可以進(jìn)入編輯、日志及輸出窗口。,是用戶和SAS系統(tǒng)交互式會(huì)話環(huán)境；點(diǎn)擊圖標(biāo)后，即可進(jìn)入的顯示管理系統(tǒng)在中有四個(gè)主要的窗口：,二、的功能鍵 用功能鍵可以代替對(duì)菜單的點(diǎn)擊，有時(shí)比較方便。 最常用的功能鍵有F1 ：顯示幫助信息

5、（）； 4 ：顯示已經(jīng)運(yùn)行的程序（）； 5 ：進(jìn)入編輯窗口（）； 6 ：進(jìn)入日志窗口（）； 7 ：進(jìn)入輸出窗口（）； 8 ：程式提交運(yùn)行（）； 9 ：顯示功能鍵（）； 以上功能鍵的用法熟悉之后，還可以定義或修改功能鍵。例如，可以定義Ctrl+E為功能鍵，用來(lái)清屏。,三、SAS語(yǔ)句,一個(gè)SAS語(yǔ)句是由SAS關(guān)鍵詞，SAS名字，特殊字符和運(yùn)算符組成的字符串，并以分號(hào)（；）結(jié)尾。,SAS關(guān)鍵詞: 如data語(yǔ)句，proc語(yǔ)句，infile語(yǔ)句等,SAS名字:最多由1-8個(gè)字符組成，首字符必須是字母或下 下劃線（字符包括：字母、數(shù)字或下劃線） 如：變量名、SAS數(shù)據(jù)集名、過(guò)程名等,SAS語(yǔ)句 一般形式

6、：,KEYWORD parameteroption;,注：SAS規(guī)定，每個(gè)語(yǔ)句以“；”作為該語(yǔ)句結(jié)束的標(biāo)志。,四、SAS程序,1.DATA語(yǔ)句：產(chǎn)生SAS數(shù)據(jù)集,2.PROC語(yǔ)句:對(duì)SAS數(shù)據(jù)集內(nèi)數(shù)據(jù)進(jìn)行分析處理,輸出結(jié)果,SAS程序：三明治,以DATA語(yǔ)句開(kāi)始(面包),告訴程序有關(guān)的數(shù)據(jù)集,最后是PROC語(yǔ)句,中間是數(shù)據(jù)(肉),一序列語(yǔ)句組成一個(gè)SAS程序，完成一定的功能.,SAS程序中語(yǔ)句分為兩種：,SAS程序構(gòu)成：,一個(gè)或多個(gè)data語(yǔ)句、一個(gè)或多個(gè)proc語(yǔ)句,一個(gè)或多個(gè)proc語(yǔ)句,一個(gè)或多個(gè)data語(yǔ)句,程序示例： data ex； input no$ sex$ age h w；

7、cards； 10 f 47 156.3 47.1 24 m 38 172.4 61.5 53 m 41 169.2 64.5 46 f 52 158.2 53.6 38 f 39 160.1 48 ； proc sort； by sex； proc means；by sex；var h w；run；,第四講 DATA步入門,一、什么是,二、創(chuàng)建SAS數(shù)據(jù)集,三、 DATA步語(yǔ)句及其類型,一、什么是DATA步,1.檢索:把輸入數(shù)據(jù)轉(zhuǎn)化為SAS數(shù)據(jù)集；,2.編輯:檢查數(shù)據(jù)中的錯(cuò)誤并修改它們，計(jì)算新變量；,3.按照用戶設(shè)計(jì)的規(guī)格輸出報(bào)告并存為磁盤文件；,4.從已經(jīng)存在的數(shù)據(jù)集中通過(guò)取子集、合并和

8、更新原數(shù)據(jù)的方法來(lái)產(chǎn)生新的數(shù)據(jù)集。,DATA步是用DATA語(yǔ)句開(kāi)始的一組或幾組SAS語(yǔ)句。其目的：,1.數(shù)據(jù)在作業(yè)流中：從作業(yè)流中輸入數(shù)據(jù)產(chǎn)生SAS數(shù)據(jù)集,二、創(chuàng)建SAS數(shù)據(jù)集,Ex.1 data d1; input x1 x2 x3; y=x1+x2+x3; cards; 1 2 3 4 5 6 ;,一般形式為：,注：數(shù)據(jù)之間以空格隔開(kāi)；數(shù)據(jù)行與行之間不用分號(hào)。,2.數(shù)據(jù)在磁盤上：從磁盤上讀入數(shù)據(jù)文件產(chǎn)生SAS數(shù)據(jù)集,二、創(chuàng)建SAS數(shù)據(jù)集,Ex.2 data d2; infile e:da1.dat; input x1 x2 x3; y=x1+x2+x3; run;,一般形式為：,3.數(shù)據(jù)來(lái)

9、自其它SAS數(shù)據(jù)集：從一個(gè)或幾個(gè)已經(jīng)存在的SAS數(shù)據(jù)集產(chǎn)生新的SAS數(shù)據(jù)集,二、創(chuàng)建SAS數(shù)據(jù)集,Ex.3 data d3; set da2; if x10; run;,一般形式為：,Ex.4 data d3; merge da1 da2; run;,4.把SAS數(shù)據(jù)集轉(zhuǎn)換為外部數(shù)據(jù)或輸出報(bào)表,二、創(chuàng)建SAS數(shù)據(jù)集,一般形式為：,三、 DATA步語(yǔ)句及其類型,只用于 DATA步的SAS語(yǔ)句共有56個(gè)，分為四類：文件操作語(yǔ)句、運(yùn)行語(yǔ)句、控制語(yǔ)句、信息語(yǔ)句。,常用的data步語(yǔ)句： 文件操作語(yǔ)句：data、input、cards、 file、 infile、by、set、merge等 運(yùn)行語(yǔ)句：賦

10、值語(yǔ)句、累加語(yǔ)句等 控制語(yǔ)句：循環(huán)語(yǔ)句、條件語(yǔ)句等 信息語(yǔ)句：（略）,（1） 語(yǔ)句 語(yǔ)句的作用是表明數(shù)據(jù)步的開(kāi)始并給出數(shù)據(jù)集的名稱。 語(yǔ)句的格式為： 數(shù)據(jù)集的名稱； 注：1）數(shù)據(jù)集的名稱必須以英文字母開(kāi)始，最長(zhǎng)不超過(guò)個(gè)字符。 2）數(shù)據(jù)集的名稱都是由符號(hào)“”分隔的兩個(gè)部分所組成。第一部分標(biāo)記數(shù)據(jù)集所存儲(chǔ)的位置，稱作一級(jí)名或庫(kù)邏輯名。第二部分標(biāo)記一個(gè)特定的數(shù)據(jù)集，稱作二級(jí)名。如果不想把數(shù)據(jù)集永久儲(chǔ)存在指定的磁盤上，可以只使用二級(jí)名。,例： ； 系統(tǒng)會(huì)自動(dòng)地把作為數(shù)據(jù)集的二級(jí)名，并添加作為數(shù)據(jù)集的一級(jí)名，因此在日志窗口顯示的信息中記該數(shù)據(jù)集為。但是這個(gè)數(shù)據(jù)集是臨時(shí)的，它僅僅在程式運(yùn)行期間有效，過(guò)后它

11、就被復(fù)蓋，不能重新調(diào)用。,（2） 語(yǔ)句 語(yǔ)句的作用是與“；”呼應(yīng)，標(biāo)志數(shù)據(jù)行的開(kāi)始與結(jié)束。 語(yǔ)句的格式為： ； 數(shù)據(jù)行 ； 注：如果使用語(yǔ)句，在的后面必須緊跟數(shù)據(jù)行，并且在一個(gè)數(shù)據(jù)步中最多只能有一個(gè)語(yǔ)句。,（3）語(yǔ)句 語(yǔ)句的作用是描述輸入記錄中的數(shù)據(jù)，并把輸入值賦給相應(yīng)的變量。 語(yǔ)句的格式為： 數(shù)據(jù)的變量名、順序及類型； 注1：用語(yǔ)句是為了讀外部文件的數(shù)據(jù)或跟在語(yǔ)句后面的數(shù)據(jù)。除非在語(yǔ)句中的變量名后有串符號(hào)或用字符的輸入格式表示、或該變量事先已被定義為字符型，否則認(rèn)定用語(yǔ)句讀入的是數(shù)值型變量的值。,注2.用語(yǔ)句時(shí)，外部文件中的數(shù)據(jù)和語(yǔ)句后面的數(shù)據(jù)都采取列表輸入的方法，各個(gè)變量的值由它們之間的空

12、格來(lái)分隔。為從一行讀入多個(gè)觀測(cè)值，應(yīng)使用行保持符限制讀數(shù)指針，使其保持在這一行上讀數(shù)，直到數(shù)據(jù)讀完為止。,Ex.1 data ex1; input x1 x2 x3; y=x1+x2+x3; cards; 1 2 3 4 5 6 ;,Ex.1 data ex1; input x1 x2 x3; y=x1+x2+x3; cards; 1 2 3 4 5 6 ;,(4)語(yǔ)句 語(yǔ)句的作用是指明外部數(shù)據(jù)文件的名稱，并從這個(gè)外部數(shù)據(jù)文件中讀取數(shù)據(jù)。 語(yǔ)句的格式為： 外部數(shù)據(jù)文件名； 例如，在驅(qū)動(dòng)器上有一個(gè)名稱為ex2的數(shù)據(jù)文件，其中數(shù)據(jù)排列為: 3.12 2.23 3.4 4.12 4.2 .5 用語(yǔ)句

13、讀取數(shù)據(jù)的程式為： ； ：ex2.； $ ； 用語(yǔ)句讀取數(shù)據(jù)時(shí)，要求外部數(shù)據(jù)文件必須是可以在操作系統(tǒng)中用命令打印出全部?jī)?nèi)容的碼文本文件。,(5) 語(yǔ)句 語(yǔ)句的作用是按照指定的條件從指定的數(shù)據(jù)集中讀取數(shù)據(jù)建立新的數(shù)據(jù)集或?qū)蓚€(gè)數(shù)據(jù)集中的觀測(cè)值縱向連接建立新的數(shù)據(jù)集。 語(yǔ)句的格式為： 數(shù)據(jù)集名表； 例如，從數(shù)據(jù)集中讀取數(shù)據(jù)建立新的數(shù)據(jù)集，可編輯語(yǔ)句： ； ； 將兩個(gè)數(shù)據(jù)集和中的觀測(cè)值縱向連接建立新的數(shù)據(jù)集，可編輯語(yǔ)句： ； ；,data d1；input g$ a b ； cards； 1 5 7 1 7 9 1 3 4 ； data d2；input g$ a c ； cards； 2 4 8

14、2 6 6 2 8 10 ； data d3；set d1 d2； proc print； run；,運(yùn)行結(jié)果 Obs g a b c 1 1 5 7 . 2 1 7 9 . 3 1 3 4 . 4 2 4 . 8 5 2 6 . 6 6 2 8 . 10,(6) 語(yǔ)句 語(yǔ)句的作用是將兩個(gè)數(shù)據(jù)集中的各個(gè)觀測(cè)值橫向合并建立新的數(shù)據(jù)集。語(yǔ)句的格式為： 數(shù)據(jù)集名表； 例如，將數(shù)據(jù)集與中的觀測(cè)值橫向合并建立新的數(shù)據(jù)集，可編輯語(yǔ)句： ； ；,data d1；input a b c ； cards； 1 5 7 1 7 9 1 3 4 1 2 6 ； data d2;input x y z ； cards

15、； 2 4 8 2 6 6 ； data d3；merge d1 d2； proc print；run； 本程式運(yùn)行的結(jié)果如下： OBS A B C X Y Z 1 1 5 7 2 4 8 2 1 7 9 2 6 6 3 1 3 4 . . . 4 1 2 6 . . .,(7) do語(yǔ)句 簡(jiǎn)單do語(yǔ)句 循環(huán)do語(yǔ)句 do while 語(yǔ)句 do until 語(yǔ)句,簡(jiǎn)單do語(yǔ)句：常常用在if-then/else語(yǔ)句里來(lái)指定一個(gè)語(yǔ)句組根據(jù)if條件成立或不成立時(shí)被執(zhí)行，格式為： do; 一些sas語(yǔ)句 end;,例： data a; input x; if x5 then do; y=x*10;

16、put x= y= ; end; cards; 12 4 11 5 ;,輸出結(jié)果： x=12 y=120 x=11 y=110,例： data ex; do x=1.645, 1.96, 2.576; p=probnorm(x); put x= p= ; end; run;,在log窗口輸出結(jié)果： x=1.645 p=0.9500150945 x=1.96 p=0.9750021049 x=2.576 p=0.9950024677,循環(huán)do語(yǔ)句： 格式為： do; 一些sas語(yǔ)句 end;,例： data a; do x=-3 to 3 by 0.25; do y=-3 to 3 by 0.2

17、5; p=x*x+y*y; put x= y= p= ; end; end; run;,第五講 PROC步入門,一、什么是PROC 步？,二、 PROC語(yǔ)句,一、什么是PROC 步？,PROC 步是用PROC語(yǔ)句開(kāi)始的一組或幾組SAS語(yǔ)句。,一般形式為：,Ex.4 proc print data=d1; run;,注： 1.不同的SAS過(guò)程可以使用的 也是不同的。,2. Proc步中常出現(xiàn)的SAS控制語(yǔ)句： （1）var 變量名; （2）by 變量名;（要求排序） （3）class 變量名; (4) ID 變量名;,二、 SAS基本輸出過(guò)程,proc print; 或 proc print d

18、ata=文件名;,2.proc plot; 語(yǔ)法：proc plot 選項(xiàng)； by 變量列表; plot 選項(xiàng)；,選項(xiàng)：如vpect=33(規(guī)定產(chǎn)生的散布圖在垂直方向占一頁(yè)的比例 hpct=50 formchar =字符串,plot 選項(xiàng)；plot vertical*horizontal 或plot vertical*horizontal=字符或變量,3.proc sort ; by variables; run;,1.proc means ; var variables; by variables; freq variables; output out=新文件名 變量名列表（統(tǒng)計(jì)關(guān)鍵詞=）

19、run;,三、 SAS基礎(chǔ)統(tǒng)計(jì)過(guò)程,Options: mean std clm alpha t prt,注1：proc means 可以計(jì)算的統(tǒng)計(jì)量由21種，分別用下述關(guān)鍵詞表示： n,nmiss,mean, std,max,min, range,sum,var,uss,css,cv,stderr,t,prt, skewness,kurtosis,clm,lclm,uclm,注2：output out=新文件名 變量名列表 （統(tǒng)計(jì)關(guān)鍵詞=）,Ex:proc means data=ex； var weight height; class sex; output out=result mean=m

20、eanweight meanheight std(weight)=stdweight std(height)=stdheight run; proc print data=result; run;,偏態(tài)與峰度分布的形狀,偏態(tài)(概念要點(diǎn)),1、數(shù)據(jù)分布偏斜程度的測(cè)度 2.偏態(tài)系數(shù)=0為對(duì)稱分布 3.偏態(tài)系數(shù) 0為右偏分布 4.偏態(tài)系數(shù) 0為左偏分布 5. 計(jì)算公式為,偏態(tài)(實(shí)例),【例】已知1997年我國(guó)農(nóng)村居民家庭按純收入分組的有關(guān)數(shù)據(jù)。試計(jì)算偏態(tài)系數(shù),農(nóng)村居民家庭村收入數(shù)據(jù)的直方圖,偏態(tài)與峰度(從直方圖上觀察),按純收入分組(元),結(jié)論：1. 為右偏分布 2. 峰度適中,偏態(tài)系數(shù)(計(jì)算結(jié)果)

21、,將計(jì)算結(jié)果代入公式得,結(jié)論：偏態(tài)系數(shù)為正值，而且數(shù)值較大，說(shuō)明農(nóng)村居民家庭純收入的分布為右偏分布，即收入較少的家庭占據(jù)多數(shù)，而收入較高的家庭則占少數(shù)，而且偏斜的程度較大,峰度(概念要點(diǎn)),1.數(shù)據(jù)分布扁平程度的測(cè)度 2.峰度系數(shù)=3扁平程度適中 3.偏態(tài)系數(shù)3為尖峰分布 5. 計(jì)算公式為,峰度系數(shù)系數(shù)(實(shí)例計(jì)算結(jié)果),代入公式得,【例】根據(jù)表4-10中的計(jì)算結(jié)果，計(jì)算農(nóng)村居民家庭純收入分布的峰度系數(shù),結(jié)論：由于=3.43，說(shuō)明我國(guó)農(nóng)村居民家庭純收入的分布為尖峰分布，說(shuō)明低收入家庭占有較大的比重,2.proc univariate var variables; by variables; fr

22、eq variables; output run;,Options:normal,注：計(jì)算分位數(shù)、頻率表、數(shù)據(jù)的正態(tài)性檢驗(yàn)等,Options: 1.Normal 2.Plot (生成一個(gè)莖葉圖或水平直方圖、 一個(gè)盒形圖、一個(gè)正態(tài)概率圖),3.proc ttest class variables; var variables; run;,4.proc corr var variables; with variables; by variables; freq variables; run;,注：計(jì)算Pearson相關(guān)系數(shù)和顯著性概率以及 單變量統(tǒng)計(jì)量,第六講 全程語(yǔ)句,一、什么是全程語(yǔ)句？,二、

23、 常用的全程語(yǔ)句,一、什么是全程語(yǔ)句？,全程語(yǔ)句： 既可以用在data步，也可以用在proc步； 既可以用在data步和proc步的外部， 也可以用在data步和proc步的內(nèi)部。 這些語(yǔ)句可以用在SAS程序的任何地方，稱為全程語(yǔ)句.,二、 常用的全程語(yǔ)句,1.注釋語(yǔ)句： * messages; 或 /* messages*/ (無(wú)需分號(hào)）,2.Title語(yǔ)句：規(guī)定同SAS輸出一起被打印的標(biāo)題行 Title (n=110),e.g. Title2 Annual costs,3.run ;,4.endsas;,應(yīng)用SAS計(jì)算二項(xiàng)分布的分位數(shù),probbnml(p,n,k)=,而PX=k= pro

24、bbnml(p,n,k)probbnml(p,n,k1),data ex; do k=1 to 10 by 1 ; p=probbnml(0.5,10,k)-probbnml(0.5,10,k-1); put p; end;,data ex; do k=1 to 10 by 1 ; do p=0.5;n=10; value=probbnml(p,n,k)-probbnml(p,n,k-1); output; end; end; proc print data=ex;run;,應(yīng)用SAS計(jì)算標(biāo)準(zhǔn)正態(tài)分布的分位數(shù) 在SAS中有probit(p)函數(shù)，用此函數(shù)可以求p分位數(shù)。SAS程序?yàn)?data

25、ex; do p=0.025,0.05,0.1,0.9,0.95,0.975; u=probit(p);put u p ; end;run;,輸出的結(jié)果如下： -1.959963985 0.025 -1.644853627 0.05 -1.281551566 0.1 1.2815515655 0.9 1.644853627 0.95 1.9599639845 0.975,當(dāng)=0.10，0.05，0.01時(shí)，應(yīng)用SAS計(jì)算雙側(cè)分位數(shù)的程序?yàn)椋?data ex;do x=0.1,0.05,0.01; p=1-x/2;u=probit(p);put x p u; end;run;,輸出的結(jié)果如下：

26、0.1 0.95 1.644853627 0.05 0.975 1.9599639845 0.01 0.995 2.5758293035,應(yīng)用SAS計(jì)算卡方分布的分位數(shù) 在SAS中有cinv(p,df)函數(shù)，用此函數(shù)可以求p分位數(shù)。SAS程序?yàn)?data ex; do df=4;do p=0.025,0.05,0.1,0.9,0.95,0.975; c=cinv(p,df);put p df c; end;end;run;,輸出的結(jié)果如下： 0.025 4 0.4844185571 0.05 4 0.7107230214 0.1 4 1.0636232168 0.9 4 7.7794403397 0.95 4 9.4877290368 0.975 4 11.143286782,應(yīng)用SAS計(jì)算t分布的分位數(shù) 在SAS中有tinv(p,df)函數(shù)，用此函數(shù)可以求