SAS软件教程90603.ppt

常用生物统计方法的SAS程序 一 SAS系统简介二 SAS系统的启动与关闭三 SAS程序结构 程序的输入 修改调试和运行四 常用生物统计方法的SAS程序 t检验 方差分析 线性回归与相关分析 一 SAS系统简介 SAS StatisticalAnalysisSystem 是美国使用最为广泛的三大著名统计分析软件 SAS SPSS SYSTAT 之一 1966年美国NorthCarolina州立大学开始开发 目前国际上最流行的一种大型统计分析系统 统计分析的标准软件 其基本部分是SAS BASE软件 在此基础上可根据需要加上SAS系统的其他模块 SAS STAT 统计分析 SAS ETS 经济预测或时间序列分析 SAS QC 质量控制管理 SAS OR 运筹学 SAS IML 矩阵运算 SAS GRAPH 绘图 等 信息存储简单 能读取任何形式的数据 语言编程能力强 对数据连续处理 统计分析方法丰富 使用简单 SAS过程傻瓜系统 二 SAS系统的启动与关闭 一 启动在WINDOWS操作系统中 直接用鼠标双击桌面上SAS系统的快捷键图标 即进入SAS系统 SAS系统软件有时间租期限制 因此只有机器时间在软件有效期内才能运行OUTPUT SAS结果输出窗口 LOG SAS日志窗口 随着SAS语句的执行 显示出SAS系统的信息和已执行的语句 PGM SAS程序编辑窗口 二 退出 四 SAS程序结构 一 程序结构在SAS系统中任何一个完整的处理过程均可分为两大步 数据步和过程步数据步 将不同来源的数据读入SAS系统建立起SAS数据集 由DATA语句开始 以RUN语句结束 数据可以由外部或过程读入 或者由CARDS语句开始 直接输入数据 数据输入完后提行单独加 结束 过程步 调用SAS系统中已编好的各种过程来处理和分析数据集中的数据 每一个过程步均以PROC语句开始 RUN语句结束 并且每个语句后均以 结束 二 程序的输入 修改调试和运行SAS程序只能在PGM窗口输入 修改 可在各行的任何位置开始SAS语句的书写 一个语句可以连续写在几行中 一行中也可以同时写上几个语句 但每个语句后面必须用 隔开 当一个程序输入完后 发送程序的命令为F10功能键或SUBMIT 当程序发送到SAS系统后 PGM的程序语句全部自动清除 LOG窗口将逐步记下程序运行的过程和出现的错误信息 用红色提示错误 如果过程步没有错误 运行完成后 通常会在OUTPUT窗口打印出结果 如果程序运行出错 则需要在PGM窗口用RECALL 或F9 命令调回已发送的程序进行修改 四 常用生物统计方法的SAS程序 一 t检验1 样本平均数与总体平均数的差异显著性检验 教材 例5 1 DATAA INPUTy y y 4 50 cards 4 334 623 894 144 784 644 524 484 554 26 PROCMEANSNMEANSTDERRTPRT RUN 程序说明 样本平均数与总体平均数的差异显著性检验可调用MEANS过程 DATA语句产生临时数据集A 表明数据步的开始 INPUT语句指明读取变量y 表示读入一条观测值后不换行 连续读入数据 减少数据输入行 CARDS语句表明以下为数据行 数据行下的 表示数据行结束 PROCMEANS语句指明调用MEANS过程对数据集A进行分析 输出样本含量N 平均数MEAN 平均数的标准误STDERR 学生氏T值和t值概率PRT RUN语句表示过程步结束 开始运行过程步 2 配对试验资料的t检验 教材 例5 5 DATAB INPUTIDx1x2 d x1 x2 CARDS 137 837 9238 239 0338 038 9437 638 4537 937 9638 139 0738 239 5837 538 6938 538 81037 939 0 PROCMEANSMEANSTDERRTPRT VARd RUN 3 非配对试验资料的t检验 教材 例5 3 DATAC INPUTbreedy CARDS 1572891120236110128211372501119239111723211042571732821532961682311118 PROCTTEST CLASSbreed VARy RUN 二 方差分析 对于一般的方差分析 平衡资料 即各处理重复数相等 可用ANOVA过程 对于非平衡资料 各处理重复数不等 的方差分析可用GLM过程 1 单因素试验 各处理重复数相等 教材 例6 1 2 单因素试验 各处理重复数不等 教材 例6 3 3 两因素交叉分组单独观察值 教材 例6 4 4 两因素交叉分组重复观察值 教材 例6 5 5 两因素系统分组次级样本含量相等 教材 例6 6 6 随机单位组试验资料 教材 例12 3 7 拉丁方试验资料 教材 例12 4 教材 例6 1 DATAE INPUTTR y cards A115A116A115A117A118A245A242A250A238A239A330A335A329A331A335A431A428A420A425A430A540A535A531A532A530 PROCanova CLASSTR modely TR MEANSTR T RUN DATAF INPUTBREED y cards B121 5B119 5B120 0B122 0B118 0B120 0B216 0B218 5B217 0B215 5B220 0B216 0B319 0B317 5B320 0B318 0B317 0B421 0B418 5B419 0B420 0B515 5B518 0B517 0B516 0 PROCGLM CLASSBREED modely BREED MEANSBREED T RUN DATAF INPUTA B y cards A1B13 5A1B22 3A1B32 0A2B12 6A2B22 0A2B31 9A3B12 0A3B21 5A3B31 2A4B11 4A4B20 8A4B30 3 PROCANOVA CLASSAB modely AB MEANSAB DUNCAN RUN DATACaP DOa 1to4 DOb 1to4 DOn 1to3 inputy output end end end DROPn CARDS 22 026 524 430 027 526 032 426 527 030 527 025 123 525 827 033 228 530 138 035 533 926 524 025 030 526 825 536 534 033 528 030 524 620 522 519 534 531 429 329 027 528 027 526 328 518 520 019 0 PROCANOVA CLASSAB MODELy ABA B RUN DATAG DOa 1to3 DOb 1to3 DOn 1to2 inputy output end end end DROPn CARDS 82 582 487 186 584 083 986 685 886 285 787 087 682 081 580 080 579 580 3 PROCANOVA CLASSAB MODELy AB A RUN DATAR INPUTT B X CARDS A1B1205A1B2168A1B3222A1B4230A2B1230A2B2198A2B3242A2B4255A3B1252A3B2248A3B3305A3B4260A4B1200A4B2158A4B3183A4B4196A5B1265A5B2275A5B3315A5B4282 PROCGLM CLASSTB MODELX TB MEANST DUNCAN RUN DATAL INPUTT ROWCOLX CARDS D1123E1221A1324B1421C1519A2122C2220E2320D2421B2522E3120A3225B3326C3422D3523B4125D4222C4325E4421A4523C5119B5220D5324A5422E5519 PROCGLM CLASSTROWCOL MODELX TROWCOLl MEANST DUNCAN RUN 三 一元线性回归与相关分析 教材 例8 1 DATAG INPUTxy CARDS 802350862400982720902500120315010226809526308324001133080105292011029601002860 PROCREGCORR MODELy x RUN 四 协方差分析 教材 例10 1 DATAK INPUTT XY CARDS CK1 5012 40CK1 8512 00CK1 3510 80CK1 4510 00CK1 4011 00CK1 4511 80CK1 5012 50CK1 5513 40CK1 4011 20CK1 5011 60CK1 6012 60CK1 7012 5011 3510 2011 209 4011 4512 2011 2010 3011 4011 3011 3011 4011 1512 8011 3010 9011 3511 6011 158 5011 3512 2011 209 3021 1510 0021 1010 6021 1010 4021 059 2021 4013 0021 4513 5021 3013 0021 7014 8021 4012 3021 4513 2021 25