SAS软件在医学数据分析中应用

SAS软件简介

SAS是美国与西欧使用最为广泛的三大著名统计分析软件（SAS，SPSS和SYSTAT）之一，是目前国际上最为流行的一种大型统计分析系统，被誉为统计分析的标准软件。是国际上公认的最权威的统计软件。以苛刻严格著称于世的美国FDA新药审批程序中,新药试验结果的统计分析规定只能用SAS进行,其他软件的计算结果一律无效!哪怕只是简单的均数和标准差也不行!由此可见SAS的权威地位。2004年开始中国新药临床试验结果规定应用SAS软件分析。第一页，共81页。SAS与SPSS比较作为统计软件大家族里的“傻瓜相机”，SPSS的功能是最弱的，但易用性排在所有统计软件之首。不过SPSS的易用性也极大的限制了它的功能和可扩展性。除了SPSS本身的菜单提供的功能，用户无法通过编程或者插件来扩展它的功能，而且SPSS提供的都只是最基本的常用的方法。第二页，共81页。SAS全球专业认证SAS专业认证是一项拥有极高国际声誉的专业认证，在欧美等国的职场上流行的一句话“IfyouhaveaSAScertification,Youwillneverloseyourjob”。销售模式：只租不售（每5年一个周期）财富500强中97%的企业都在使用SAS商业智能解决方案

第三页，共81页。第四页，共81页。2003年(301医院与军事医科院博士生）各统计软件使用情况第五页，共81页。Pubmed检索SCI英文文献各统计软件使用情况第六页，共81页。SAS软件功能进行数据管理数据输入、建库、保存……进行统计分析t检验方差分析卡方检验相关、回归秩和检验……第七页，共81页。SAS窗口PGM窗口（ProgramEditor)：提供编写SAS程序的文本编辑器LOG窗口：显示执行程序过程中的相关信息OUTPUT窗口：显示程序运行的结果命令窗口与探索者窗口第八页，共81页。程序编辑窗口日志窗口结果输出窗口探索者窗口命令窗口第九页，共81页。第十页，共81页。第十一页，共81页。第十二页，共81页。SAS程序的编写一个完整的SAS程序一般由数据步（DATA步）和过程步（PROC步）两部分组成。DATA步以DATA语句开始，通过指定数据集名称、定义变量名及类型和读入数据等建立SAS数据集。过程步以PROC语句开始，通过调用不同的分析模块来分析SAS数据集。第十三页，共81页。其它格式的数据库直接导入转换：excel(.xls);dbase(.dbf);lotus(.wk)等将数据读入系统，建立数据集（直接输入Editor窗口）。SAS数据集建立第十四页，共81页。外部数据导入例：有一批体检数据，包括受检者的姓名、性别、年龄、身高、体重信息。要求按性别分别计算身高、体重的描述性统计量：均值、标准差、最大值、最小值。该数据以Excel.xls文件形式储存。文件名称：tijiandata.xls第十五页，共81页。第十六页，共81页。ImportData第十七页，共81页。第十八页，共81页。第十九页，共81页。第二十页，共81页。库标记导入的SAS数据库文件名第二十一页，共81页。第二十二页，共81页。探索者窗口第二十三页，共81页。第二十四页，共81页。第二十五页，共81页。在Editor程序编辑窗口输入源程序数据步：将数据读入系统，建立数据集:DATA语句开始,后面加上数据库的名字。INPUT语句描述数据变量,用$表示变量是字符型CARDS语句指示后面是数据,每行为一个观测数据，数据间无分号。数据最后以“；”开头的空行结束。每个语句都以“；”结束。第二十六页，共81页。SAS程序数据步格式Data<数据库名称>；Input<变量名1>$<变量名2>….<变量名K>@@;Cards;..............................;第二十七页，共81页。数据步

过程步第二十八页，共81页。SAS过程步对数据集中数据进行分析、处理。以PROC语句开始，调用SAS过程。每个语句都以“；”结束。整个程序最后要有RUN语句。第二十九页，共81页。过程步第三十页，共81页。SAS程序的书写格式不分大小写单个名称中不能空格，不同名称间必须空格一行可写多条语句，用分号分割一条语句可写多行数据行后不加分号，数据最后单独加分号行程序用“RUN；”语句结束第三十一页，共81页。运行SAS程序方式一：鼠标点击执行按钮方式二：按功能键“F8”方式三：执行菜单命令Run→Submit方式四：在命令框输入SUBMIT命令并执行第三十二页，共81页。第三十三页，共81页。SAS程序实例例子

某学生上学期8门课程的分数分别为88，90，87，76，92，78，80，86。试用SAS编程求其平均分数。（1）

在PGM窗口键入下列程序

datafenshu;inputchengji@@;cards;88

90

87

76

92

78

80

86;procmeans;varchengji;run;

过程步数据步第三十四页，共81页。查看结果（output窗口）AnalysisVariable:chengji

N

Mean

StdDev

Minimum

Maximum

-----------------------------------------------------------------------

8

84.6250000

5.8782286

76.0000000

92.0000000

-----------------------------------------------------------------------

从结果看出，该学生上学期8门课程的平均成绩为84.625。第三十五页，共81页。查看运行过程及出错信息激活LOG窗口(F6)

黑色文字：源程序及运行情况蓝色文字：程序运行情况的说明信息红色文字：出错信息绿色文字：警告信息第三十六页，共81页。第三十七页，共81页。保存SAS程序及其他窗口内容激活相应窗口，然后选择：方式一：鼠标点击“存盘”工具按钮方式二：执行菜单命令【File】→【Save】/【Saveas】可保存程序、数据文件、运行结果、日志等。SAS保存成纯文本格式文件（.txt)。第三十八页，共81页。保存的文件的名第三十九页，共81页。在SAS系统内打开已保存的文件只能在ProgramEdit窗口打开文件，选择：方式一：鼠标点击“打开”工具按钮方式二：执行菜单命令【File】→【Open】可打开以前保存的程序、数据文件、运行结果、日志等文件及其他文本格式的文件。第四十页，共81页。SAS程序主要关键词means;(n,mean,std,Q1,Q3,CV,median,clm…)

univariatenormal;(正态性检验，均数是否为0的假设检验：参数与非参数检验）ttest;(两个样本均数比较的t检验）anova;(各样本含量相等的方差分析）

glm;(各样本含量不等的方差分析）npar1way;(非参数检验）

reg;(线性回归方程：单元或多元）corr;(线性相关：线性相关或秩相关）freq;(卡方、Fisher检验）logistic;(非线性回归）

lifetest;(生存分析）proc第四十一页，共81页。定量资料的描述性分析第四十二页，共81页。分析结果第四十三页，共81页。定量资料的正态性检验第四十四页，共81页。正态性检验结果第四十五页，共81页。配对设计的t-test为了比较两种方法对乳酸饮料中脂肪的含量测定结果是否不同，随机抽取了10份乳酸饮料制品，分别采用哥特里-罗紫法和脂肪酸水解法测定其结果，问两种方法检测结果是否不同？编号哥特里-罗紫法脂肪酸水解法10.8400.58020.5910.50930.6740.50040.6320.31650.6870.33760.9780.51770.7500.45480.7300.51291.2000.997100.8700.506第四十六页，共81页。第四十七页，共81页。第四十八页，共81页。成组设计资料的t-test国产药进口药-5-176448634777607458-2-14897277486361705436228266-14采用国产药物和进口药物治疗绝经后妇女骨质疏松症，试比较两种药物的疗效是否有差别第四十九页，共81页。第五十页，共81页。第五十一页，共81页。完全随机设计资料的方差分析安慰剂组2.4g组4.8g组7.2g组3.532.422.860.894.593.362.281.064.344.322.391.082.662.342.281.273.592.682.481.633.132.952.281.894.042.632.321.743.532.862.612.163.652.933.643.373.852.172.582.974.072.723.651.693.932.223.682.114.192.902.652.812.962.973.022.52第五十二页，共81页。第五十三页，共81页。方差齐性检验第五十四页，共81页。第五十五页，共81页。第五十六页，共81页。第五十七页，共81页。第五十八页，共81页。第五十九页，共81页。第六十页，共81页。四格表资料的χ2-test组别有效无效合计有效率胞磷胆碱4665288.5%神经节苷脂1882669.2%合计64147882.1%某医师欲比较胞磷胆碱与神经节苷脂治疗脑血管疾病的疗效，将78例脑血管疾病患者随机分为两组，试比较两种药物的疗效第六十一页，共81页。第六十二页，共81页。第六十三页，共81页。第六十四页，共81页。配对四格表资料的χ2-test配对设计的卡方检验第六十五页，共81页。第六十六页，共81页。第六十七页，共81页。第六十八页，共81页。非参数检验（配对）第六十九页，共81页。第七十页，共81页。第七十一页，共81页。非参数检验9D11CDSC253255266366466477487510971210711