1、醫(yī)學(xué)數(shù)據(jù)統(tǒng)計(jì)處理及SAS軟件的應(yīng)用1SAS軟件簡(jiǎn)介SAS是美國(guó)與西歐使用最為廣泛的三大著名統(tǒng)計(jì)分析軟件（SAS，SPSS和SYSTAT）之一，是目前國(guó)際上最為流行的一種大型統(tǒng)計(jì)分析系統(tǒng)，被譽(yù)為統(tǒng)計(jì)分析的標(biāo)準(zhǔn)軟件。是國(guó)際上公認(rèn)的最權(quán)威的統(tǒng)計(jì)軟件。以苛刻嚴(yán)格著稱于世的美國(guó)FDA新藥審批程序中,新藥試驗(yàn)結(jié)果的統(tǒng)計(jì)分析規(guī)定只能用SAS進(jìn)行,其他軟件的計(jì)算結(jié)果一律無(wú)效!哪怕只是簡(jiǎn)單的均數(shù)和標(biāo)準(zhǔn)差也不行!由此可見(jiàn)SAS的權(quán)威地位。 2004年開(kāi)始中國(guó)新藥臨床試驗(yàn)結(jié)果規(guī)定應(yīng)用SAS軟件分析。 2SAS與SPSS比較作為統(tǒng)計(jì)軟件大家族里的“傻瓜相機(jī)” ，SPSS的功能是最弱的，但易用性排在所有統(tǒng)計(jì)軟件之首。不

2、過(guò)SPSS的易用性也極大的限制了它的功能和可擴(kuò)展性。除了SPSS本身的菜單提供的功能，用戶無(wú)法通過(guò)編程或者插件來(lái)擴(kuò)展它的功能，而且SPSS提供的都只是最基本的常用的方法。3SAS全球?qū)I(yè)認(rèn)證SAS專業(yè)認(rèn)證是一項(xiàng)擁有極高國(guó)際聲譽(yù)的專業(yè)認(rèn)證，在歐美等國(guó)的職場(chǎng)上流行的一句話“If you have a SAS certification, You will never lose your job ”。銷售模式：只租不售 （每5年一個(gè)周期）財(cái)富500強(qiáng)中97%的企業(yè)都在使用SAS商業(yè)智能解決方案 452003年(301醫(yī)院與軍事醫(yī)科院博士生）各統(tǒng)計(jì)軟件使用情況6Pubmed檢索SCI英文文獻(xiàn)各統(tǒng)計(jì)軟件

3、使用情況7SAS軟件功能進(jìn)行數(shù)據(jù)管理數(shù)據(jù)輸入、建庫(kù)、保存 進(jìn)行統(tǒng)計(jì)分析t檢驗(yàn)方差分析卡方檢驗(yàn)相關(guān)、回歸秩和檢驗(yàn)8SAS的安裝 SAS8.12版本是windows操作系統(tǒng)下的32位軟件，和其他windows系統(tǒng)下的軟件安裝方法類似。用鼠標(biāo)左鍵雙擊光盤(pán)根目錄下的setup圖標(biāo)，然后按照提示安裝。9SAS的啟動(dòng) 單擊開(kāi)始菜單“程序”項(xiàng)。桌面上雙擊此快捷方式也可啟動(dòng)SAS。 10執(zhí)行開(kāi)始菜單中程序單項(xiàng)【程序】【The SAS System】【 The SAS System for Windows】 11SAS退出方式1：選擇關(guān)閉按鈕方式2：執(zhí)行菜單命令【File】【Exit】方式3：在命令框執(zhí)行BYE

4、或ENDSAS命令12SAS窗口PGM窗口（Program Editor)： 提供編寫(xiě)SAS程序的文本編輯器LOG窗口： 顯示執(zhí)行程序過(guò)程中的相關(guān)信息OUTPUT窗口： 顯示程序運(yùn)行的結(jié)果 命令窗口與探索者窗口13程序編輯窗口日志窗口結(jié)果輸出窗口探索者窗口命令窗口14151617SAS程序的編寫(xiě) 一個(gè)完整的SAS程序一般由數(shù)據(jù)步（DATA步）和過(guò)程步（PROC步）兩部分組成。DATA步以DATA語(yǔ)句開(kāi)始，通過(guò)指定數(shù)據(jù)集名稱、定義變量名及類型和讀入數(shù)據(jù)等建立SAS數(shù)據(jù)集。過(guò)程步以PROC語(yǔ)句開(kāi)始，通過(guò)調(diào)用不同的分析模塊來(lái)分析SAS數(shù)據(jù)集。18其它格式的數(shù)據(jù)庫(kù)直接導(dǎo)入轉(zhuǎn)換： excel(.xls)

5、; dbase(.dbf); lotus(.wk)等將數(shù)據(jù)讀入系統(tǒng)，建立數(shù)據(jù)集（直接輸入Editor窗口）。SAS數(shù)據(jù)集建立19外部數(shù)據(jù)導(dǎo)入例：有一批體檢數(shù)據(jù)，包括受檢者的姓名、性別、年齡、身高、體重信息。要求按性別分別計(jì)算身高、體重的描述性統(tǒng)計(jì)量：均值、標(biāo)準(zhǔn)差、最大值、最小值。該數(shù)據(jù)以Excel.xls文件形式儲(chǔ)存。文件名稱：tijiandata.xls2021Import Data22232425庫(kù)標(biāo)記導(dǎo)入的SAS數(shù)據(jù)庫(kù)文件名2627探索者窗口282930在Editor 程序編輯窗口輸入源程序數(shù)據(jù)步：將數(shù)據(jù)讀入系統(tǒng)，建立數(shù)據(jù)集 :DATA語(yǔ)句開(kāi)始,后面加上數(shù)據(jù)庫(kù)的名字。INPUT語(yǔ)句描述數(shù)

6、據(jù)變量 , 用$表示變量是字符型CARDS語(yǔ)句指示后面是數(shù)據(jù),每行為一個(gè)觀測(cè)數(shù)據(jù)，數(shù)據(jù)間無(wú)分號(hào)。 數(shù)據(jù)最后以“；”開(kāi)頭的空行結(jié)束 。 每個(gè)語(yǔ)句都以“；”結(jié)束。31SAS程序數(shù)據(jù)步格式Data ；Input $ .;Cards;. . . . . . . . . . . . . . . . . . . . . . . . . . . .;32數(shù)據(jù)步 過(guò)程步33SAS過(guò)程步對(duì)數(shù)據(jù)集中數(shù)據(jù)進(jìn)行分析、處理。以PROC語(yǔ)句開(kāi)始，調(diào)用SAS過(guò)程。每個(gè)語(yǔ)句都以“；”結(jié)束。整個(gè)程序最后要有RUN語(yǔ)句。34過(guò)程步35SAS程序的書(shū)寫(xiě)格式不分大小寫(xiě)單個(gè)名稱中不能空格，不同名稱間必須空格一行可寫(xiě)多條語(yǔ)句，用分號(hào)分割

7、一條語(yǔ)句可寫(xiě)多行數(shù)據(jù)行后不加分號(hào)，數(shù)據(jù)最后單獨(dú)加分號(hào)行程序用“RUN；”語(yǔ)句結(jié)束36運(yùn)行SAS程序方式一：鼠標(biāo)點(diǎn)擊執(zhí)行按鈕 方式二：按功能鍵“F8”方式三：執(zhí)行菜單命令 RunSubmit方式四：在命令框輸入SUBMIT 命令并執(zhí)行3738SAS程序?qū)嵗?例子 某學(xué)生上學(xué)期8門(mén)課程的分?jǐn)?shù)分別為88，90，87，76，92，78，80，86。試用SAS編程求其平均分?jǐn)?shù)。（1） 在PGM窗口鍵入下列程序data fenshu;input chengji;cards;88 90 87 76 92 78 80 86;proc means;var chengji;run; 過(guò)程步 數(shù)據(jù)步39查看結(jié)果 （

8、output 窗口）Analysis Variable : chengji N Mean Std Dev Minimum Maximum -8 84.6250000 5.8782286 76.0000000 92.0000000- 從結(jié)果看出，該學(xué)生上學(xué)期8門(mén)課程的平均成績(jī)?yōu)?4.625。40查看運(yùn)行過(guò)程及出錯(cuò)信息激活LOG窗口( F6 ) 黑色文字：源程序及運(yùn)行情況藍(lán)色文字：程序運(yùn)行情況的說(shuō)明信息紅色文字：出錯(cuò)信息綠色文字：警告信息4142保存SAS程序及其他窗口內(nèi)容激活相應(yīng)窗口，然后選擇： 方式一：鼠標(biāo)點(diǎn)擊“存盤(pán)”工具按鈕 方式二：執(zhí)行菜單命令【File】【Save】/ 【Save as】

9、可保存程序、數(shù)據(jù)文件、運(yùn)行結(jié)果、日志等。SAS保存成純文本格式文件（.txt)。43保存的文件的名44在SAS系統(tǒng)內(nèi)打開(kāi)已保存的文件只能在Program Edit窗口打開(kāi)文件，選擇：方式一：鼠標(biāo)點(diǎn)擊“打開(kāi)”工具按鈕 方式二：執(zhí)行菜單命令【File】【Open】 可打開(kāi)以前保存的程序、數(shù)據(jù)文件、運(yùn)行結(jié)果、日志等文件及其他文本格式的文件。45SAS程序主要關(guān)鍵詞 means; (n, mean, std, Q1, Q3, CV, median, clm) univariate normal; (正態(tài)性檢驗(yàn)，均數(shù)是否為0的 假設(shè)檢驗(yàn)：參數(shù)與非參數(shù)檢驗(yàn)） ttest; (兩個(gè)樣本均數(shù)比較的t檢驗(yàn)） a

10、nova; (各樣本含量相等的方差分析） glm; (各樣本含量不等的方差分析） npar1way; (非參數(shù)檢驗(yàn)） reg; (線性回歸方程：?jiǎn)卧蚨嘣?corr; (線性相關(guān)：線性相關(guān)或秩相關(guān)） freq; (卡方、Fisher檢驗(yàn)） logistic; (非線性回歸） lifetest; (生存分析） proc46定量資料的描述性分析47分析結(jié)果48定量資料的正態(tài)性檢驗(yàn)49正態(tài)性檢驗(yàn)結(jié)果50配對(duì)設(shè)計(jì)的t-test 為了比較兩種方法對(duì)乳酸飲料中脂肪的含量測(cè)定結(jié)果是否不同，隨機(jī)抽取了10份乳酸飲料制品，分別采用哥特里-羅紫法和脂肪酸水解法測(cè)定其結(jié)果，問(wèn)兩種方法檢測(cè)結(jié)果是否不同？編號(hào)哥特里-

11、羅紫法脂肪酸水解法10.840 0.580 20.591 0.509 30.674 0.500 40.632 0.316 50.687 0.337 60.978 0.517 70.750 0.454 80.730 0.512 91.200 0.997 100.870 0.506 515253成組設(shè)計(jì)資料的t-test采用國(guó)產(chǎn)藥物和進(jìn)口藥物治療絕經(jīng)后婦女骨質(zhì)疏松癥，試比較兩種藥物的療效是否有差別545556完全隨機(jī)設(shè)計(jì)資料的方差分析（教科書(shū)73頁(yè)）安慰劑組2.4g組4.8g組7.2g組3.53 2.42 2.86 0.89 4.59 3.36 2.28 1.06 4.34 4.32 2.39 1

12、.08 2.66 2.34 2.28 1.27 3.59 2.68 2.48 1.63 3.13 2.95 2.28 1.89 4.04 2.63 2.32 1.74 3.53 2.86 2.61 2.16 3.65 2.93 3.64 3.37 3.85 2.17 2.58 2.97 4.07 2.72 3.65 1.69 3.93 2.22 3.68 2.11 4.19 2.90 2.65 2.81 2.96 2.97 3.02 2.52 5758方差齊性檢驗(yàn)59606162636465四格表資料的2-test組別有效無(wú)效合計(jì)有效率胞磷膽堿4665288.5%神經(jīng)節(jié)苷脂1882669.2%合計(jì)64147882.1% 某醫(yī)師欲比較胞磷膽堿與神經(jīng)節(jié)苷脂治療腦血管疾病的療效，將78例腦血管疾病患者隨機(jī)分為兩組，試比較兩種藥物的療效66676869配對(duì)四格表資料的2-test配對(duì)