常用生物统计方法的sas 程

常用生物统计方法的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.2512.0021.3012.8031.2012.4031.009.8031.1511.6031.1010.6031.009.2031.4513.9031.3512.8031.159.3031.109.6031.2012.4031.0511.2031.1011.00;PROCGLM;CLASST;MODELY=TX/SOLUTION;MEANST/DUNCAN;LSMEANST/STDERRPDIFFTDIFF;RUN;