SAS统计分析系统操作手册idl.pptx

《SAS统计分析系统操作手册idl.pptx》由会员分享，可在线阅读，更多相关《SAS统计分析系统操作手册idl.pptx（124页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、SAS 统计分析系统统计分析系统Statistical Analysis System公共卫生学院信息数据处理教学实验室SAS软件简介软件简介1966年美国年美国North Carolina大学的两位大学的两位生物统计学研究生开始研制生物统计学研究生开始研制SAS系统，直系统，直至至1976年成立了年成立了SAS软件研究所，正式推软件研究所，正式推出了出了SAS软件，软件，1985年推出微机版年推出微机版SAS。SAS是用于决策支持的大型集成信息系是用于决策支持的大型集成信息系统，软件系统最早的功能限于统计分析，统，软件系统最早的功能限于统计分析，现在统计分析功能也仍是它的重要组成部现在统计分

2、析功能也仍是它的重要组成部分和核心功能。分和核心功能。SAS现在的最新版本为现在的最新版本为9.0版，大小约为版，大小约为1GB。在数据处理和统计分析领域，在数据处理和统计分析领域，SAS系系统被誉为国际上的标准软件系统，并在统被誉为国际上的标准软件系统，并在9697年度被评选为建立数据库的首选年度被评选为建立数据库的首选产品，堪称统计软件界的巨无霸。产品，堪称统计软件界的巨无霸。经过多年的发展，经过多年的发展，SAS已经被全世界已经被全世界120多个国家和地区的近三万家机构所采多个国家和地区的近三万家机构所采用，直接用户超过三百万人，遍及金融、用，直接用户超过三百万人，遍及金融、医药卫生、生

3、产、运输、通讯、政府和医药卫生、生产、运输、通讯、政府和教育科研等领域。教育科研等领域。由于由于SAS系统是从大型机上的系统发展系统是从大型机上的系统发展而来，在设计上也是完全针对专业用户进行而来，在设计上也是完全针对专业用户进行设计，因此操作以设计，因此操作以SAS语言编程为主，另外语言编程为主，另外也有非编程方式的菜单系统，但人机对话的也有非编程方式的菜单系统，但人机对话的界面不太友好。界面不太友好。学习学习SAS软件时，需要用户有基本的统软件时，需要用户有基本的统计学知识，对所选用的统计方法有较清楚的计学知识，对所选用的统计方法有较清楚的了解，没有统计学基础知识要掌握了解，没有统计学基础

4、知识要掌握SAS软件软件较为困难。较为困难。第一节第一节 SAS软件的使用基础软件的使用基础一、一、SAS的启动和退出的启动和退出 vv启动SASvv双击桌面SAS图标vv开始菜单 程序 The Sas System The Sas System for Windows V*.*vv直接双击SAS应用程序vv退出SASvv标题栏关闭按钮vv菜单栏 File Exitvv在命令框键入“Bye”或“Endsas”标题栏、菜单栏、工具栏标题栏、菜单栏、工具栏vPROGRAM EDITOR 程序编辑窗口vLOG 运行过程日志窗口vOUTPUT 结果输出窗口 SAS常用功能键常用功能键三个窗口间切换三个

5、窗口间切换F4 F4 RECALL (RECALL (调回程序调回程序)F5F5 PGM (PGM (程序程序)F8F8 SUBMIT (SUBMIT (发送程序发送程序)F6F6 LOG (LOG (日志日志)F1F1 HELP (HELP (帮助帮助)F7F7 OUTPUT (OUTPUT (结果结果)F9F9 KEYS (KEYS (功能键列表功能键列表)二、二、SAS的工作窗口的工作窗口例：一组体检资料，包括姓名、性别、年龄、身高、体重，要例：一组体检资料，包括姓名、性别、年龄、身高、体重，要求分性别计算身高、体重的均数、标准差、最大值、最小值。求分性别计算身高、体重的均数、标准差、最

6、大值、最小值。姓名姓名 性别性别 年龄年龄 身高身高 体重体重 ZHANGLIN F 47 156.3 47.1 ZHAOHUA M 38 172.4 61.5 WANGQANG M 41 169.2 64.5 LIULI F 52 158.2 53.6 SHIDONG F 39 160.1 48.0 KONGYING M 29 174.1 64.6 LILING F 37 152.3 42.2 GUANFEN M 32 166.2 60.2 MIAOQING M 26 180.3 74.3 NIUHUA M 31 178.2 80.2 TANGJIE F 35 169.2 60.8 GAOJ

7、UN M 24 176.0 73.3 SUNHONG F 27 158.3 49.9 三、三、SAS程序构成程序构成DATA AA1;INPUT NAME$SEX AGE HEIGHT WEIGHT;CARDS;ZHANGLIN 2 47 156.3 47.1 ZHAOHUA 1 38 172.4 61.5 WANGQANG 1 41 169.2 64.5 LIULI 2 52 158.2 53.6 SHIDONG 2 39 160.1 48.0 KONGYING 1 29 174.1 64.6 LILING 2 37 152.3 42.2 GUANFEN 1 32 166.2 60.2 MI

8、AOQING 1 26 180.3 74.3 NIUHUA 1 31 178.2 80.2 TANGJIE 2 35 169.2 60.8 GAOJUN 2 24 176.0 73.3 SUNHONG 2 27 158.3 49.9 ;PROC MEANS;CLASS SEX;VAR HEIGHT WEIGHT;RUN;数据步数据步数据步数据步过程步过程步过程步过程步v DATA STEP(数据步数据步)以以DATA语句开始语句开始将数据读入将数据读入SAS系统，建立系统，建立SAS数据集数据集v PROCEDURE STEP(过程步过程步)以以PROC语句语句开始开始调用各种已编好的调用各种

9、已编好的过程过程处理和分析数据处理和分析数据v RUN(结束结束)整个程序的最后必须要有整个程序的最后必须要有RUN语句结束，语句结束，表示要执行以上全部任务。表示要执行以上全部任务。1.必须在必须在PGM窗口。窗口。SAS语句书写格式自由，语句书写格式自由，每个语句以每个语句以“；”结束，结束，可以在一行中写几个语句或一个语句写在几行中，不区分字母大小写。2.运行运行SAS程序程序v工具栏执行按钮v功能键“F8”（Submit）v菜单栏 Locals Submit v命令框键入“Submit”命令自动激活OUTPUT窗口，出现程序运行结果四、四、SAS程序编写、运行、调试、保存程序编写、运行

10、、调试、保存3.调试调试SAS程序程序 没出现结果窗口，或结果不理想时。v激活LOG窗口（F6）查看运行过程及出错信息黑色文字：源程序及运行情况蓝色文字：程序运行情况说明信息红色文字：出错信息绿色文字：警告信息v激活PGM窗口（F5）v调回SAS程序（F4）修改错误后重新发送4.保存保存SAS程序或其他窗口内容程序或其他窗口内容 激活相应窗口激活相应窗口v常用工具栏 存盘按钮v菜单栏 File Save（或Save as）v命令框输入File命令可以分别保存程序、运行结果、日志5.打开已保存的程序、结果等打开已保存的程序、结果等 只能在只能在PGM窗口窗口v常用工具栏 打开按钮v菜单栏 Fil

11、e Openv命令框输入Include命令可打开以前保存的程序、运行结果、日志等文件及其他文本格式的文件。第二节第二节 SAS数据集数据集SAS数据集数据集(SAS Dataset)数据值被安排在一个矩阵式的表格结构中数据值被安排在一个矩阵式的表格结构中 每一行数据表示一个每一行数据表示一个观测观测（Observation）每一列数据表示一个每一列数据表示一个变量变量（Variable）变量命名：变量命名：18个字符，第个字符，第1个字符必须是字个字符必须是字母，后面母，后面 可以是字母、数字和下划线可以是字母、数字和下划线 变量类型：数值型、字符型和日期时间型变量类型：数值型、字符型和日期时

12、间型 字符型变量名后加字符型变量名后加“$”表示，表示，日期型日期型:16OCT9816OCT98 D，时间型，时间型:9:25:30 T 变量缺失值：用变量缺失值：用“”表示表示 一、一、SAS数据集的种类数据集的种类v临时SAS数据集不保存，一旦退出SAS就消失。数据集名称为1个单词，如DATA AA；v永久SAS数据集保存在用户指定的文件夹里。数据集名称为2个单词，中间用“.”分开，也称二级名，如DATA NEW.AA；LIBNAME NEW D:MY DATA;DATA NEW.AA;在数据集的二级名称 NEW.AA中；n n第1个词NEW是库标记，代表用户指定的保存SAS数据文件的位

13、置，前面要用LIBNAME语句指定它代表的文件夹和位置。n n第2个词AA是保存的SAS数据集的文件名。执行这个过程的结果将会在D:MY DATA文件夹中产生1个永久SAS数据文件 AA.SD2 SD2是自动产生的SAS数据文件默认扩展名二、创建二、创建SAS数据集数据集 1.在在PGM窗口创建数据集窗口创建数据集 DATA AA;命名将要创建的数据集名命名将要创建的数据集名 INPUT X1-X4；描述数据集中的变量名描述数据集中的变量名 CARDS;以下开始是数据以下开始是数据 2 47 156.3 47.1 1 38 172.4 61.5 ；数据结束数据结束 INPUT X YINPUT

14、 X Y；连续读数据不换行；连续读数据不换行 INPUT NAME$AGE;INPUT NAME$AGE;字符型变量名加字符型变量名加2.在在VIEWTABLE窗口创建数据集窗口创建数据集菜单栏Globals Manage Open tables打开VIEWTABLE窗口直接输入数据保存数据：File Save(或Save as)在Tables右边框中输入数据集名称在Path右边框中n n选默认位置*:sassaswork 建立临时数据集n n用户指定位置可以建立永久数据集3.用分析员模块创建数据集用分析员模块创建数据集 菜单栏 Globals Analyze Sta-tistical Ana

15、lysis打开Analyze(分析员)窗口直接输入数据保存数据：File Save Library：库标记名（代表要保存文件的 位置，选Work可以建临时数据集）Table Name：数据集文件名 Description：可以对文件进行描述4.由其它类型数据文件转换成由其它类型数据文件转换成SAS数据集数据集 菜单栏菜单栏 File importSAS可以直接与下列格式文件互相转换l ldBASE/FoxBASE File(*.dbf)l lExcel 4.0/5.0 File(*.xls)l lLotus 1/3/4 File(*.wk1/3/4)l lDelimited File(*.*)

16、空格分隔l lTab Delimited File(*.txt)制表符分隔l lComma Separated Velues(*.csv)逗号分隔 注意：注意：n nVisual Foxpro的数据表（*.dbf)在VFP系统中通过菜单栏 文件 导出 类型选 dBASE、FoxBASE、Excel 5.0、Delimited File 均可以n nExcel 97/2000/2002的文件（*.xls)在Excel系统中通过菜单栏 文件 另存为 保存类型选 Excel5.0/4.0(*.xls)或逗号分隔(*.csv)制表符分隔(*.txt)均可以n nExcel文件中要求第1行是变量名另外：

17、在写字板、记事簿、Word 或 Excel 中的数据也可以直接采用复制、粘贴的方式，通过PGM窗口创建SAS数据集。通过 菜单栏菜单栏 File Export也可以把SAS数据集转换成相应的其它类型数据文件 5.从已有的从已有的SAS数据集中创建新数据集数据集中创建新数据集 DATA A1(子集)；SET AA(原数据集)；(1)选择部分变量（列）选择部分变量（列）l l保留部分变量保留部分变量：KEEP DATA A1；SET AA；KEEP X1 X2；去掉部分变量去掉部分变量：DROP DATA A1；SET AA；DROP X3 X4；赋值新的变量赋值新的变量 Y=X3+X4；LY=L

18、OG(Y)；常用常用SAS函数：函数：LOG(x)以e为底的自然对数LOG10(x)以10为底的常用对数EXP(x)以e为底的反对数SQRT(x)平方根ABS(x)绝对值 INT(x)取整数(不四舍五入)MOD(x,y)取x/y的余数SIN(x)正弦COS(x)余弦 三角函数(x是弧度)TAN(x)正切(2)选择部分观测（行）选择部分观测（行）建子集语句：建子集语句：IF DATA A1；SET AA；IF SEX=1；DATA A1；SET AA；IF AGE40；IF X1=3 (不等于不等于)IF SEX=1 AND AGE60 IF AGE20 AND AGE=40 IF AGE40(

19、3)数据集连接数据集连接n n纵向连接 A1 DATA AA；SET A1 A2；A2 要求 A1和 A2中的变量相同 A1中有10行，A2中有15行，则AA有25行。n n横向合并 A1 A2 DATA AA；MERGE A1 A2；BY Num；匹配合并时，要求 A1和 A2中必须有1个相同的变量，且排列顺序也相同，作为标识变量用 BY 语句指定。(4)数据集转置（行列互换）数据集转置（行列互换）PROC TRANSPOSE；BY G；VAR X1-X4；OBS X1 X2 X3OBS X1 X2 X31 1 1 2 3 1 2 32 2 4 5 6 4 5 63 3 7 8 9 7 8

20、9OBS _Name_ Col1 Col2 Col3OBS _Name_ Col1 Col2 Col31 1X1 1 4 7X1 1 4 72 2X2 2 5 8X2 2 5 83 3X3 3 6 9X3 3 6 9选择项有：OUT=新数据集名（缺省按DATA n自动命名）PREFIX=新变量名前缀（缺省按Col n自动命名）NAME=为原变量名列命名（缺省自动命名为_Name_）第三节第三节 统计描述统计描述一、定量资料的统计描述一、定量资料的统计描述n n简单描述：MEANS过程 输出样本含量、均数、标准差、最小值、最大值输出样本含量、均数、标准差、最小值、最大值 (N (N、MeanMe

21、an、Std DevStd Dev、MinimumMinimum、Maximum)Maximum)l l详细描述：UNIVARIATE过程输出更多指标：样本含量、均数、标准差、最小输出更多指标：样本含量、均数、标准差、最小值、最大值、方差、标准误值、最大值、方差、标准误 、中位数、四分位数、中位数、四分位数间距、常用百分位数、众数、平方和、离均差平间距、常用百分位数、众数、平方和、离均差平方和等，可以作正态性检验方和等，可以作正态性检验另外还显示原数据中的部分极端值另外还显示原数据中的部分极端值PROC MEANS;CLASS SEX;VAR HEIGHT WEIGHT;RUN;CLASSCL

22、ASS 指定分组变量（可以不写）指定分组变量（可以不写）指定分组变量（可以不写）指定分组变量（可以不写）VAR VAR 指定需要分析的变量指定需要分析的变量指定需要分析的变量指定需要分析的变量(不写即分析全部变量）不写即分析全部变量）不写即分析全部变量）不写即分析全部变量）MAXDEC=nMAXDEC=n 指定计算结果保留小数位数指定计算结果保留小数位数指定计算结果保留小数位数指定计算结果保留小数位数(默认默认默认默认2 2位位位位)可加的选项有：标准误可加的选项有：标准误(StdErr)(StdErr)、方差、方差(Var)(Var)、总和、总和(Sum)(Sum)、平方和、平方和(USS)

23、(USS)、离均差平方和、离均差平方和(CSS)(CSS)、全距、全距(Range)(Range)、变异系数、变异系数(CV)(CV)、偏度系数、偏度系数(Skewness)(Skewness)、峰度、峰度系数系数(Kurtosis)(Kurtosis)、可信区间、可信区间(CLM)(CLM)等等还可以计算总体均数还可以计算总体均数=0=0的的 t t 检验检验(T(T、Prt)Prt)不加任何选项时，自动输出不加任何选项时，自动输出5 5项结果：项结果：N N、MeanMean、Std DevStd Dev、MinimumMinimum、MaximumMaximumMEANS过程过程频数表资

24、料：频数表资料：X-X-f f2.4-2.4-5 52.6-2.6-10102.8-2.8-12123.0-3.0-23233.2-3.2-24243.4-3.4-22223.6-3.6-17173.8-3.8-7 74.0-4.0-3 3DATA AA2;INPUT X F;CARDS;2.5 5 2.7 10 2.9 12 3.1 23 3.3 24 3.5 22 3.7 17 3.9 7 4.1 3 ；PROC MEANS N MEAN STD STDERR CLM；VAR X；FREQ F；RUN；数据步中：X 是组中值，F是频数过程步中：过程步中：VAR XVAR X；指定分析变量是

25、；指定分析变量是 X X FREQ F FREQ F；指定；指定 F F 是频数是频数计算几何均数计算几何均数滴度滴度人数人数1:41:4 1 11:81:8 5 51:161:16 6 61:321:32 2 21:641:64 7 71:1281:12810101:2561:256 4 41:5121:512 5 5合计合计4040DATA AA3;INPUT X F;Y=LOG10(X);CARDS;4 1 8 5 16 6 32 2 64 7 128 10 256 4 512 5 ;PROC MEANS;VAR Y;FREQ F;OUTPUT OUT=C MEAN=M;&DATA B2

26、;SET C;G=10*M;#PROC PRINT DATA=B2;VAR G;*RUN;&输出临时数据集C，变量M=均数，#新建数据集B2，取自C，求反对数G=10M*输出过程步，输出几何均数 GPROC SORT;BY SEX;PROC UNIVARIATE;BY SEX;VAR HEIGHT WEIGHT;RUN;BY 指定分组变量（可以不写）指定分组变量（可以不写）SORT 用用BY之前必须先按分组变量排序之前必须先按分组变量排序加选项 NORMAL 可以作正态性检验（W）在已知资料不服从正态分布时或对资料的分布类型不确定时可以用UNIVARIATE过程作探索性分析。UNIVARIAT

27、E过程过程MomentsMoments样本含量样本含量 N 7 N 7 频数之和频数之和 Sum Wgts 7 Sum Wgts 7均数均数 Mean 53.55714 Mean 53.55714 X X Sum 374.9Sum 374.9标准差标准差Std Dev 10.46564 Std Dev 10.46564 方差方差Variance 109.5295 Variance 109.5295 偏度系数偏度系数Skewness 1.246 Skewness 1.246 峰度系数峰度系数Kurtosis 1.376 Kurtosis 1.376 XX2 2 USS 20735.75 CSS

28、657.1771 USS 20735.75 CSS 657.1771变异系数变异系数 CV 19.54106 CV 19.54106 标准误标准误标准误标准误 Std Mean Std Mean 3.955638 3.955638 =0=0的的t t检验检验 T:Mean=0 13.53944 Pr|T|0.0001 T:Mean=0 13.53944 Pr|T|0.0001非非0 0变量个数变量个数 Num=0 7 Num=0 7 00变量个数变量个数 Num 0 7 Num 0 7符号检验符号检验 M(Sign)3.5 Pr=|M|0.0156 M(Sign)3.5 Pr=|M|0.015

29、6 符号秩和检验符号秩和检验 Sgn Rank 14 Pr=|S|0.0156 Sgn Rank 14 Pr=|S|0.0156正态性检验正态性检验 W:Normal 0.900687 PrW 0.3488 W:Normal 0.900687 PrW 0.3488Quantiles(Def=5)Quantiles(Def=5)P100P100 100%Max 73.3 100%Max 73.3 P99P99 99%73.3 99%73.3P75P75 75%Q3 60.8 75%Q3 60.8 P95 P95 95%73.3 95%73.3MM 50%Med 49.9 50%Med 49.9

30、P90P90 90%73.3 90%73.3P25P25 25%Q1 47.1 25%Q1 47.1 P10 P10 10%42.2 10%42.2P0P0 0%Min 42.2 0%Min 42.2 P5 P5 5%42.2 5%42.2 P1 P1 1%42.2 1%42.2全距全距 Range 31.1 Range 31.1四分位数间距四分位数间距 Q3-Q1 13.7 Q3-Q1 13.7众数众数 Mode 42.2 Mode 42.2 极端值Extremes Lowest Obs Highest Obs 最小 序号 最大 序号 42.2 (4)48 (3)47.1 (1)49.9 (

31、7)48 (3)53.6 (2)49.9 (7)60.8 (5)53.6 (2)73.3 (6)定量资料统计描述定量资料统计描述定量资料统计描述定量资料统计描述非编程方式（分析员模块）非编程方式（分析员模块）非编程方式（分析员模块）非编程方式（分析员模块）(1)打开分析员窗口菜单栏 Globals Analyze Statistical Analysis(2)打开SAS数据集菜单栏 File Open by SAS name(3)作简单统计分析菜单栏 Statistics Descriptive Summary Statistics在Analysis窗口选入要分析的变量在Classificat

32、ion窗口选入分组变量用用StatisticsStatistics可以选更多的统计分析结果可以选更多的统计分析结果用用VariablesVariables可以对变量有更多的选择方式可以对变量有更多的选择方式选选OKOK完成完成二、定性资料统计描述二、定性资料统计描述FREQ 过程过程 可以产生频数表和 n 维列联表PROC FREQ；TABLES A；RUN；TABLES TABLES 指定要分析的定性变量指定要分析的定性变量TABLES A TABLES A 产生一维频数表产生一维频数表TABLES A*B TABLES A*B 产生二维列联表产生二维列联表(A(A是行是行B B是列是列)T

33、ABLES A*B*C TABLES A*B*C 产生三维列联表产生三维列联表 (A (A是分层，是分层，B B是行，是行，C C是列是列)TABLES A*TABLES A*（B CB C）产生产生A*BA*B、A*CA*C两两个二维表个二维表TABLES X1-X3 TABLES X1-X3 产生产生X1X1、X2X2、X3X3三个一维表三个一维表TABLES (X1-X4)*Y TABLES (X1-X4)*Y 产生产生X1*Y X2*Y X3*Y X4*YX1*Y X2*Y X3*Y X4*Y 四个二维表四个二维表如描述以下三个定性变量（258人）A文化程度（1=文盲2=初中以下3=高

34、中4=大学）B工作强度（1=轻 2=中 3=重）C月收入 （1=|T|Mean Std Dev Std Error T Prob|T|236.67 586.60 151.46 1.56 0.1405 236.67 586.60 151.46 1.56 0.1405P0.05P0.05，不拒绝，不拒绝H H0 0，尚不能认为开展孕期保健指导对，尚不能认为开展孕期保健指导对新生儿体重有影响新生儿体重有影响 。二、配对设计二、配对设计 t 检验检验为了比较新旧两个测声计的测定效果，随机选定10个场地，在每地同时用新旧两个测声计测量数据如下，问新旧两个测声计的测量结果是否不同？场地 1 2 3 4 5

35、 6 7 8 9 10旧：87 65 74 95 65 55 63 88 61 54新：86 66 77 95 60 53 62 85 59 55 配对设计配对设计 t 检验用检验用MEANS过程过程DATA AA7;INPUT X1 X2;DATA AA7;INPUT X1 X2;Y=X1-X2Y=X1-X2;CARDS;CARDS;87 86 65 66 74 77 95 95 65 60 87 86 65 66 74 77 95 95 65 60 55 53 63 62 88 85 61 59 54 55 55 53 63 62 88 85 61 59 54 55 ;PROC MEANS

36、 MEAN STD STDERR T PRT;VAR Y;PROC MEANS MEAN STD STDERR T PRT;VAR Y;RUN;RUN;结果结果 Mean Std Dev Std Error T Prob|T|Mean Std Dev Std Error T Prob|T|0.90 2.28 0.72 1.25 0.2440 0.90 2.28 0.72 1.25 0.2440 P0.05 P0.05，不拒绝，不拒绝H H0 0，尚不能认为新旧两个测声计的测定，尚不能认为新旧两个测声计的测定结果有差别。结果有差别。三、两样本均数比较三、两样本均数比较试比较克山病人与健康人的血磷

37、值(mmol/L)是否不同？克山病人 0.84 1.05 1.20 1.20 1.39 1.53 1.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两样本均数比较用两样本均数比较用TTEST过程过程DATA AA8;INPUT X;IF _N_11 THEN G=1;ELSE G=0;CARDS;0.84 1.05 1.20 1.20 1.39 1.53 1.67 1.80 1.87 2.07 2.11 0.54 0.64 0.64 0.75 0.76 0.81

38、 1.16 1.20 1.34 1.35 1.48 1.56 1.87 ;PROC TTEST;CLASS G;VAR X;RUN;在TTEST过程中，分类变量G只能是两分类。条件语句 IF THEN/ELSEIF 后面是条件的表达式，如果满足条件就执行 THEN 后面操作，否则(不满足条件)，就执行 ELSE 后面的操作(ELSE可选)。IF X4=0 THEN X4=.;IF CJ2 AND X5=1 THEN Y=1;ELSE Y=2;_N_ 是SAS数据步中自定义的变量，它的值表示观测的序号。IF _N_11 THEN G=1;ELSE G=0;IF _N_11 THEN G=1;EL

39、SE G=0;也可用也可用G=(_N_11);G=(_N_11);赋值分组变量赋值分组变量G G：前：前1111个观测个观测G=0G=0，超过，超过1111时时G=1G=1。t 检验结果(1)G N Mean Std Dev Std Error 1 11 1.52 0.4218 0.1272 2 13 1.08 0.4221 0.1171(2)Variances方差 T DF Prob|T|Unequal 不齐 2.5239 21.4 0.0196 Equal 齐 2.5237 22.0 0.0193(3)For H0:Variances are equal,方差齐性检验 F=1.00 DF=

40、(12,10)ProbF=1.00结论：结论：P0.05P0.05，拒绝，拒绝H H0 0，接受，接受H H1 1，可以认为两组人群，可以认为两组人群血磷值不同，克山病人血磷值较高。血磷值不同，克山病人血磷值较高。两大样本均数比较时，如只有均数标准差，而无原始数据时，无法直接调用SAS软件中的某个过程进行分析，只能在数据步中按照计算公式编写程序完成。调查健康人群红细胞数(1012/L)资料如下：男性360人，均数466.02，标准差57.46；女性255人，均数417.80，标准差29.10。试比较男女红细胞数有无差别？DATA BB;INPUT N1 X1 S1 N2 X2 S2;SE=SQ

41、RT(S1*S1/N1+S2*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 P0.0001 THEN P=0.0001;FILE PRINT;PUT#2 10 U=25 P=;CARDS;360 466.02 57.46 255 417.80 29.10 ;RUN;定量资料假设检验定量资料假设检验定量资料假设检验定量资料假设检验非编程方式（分析员模块）非编程方式（分析员模块）非编程方式（分析员模块）非编程方式（分析员模块）菜单栏 Globals Analyze Stat

42、istical Analysis 菜单栏 File Open by SAS name(1)样本均数与总体均数比较样本均数与总体均数比较 菜单栏 Statistics Hypothesis Tasts One-Sample t-test for Means 往往VariableVariable窗口选入变量，窗口选入变量，Mean Mean 窗口输入窗口输入,OK,OKn n和己知时可选 One-Sample z-test for Means Mean Mean 窗口输入窗口输入,Std DevStd Dev 窗口输入窗口输入。(2)配对配对 t检验检验菜单栏 Statistics Hypothe

43、sis Tasts Two-Sample Paired t-test for Means 分别选入第一组和第二组变量分别选入第一组和第二组变量 ，OKOK(3)两样本均数比较两样本均数比较 t检验检验菜单栏 Statistics Hypothesis Tasts Two-Sample t-test for Means 往往DependentDependent窗口选入要分析的变量窗口选入要分析的变量 往往GroupGroup窗口选入分类变量，窗口选入分类变量，OKOK四、单因素方差分析四、单因素方差分析用于完全随机设用于完全随机设计（成组设计）计（成组设计）资料资料例：比较三组人群例：比较三组人

44、群载脂蛋白测定值有载脂蛋白测定值有无差别？无差别？糖尿病糖尿病 IGTIGT异常异常健康人健康人85.785.796.096.0144.0144.0105.2105.2124.5124.5117.0117.0109.5109.5105.1105.1110.0110.096.096.076.476.4109.0109.0115.2115.295.395.3103.0103.095.395.3110.0110.0123.0123.0110.0110.095.295.2127.0127.0106.5106.599.099.0121.0121.0125.6125.6120.0120.0115.0115

45、.0100.0100.0159.0159.0111.0111.0单因素方差分析用单因素方差分析用ANOVA过程过程DATA AA9;INPUT X;IF _N_20 THEN G=3;ELSE G=2;CARDS;85.7 105.2 109.5 96.0 115.2 95.3 ；PROC MEANS;CLASS G;VAR X;PROC ANOVA;CLASS G;MODEL X=G;MEANS G/SNK ;RUN;在SNK后面加HOVTEST可以作方差齐性检验Analysis of Variance ProcedureDependent Variable:XSource DF Sum o

46、f Squares Mean Square F Value Pr F变异来源 自由度 SS MS F PModel模型 2 2384.02550505 1192.01275253 5.85 0.0077Error误差 27 5497.83616162 203.62356154Total总 29 7881.86166667 R-Square R2 C.V.Root MSE Sy.x X Mean 0.302470 12.93519 14.26967279 110.31666667Source DF Anova SS Mean Square F Value Pr F变异来源 自由度 SS MS F

47、 P G组间 2 2384.02550505 1192.01275253 5.85 0.0077 P0.01，拒绝，拒绝H0，接受，接受H1，可以认为三组人群载脂蛋白不同。，可以认为三组人群载脂蛋白不同。Student-Newman-Keuls test for variable:X 均数间两两比较均数间两两比较q检验（检验（SNK法）法）Means with the same letter are not significantly different.相同字母表示的均数无显著性差别 SNK Grouping Mean N G A 122.800 10 3 B 105.455 11 1 B

48、B 102.389 9 2第三组高于前两组，即健康人的载脂蛋白高于两组病人。单因素方差分析时，如只有各组的均数、标准差，而无原始数据时，无法直接调用SAS软件中的某个过程进行分析，只能在数据步中按照计算公式编写程序完成。例：研究不同季节人体肾上腺皮质3-HSD活性变化季节季节N NMeansMeansStd Std 春季春季42420.780.780.130.13夏季夏季40400.690.690.220.22秋季秋季32320.680.680.140.14冬季冬季36360.580.580.200.20DATA AA;INPUT K XX;DO I=1 TO K;INPUT N X S;OU

49、TPUT;SS1+(XX-X)*(XX-X)*N;SS2+S*S*(N-1);NN+N;END;SS1=SS1/(K-1);SS2=SS2/(NN-K);F=SS1/SS2;P=1-PROBF(F,K-1,NN-K);F=ROUND(F,0.001);P=ROUND(p,0.00001);IF P FModel模型 9 2521.29375 280.14375 11.52 0.0001Error误差 14 340.54250 24.32446Total总 23 2861.83625 R-Square C.V.Root MSE X Mean 0.881006 8.860514 4.9319838

50、1 55.6625 Source DF Anova SS Mean Square F Value Pr F I 处理间 2 144.9175 72.45875 2.98 0.0836 J 区组间 7 2376.37625 339.48232 13.96 0.0001处理间处理间P0.05，不拒绝，不拒绝H0，尚不能认为三种营养素喂养小鼠增重有，尚不能认为三种营养素喂养小鼠增重有差别；区组间差别；区组间P FModel 9 6396.5625 710.7292 4.02 0.0521Error 6 1059.8750 176.6458Total 15 7456.4375 R-Square C.V