常用生物统计方法的SAS程序

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,常用生物统计方法的,SAS,程序,一、,SAS,系统简介,二、,SAS,系统的启动与关闭,三、,SAS,程序结构、程序的输入、修改 调试和运行,四、常用生物统计方法的,SAS,程序,t,检验,方差分析,线性回归与相关分析,一、,SAS,系统简介,SAS,（,Statistical Analysis System,）,是美国使用最为广泛的,三大著名统计分析软件（,SAS,，,SPSS,，,SYSTAT,）,之一；,1966,年美国,North Carolina,州立大学开始开发,目前国际上最流行的一种,大型统计分析系统,统计分析的,标准软件,其基本部分是,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】,),DATA,A;,INPUT y;,y=y-4.50;,cards,;,4.33 4.62 3.89 4.14 4.78 4.64 4.52 4.48 4.55 4.26,;,PROC MEANS,N MEAN STDERR T PRT;,RUN;,程序说明：,样本平均数与总体平均数的差异显著性检验可调用,MEANS,过程。,DATA,语句,产生临时数据集,A,，,表明数据步的开始；,INPUT,语句,指明读取变量,y,，,表示读入一条观测值后不换行，连续读入数据，减少数据输入行；,CARDS,语句,表明以下为数据行，数据行下的“；”表示数据行结束；,PROC MEANS,语句指明调用,MEANS,过程对数据集,A,进行分析，输出样本含量,N,、,平均数,MEAN,、,平均数的标准误,STDERR,、,学生氏,T,值和,t,值概率,PRT,；,RUN,语句,表示过程步结束，开始运行过程步。,2,、,配对试验资料的,t,检验,（教材,【,例,5.5】,）,DATA B;,INPUT ID x1 x2;,d=x1-x2;,CARDS;,1 37.8 37.9,2 38.2 39.0,3 38.0 38.9,4 37.6 38.4,5 37.9 37.9,6 38.1 39.0,7 38.2 39.5,8 37.5 38.6,9 38.5 38.8,10 37.9 39.0,;,PROC MEANS,MEAN STDERR T PRT;,VAR d;,RUN;,3,、,非配对试验资料的,t,检验,（教材,【,例,5.3】,）,DATA C;,INPUT breed y;,CARDS,;,1 57 2 89 1 120 2 36 1 101 2 82,1 137 2 50 1 119 2 39 1 117 2 32,1 104 2 57 1 73 2 82 1 53 2 96,1 68 2 31 1 118,;,PROC TTEST,;,CLASS,breed;,VAR y;,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】,DATA E;,INPUT TR$y;,cards;,A1 15 A1 16 A1 15 A1 17 A1 18,A2 45 A2 42 A2 50 A2 38 A2 39,A3 30 A3 35 A3 29 A3 31 A3 35,A4 31 A4 28 A4 20 A4 25 A4 30,A5 40 A5 35 A5 31 A5 32 A5 30,;,PROC,anova,;,CLASS,TR;,model,y=TR;,MEANS,TR/T;,RUN;,DATA F;,INPUT BREED$y;,cards;,B1 21.5 B1 19.5 B1 20.0 B1 22.0 B1 18.0 B1 20.0,B2 16.0 B2 18.5 B2 17.0 B2 15.5 B2 20.0 B2 16.0,B3 19.0 B3 17.5 B3 20.0 B3 18.0 B3 17.0,B4 21.0 B4 18.5 B4 19.0 B4 20.0,B5 15.5 B5 18.0 B5 17.0 B5 16.0,;,PROC GLM,;,CLASS,BREED;,model,y=BREED;,MEANS BREED/T;,RUN;,DATA F;,INPUT A$B$y;,cards;,A1 B1 3.5 A1 B2 2.3 A1 B3 2.0,A2 B1 2.6 A2 B2 2.0 A2 B3 1.9,A3 B1 2.0 A3 B2 1.5 A3 B3 1.2,A4 B1 1.4 A4 B2 0.8 A4 B3 0.3,;,PROC ANOVA,;,CLASS A B;,model y=A B,;,MEANS A B/DUNCAN;,RUN;,DATA,CaP,;,DO a=1 to 4;,DO b=1 to 4;,DO n=1 to 3;,input y;,output;,end;,end;,end;,DROP,n;,CARDS;,22.0 26.5 24.4 30.0 27.5 26.0 32.4 26.5 27.0,30.5 27.0 25.1 23.5 25.8 27.0 33.2 28.5 30.1,38.0 35.5 33.9 26.5 24.0 25.0 30.5 26.8 25.5,36.5 34.0 33.5 28.0 30.5 24.6 20.5 22.5 19.5,34.5 31.4 29.3 29.0 27.5 28.0 27.5 26.3 28.5,18.5 20.0 19.0,;,PROC ANOVA,;CLASS A B;,MODEL y=A B A*B;,RUN;,DATA G;,DO a=1 to 3;,DO b=1 to 3;,DO n=1 to 2;,input y;,output;,end;,end;,end;,DROP,n;,CARDS;,82.5 82.4 87.1 86.5 84.0 83.9 86.6 85.8,86.2 85.7 87.0 87.6 82.0 81.5 80.0 80.5,79.5 80.3,;,PROC ANOVA,;,CLASS A B;,MODEL y=A B(A),;,RUN;,DATA R;,INPUT T$B$X;,CARDS;,A1 B1 205 A1 B2 168 A1 B3 222 A1 B4 230,A2 B1 230 A2 B2 198 A2 B3 242 A2 B4 255,A3 B1 252 A3 B2 248 A3 B3 305 A3 B4 260,A4 B1 200 A4 B2 158 A4 B3 183 A4 B4 196,A5 B1 265 A5 B2 275 A5 B3 315 A5 B4 282,;,PROC GLM,;,CLASS,T B;,MODEL,X=T B;,MEANS,T/,DUNCAN,;,RUN;,DATA L;,INPUT T$ROW COL X;,CARDS;,D 1 1 23 E 1 2 21 A 1 3 24 B 1 4 21 C 1 5 19,A 2 1 22 C 2 2 20 E 2 3 20 D 2 4 21 B 2 5 22,E 3 1 20 A 3 2 25 B 3 3 26 C 3 4 22 D 3 5 23,B 4 1 25 D 4 2 22 C 4 3 25 E 4 4 21 A 4 5 23,C 5 1 19 B 5 2 20 D 5 3 24 A 5 4 22 E 5 5 19,;,PROC GLM,;,CLASS,T ROW COL;,MODEL,X=T ROW,COLl,;,MEANS,T/,DUNCAN,;,RUN;,（三）一元线性回归与相关分析,（教材,【,例,8.1】,）,DATA G；,INPUT,x,y,；,CARDS；,80 2350 86 2400 98 2720 90 2500 120 3150,102 2680 95 2630 83 2400 113 3080 105 2920,110 2960 100 2860,；,PROC,REG,CORR,；,MODEL,y,=,x,；,RUN；,（四）,协方差分析,（教材,【,例,10.1】,）,DATA K;,INPUT T$X Y;,CARDS;,CK 1.50 12.40 CK 1.85 12.00 CK 1.35 10.80,CK 1.45 10.00 CK 1.40 11.00 CK 1.45 11.80,CK 1.50 12.50 CK 1.55 13.40 CK 1.40 11.20,CK 1.50 11.60 CK 1.60 12.60 CK 1.70 12.50,1 1.35 10.20 1 1.20 9.40 1 1.45 12.20,1 1.20 10.30 1 1.40 11.30 1 1.30 11.40,1 1.