常用生物统计方法的程

常用生物统计方法的程第一页，共二十三页，编辑于2023年，星期六一、SAS系统简介

—SAS（StatisticalAnalysisSystem）是美国使用最为广泛的三大著名统计分析软件（SAS，SPSS，SYSTAT）之一；

—1966年美国NorthCarolina州立大学开始开发

—目前国际上最流行的一种大型统计分析系统

—统计分析的标准软件

—其基本部分是SAS/BASE软件，在此基础上可根据需要加上SAS系统的其他模块第二页，共二十三页，编辑于2023年，星期六

SAS/STAT（统计分析）；SAS/ETS（经济预测或时间序列分析）；SAS/QC（质量控制管理）；SAS/OR（运筹学）；SAS/IML（矩阵运算）；SAS/GRAPH（绘图）等

—信息存储简单（能读取任何形式的数据）

—语言编程能力强

—对数据连续处理

—统计分析方法丰富，使用简单

—SAS过程傻瓜系统……第三页，共二十三页，编辑于2023年，星期六二、SAS系统的启动与关闭

（一）启动

在WINDOWS操作系统中，直接用鼠标双击桌面上SAS系统的快捷键图标，即进入SAS系统。☼SAS系统软件有时间租期限制，因此只有机器时间在软件有效期内才能运行

OUTPUT（SAS结果输出窗口）

LOG（SAS日志窗口，随着SAS语句的执行，显示出SAS系统的信息和已执行的语句）

PGM（SAS程序编辑窗口）（二）退出

第四页，共二十三页，编辑于2023年，星期六四、SAS程序结构（一）程序结构在SAS系统中任何一个完整的处理过程均可分为两大步——数据步和过程步数据步——将不同来源的数据读入SAS系统建立起SAS数据集。由DATA语句开始，以RUN语句结束。

—数据可以由外部或过程读入；或者由CARDS语句开始，直接输入数据，数据输入完后提行单独加“;”结束第五页，共二十三页，编辑于2023年，星期六

过程步——调用SAS系统中已编好的各种过程来处理和分析数据集中的数据。

—每一个过程步均以PROC语句开始，RUN语句结束，并且每个语句后均以“;”结束。（二）程序的输入、修改调试和运行

SAS程序只能在PGM窗口输入、修改，可在各行的任何位置开始SAS语句的书写。一个语句可以连续写在几行中，一行中也可以同时写上几个语句，但每个语句后面必须用“；”隔开。第六页，共二十三页，编辑于2023年，星期六

当一个程序输入完后，发送程序的命令为F10功能键或SUBMIT。当程序发送到SAS系统后，PGM的程序语句全部自动清除，LOG窗口将逐步记下程序运行的过程和出现的错误信息（用红色提示错误）。

—如果过程步没有错误，运行完成后，通常会在OUTPUT窗口打印出结果

—如果程序运行出错，则需要在PGM窗口用RECALL（或F9）命令调回已发送的程序进行修改。第七页，共二十三页，编辑于2023年，星期六四、常用生物统计方法的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;第八页，共二十三页，编辑于2023年，星期六程序说明：样本平均数与总体平均数的差异显著性检验可调用MEANS过程。DATA语句产生临时数据集A，表明数据步的开始；INPUT语句指明读取变量y，@@表示读入一条观测值后不换行，连续读入数据，减少数据输入行；CARDS语句表明以下为数据行，数据行下的“；”表示数据行结束；PROCMEANS

语句指明调用MEANS过程对数据集A进行分析，输出样本含量N、平均数MEAN、平均数的标准误STDERR、学生氏T值和t值概率PRT；RUN语句表示过程步结束，开始运行过程步。第九页，共二十三页，编辑于2023年，星期六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;第十页，共二十三页，编辑于2023年，星期六3、非配对试验资料的t检验（教材【例5.3】）

DATAC;INPUTbreedy@@;CARDS;1572891120236110128211372501119239111723211042571732821532961682311118;

PROCTTEST;

CLASSbreed;VARy;RUN;第十一页，共二十三页，编辑于2023年，星期六（二）方差分析

对于一般的方差分析（平衡资料，即各处理重复数相等）可用ANOVA过程；对于非平衡资料（各处理重复数不等）的方差分析可用GLM过程。

第十二页，共二十三页，编辑于2023年，星期六1.单因素试验（各处理重复数相等，教材【例6.1】

）2.单因素试验（各处理重复数不等，教材【例6.3】

）3.两因素交叉分组单独观察值（教材【例6.4】

）4.两因素交叉分组重复观察值（教材【例6.5】

）5.两因素系统分组次级样本含量相等（教材【例6.6】

）6.随机单位组试验资料（教材【例12.3】）7.拉丁方试验资料（教材【例12.4】）第十三页，共二十三页，编辑于2023年，星期六教材【例6-1】

DATAE;INPUTTR$y@@;cards;A115A116A115A117A118A245A242A250A238A239A330A335A329A331A335A431A428A420A425A430A540A535A531A532A530;

PROCanova;

CLASSTR;

modely=TR;

MEANSTR/T;RUN;第十四页，共二十三页，编辑于2023年，星期六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;第十五页，共二十三页，编辑于2023年，星期六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;第十六页，共二十三页，编辑于2023年，星期六

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;第十七页，共二十三页，编辑于2023年，星期六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;第十八页，共二十三页，编辑于2023年，星期六DATAR;INPUTT$B$X@@;CARDS;A1B1205A1B2168A1B3222A1B4230A2B1230A2B2198A2B3242A2B4255A3B1252A3B2248A3B3305A3B4260A4B1200A4B2158A4B3183A4B4196A5B1265A5B2275A5B3315A5B4282;PROCGLM;CLASSTB;MODELX=TB;MEANST/DUNCAN;RUN;第十九页，共二十三页，编辑于2023年，星期六DATAL;INPUTT$ROWCOLX@@;CARDS;D1123E1221A1324B1421C1519A2122C2220E2320D2421B2522E3120A3225B3326C3422D3523B4125D4222C4325E4421A4523C5119B5220D5324A5422E5519;PROCGLM;CLASSTROWCOL;MODELX=TROWCOLl;MEANST/DUNCAN;RUN;第二十页，共二十三页，编辑于2023年，星期六（三）一元线性回归与相关分析（教材【例8.1】）

DATAG；INPUTx

y@@；CARDS；

802350862400982720902500120315010226809526308324001133080105292011029601002860；PROC

REG

CORR；MODELy=x

；RUN；第二十一页，共二十三页，编辑于2023年，星期六（四）协方差分析（教材【例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.35