SAS統(tǒng)計(jì)分析軟件應(yīng)用

1、1 tatistical nalysis ystem 鹽城師范學(xué)院數(shù)學(xué)科學(xué)學(xué)院鹽城師范學(xué)院數(shù)學(xué)科學(xué)學(xué)院 統(tǒng)計(jì)系統(tǒng)計(jì)系 李萬斌李萬斌 2014.9 2 第一節(jié)第一節(jié) SASSAS軟件簡(jiǎn)介軟件簡(jiǎn)介 1966年美國(guó)年美國(guó)North Carolina州立大學(xué)的兩位生物州立大學(xué)的兩位生物 統(tǒng)計(jì)學(xué)研究生開始研制統(tǒng)計(jì)學(xué)研究生開始研制SAS系統(tǒng)，直至系統(tǒng)，直至1976年成立了年成立了 SAS軟件研究所，正式推出了軟件研究所，正式推出了SAS軟件，軟件，1985年推出年推出 微機(jī)版微機(jī)版SAS。 SAS是用于決策支持的大型集成信息系統(tǒng)，軟件是用于決策支持的大型集成信息系統(tǒng)，軟件 系統(tǒng)最早的功能限于統(tǒng)計(jì)分析，逐漸

2、成為一個(gè)用來管系統(tǒng)最早的功能限于統(tǒng)計(jì)分析，逐漸成為一個(gè)用來管 理、分析數(shù)據(jù)和編寫報(bào)告的大型集成應(yīng)用軟件系統(tǒng)，理、分析數(shù)據(jù)和編寫報(bào)告的大型集成應(yīng)用軟件系統(tǒng)， 具有完備的數(shù)據(jù)訪問、管理、分析、呈現(xiàn)及應(yīng)用開發(fā)具有完備的數(shù)據(jù)訪問、管理、分析、呈現(xiàn)及應(yīng)用開發(fā) 等功能，完全超出了單純統(tǒng)計(jì)應(yīng)用的功能。屬于世界等功能，完全超出了單純統(tǒng)計(jì)應(yīng)用的功能。屬于世界 領(lǐng)先，使用最為廣泛的統(tǒng)計(jì)軟件之一。領(lǐng)先，使用最為廣泛的統(tǒng)計(jì)軟件之一。 3 由于由于SAS系統(tǒng)是從大型機(jī)上的系統(tǒng)發(fā)展系統(tǒng)是從大型機(jī)上的系統(tǒng)發(fā)展 而來，在設(shè)計(jì)上也是完全針對(duì)專業(yè)用戶進(jìn)行而來，在設(shè)計(jì)上也是完全針對(duì)專業(yè)用戶進(jìn)行 設(shè)計(jì)，因此操作以設(shè)計(jì)，因此操作以SA

3、S語言編程為主，另外語言編程為主，另外 也有非編程方式的菜單系統(tǒng)，但人機(jī)對(duì)話的也有非編程方式的菜單系統(tǒng)，但人機(jī)對(duì)話的 界面不太友好。界面不太友好。 學(xué)習(xí)學(xué)習(xí)SAS軟件時(shí)，需要用戶有基本的統(tǒng)軟件時(shí)，需要用戶有基本的統(tǒng) 計(jì)學(xué)知識(shí)，對(duì)所選用的統(tǒng)計(jì)方法有較清楚的計(jì)學(xué)知識(shí)，對(duì)所選用的統(tǒng)計(jì)方法有較清楚的 了解，沒有統(tǒng)計(jì)學(xué)基礎(chǔ)知識(shí)要掌握了解，沒有統(tǒng)計(jì)學(xué)基礎(chǔ)知識(shí)要掌握SAS軟件軟件 較為困難。較為困難。 4 SAS系統(tǒng)是一個(gè)模塊化的組合軟件系統(tǒng)，它系統(tǒng)是一個(gè)模塊化的組合軟件系統(tǒng)，它 提供了約提供了約20多個(gè)模塊，各個(gè)模塊之間既相互獨(dú)立多個(gè)模塊，各個(gè)模塊之間既相互獨(dú)立 又相互交融補(bǔ)充。本課程用得最多的是又相互交

4、融補(bǔ)充。本課程用得最多的是Base SAS 模塊和模塊和SAS/STAT模塊。模塊。 Base SAS是是SAS系統(tǒng)的基礎(chǔ)。系統(tǒng)的基礎(chǔ)。它既可以單獨(dú)它既可以單獨(dú) 使用，也可以與其他模塊組成一個(gè)用戶化的使用，也可以與其他模塊組成一個(gè)用戶化的SAS系系 統(tǒng)，但是其他模塊必須與之結(jié)合起來才能使用。統(tǒng)，但是其他模塊必須與之結(jié)合起來才能使用。 Base SAS主要承擔(dān)著數(shù)據(jù)及用戶使用環(huán)境的管理、主要承擔(dān)著數(shù)據(jù)及用戶使用環(huán)境的管理、 SAS語言程序的處理，并具有基本的數(shù)據(jù)分析和報(bào)語言程序的處理，并具有基本的數(shù)據(jù)分析和報(bào) 告等統(tǒng)計(jì)功能。告等統(tǒng)計(jì)功能。 5 SAS/STAT提供了當(dāng)今流行的主要統(tǒng)計(jì)提供了當(dāng)今流

5、行的主要統(tǒng)計(jì) 分析方法分析方法，是國(guó)際上統(tǒng)計(jì)分析領(lǐng)域的標(biāo)準(zhǔn)權(quán)，是國(guó)際上統(tǒng)計(jì)分析領(lǐng)域的標(biāo)準(zhǔn)權(quán) 威軟件。它具有回歸分析、方差分析、屬性威軟件。它具有回歸分析、方差分析、屬性 數(shù)據(jù)分析、多元分析、聚類分析、判別分析、數(shù)據(jù)分析、多元分析、聚類分析、判別分析、 非參數(shù)分析、生存分析和心理測(cè)量分析等統(tǒng)非參數(shù)分析、生存分析和心理測(cè)量分析等統(tǒng) 計(jì)功能。計(jì)功能。 6 第二節(jié)第二節(jié) SASSAS軟件的使用基礎(chǔ)軟件的使用基礎(chǔ) 一、一、SAS的啟動(dòng)和退出的啟動(dòng)和退出 v啟動(dòng)啟動(dòng)SAS v雙擊桌面SAS圖標(biāo) v開始菜單 程序 The Sas System The Sas System for Windows V*.*

6、v直接雙擊SAS應(yīng)用程序 v退出退出SAS v標(biāo)題欄關(guān)閉按鈕 v菜單欄 File Exit 7 標(biāo)題欄、菜單欄、工具欄標(biāo)題欄、菜單欄、工具欄 vEditor 程序編輯器窗口 vLOG 運(yùn)行過程日志窗口 vOutput 結(jié)果輸出窗口 vExplorer 資源管理器窗口 vResults 結(jié)果索引窗口 三個(gè)窗口間切換F4 Recall (調(diào)回程序) F5 Editor (程序)F8 Submit (發(fā)送程序) F6 LOG (日志)F1 Help (幫助) F7 Output (結(jié)果)F9 Keys (功能鍵列表) 二、二、SAS的工作窗口的工作窗口 SAS常用功能鍵常用功能鍵 8 第二節(jié)第二節(jié)

7、SASSAS軟件的使用基礎(chǔ)軟件的使用基礎(chǔ) 9 vEditor窗口窗口 主要用于打開SAS程序文件（*.sas）、編輯 和修改SAS程序、并提交全部或部分SAS程序。 vLOG窗口窗口 顯示有關(guān)的SAS會(huì)話和提交SAS程序的信息， 包括程序的出錯(cuò)信息等（*.log) vOutput窗口窗口 主要顯示提交SAS程序后的運(yùn)行結(jié)果（*.lst） vExplorer窗口窗口 主要用于查看和管理所有SAS文件，而且 可以對(duì)非SAS文件創(chuàng)建快捷方式。 vResults 窗口窗口 主要用于操作和管理提交SAS程序后的輸 出結(jié)果。他的內(nèi)容與Output窗口的內(nèi)容一一對(duì)應(yīng)，可以 看成是Output窗口內(nèi)容的名稱。

8、 二、二、SAS的工作窗口的工作窗口 10 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 一、一、SAS 語句語句 一個(gè)SAS語句就是要求SAS系統(tǒng)執(zhí)行某種操作 或給SAS系統(tǒng)提供一些信息的命令。 SAS語句通常由語句通常由SAS關(guān)鍵字、關(guān)鍵字、SAS名稱、運(yùn)算符及特名稱、運(yùn)算符及特 殊字符組成殊字符組成。 規(guī)則： （1）以 ；結(jié)尾 （2）以SAS關(guān)鍵字開始。如： data,proc,input,cards,model,class,if,keep,set,run 11 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 一、一、SAS 語句語句 （3）注釋語句可以用 *

9、開始，或者用/*/表示中間 內(nèi)容是注釋語句。 例1.1 Data sas1_1; /*sas1_1是數(shù)據(jù)集名*/ Input id name$ height weight; Bmi=weight/height*2; Cards; /*數(shù)據(jù)流開始*/ 1 judy 1.56 45 2 lucy 1.67 53 3 lili 1.58 50 ; /*數(shù)據(jù)流結(jié)束*/ Proc print data=sas1_1; Run; 12 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 一、一、SAS 語句語句 1.SAS關(guān)鍵字關(guān)鍵字 幾乎所有的SAS語句都是由SAS關(guān)鍵字開始的， 說明SAS語句

10、的類型。data語句，input語句，proc語句， run語句 2.SAS數(shù)據(jù)集名和變量名數(shù)據(jù)集名和變量名 （1）32個(gè)字符之內(nèi)， （2）第一個(gè)字母必須為字母或_；第二個(gè)以后可以為字母或 數(shù)字。 （3）字母不區(qū)分大小寫。 （4）不能使用空格和%￥#$等特殊字符 13 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 一、一、SAS 語句語句 3.SAS運(yùn)算符運(yùn)算符 包括比較運(yùn)算、算術(shù)運(yùn)算、邏輯運(yùn)算 （1）算數(shù)運(yùn)算符 + - * / * （2）比較運(yùn)算符 = = = 2 and y3 or/| 邏輯或 x2 or y3 Not/ 邏輯非 最大 35；結(jié)果為5 5；結(jié)果為3 | 連接

11、A=my name is; B=SAS; C=A|B; 那么C=“my name is SAS” 14 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 一、一、SAS 語句語句 4.SAS函數(shù)函數(shù) 函數(shù)的一般形式：函數(shù)的一般形式： 函數(shù)名（自變量，自變量，函數(shù)名（自變量，自變量，） (1)SAS常用的概率和密度函數(shù)常用的概率和密度函數(shù) 標(biāo)準(zhǔn)正態(tài)分布密度函數(shù)：標(biāo)準(zhǔn)正態(tài)分布密度函數(shù)：probnorm(x) 如：如：y=probnorm(-1.96); 結(jié)果為結(jié)果為0.025 t分布密度函數(shù)：分布密度函數(shù)：probt(x,df) 如：如：y1= probt(0.95,100);結(jié)果為結(jié)果

12、為0.82780 F分布密度函數(shù)：分布密度函數(shù)：probF(x,df1,df2) 2分布密度函數(shù)：分布密度函數(shù)：probCHI(x,df) 二項(xiàng)分布概率函數(shù)：二項(xiàng)分布概率函數(shù)：probBNML(p,n,m) 泊松分布概率函數(shù)：泊松分布概率函數(shù)：poisson(m,n) 15 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 正態(tài)分布分位數(shù)函數(shù)：正態(tài)分布分位數(shù)函數(shù)：probit(p) 0p1 如：如：y=probit (0.025), 結(jié)果為結(jié)果為-1.96 t分布分位數(shù)函數(shù)：分布分位數(shù)函數(shù)：TINV(p,df) F分布分位數(shù)函數(shù)：分布分位數(shù)函數(shù)：FINV(p,df1,df2) 2分布

13、分位數(shù)函數(shù)：分布分位數(shù)函數(shù)：CINV(p,df) 其它常用函數(shù)其它常用函數(shù) ABS(x)-取取X的絕對(duì)值的絕對(duì)值 SQRT(x)-計(jì)算計(jì)算X的平方根的平方根 MOD(x,y)-求求x/y的余數(shù)的余數(shù) 16 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 EXP(x)-計(jì)算計(jì)算e的的x次冪次冪 LOG(x)-對(duì)自變量對(duì)自變量x求以求以e為底的自然對(duì)數(shù)為底的自然對(duì)數(shù) LOG10(x)-對(duì)自變量對(duì)自變量x求以求以10為底的對(duì)數(shù)為底的對(duì)數(shù) CEIL(x)-取取自變量自變量x的最小整數(shù)。的最小整數(shù)。 如如 CEIL(5.7)=6 FLOOR(x)-取取自變量自變量x的最大整數(shù)。的最大整數(shù)。

14、如如 CEIL(5.7)=5 INT(x)-取取X的整數(shù)部分的整數(shù)部分 ROUND(x,n)-X按按n的精度取舍入值的精度取舍入值 17 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 data sas1_2; y1=probnorm(-1.96); /*標(biāo)準(zhǔn)正態(tài)分布密度函數(shù)標(biāo)準(zhǔn)正態(tài)分布密度函數(shù)*/ y2=probt(0.95,100); /* t分布密度函數(shù)分布密度函數(shù)*/ y3=probit(0.025); /*正態(tài)分布分位數(shù)函數(shù)正態(tài)分布分位數(shù)函數(shù)*/ x1=log10(100); /*對(duì)自變量對(duì)自變量x求以求以10為底的對(duì)數(shù)為底的對(duì)數(shù)*/ x2=exp(2); /*計(jì)算計(jì)算e

15、的的x次冪次冪*/ x3=CEIL(5.7); /*取取自變量自變量x的最小整數(shù)的最小整數(shù)*/ x4=ABS(-8); /*取取X的絕對(duì)值的絕對(duì)值*/ x5=SQRT(36); /*計(jì)算計(jì)算X的平方根的平方根*/ x6=MOD(10,3); /*求求x/y的余數(shù)的余數(shù)*/ x7=FLOOR(5.7); /*取取自變量自變量x的最大整數(shù)的最大整數(shù)*/ x8=INT(5.6874); /*取取X的整數(shù)部分的整數(shù)部分*/ x9=ROUND(5.4874,0.001); /* X按按n的精度取舍入值的精度取舍入值*/ A=my name is ; B=SAS; C=A|B; proc print da

16、ta=sas1_2; run; 18 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 二、二、SAS 程序程序 將一系列將一系列SAS語句按邏輯順序排列起來，構(gòu)成語句按邏輯順序排列起來，構(gòu)成SAS程程 序。通常包括數(shù)據(jù)步和過程步兩部分。序。通常包括數(shù)據(jù)步和過程步兩部分。 數(shù)據(jù)步數(shù)據(jù)步以以data語句開頭，以語句開頭，以run語句結(jié)束。主要作用是語句結(jié)束。主要作用是 建立數(shù)據(jù)集。建立數(shù)據(jù)集。 過程步過程步以以proc開頭，以開頭，以run語句結(jié)束。主要作用是激活語句結(jié)束。主要作用是激活 SAS過程對(duì)數(shù)據(jù)進(jìn)行處理和分析。過程對(duì)數(shù)據(jù)進(jìn)行處理和分析。 提交提交SAS程序：（程序：（1）工具

17、欄提交圖標(biāo)，）工具欄提交圖標(biāo)， （2）run/submit， （3）使用）使用F3功能鍵功能鍵 19 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 二、二、SAS 程序程序 程序的儲(chǔ)存：（程序的儲(chǔ)存：（1）鍵入）鍵入“file 路徑路徑+文件名文件名.擴(kuò)展名擴(kuò)展名” 如：如：file d:usersas1-1.sas; （2）“文件文件”菜單菜單“保存保存” 程序的調(diào)用：（程序的調(diào)用：（1）鍵入）鍵入infile路徑路徑+文件名文件名+擴(kuò)展名擴(kuò)展名 如：如：infile d:usersas1-1.sas; （2）“文件文件”菜單菜單“打開打開” 20 第三節(jié)第三節(jié) SASSAS語

18、言的語句和程序語言的語句和程序 例例1.3 SAS程序舉例程序舉例 data sas1_3; /*創(chuàng)建名為sas1_3的SAS數(shù)據(jù)集*/ Input id nam$ height weight; bmi=weight/height*2; cards; 1 Judy 1.56 45 /*數(shù)據(jù)流中每個(gè)數(shù)據(jù)值之間最少有1各空格*/ 2 Lucy 1.67 53 ; run; /*數(shù)據(jù)步結(jié)束*/ proc print data=sas1_3; /*顯示輸出sas1_3數(shù)據(jù)集到output窗口*/ run; /*過程步結(jié)束*/ 21 第二章第二章 建立建立SASSAS數(shù)據(jù)集數(shù)據(jù)集 第一節(jié)第一節(jié) SAS數(shù)

19、據(jù)集概述數(shù)據(jù)集概述 一、一、SAS數(shù)據(jù)集基本格式數(shù)據(jù)集基本格式 SAS數(shù)據(jù)集可以看做由若干行和若干列組成的表格，數(shù)據(jù)集可以看做由若干行和若干列組成的表格， 一個(gè)數(shù)據(jù)集等價(jià)于關(guān)系型數(shù)據(jù)庫中的一個(gè)表。一個(gè)數(shù)據(jù)集等價(jià)于關(guān)系型數(shù)據(jù)庫中的一個(gè)表。 idnameheightweight 1Judy15645 2Lucy16753 3Susan15655 4tony17870 變量 觀觀 測(cè)測(cè) 變量有兩種類型：變量有兩種類型： 數(shù)值型變量數(shù)值型變量 字符型變量字符型變量 22 第二章第二章 建立建立SASSAS數(shù)據(jù)集數(shù)據(jù)集 第一節(jié)第一節(jié) SAS數(shù)據(jù)集概述數(shù)據(jù)集概述 一、一、SAS數(shù)據(jù)集基本格式數(shù)據(jù)集基本格式

20、 例例2_1 缺失值的例子缺失值的例子 Data sas2_1; Input id name$ height weight; If name=999 then name=.; If weight=999 then weight=.; Cards; 1 Judy 156 999 2 Lucy 167 58 3 999 165 55 4 Tomy 178 62 ; Proc print data=sas2_1; run; 23 第二章第二章 建立建立SASSAS數(shù)據(jù)集數(shù)據(jù)集 二、臨時(shí)二、臨時(shí)SAS數(shù)據(jù)集和永久數(shù)據(jù)集和永久SAS數(shù)據(jù)集數(shù)據(jù)集 臨時(shí)數(shù)據(jù)集臨時(shí)數(shù)據(jù)集只能在只能在SAS會(huì)話過程中創(chuàng)建使用，

21、一旦退會(huì)話過程中創(chuàng)建使用，一旦退 出出SAS，數(shù)據(jù)集就不存在了。，數(shù)據(jù)集就不存在了。 在在Explorer窗口中可以看到一個(gè)名為窗口中可以看到一個(gè)名為work的的SAS數(shù)據(jù)數(shù)據(jù) 庫，主要存放庫，主要存放SAS臨時(shí)數(shù)據(jù)集。創(chuàng)建時(shí)可以用兩水平臨時(shí)數(shù)據(jù)集。創(chuàng)建時(shí)可以用兩水平 命名方式：如命名方式：如work.數(shù)據(jù)集名數(shù)據(jù)集名，也可以用，也可以用一水平的命一水平的命 名方式名方式：如：如數(shù)據(jù)集名。數(shù)據(jù)集名。 24 二、臨時(shí)二、臨時(shí)SAS數(shù)據(jù)集和永久數(shù)據(jù)集和永久SAS數(shù)據(jù)集數(shù)據(jù)集 永久性數(shù)據(jù)集永久性數(shù)據(jù)集儲(chǔ)存在硬盤里，因此，在以后的每一次儲(chǔ)存在硬盤里，因此，在以后的每一次 SAS會(huì)話中都可以在打開。在創(chuàng)

22、建和使用會(huì)話中都可以在打開。在創(chuàng)建和使用SAS永久性永久性 數(shù)據(jù)集之前，需要先建立一個(gè)數(shù)據(jù)集之前，需要先建立一個(gè)SAS數(shù)據(jù)庫來制定永久數(shù)據(jù)庫來制定永久 數(shù)據(jù)集存放的路徑。一個(gè)數(shù)據(jù)集存放的路徑。一個(gè)SAS數(shù)據(jù)庫相當(dāng)于硬盤上的數(shù)據(jù)庫相當(dāng)于硬盤上的 一個(gè)文件夾，可以將一個(gè)文件夾，可以將SAS數(shù)據(jù)集寫入或讀出。數(shù)據(jù)集寫入或讀出。 命名永久性數(shù)據(jù)集時(shí)必須用命名永久性數(shù)據(jù)集時(shí)必須用兩水平的命名方式兩水平的命名方式，即，即數(shù)數(shù) 據(jù)庫名據(jù)庫名.數(shù)據(jù)集名數(shù)據(jù)集名，數(shù)據(jù)庫名實(shí)際是一個(gè)，數(shù)據(jù)庫名實(shí)際是一個(gè)SAS數(shù)據(jù)庫的數(shù)據(jù)庫的 邏輯名。邏輯名。 第二章第二章 建立建立SASSAS數(shù)據(jù)集數(shù)據(jù)集 25 例例2_2 SA

23、S程序舉例程序舉例 libname tj d:user; /*創(chuàng)建名為tj的sas數(shù)據(jù)庫*/ data tj.sas2_2; /*創(chuàng)建永久性數(shù)據(jù)集sas2_2，保存在d:user 下*/ Input id name$ height weight; bmi=weight/height*2; cards; 1 Judy 1.56 45 2 Lucy 1.67 53 ; run; proc print data=tj.sas2_2; /*將保存在d:user的數(shù)據(jù) 集sas2_2顯示輸出到output窗口*/ run; 第二章第二章 建立建立SASSAS數(shù)據(jù)集數(shù)據(jù)集 26 第二章第二章 建立建立SA

24、SSAS數(shù)據(jù)集數(shù)據(jù)集 第二節(jié)第二節(jié) SAS數(shù)據(jù)集的建立數(shù)據(jù)集的建立 一、用數(shù)據(jù)步創(chuàng)建一、用數(shù)據(jù)步創(chuàng)建SAS數(shù)據(jù)集數(shù)據(jù)集 創(chuàng)建創(chuàng)建SAS數(shù)據(jù)集最主要的方式是使用數(shù)據(jù)集最主要的方式是使用Base SAS中的數(shù)據(jù)步。中的數(shù)據(jù)步。 一個(gè)一個(gè)SAS數(shù)據(jù)步由數(shù)據(jù)步由data語句開始，后面接著其他程序語句，語句開始，后面接著其他程序語句， 通過從原始數(shù)據(jù)文件中讀入數(shù)據(jù)或?qū)σ汛嬖诘耐ㄟ^從原始數(shù)據(jù)文件中讀入數(shù)據(jù)或?qū)σ汛嬖诘腟AS數(shù)據(jù)集進(jìn)數(shù)據(jù)集進(jìn) 行操作來建立新的行操作來建立新的SAS數(shù)據(jù)集。在數(shù)據(jù)集。在SAS程序中，用戶可以使程序中，用戶可以使 用不同的來源輸入數(shù)據(jù)。常用的數(shù)據(jù)輸入來源包括用不同的來源輸入數(shù)據(jù)。

25、常用的數(shù)據(jù)輸入來源包括： （1）原始數(shù)據(jù)以數(shù)據(jù)流的方式；）原始數(shù)據(jù)以數(shù)據(jù)流的方式； （2）原始數(shù)據(jù)以外部文件的方式；）原始數(shù)據(jù)以外部文件的方式； （3）已經(jīng)存在的）已經(jīng)存在的SAS數(shù)據(jù)集；數(shù)據(jù)集； （4）通過程序語句來建立數(shù)據(jù)集。）通過程序語句來建立數(shù)據(jù)集。 27 1、從數(shù)據(jù)流中讀入數(shù)據(jù)、從數(shù)據(jù)流中讀入數(shù)據(jù) 在在SAS程序中，用程序中，用cards語句定義數(shù)據(jù)流的開始。語句定義數(shù)據(jù)流的開始。 Input語句用于從數(shù)據(jù)流中讀入數(shù)據(jù)。語句用于從數(shù)據(jù)流中讀入數(shù)據(jù)。 28 例例2.3 從數(shù)據(jù)流讀入數(shù)據(jù)建立數(shù)據(jù)集從數(shù)據(jù)流讀入數(shù)據(jù)建立數(shù)據(jù)集 data sas2_3; /*創(chuàng)建名為sas2_3的SAS數(shù)據(jù)集

26、*/ Input id name$ height weight; /*從數(shù)據(jù)流中讀入數(shù)據(jù)*/ cards; 1 Judy 156 45 2 Lucy 167 53 3 Susan 165 55 4 tony 178 70 ; run; /*數(shù)據(jù)步結(jié)束*/ proc print data=sas2_3; /*顯示輸出sas2_3數(shù)據(jù)集到output 窗口*/ run; /*過程步結(jié)束*/ 29 2、從外部文件中讀入數(shù)據(jù)、從外部文件中讀入數(shù)據(jù) Infile語句用來指定外部文件，通常外部文件并不是語句用來指定外部文件，通常外部文件并不是 SAS文件。用文件。用input語句將外部文件中的數(shù)據(jù)讀入語句

27、將外部文件中的數(shù)據(jù)讀入 SAS中建立中建立SAS數(shù)據(jù)集。外部文件中的記錄可以按數(shù)據(jù)集。外部文件中的記錄可以按 一個(gè)變量一列的方式排列，也可以沒用格式，不安一個(gè)變量一列的方式排列，也可以沒用格式，不安 列方式排列。列方式排列。 30 例例2.4 從外部文件讀入數(shù)據(jù)從外部文件讀入數(shù)據(jù) 外部文本文件外部文本文件k:userstudent.txt 1 Judy 1.56 45 2 Lucy 1.67 53 3 Susan 1.65 55 4 tony 1.78 70 SAS程序：程序： Data sas2_4; Infile d:userstudent.txt; /*指定外部文件*/ Input id

28、 name$ height weight; /*從外部文件中讀取記錄*/ Run; Proc print data=sas2_4; Run; 31 3、從、從SAS數(shù)據(jù)集中讀入數(shù)據(jù)數(shù)據(jù)集中讀入數(shù)據(jù) 用用set語句可以從已有的語句可以從已有的SAS數(shù)據(jù)集中創(chuàng)建新的數(shù)據(jù)集。數(shù)據(jù)集中創(chuàng)建新的數(shù)據(jù)集。 運(yùn)行運(yùn)行set語句，語句，SAS會(huì)從已有的會(huì)從已有的SAS數(shù)據(jù)集中讀入觀數(shù)據(jù)集中讀入觀 測(cè)，包括原始數(shù)據(jù)集中的所有變量，還可以在原數(shù)測(cè)，包括原始數(shù)據(jù)集中的所有變量，還可以在原數(shù) 據(jù)集變量基礎(chǔ)上創(chuàng)建新的變量。據(jù)集變量基礎(chǔ)上創(chuàng)建新的變量。 32 例例2.5 從已有的從已有的SAS數(shù)據(jù)集中讀入數(shù)據(jù)數(shù)據(jù)集中讀入

29、數(shù)據(jù) Data sas2_5; set sas2_4; /*從已有的數(shù)據(jù)集sas2_4中讀入數(shù)據(jù)*/ Bmi=weight/height*2; /*生成一個(gè)新變量*/ Run; Proc print data=sas2_5; Run; 33 例例2.6 利用程序語句建立利用程序語句建立SAS數(shù)據(jù)集數(shù)據(jù)集 Data sas2_6; Do i=1 to 10; /*do與end語句之間的語句被重復(fù)執(zhí)行10次*/ Y=rannor(2008); /*產(chǎn)生來自標(biāo)準(zhǔn)正態(tài)分布的隨機(jī)數(shù)*/ Output; End; Proc print data=sas2_6; Run; 4、利用程序語句建立、利用程序語句

30、建立SAS數(shù)據(jù)集數(shù)據(jù)集 34 SAS中中Import/Export向?qū)Э捎糜趯⑵渌愊驅(qū)Э捎糜趯⑵渌?型的文件，如型的文件，如Excel、Access、dBase、文、文 本文件等轉(zhuǎn)為本文件等轉(zhuǎn)為SAS數(shù)據(jù)集。在數(shù)據(jù)集。在file菜單下點(diǎn)菜單下點(diǎn) 擊擊import/export，按向?qū)е甘静僮骷纯伞?，按向?qū)е甘静僮骷纯伞?二、二、Import/Export向?qū)?chuàng)建向?qū)?chuàng)建SAS數(shù)據(jù)集數(shù)據(jù)集 35 第三節(jié)第三節(jié) SASSAS數(shù)據(jù)集的修改數(shù)據(jù)集的修改 例例2.7 按觀測(cè)連接按觀測(cè)連接SAS數(shù)據(jù)集數(shù)據(jù)集 IdnameSexheightweight 4Tonymale17870 5Tommale17

31、065 IdnameSexheigh t weight 1Judyfemale15645 2Lucyfemale16753 3Susanfemale16555 female male 需要連接的SAS數(shù)據(jù)集 一、一、SAS數(shù)據(jù)集的合并數(shù)據(jù)集的合并 1、按觀測(cè)連接、按觀測(cè)連接SAS數(shù)據(jù)集數(shù)據(jù)集 36 第三節(jié)第三節(jié) SASSAS數(shù)據(jù)集的修改數(shù)據(jù)集的修改 例例2.7 連接連接SAS數(shù)據(jù)集的例子數(shù)據(jù)集的例子 Data sas2_7; Set tj.female tj.male; Run; Proc print data=sas2_7; Run; 使用使用set語句可以將兩個(gè)或多個(gè)數(shù)據(jù)集按觀測(cè)連成一語句

32、可以將兩個(gè)或多個(gè)數(shù)據(jù)集按觀測(cè)連成一 個(gè)大的數(shù)據(jù)集，這也是將新變量追加已有數(shù)據(jù)集上個(gè)大的數(shù)據(jù)集，這也是將新變量追加已有數(shù)據(jù)集上 的一種方式。的一種方式。 37 第三節(jié)第三節(jié) SASSAS數(shù)據(jù)集的修改數(shù)據(jù)集的修改 例例2.8 一對(duì)一匹配合并一對(duì)一匹配合并SAS數(shù)據(jù)集數(shù)據(jù)集 Idheightweight 115645 216753 316555 IdnameSexAge 1Judyfemale21 2Lucyfemale22 3Susanfemale19 student1student2 需要連接的SAS數(shù)據(jù)集 語句語句merge可以合并兩個(gè)或多個(gè)數(shù)據(jù)集中的變量，這可以合并兩個(gè)或多個(gè)數(shù)據(jù)集中的變量，

33、這 些數(shù)據(jù)集中的觀測(cè)是一一對(duì)應(yīng)的。些數(shù)據(jù)集中的觀測(cè)是一一對(duì)應(yīng)的。 2、按變量合并、按變量合并SAS數(shù)據(jù)集數(shù)據(jù)集 38 第三節(jié)第三節(jié) SASSAS數(shù)據(jù)集的修改數(shù)據(jù)集的修改 例例2.8 按變量合并按變量合并SAS數(shù)據(jù)集數(shù)據(jù)集 Data sas2_8; Merge student1 student2; Run; Proc print data=sas2_8; Run; 39 第三節(jié)第三節(jié) SASSAS數(shù)據(jù)集的修改數(shù)據(jù)集的修改 例例2.9 利用程序語句建立利用程序語句建立SAS數(shù)據(jù)集數(shù)據(jù)集 Idheightweight 115645 216753 417870 517765 IdnameSexAge

34、1Judyfemale21 2Lucyfemale22 3Susanfemale19 student3 student4 需要連接的SAS數(shù)據(jù)集 如果需要合并的數(shù)據(jù)集中的觀測(cè)并不是一一對(duì)應(yīng)的關(guān)系，則需如果需要合并的數(shù)據(jù)集中的觀測(cè)并不是一一對(duì)應(yīng)的關(guān)系，則需 要用要用by語句指定各數(shù)據(jù)集中用于匹配觀測(cè)的變量。（使用語句指定各數(shù)據(jù)集中用于匹配觀測(cè)的變量。（使用by語語 句之前，需要先對(duì)關(guān)鍵字段進(jìn)行排序。句之前，需要先對(duì)關(guān)鍵字段進(jìn)行排序。 40 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 例例2.9 按變量合并按變量合并SAS數(shù)據(jù)集數(shù)據(jù)集 Proc sort data student3

35、;by id; Proc sort data student4;by id; Data sas2_9; Merge student3 student4; Run; Proc print data=sas2_9; Run; 41 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 如果只想保留數(shù)據(jù)集中的部分觀測(cè)，有兩種如果只想保留數(shù)據(jù)集中的部分觀測(cè)，有兩種 方式建立啞集。方式建立啞集。 第一種可以用第一種可以用if 條件語句告訴條件語句告訴SAS需要保留需要保留 符合條件的觀測(cè)；另一種用符合條件的觀測(cè)；另一種用ifthan delete 語句告訴語句告訴SAS將不符合條件的觀測(cè)刪除。將不

36、符合條件的觀測(cè)刪除。 二、建立二、建立SAS數(shù)據(jù)集的啞集數(shù)據(jù)集的啞集 1、刪除觀測(cè)、刪除觀測(cè) 42 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 例例2.10 刪除觀測(cè)的例子刪除觀測(cè)的例子 IdnameSexAgeheightweight 1Judyfemale2115645 2Lucyfemale2216753 3Susanfemale1916555 4TonyMale2317870 5Tommale2517765 student 43 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 例例2.10 程序程序1 Data sas2_10; Set student; I

37、f sex=female; Run; Proc print data=sas2_10; Run; 44 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 例例2.10 程序程序2 Data sas2_10; Set student; If sex=female then delete; Run; Proc print data=sas2_10; Run; 45 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 有時(shí)并不需要數(shù)據(jù)集中的所有變量，可以有時(shí)并不需要數(shù)據(jù)集中的所有變量，可以 用用drop或或keep語句來刪除或保留數(shù)據(jù)集語句來刪除或保留數(shù)據(jù)集 中的變量，注意同一個(gè)數(shù)

38、據(jù)步中中的變量，注意同一個(gè)數(shù)據(jù)步中drop和和 keep語句不能同時(shí)出現(xiàn)。語句不能同時(shí)出現(xiàn)。 2、刪除變量、刪除變量 46 第三節(jié)第三節(jié) SASSAS語言的語句和程序語言的語句和程序 例例2.11 程序程序1 Data sas2_11; Set student; Keep id name sex age; Run; Proc print data=sas2_11; Run; 刪除變量的例子，僅保留數(shù)據(jù)集刪除變量的例子，僅保留數(shù)據(jù)集student中的部分變量中的部分變量 （id,name,sex,age）?？梢杂孟旅鎯蓚€(gè)）?？梢杂孟旅鎯蓚€(gè)SAS程序。程序。 47 第三節(jié)第三節(jié) SASSAS語言的

39、語句和程序語言的語句和程序 例例2.11 程序程序2 Data sas2_11; Set student; Drop height weight; Run; Proc print data=sas2_11; Run; 48 例例3.1 抽樣調(diào)查某地抽樣調(diào)查某地120名名1835歲健康男性居民血?dú)q健康男性居民血 清鐵含量（清鐵含量（mol/L），數(shù)據(jù)如下。試編制此血清鐵資），數(shù)據(jù)如下。試編制此血清鐵資 料的頻數(shù)分布表并描述該資料。料的頻數(shù)分布表并描述該資料。 第三章第三章 常用統(tǒng)計(jì)描述常用統(tǒng)計(jì)描述 7.42 8.65 23.02 21.61 21.31 21.46 9.97 22.73 14.9

40、4 20.18 21.62 23.07 20.38 8.40 17.32 29.64 19.69 21.69 23.90 17.45 19.08 20.52 24.14 23.77 18.36 23.04 24.22 24.13 21.53 11.09 18.89 18.26 23.29 17.67 15.38 18.61 14.27 17.40 22.55 17.55 16.10 17.98 20.13 21.00 14.56 19.89 19.82 17.48 14.89 18.37 19.50 17.08 18.12 26.02 11.34 13.81 10.25 15.94 15.83

41、 18.54 24.52 19.26 26.13 16.99 18.89 18.46 20.87 17.51 13.12 11.75 17.40 21.36 17.14 13.77 12.50 20.40 20.30 19.38 23.11 12.67 23.02 24.36 25.61 19.53 14.77 14.37 24.75 12.73 17.25 19.09 16.79 17.19 19.32 19.59 19.12 15.31 21.75 19.47 15.51 10.86 27.81 21.65 16.32 20.75 22.11 13.17 17.55 19.26 12.65

42、 18.48 19.83 23.12 19.22 19.22 16.72 27.90 11.74 24.66 14.18 16.52 49 第三章第三章 常用統(tǒng)計(jì)描述常用統(tǒng)計(jì)描述 一、定量資料的統(tǒng)計(jì)描述一、定量資料的統(tǒng)計(jì)描述 1、定量資料的頻數(shù)表、定量資料的頻數(shù)表 l先用先用MEANS過程找出最大值和最小值過程找出最大值和最小值 MIN=7.42 MAX=29.64 RANGE=22.22 l確定組距確定組距(I=2)和第一組段的下限和第一組段的下限(L=7) l然后建立新變量（然后建立新變量（x0），將變量），將變量X轉(zhuǎn)化為該數(shù)據(jù)所在組段的下限值轉(zhuǎn)化為該數(shù)據(jù)所在組段的下限值 賦值新的變量：賦

43、值新的變量：x0=INT(X-L)/I)*I+L 或：或：x0=x-mod(x-L,I) 或：x0=I*INT(x/I)+I/2 l用用FREQ過程作頻數(shù)表：過程作頻數(shù)表： proc freq; tables x0; Run; 50 第三章第三章 常用統(tǒng)計(jì)描述常用統(tǒng)計(jì)描述 2、簡(jiǎn)單描述簡(jiǎn)單描述：MEANS過程 輸出樣本含量、均數(shù)、標(biāo)準(zhǔn)差、最小值、最大值 ( N、Mean、Std Dev、Minimum、Maximum) 3、詳細(xì)描述詳細(xì)描述：UNIVARIATE過程 輸出更多指標(biāo)：樣本含量、均數(shù)、標(biāo)準(zhǔn)差、最小值、最大值、 方差、標(biāo)準(zhǔn)誤 、中位數(shù)、四分位數(shù)間距、常用百分位數(shù)、 眾數(shù)、平方和、離均

44、差平方和等，可以作正態(tài)性檢驗(yàn) 另外還顯示原數(shù)據(jù)中的部分極端值 51 第三章第三章 常用統(tǒng)計(jì)描述常用統(tǒng)計(jì)描述 PROC MEANS; CLASS SEX; VAR HEIGHT WEIGHT; RUN; CLASS 指定分組變量（可以不寫）指定分組變量（可以不寫） VAR 指定需要分析的變量指定需要分析的變量(不寫即分析全部變量）不寫即分析全部變量） MAXDEC=n 指定計(jì)算結(jié)果保留小數(shù)位數(shù)指定計(jì)算結(jié)果保留小數(shù)位數(shù)(默認(rèn)默認(rèn)2位位) 可加的選項(xiàng)有：標(biāo)準(zhǔn)誤(StdErr)、方差(Var)、總和(Sum)、平方 和(USS)、離均差平方和(CSS)、全距(Range)、變異系數(shù) (CV)、偏度系數(shù)

45、(Skewness)、峰度系數(shù)(Kurtosis)、可信區(qū)間 (CLM)等 還可以計(jì)算總體均數(shù)=0的 t 檢驗(yàn)(T、Prt) 不加任何選項(xiàng)時(shí)，自動(dòng)輸出5項(xiàng)結(jié)果： N、Mean、Std Dev、Minimum、Maximum 52 第三章第三章 常用統(tǒng)計(jì)描述常用統(tǒng)計(jì)描述 Data sas3_1;Input x; Gr=2; x0=gr*int(x/gr)+gr/2; /*計(jì)算組中值計(jì)算組中值X0*/ Cards; 7.42 8.65 23.02 21.6124.66 14.18 16.52 ; proc freq data=sas3_1; Tables x0/out=iron; /*out=ir

46、on表示輸出變量為X0的一維頻數(shù)分布表*/ Proc univariate data =sas3_1;Var x; Histogram x/midpoints=7 to 29 by 2; /*組中值的取值范圍從729,組距為2*/ Run; Proc means data=sas3_1;Var x; Run; 53 第三章第三章 常用統(tǒng)計(jì)描述常用統(tǒng)計(jì)描述 例例3.2隨機(jī)抽取隨機(jī)抽取15名學(xué)生，記錄他們的性別（名學(xué)生，記錄他們的性別（sex)、年齡、年齡 （age）、體重（）、體重（w,kg）、和身高（）、和身高（h,cm），求：），求： （1）按性別分組分別計(jì)算學(xué)生們身高和體重的均數(shù)、標(biāo)準(zhǔn)）按

47、性別分組分別計(jì)算學(xué)生們身高和體重的均數(shù)、標(biāo)準(zhǔn) 差、變異系數(shù)、中位數(shù)和差、變異系數(shù)、中位數(shù)和P50。 （2）計(jì)算學(xué)生身高）計(jì)算學(xué)生身高95%可信區(qū)間。要求結(jié)果保留小數(shù)點(diǎn)后可信區(qū)間。要求結(jié)果保留小數(shù)點(diǎn)后3 位。位。 F 15 46 156 f 14 41 149 f 18 65 165 m 15 50 160 m 13 48 155 m 18 70 180 m 14 38 150 m 16 55 165 m 17 68 176 m 16 60 170 f 17 50 160 f 17 58 160 f 16 60 165 m 17 65 175 f 18 61 162 54 Data sas3_2

48、; Input sex$ age w h ; Cards; f 15 46 156 f 14 41 149 f 18 65 165 m 15 50 160 m 13 48 155 m 18 70 180 m 14 38 150 m 16 55 165 m 17 68 176 m 16 60 170 f 17 50 160 f 17 58 160 f 16 60 165 m 17 65 175 f 18 61 162 ; Proc means data=sas3_2 mean std cv median p50 clm maxdec=3; /*選項(xiàng)maxdec指小數(shù)點(diǎn)后保留的位數(shù)*/ Var h

49、 w; Class sex; Run; 55 第三章第三章 常用統(tǒng)計(jì)描述常用統(tǒng)計(jì)描述 例例3.3 對(duì)對(duì)40名麻疹易感兒童接種麻疹疫苗一個(gè)月后，名麻疹易感兒童接種麻疹疫苗一個(gè)月后， 測(cè)得其血凝抑制抗體滴度如下，求其平均滴度。測(cè)得其血凝抑制抗體滴度如下，求其平均滴度。 56 Data sas3_3; /*計(jì)算幾何均數(shù)計(jì)算幾何均數(shù)*/ Input f x ; Y=log10(x); Cards; 1 4 5 8 6 16 2 32 7 64 10 128 4 256 5 512 ; Proc means noprint; /*選項(xiàng)noprint表示不打印程序運(yùn)行結(jié)果*/ Var y; Freq f;

50、 Output out=b mean=my; /*選項(xiàng)out=b和mean=my表示把該過程步運(yùn)行所得均數(shù)叫my,并 放入數(shù)據(jù)集b中*/ Data c; Set b; mx=10*my; Proc print data=c; Var mx; Run; 57 二、定性資料統(tǒng)計(jì)描述二、定性資料統(tǒng)計(jì)描述 FREQ 過程過程 可以產(chǎn)生頻數(shù)表和 n 維列聯(lián)表 PROC FREQ；TABLES A ；RUN； TABLES 指定要分析的定性變量 TABLES A 產(chǎn)生一維頻數(shù)表 TABLES A*B 產(chǎn)生二維列聯(lián)表(A是行B是列) TABLES A*B*C 產(chǎn)生三維列聯(lián)表 (A是分層，B是行，C是列) T

51、ABLES A*（B C） 產(chǎn)生A*B、A*C兩個(gè)二維表 TABLES X1-X3 產(chǎn)生X1、X2、X3三個(gè)一維表 TABLES (X1-X4)*Y 產(chǎn)生X1*Y X2*Y X3*Y X4*Y四個(gè)二維表 第三章第三章 常用統(tǒng)計(jì)描述常用統(tǒng)計(jì)描述 58 如描述以下三個(gè)定性變量（如描述以下三個(gè)定性變量（258人）人） A文化程度（文化程度（1=文盲文盲2=初中以下初中以下3=高中高中4=大學(xué)）大學(xué)） B工作強(qiáng)度（工作強(qiáng)度（1=輕輕 2=中中 3=重）重） C月收入月收入 （1=11 THEN G=2; ELSE G=1; CARDS; 0.84 1.05 1.20 1.20 1.39 1.53 1.

52、67 1.80 1.87 2.07 2.11 0.54 0.64 0.64 0.75 0.76 0.81 1.16 1.20 1.34 1.35 1.48 1.56 1.87 ; PROC TTEST; CLASS G; VAR X; RUN; 70 兩樣本均數(shù)比較用兩樣本均數(shù)比較用TTEST過程（方式二）過程（方式二） DATA sas4_4; INPUT X g; CARDS; 0.84 1 1.05 1 1.20 1 1.20 1 1.39 1 1.53 1 1.67 1 1.80 1 1.87 1 2.07 1 2.11 1 0.54 2 0.64 2 0.64 2 0.75 2 0.

53、76 2 0.81 2 1.16 2 1.20 2 1.34 2 1.35 2 1.48 2 1.56 2 1.87 2 ; PROC TTEST; CLASS g; VAR X; RUN; 71 結(jié)果：結(jié)果： Statistics Lower CL Upper CL Lower CL Upper CL Variable G N Mean Mean Mean Std Dev Std Dev Std Dev Std Err X 13 0.8295 1.0846 1.3397 0.3027 0.4221 0.6969 0.1171 1 X 11 1.2375 1.5209 1.8043 0.294

54、7 0.4218 0.7402 0.1272 2 X Diff (1-2) -0.795 -0.436 -0.078 0.3264 0.422 0.5973 0.1729 T-Tests Variable Method Variances DF t Value Pr |t| X Pooled Equal 22 -2.52 0.0193 X Satterthwaite Unequal 21.4 -2.52 0.0196 Equality of Variances Variable Method Num DF Den DF F Value Pr F X Folded F 12 10 1.00 1.

55、0000 72 兩大樣本均數(shù)比較時(shí)，如只有均數(shù)標(biāo)準(zhǔn)差，而無原始兩大樣本均數(shù)比較時(shí)，如只有均數(shù)標(biāo)準(zhǔn)差，而無原始 數(shù)據(jù)時(shí)，無法直接調(diào)用數(shù)據(jù)時(shí)，無法直接調(diào)用SAS軟件中的某個(gè)過程進(jìn)行分軟件中的某個(gè)過程進(jìn)行分 析，只能在數(shù)據(jù)步中按照計(jì)算公式編寫程序完成。析，只能在數(shù)據(jù)步中按照計(jì)算公式編寫程序完成。 調(diào)查健康人群紅細(xì)胞數(shù) (1012/L)資料如下： 男性360人，均數(shù)466.02，標(biāo)準(zhǔn)差57.46； 女性255人，均數(shù)417.80，標(biāo)準(zhǔn)差29.10。 試比較男女紅細(xì)胞數(shù)有無差別？ （方式三） 73 DATA BB; INPUT N1 X1 S1 N2 X2 S2; SE=SQRT(S1*S1/N1+S2

56、*S2/N2); U=ABS(X1-X2)/SE; P1=1-PROBNORM(U); P=2*P1; U=ROUND(U,0.001); P=ROUND(P,0.00001); IF P22表作表作Fisher確切概率確切概率 EXPECTED可以顯示每個(gè)格子的理論值可以顯示每個(gè)格子的理論值 CMH可以計(jì)算隊(duì)列研究和病例對(duì)照研究的可以計(jì)算隊(duì)列研究和病例對(duì)照研究的RR值值 NOROW、NOCOL、NOPERCENT 不顯示相應(yīng)百分不顯示相應(yīng)百分 數(shù)數(shù) NOPRINT 不顯示表格，只顯示計(jì)算的統(tǒng)計(jì)量不顯示表格，只顯示計(jì)算的統(tǒng)計(jì)量 97 例例6-2 某醫(yī)師為研究乙肝免疫球蛋白預(yù)防胎兒宮內(nèi)某醫(yī)師為研

57、究乙肝免疫球蛋白預(yù)防胎兒宮內(nèi) 感染感染HBV的效果，將的效果，將33例例HBsAg陽性孕婦隨機(jī)分為陽性孕婦隨機(jī)分為 預(yù)防注射組和非預(yù)防組，結(jié)果見表預(yù)防注射組和非預(yù)防組，結(jié)果見表7-4。問兩組新生。問兩組新生 兒的兒的HBV總體感染率有無差別？總體感染率有無差別？ 組別組別陽性陽性陰性陰性合計(jì)合計(jì) 感染率感染率 （%） 預(yù)防注射組預(yù)防注射組4182218.18 非預(yù)防組非預(yù)防組5(3)61145.45 合計(jì)合計(jì)9243327.27 四格表確切概率法四格表確切概率法 98 第六章第六章 2 2 檢驗(yàn) 檢驗(yàn) data sas6_2; input r c f ; cards; 1 1 4 1 2 18

58、 2 1 5 2 2 6 ; proc freq; weight f; tables r*c/chisq expected; run; 99 例例6-3 某實(shí)驗(yàn)室分別用乳膠凝集法和免疫熒光法對(duì)某實(shí)驗(yàn)室分別用乳膠凝集法和免疫熒光法對(duì)58名可疑系統(tǒng)名可疑系統(tǒng) 紅斑狼瘡患者血清中抗核抗體進(jìn)行測(cè)定，結(jié)果見表紅斑狼瘡患者血清中抗核抗體進(jìn)行測(cè)定，結(jié)果見表7-3。問兩種。問兩種 方法的檢測(cè)結(jié)果有無差別？方法的檢測(cè)結(jié)果有無差別？ 第六章第六章 2 2 檢驗(yàn) 檢驗(yàn) 二、配對(duì)四格表的二、配對(duì)四格表的2 檢驗(yàn)檢驗(yàn) 免疫熒光法免疫熒光法 乳膠凝集法乳膠凝集法 合計(jì)合計(jì) 11（a）12（b）23 2（c）33（d）35

59、 合計(jì)合計(jì)134558 100 data sas6_3; input r c f ; cards; 1 1 11 1 2 12 2 1 2 2 2 33 ; proc freq; weight f; tables r*c/chisq agree; run; 第六章第六章 2 2 檢驗(yàn) 檢驗(yàn) 二、配對(duì)四格表的二、配對(duì)四格表的2 檢驗(yàn)（方法一：適合檢驗(yàn)（方法一：適合b+c40） Tables r*c/agree 對(duì)列聯(lián)表作 MCNemar檢驗(yàn)和 一致性檢驗(yàn) 101 DATA b; INPUT f11 f12 f21 f22; IF f12+f2140 THEN chisq=(ABS(f12-f21

60、)-1)*2/(f12+f21); ELSE chisq=ABS(f12-f21)*2/(f12+f21); p=1-PROBCHI(chisq,1); CARDS; 11 12 2 33 ; PROC PRINT data=b; RUN; 第六章第六章 2 2 檢驗(yàn) 檢驗(yàn) 二、配對(duì)四格表的二、配對(duì)四格表的2 檢驗(yàn)（方法二）檢驗(yàn)（方法二） 102 療法 有效 無效 合計(jì) 有效率（%） 物理療法組 199 7 206 96.60 藥物治療組 164 18 182 90.11 外用膏藥組 118 26 144 81.94 合計(jì) 481 51 532 90.41 例6-4 某醫(yī)師研究物理療法、藥物治