醫(yī)學(xué)數(shù)據(jù)統(tǒng)計處理及SAS軟件的應(yīng)用

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

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

3、使用情況7SAS軟件功能進行數(shù)據(jù)管理數(shù)據(jù)輸入、建庫、保存 進行統(tǒng)計分析t檢驗方差分析卡方檢驗相關(guān)、回歸秩和檢驗8SAS的安裝 SAS8.12版本是windows操作系統(tǒng)下的32位軟件，和其他windows系統(tǒng)下的軟件安裝方法類似。用鼠標(biāo)左鍵雙擊光盤根目錄下的setup圖標(biāo)，然后按照提示安裝。9SAS的啟動 單擊開始菜單“程序”項。桌面上雙擊此快捷方式也可啟動SAS。 10執(zhí)行開始菜單中程序單項【程序】【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)： 提供編寫SAS程序的文本編輯器LOG窗口： 顯示執(zhí)行程序過程中的相關(guān)信息OUTPUT窗口： 顯示程序運行的結(jié)果 命令窗口與探索者窗口13程序編輯窗口日志窗口結(jié)果輸出窗口探索者窗口命令窗口14151617SAS程序的編寫 一個完整的SAS程序一般由數(shù)據(jù)步（DATA步）和過程步（PROC步）兩部分組成。DATA步以DATA語句開始，通過指定數(shù)據(jù)集名稱、定義變量名及類型和讀入數(shù)據(jù)等建立SAS數(shù)據(jù)集。過程步以PROC語句開始，通過調(diào)用不同的分析模塊來分析SAS數(shù)據(jù)集。18其它格式的數(shù)據(jù)庫直接導(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ù)，包括受檢者的姓名、性別、年齡、身高、體重信息。要求按性別分別計算身高、體重的描述性統(tǒng)計量：均值、標(biāo)準(zhǔn)差、最大值、最小值。該數(shù)據(jù)以Excel.xls文件形式儲存。文件名稱：tijiandata.xls2021Import Data22232425庫標(biāo)記導(dǎo)入的SAS數(shù)據(jù)庫文件名2627探索者窗口282930在Editor 程序編輯窗口輸入源程序數(shù)據(jù)步：將數(shù)據(jù)讀入系統(tǒng)，建立數(shù)據(jù)集 :DATA語句開始,后面加上數(shù)據(jù)庫的名字。INPUT語句描述數(shù)

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

7、一條語句可寫多行數(shù)據(jù)行后不加分號，數(shù)據(jù)最后單獨加分號行程序用“RUN；”語句結(jié)束36運行SAS程序方式一：鼠標(biāo)點擊執(zhí)行按鈕 方式二：按功能鍵“F8”方式三：執(zhí)行菜單命令 RunSubmit方式四：在命令框輸入SUBMIT 命令并執(zhí)行3738SAS程序?qū)嵗?例子 某學(xué)生上學(xué)期8門課程的分?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; 過程步 數(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門課程的平均成績?yōu)?4.625。40查看運行過程及出錯信息激活LOG窗口( F6 ) 黑色文字：源程序及運行情況藍(lán)色文字：程序運行情況的說明信息紅色文字：出錯信息綠色文字：警告信息4142保存SAS程序及其他窗口內(nèi)容激活相應(yīng)窗口，然后選擇： 方式一：鼠標(biāo)點擊“存盤”工具按鈕 方式二：執(zhí)行菜單命令【File】【Save】/ 【Save as】

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

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

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è)計資料的t-test采用國產(chǎn)藥物和進口藥物治療絕經(jīng)后婦女骨質(zhì)疏松癥，試比較兩種藥物的療效是否有差別545556完全隨機設(shè)計資料的方差分析（教科書73頁）安慰劑組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方差齊性檢驗59606162636465四格表資料的2-test組別有效無效合計有效率胞磷膽堿4665288.5%神經(jīng)節(jié)苷脂1882669.2%合計64147882.1% 某醫(yī)師欲比較胞磷膽堿與神經(jīng)節(jié)苷脂治療腦血管疾病的療效，將78例腦血管疾病患者隨機分為兩組，試比較兩種藥物的療效66676869配對四格表資料的2-test配對