生物统计实验saPPT演示文稿

《生物统计与试验设计》《生物统计与试验设计》实验课（实验课（SASSAS））山东农业大学信息科学与工程学院数学系山东农业大学信息科学与工程学院数学系陈陈 茂茂 学学1 SASSAS系统是用于数据分析与决策支持的大型集系统是用于数据分析与决策支持的大型集成式模块化软件系统成式模块化软件系统 其早期的名称为统计分析系统其早期的名称为统计分析系统-------------- Statistical Analysis System Statistical Analysis System（缩写为（缩写为SASSAS）） 在数据处理和统计分析领域，在数据处理和统计分析领域，SASSAS系统被誉为系统被誉为国际上的标准软件系统国际上的标准软件系统第一章第一章 SAS SAS系统简介系统简介21.1 SAS1.1 SAS软件简史软件简史 早在早在6060年代，美国北卡罗来纳州年代，美国北卡罗来纳州 Raleigh Raleigh的的SASSAS软件有限公司开始研制统计分析系统软件有限公司开始研制统计分析系统------ ------ SASSAS软件包。

软件包19761976年研制成功，并作为产品推向市年研制成功，并作为产品推向市场目前已被场目前已被120120多个国家和地区的多个国家和地区的2900029000多个机多个机构所采用，直接用户构所采用，直接用户300300多万人起初，多万人起初， SAS SAS是建是建立在立在DOSDOS环境下的有适于不同机型的各种版本环境下的有适于不同机型的各种版本在在WINDOWSWINDOWS出现不久，出现不久， SAS SAS公司就研制出公司就研制出WINDOWSWINDOWS环境下运行的环境下运行的SASSAS系统31.2 SAS1.2 SAS软件的结构软件的结构SASSAS软件是模块式结构，整个系统有软件是模块式结构，整个系统有3030多个模块组成多个模块组成SAS/BASE----SAS/BASE----基础核心模块基础核心模块SAS/STAT----SAS/STAT----常用和多元统计分析模块常用和多元统计分析模块SAS/IML-----SAS/IML-----矩阵运算模块矩阵运算模块SAS/ETS-----SAS/ETS-----时间序列预测模块时间序列预测模块SAS/QC------SAS/QC------工业生产质量控制和实验设计模块。

工业生产质量控制和实验设计模块SAS/ASSIST--SAS/ASSIST--处理一般简单问题模块处理一般简单问题模块SAS/INSIGHT-SAS/INSIGHT-交互式统计分析模块交互式统计分析模块SAS/OR------SAS/OR------运筹学模块运筹学模块SAS/AA------Analyst applicationSAS/AA------Analyst application分析员应用模块分析员应用模块41.3 SAS1.3 SAS软件的功能软件的功能 功能十分强大，包括客户机与服务器之间的功能十分强大，包括客户机与服务器之间的信息交换和计算、数据访问、数据存储和管理、信息交换和计算、数据访问、数据存储和管理、数据报告和分析、质量控制和项目管理、图形处数据报告和分析、质量控制和项目管理、图形处理和实验设计、应用开发理和实验设计、应用开发51.4 1.4 进入进入SASSAS系统的方法和标志系统的方法和标志 在在WIN95WIN95或或WIN98WIN98下（假设已经装好下（假设已经装好SASSAS），直），直接调用接调用SASSAS快捷键或选择快捷键或选择““开始开始”--“”--“程序程序”--”--“The SAS System”“The SAS System”即可进入即可进入SASSAS系统。

系统 进入进入SASSAS后，屏幕出现两个窗口，即后，屏幕出现两个窗口，即LOGLOG（记（记录）、录）、PROGRAM EDITERPROGRAM EDITER（程序编辑器，简写（程序编辑器，简写PGMPGM））窗口，通过菜单条窗口，通过菜单条“Window”“Window”或或“Globals”“Globals”选项选项可以进入可以进入 OUTPUT OUTPUT（输出）窗口输出）窗口6SAS6.12界面界面7SAS8.2界面界面8SAS9.0界面界面9菜单中的几个子菜单菜单中的几个子菜单1 1 文件子菜单文件子菜单102 2 编辑子菜单编辑子菜单113 3 查看子菜单查看子菜单124 4 工具子菜单工具子菜单135 5 运行子菜单运行子菜单145 5 解决方案子菜单解决方案子菜单156 6 帮助子菜单帮助子菜单162.1 SAS2.1 SAS过程与过程与SASSAS程序程序 SAS SAS过程过程(SAS PROCEDURE)(SAS PROCEDURE)是是SASSAS研制者为解决特定研制者为解决特定问题而编写的、经过编译的问题而编写的、经过编译的SASSAS程序，以程序，以“.DLL”“.DLL”为为其扩展名，存储在其扩展名，存储在SASSAS各子目录下，供用户调用。

各子目录下，供用户调用 SAS SAS程序程序(SAS PROGRAM)(SAS PROGRAM)是是SASSAS用户运用用户运用SASSAS语言编语言编写的一段程序其目的是为了将用户的数据告诉写的一段程序其目的是为了将用户的数据告诉SASSAS系统调用特定的系统调用特定的SASSAS过程完成某项任务过程完成某项任务第二章第二章 用用SASSAS需了解的几个基本概念需了解的几个基本概念172.2 SAS2.2 SAS文件与文件与SASSAS数据集数据集 SAS SAS文件是由文件是由SASSAS系统创建和管理的有特殊系统创建和管理的有特殊结构的文件如：在程序编辑窗口，输入结构的文件如：在程序编辑窗口，输入n n行数据存盘得行数据存盘得到的文件到的文件 SAS SAS数据集是数据集是SASSAS文件的一种，其扩展名为文件的一种，其扩展名为. .SD2SD2产生SASSAS数据集的主要方法有：数据集的主要方法有： 1 PGM 1 PGM窗口运行一段窗口运行一段SASSAS程序（程序（ SAS SAS数据步）来产生数据步）来产生 2 2 用用VIEWTABLEVIEWTABLE界面输入并存储的方式来产生。

界面输入并存储的方式来产生 3 3 其他方式（其他方式（SAS/INSIGHT, SAS/AA,SAS/INSIGHT, SAS/AA,转换等）转换等）182.3 SAS2.3 SAS程序结构程序结构 SAS SAS程序包括：程序包括： 1 SAS 1 SAS数据步（数据步（SAS DATA STEP)SAS DATA STEP) 作用是建立作用是建立SASSAS数据集 2 SAS 2 SAS过程步（过程步（SAS PROCEDURE STEP)SAS PROCEDURE STEP) 作用是激活作用是激活SASSAS过程 例如：调查并测得例如：调查并测得8 8名名14--1614--16岁中学生有关指标岁中学生有关指标的数据如下，写一段的数据如下，写一段SASSAS程序，分析男女生体重均程序，分析男女生体重均值之间的差异是否显著值之间的差异是否显著19DATA students;DATA students;INPUT name$ sex$ age weight height;INPUT name$ sex$ age weight height;CARDS;CARDS;Wang M 14 40 150Wang M 14 40 150Zhang M 16 42 170Zhang M 16 42 170Li F 15 44 149Li F 15 44 149Tang M 15 38 162Tang M 15 38 162Liu F 14 47 162Liu F 14 47 162Chen F 16 52 168Chen F 16 52 168Zhu M 14 45 158Zhu M 14 45 158Jia F 16 45 167Jia F 16 45 167; ;PROC TTEST;PROC TTEST;CLASS sex;CLASS sex;VAR weight;VAR weight;run;run;从从“DADA”“DADA”到独占一行的到独占一行的“;”“;”为为SASSAS数据步；数据步；从从“PROC”“PROC”到到“RUN“RUN；；””为为SASSAS过程步。

过程步students--students--产生产生SASSAS数据集的名字；数据集的名字；INPUT--INPUT--产生产生5 5个变量；（个变量；（$--$--字符型）字符型）CARDS--CARDS--表明各变量的具体取值；表明各变量的具体取值；TTEST--tTTEST--t检验；检验；CLASS--CLASS--以性别以性别sexsex变量为分组变量；变量为分组变量；VAR--VAR--指明分析的变量（如指明分析的变量（如weight)weight)；；RUN--RUN--表明程序的结束表明程序的结束20程序运行结果程序运行结果(6.12)(6.12)21程序运行结果程序运行结果(SAS 9.0)(SAS 9.0)：：222.4 SAS2.4 SAS程序的产生办法程序的产生办法 1 1 在在PGMPGM窗口直接编写窗口直接编写SASSAS程序 2 2 利用利用VIEWTABLEVIEWTABLE建立好建立好SASSAS数据集后，借助部分非数据集后，借助部分非编程编程SASSAS模块自动产生模块自动产生SASSAS程序程序( (如如SAS/AA)SAS/AA)。

3 3 直接调用由直接调用由SASSAS软件应用开发者编好的软件应用开发者编好的SASSAS程序2.5 SAS2.5 SAS程序的运行程序的运行 在在PGMPGM窗口内有了窗口内有了SASSAS程序后，只需按一下工具拦程序后，只需按一下工具拦中的小人像图标，中的小人像图标，SASSAS系统就开始执行操作，完成指系统就开始执行操作，完成指定的任务如果程序有错，错误信息显示在定的任务如果程序有错，错误信息显示在LOGLOG窗口窗口内如果程序正确，计算结果显示在内如果程序正确，计算结果显示在OUTPUTOUTPUT窗口内23第三章第三章 SAS SAS语言语言3.1 SAS程序设计概述程序设计概述 SAS程序分：数据（程序分：数据（DATA)步和过程（步和过程（PROC)步 DATA步的基本语句有步的基本语句有DATA,INPUT,CARDS,….3.2 SAS常量、变量、函数、表达式和数组常量、变量、函数、表达式和数组3.2.1 常量常量 数值、字符、日期、时间常量数值、字符、日期、时间常量 如：数值常量，如：数值常量，1024，，-10，，1.5E23 字符常量，字符常量，“JAC”，，“1232”3.2.2 变量变量 如如 x,xy,x1,abc;x$,name$,…243.2.3 函数函数 (1) 算术函数：绝对值算术函数：绝对值ABS( ),平方根平方根SQRT( ),... (2) 截取函数：取整截取函数：取整INT( ),舍入舍入ROUND( ),… (3) 数学函数：数学函数：EXP( ),LOG( ),LOG2( ),LOG10( )… (4) 三角函数：三角函数：SIN( ),COS( ),TAN( ),ARSIN( )… (5) 字符函数：字符串长度字符函数：字符串长度LENGTH( ),取子字符串取子字符串SUBSTR( ),如如SUBSTR(“ABCDEF”,2,3)=“BCD” (6) 日期与时间函数：日期与时间函数：DATE( ),TIME( ),YEAR( ),…25(7) 随机函数：随机函数： UNIFORM( )、、RANUNI( )均匀分布均匀分布 NORMAL( )、、 RANNOR( )正态分布正态分布 RANEXP( )指数分布指数分布 RANBIN( )二项分布，二项分布，RANPOI( ) 普松分布，普松分布，...(8) 概率分布函数：概率分布函数：P(ξ<=x) PROBBNML(p, n, x)，二项分布函数；，二项分布函数； POISSON（（lambda, x), 普松分布函数；普松分布函数； PROBNORM(x), 标准正态分布函数；标准正态分布函数；……26(9) 分位数函数：分位数函数： PROBIT(p), 标准正态分布的分位数，标准正态分布的分位数， 如如PROBIT(0.025)=-1.96;…..(10) 样本统计函数：样本统计函数： MEAN( ),均值；如均值；如MEAN(2,4,9)=5 MAX( ),最大值；最大值；MIN( ),最小值；最小值；SUM( ),求和；求和； VAR( ),方差；方差；STD( ),标准差；标准差；…… 273.2.4 SAS运算符运算符 (1) 算术运算符：算术运算符：+，，-，，*，，/，，** (2) 比较运算符：比较运算符：=或或EQ,^=或或NE(不等于），不等于），>或或GT, <或或LT, >=或或GE, <=或或LE (3) 逻辑运算符：逻辑运算符：&或或AND, |或或OR, ^或或NOT3.2.5 表达式表达式 由运算符将常量、变量、函数连接在一起的式子。

由运算符将常量、变量、函数连接在一起的式子3.2.6 数组数组 ARRAY 数组名（数组名（{n}）） [$] 变量表变量表 [初始值初始值]283.3 SAS的基本输入的基本输入 INPUT 语句语句 (1) 自由格式自由格式 INPUT 变量变量[$][&];$表示字符型变量；当字符常量（表示字符型变量；当字符常量（<=8个字符）中有空个字符）中有空格时变量后加格时变量后加&,若字符中含有一个空格，则下一常量若字符中含有一个空格，则下一常量与该字符之间至少空两个空格例如：与该字符之间至少空两个空格例如：DATA new;INPUT name$& age hgt;CARDS;Carl Hui 32 170.0C. Jiang 46 180.8;PROC PRINT;RUN;运行结果运行结果29(2) 列方式列方式 INPUT 变量变量[$]起始列号起始列号[-终止列号终止列号][.小数位数小数位数]; 字符长度可以为字符长度可以为1-200 例如：例如：DATA one;INPUT name$1-10 hgt 13-16 .1 age 11-12;CARDS;Carl Hui 301700C.J. Jiang401808;PROC PRINT;RUN;30 另外还有命令另外还有命令方式、指针方式等。

方式、指针方式等 仅对行固定表示符的使用进行说明仅对行固定表示符的使用进行说明 行表示符：行表示符：@@@@例：例：DATA new;INPUT x1-x3 @@;CARDS;12 34 16 25 27 2323 56 24 76 98 55;PROC PRINT;RUN;313.4 SAS3.4 SAS程序结构程序结构3.4.1 3.4.1 顺序结构顺序结构3.4.2 3.4.2 转移结构转移结构SELECT [SELECT [表达式表达式] ]；； WHEN < WHEN <条件表达式条件表达式1> 1> 语句语句1 1；； …… …… [OTHERWISE [OTHERWISE 语句语句] ]；；END;END;多分支多分支结构：结构：IF THEN > THEN 语句语句 [ [；；ELSE ELSE 语句语句] ]；；3.4.3 3.4.3 分支结构分支结构GOTO GOTO 标号；（标号标号；（标号 OK:y=x+20;) OK:y=x+20;)IF THEN GOTO > THEN GOTO 标号；标号；323.4.4 3.4.4 DODO循环结构循环结构DATA test;DATA test;ARRAY x(10) t1-t10;ARRAY x(10) t1-t10;INPUT t1-t10; SUM=0;INPUT t1-t10; SUM=0;DO I=1 to 10;DO I=1 to 10; sum=sum+x(I); sum=sum+x(I);END;END;CARDS;CARDS;1 2 3 4 5 6 7 8 9 101 2 3 4 5 6 7 8 9 10; ;PROC PRINT; RUN;PROC PRINT; RUN;DODO [ [循环变量循环变量= =初值初值[TO [TO 终值终值[BY [BY 步长步长]]]]；； SAS SAS语句；语句；END;END;例：例：333.5 3.5 数据预处理数据预处理3.5.1 3.5.1 产生新变量产生新变量 DATA one;DATA one;INPUT x y@@;INPUT x y@@;x2=x*x;x2=x*x;x3=x*x*x;x3=x*x*x;z=x/y*100;z=x/y*100;yp=LOG(y);yp=LOG(y);CARDS;CARDS;1 232 3 234 2 223 4 3541 232 3 234 2 223 4 354; ;PROC PRINT;RUN;PROC PRINT;RUN;343.5.2 3.5.2 修改变量修改变量 例如：例如：x=x-100;y=LOG(y);x=x-100;y=LOG(y);3.5.3 3.5.3 删除变量删除变量 例如：例如：DATA one;DATA one;INPUT x1 x2 x3 y;INPUT x1 x2 x3 y;DROP x2 x3;DROP x2 x3;CARDS;CARDS;1 2 3 401 2 3 402 4 5 802 4 5 803 5 7 833 5 7 83; ;PROC PRINT;PROC PRINT;RUN;RUN;354.1 4.1 完全随机化设计完全随机化设计例例1 1 夏玉米硼肥肥效试验。

处理夏玉米硼肥肥效试验处理A—A—施硼肥施硼肥; ;处理处理B—B—不施不施硼肥硼肥; ;完全随机化设计完全随机化设计, ,各重复各重复5 5次次, ,小区产量结果如下小区产量结果如下: : A—39, 37.6, 44.3, 42.4, 44.1 A—39, 37.6, 44.3, 42.4, 44.1 B—33.9, 35.9, 39.1, 33.5, 38.1 B—33.9, 35.9, 39.1, 33.5, 38.1试分析两处理平均产量间有无显著差异试分析两处理平均产量间有无显著差异第四章第四章 两均值检验及统计量计算两均值检验及统计量计算DATA tem; INPUT g$ n;DATA tem; INPUT g$ n;DO i=1 to n; INPUT x@@; OUTPUT; END;DO i=1 to n; INPUT x@@; OUTPUT; END;CARDS;CARDS;A 5A 539.0 37.6 44.3 42.4 44.139.0 37.6 44.3 42.4 44.1B 5B 533.9 35.9 39.1 33.5 38.133.9 35.9 39.1 33.5 38.1; ;PROC TTEST; CLASS g; VAR x; RUN;PROC TTEST; CLASS g; VAR x; RUN;3637例例2 2 为鉴定某小麦新品系为鉴定某小麦新品系A A的生产能力，以当地的优良的生产能力，以当地的优良小麦品种小麦品种B B为对照，两两组成一对分别种植于相邻的两为对照，两两组成一对分别种植于相邻的两个小区上，重复个小区上，重复8 8次，试验结果如下：次，试验结果如下：A—21 20 18 17 16 17 19 20A—21 20 18 17 16 17 19 20B—15 17 15 16 16 16 17 19B—15 17 15 16 16 16 17 19试分析两种子产量之间有无显著差异。

试分析两种子产量之间有无显著差异程序如下：程序如下：4.2 4.2 配对设计配对设计38DATA tem; DATA tem; INPUT x1 x2@@; INPUT x1 x2@@; d=x1-x2;d=x1-x2;CARDS;CARDS;21 15 20 17 18 15 17 16 21 15 20 17 18 15 17 16 16 16 17 16 19 17 20 1916 16 17 16 19 17 20 19; ;PROC MEANS MEAN STDERR T PRT MAXDEC=4;PROC MEANS MEAN STDERR T PRT MAXDEC=4;VAR d;VAR d;RUN;RUN;39例例3 某班某班31个学生，某门课程的考试成绩如下：个学生，某门课程的考试成绩如下： 25 45 50 54 55 61 64 68 72 75 75 78 79 81 83 84 84 84 85 86 86 86 87 89 89 89 90 91 91 92 1004.3 4.3 统计量的计算统计量的计算计算均值、方差、标准差、变异系数、中位数、极差。

计算均值、方差、标准差、变异系数、中位数、极差data ex3;input x @@;cards;25 45 50 54 55 61 64 68 72 75 75 78 79 81 83 8484 84 85 86 86 86 87 89 89 89 90 91 91 92 100;proc means mean var std cv median range maxdec=2;var x; run; quit;解解 SAS程序：程序：40data ex3;input x @@;cards;25 45 50 54 55 61 64 68 72 75 75 78 79 81 83 8484 84 85 86 86 86 87 89 89 89 90 91 91 92 100;proc means mean var std cv median range maxdec=2;var x; run; quit;输出结果输出结果:41 例例1 1设有设有5 5个品种的单因素完全随机化设计，各重复个品种的单因素完全随机化设计，各重复4 4次，小区产量结果如下，试进行方差分析。

次，小区产量结果如下，试进行方差分析品种品种 重重 复复 1 2 3 4A1A2A3A4A524 30 28 2627 24 21 2631 28 25 3032 33 33 2821 22 16 215.1 5.1 单因素完全随机化设计的方差分析单因素完全随机化设计的方差分析第五章第五章 方差分析方差分析42DATA one; INPUT a$ n;DATA one; INPUT a$ n;DO i=1 TO n; INPUT y@@; OUTPUT; END;DO i=1 TO n; INPUT y@@; OUTPUT; END;CARDS;CARDS;A1 4 24 30 28 26A2 4 27 24 21 26A3 4 31 28 25 30A4 4 32 33 33 28A5 4 21 22 16 21; ;PROC ANOVA; CLASS a; MODEL y=a; PROC ANOVA; CLASS a; MODEL y=a; MEANS a /LSD ALPHA=0.05; RUN; QUIT;MEANS a /LSD ALPHA=0.05; RUN; QUIT;43方差分析结果方差分析结果44多重比较结果多重比较结果45 例例2 2 有有6 6个大豆品种的比较试验，随机区组设计，个大豆品种的比较试验，随机区组设计，4 4个区组，试验结果如下表。

试进行方差分析个区组，试验结果如下表试进行方差分析品品种种 区区 组组 1 2 3 4A1A2A3A4A5A658 54 50 4942 38 41 3632 36 29 3546 45 43 4635 31 34 3444 42 36 385.2 5.2 单因素随机区组设计的方差分析单因素随机区组设计的方差分析46DATA two;DATA two;DO DO a=‘A1’,’A2’,’A3’,’A4’,’A5’,’A6a=‘A1’,’A2’,’A3’,’A4’,’A5’,’A6’;’;DO b=1 to 4;INPUT y@@; OUTPUT; END; END;DO b=1 to 4;INPUT y@@; OUTPUT; END; END;CARDS;CARDS;58 54 50 49 42 38 41 3632 36 29 35 46 45 43 4635 31 34 34 44 42 36 38; ;PROC ANOVA; CLASS a b;PROC ANOVA; CLASS a b;MODEL y=a b;MODEL y=a b;MEANS a /TUKEY;MEANS a /TUKEY;RUN;RUN;47方差分析结果方差分析结果48多重比较结果多重比较结果49例例3 3：玉米试验。

玉米试验A A施肥期：前期施肥期：前期A1 A1 ，后期，后期A2 A2 B B施肥量：施肥量：B1—N5B1—N5；；B2—N10B2—N10；；B3—N15B3—N15设置设置5 5个区组，试验结果如下，进行方差分析个区组，试验结果如下，进行方差分析 处理处理 区区 组组 1 2 3 4 5A1B1A1B2A1B3A2B1A2B2A2B350 59 45 52 5860 53 58 62 5564 66 67 71 6355 60 56 62 4571 65 75 78 8055 59 58 63 605.3 5.3 多因素随机区组设计的方差分析多因素随机区组设计的方差分析50DATA ex3;DATA ex3;DO a=‘A1’,’A2’; DO DO a=‘A1’,’A2’; DO b=‘B1’,’B2’,’B3’;b=‘B1’,’B2’,’B3’;DO c=1 to 5’;DO c=1 to 5’;INPUT y@@; OUTPUT; END; END; END;INPUT y@@; OUTPUT; END; END; END;CARDS;CARDS;50 59 45 52 58 60 53 58 62 5564 66 67 71 63 55 60 56 62 4571 65 75 78 80 55 59 58 63 60; ;PROC ANOVA;CLASS a b c;PROC ANOVA;CLASS a b c;MODEL y=a b c a*b;MODEL y=a b c a*b;MEANS a b a*b/LSD;MEANS a b a*b/LSD;RUN;RUN;51方差分析结果方差分析结果52A因素多重比较结果因素多重比较结果53B因素多重比较结果因素多重比较结果54 例例4 4 A A、、B B、、C C、、D D、、E E 五个水稻品种试验的拉丁五个水稻品种试验的拉丁方设计。

试验设计及试验结果如下表试进行方差方设计试验设计及试验结果如下表试进行方差分析 行行 列列 1 2 3 4 5 1 2 3 4 5D 37 A 38 C 38 B 44 E 38B 48 E 40 D 36 C 32 A 35C 27 B 32 A 32 E 30 D 26E 28 D 37 B 43 A 38 C 41A 34 C 30 E 27 D 30 B 415.4 5.4 拉丁方设计的方差分析拉丁方设计的方差分析55DATA ex4;DATA ex4;DO r=1 TO 5; DO c=1 TO 5;DO r=1 TO 5; DO c=1 TO 5;INPUT t$ y@@; OUTPUT; END; END;INPUT t$ y@@; OUTPUT; END; END;CARDS;CARDS;D 37 A 38 C 38 B 44 E 38B 48 E 40 D 36 C 32 A 35C 27 B 32 A 32 E 30 D 26E 28 D 37 B 43 A 38 C 41A 34 C 30 E 27 D 30 B 41; ;PROC ANOVA;CLASS r c t;PROC ANOVA;CLASS r c t;MODEL y=r c t;MODEL y=r c t;MEANS t/LSD; RUN;MEANS t/LSD; RUN;56方差分析结果方差分析结果57处理的多重比较结果处理的多重比较结果58第六章第六章 相关与回归分析相关与回归分析6.1 相关系数的计算相关系数的计算 例例1 1 某物质在凝固时放出的热量某物质在凝固时放出的热量y y（卡（卡/ /克）与此物克）与此物质中质中4 4种化学成分（种化学成分（% %））x1,x2,x3,x4x1,x2,x3,x4有关有关, ,求这求这5 5个变个变量间的相关系数。

数据与程序如下：量间的相关系数数据与程序如下：59DATA hald;INPUT x1 x2 x3 x4 y;CARDS; 7 26 6 60 78.5 1 29 15 52 74.311 56 8 20 104.311 31 8 47 87.6 7 52 6 33 95.911 55 9 22 109.2 3 71 17 6 102.7 1 31 22 44 72.5 2 54 18 22 93.121 47 4 26 115.9 1 40 23 34 83.811 66 9 12 113.310 68 8 12 109.4;PROC CORR;VAR x1-x4 y;RUN;60相关分析结果相关分析结果61例例2 某植物的苗龄某植物的苗龄x(d)与株高的试验数据如下，求与株高的试验数据如下，求x与与y的直线回归方程的直线回归方程6.2 简单线性回归简单线性回归苗龄苗龄x(d) 5 10 15 20 25 30 35株高株高y(cm) 2 5 9 14 19 25 33其其SASSAS程序如下：程序如下：DATA two; INPUT x yDATA two; INPUT x y@@; ;CARDS;CARDS;5 2 10 5 15 9 20 14 25 19 30 25 35 33;PROC REG; MODEL y=x; RUN;PROC REG; MODEL y=x; RUN;62方差分析与参数估计方差分析与参数估计63如果输出预测值，如果输出预测值，在在 MODEL MODEL 语句后加语句后加/p/pPROC REG;PROC REG; MODEL y=x/P; MODEL y=x/P;RUN;RUN;PROC GPLOT PROC GPLOT ；；PLOT PLOT 纵轴变量名纵轴变量名* *横轴变量名；横轴变量名；SYMBOL V=SYMBOL V=符号符号 CV= CV=颜色颜色 I=NONE I=NONE 其他图形选项；其他图形选项；RUN;RUN;用用GPLOTGPLOT过程作图过程作图其中：符号其中：符号——散点符号；散点符号；I=NONE—I=NONE—散点间不连线；散点间不连线； I=RL I=RL表示线性回归（表示线性回归（RQRQ、、RCRC是二、三次回归）是二、三次回归） 64PROC GPLOT;PROC GPLOT;PLOT y*x;PLOT y*x; SYMBOL V=star I=RL CV=orange CI=blue; SYMBOL V=star I=RL CV=orange CI=blue;RUN;RUN;其中：其中：CVCV、、CL—CL—分别表示点的符号和回归线的颜色分别表示点的符号和回归线的颜色上例作上例作y关于关于x的回归和散点图。

程序如下：的回归和散点图程序如下：65666.3 曲线回归曲线回归 例例3 苹果叶面积苹果叶面积y(cm2)生长量与前期间隔天数生长量与前期间隔天数x(d)的数的数据如下表，求据如下表，求x与与y的指数曲线回归方程的指数曲线回归方程x 0 5 10 15 20 25 30 0 5 10 15 20 25 30 y 5.7 43.7 76.7 102.3 183.4 225.1 344.25.7 43.7 76.7 102.3 183.4 225.1 344.2其其SASSAS程序如下：程序如下：DATA three;INPUT x y@@; yp=LOG(y);DATA three;INPUT x y@@; yp=LOG(y);CARDS;CARDS;0 5.7 5 43.7 10 76.7 15 102.3 20 183.4 25 225.1 30 344.2;PROC REG; MODEL yp=x; RUN;PROC REG; MODEL yp=x; RUN;67输出结果：输出结果：指数方程：指数方程：y=e2.6594*e0.1175 x=14.2877e0.1175 x68在例在例1中，建立中，建立y与与x1,x2,x3,x4的多元线性回归方程。

的多元线性回归方程数据步）数据步）PROC REG; MODEL y=x1 x2 x3 x4;RUN;6.4 多元线性回归多元线性回归69输出结果输出结果70 采用逐步回归方法筛选自变量采用逐步回归方法筛选自变量, ,建立建立y y与与x1,x2,x3,x4x1,x2,x3,x4的多元线性回归方程程序如下：的多元线性回归方程程序如下：PROC REG;PROC REG;MODEL y=x1 x2 x3 x4/MODEL y=x1 x2 x3 x4/SELECTION=STEPWISESELECTION=STEPWISE; ;RUN;RUN;注意：在注意：在REGREG过程中，可以有多个过程中，可以有多个MODELMODEL语句71 引入、剔除变量显著性水平的选择引入、剔除变量显著性水平的选择 在在 MODEL MODEL 语句中的选项语句中的选项 “SELECTION= “SELECTION=方法方法” ” 后后面加入选项面加入选项 1 1 引入变量显著性水平引入变量显著性水平 STEPWISE STEPWISE法用：法用：SLE=SLE=水平值，缺省值为水平值，缺省值为0.15.0.15. 2 2 剔除变量显著性水平剔除变量显著性水平 STEPWISE STEPWISE法用：法用：SLS=SLS=水平值，缺省值为水平值，缺省值为0.15.0.15.72DATA hald;INPUT x1 x2 x3 x4 y;CARDS; 7 26 6 60 78.5 1 29 15 52 74.311 56 8 20 104.311 31 8 47 87.6 7 52 6 33 95.911 55 9 22 109.2 3 71 17 6 102.7 1 31 22 44 72.5 2 54 18 22 93.121 47 4 26 115.9 1 40 23 34 83.811 66 9 12 113.310 68 8 12 109.4;PROC REG;MODEL y=x1 x2 x3 x4/SELECTION=STEPWISE SLE=0.1 SLS=0.1;RUN;73逐步回归分析结果逐步回归分析结果747576 7.1.1 离散型分布的检验离散型分布的检验 例例1 研究大豆蛋白质型类的遗传，蛋白质型类研究大豆蛋白质型类的遗传，蛋白质型类A与蛋白与蛋白质型类质型类B杂交第一代（杂交第一代（F1）为）为AB，当，当F1自交得自交得F2共分离共分离成成A型类为型类为24，，AB型类为型类为35，，B型类为型类为21，问这种现象，问这种现象是否符合孟德尔是否符合孟德尔1：：2：：1的遗传分离定律。

的遗传分离定律7.1 适合性检验适合性检验第七章第七章 非参数统计非参数统计解：解：SAS程序如下：程序如下：77data ex1; input f$ y @@;cards;A 24 AB 35 B 21;proc freq; weight y;tables f/nocum testp=(25 50 25); run; quit;输出结果输出结果:787.1.2 连续型分布的检验连续型分布的检验检验是否服从正态分布检验是否服从正态分布例例2 对对100名女生测定血清总蛋白含量名女生测定血清总蛋白含量(g/L), 数据如下数据如下:74.3 78.8 68.8 78.0 70.4 80.5 80.5 69.7 71.2 73.579.5 75.6 75.0 78.8 72.0 72.0 72.0 74.3 71.2 72.075.0 73.5 78.8 74.3 75.8 65.0 74.3 71.2 69.7 68.073.5 75.0 72.0 64.3 75.8 80.3 69.7 74.3 73.5 73.575.8 75.8 68.8 76.5 70.4 71.2 81.2 75.0 70.4 68.070.4 72.0 76.5 74.3 76.5 77.6 67.3 72.0 75.0 74.373.5 79.5 73.5 74.7 65.0 76.5 81.6 75.4 72.7 72.767.2 76.5 72.7 70.4 77.2 68.8 67.3 67.3 67.3 72.775.8 73.5 75.0 73.5 73.5 73.5 72.7 81.6 70.3 74.373.5 79.5 70.4 76.5 72.7 77.2 84.3 75.0 76.5 70.4解：解：SAS程序如下：程序如下：79解：解：SAS程序如下：程序如下：data ex1; input x@@;cards;74.3 78.8 68.8 78.0 70.4 80.5 80.5 69.7 71.2 73.579.5 75.6 75.0 78.8 72.0 72.0 72.0 74.3 71.2 72.075.0 73.5 78.8 74.3 75.8 65.0 74.3 71.2 69.7 68.073.5 75.0 72.0 64.3 75.8 80.3 69.7 74.3 73.5 73.575.8 75.8 68.8 76.5 70.4 71.2 81.2 75.0 70.4 68.070.4 72.0 76.5 74.3 76.5 77.6 67.3 72.0 75.0 74.373.5 79.5 73.5 74.7 65.0 76.5 81.6 75.4 72.7 72.767.2 76.5 72.7 70.4 77.2 68.8 67.3 67.3 67.3 72.775.8 73.5 75.0 73.5 73.5 73.5 72.7 81.6 70.3 74.373.5 79.5 70.4 76.5 72.7 77.2 84.3 75.0 76.5 70.4;proc capability graphics ;histogram x / normal; run;8081卡平方检验：卡平方检验： 2=4.06491174；自由度；自由度=5；；p=0.540。

827.2 独立性检验独立性检验(列联表分析列联表分析)解：解：SAS程序如下程序如下例例2 考察不同灌溉方式对水稻叶子衰老情况的影响考察不同灌溉方式对水稻叶子衰老情况的影响, 观观测结果如下表，试检验叶子衰老与灌溉方式是否有关测结果如下表，试检验叶子衰老与灌溉方式是否有关灌溉方式灌溉方式绿叶绿叶1黄叶黄叶2枯叶枯叶3Ti.深水深水A114677160浅水浅水A2183913205湿润湿润A31521416182T.j481303654783data ex2;do a=1 to 3; do b=1 to 3;input f@@; output; end; end;cards;146 7 7183 9 13152 14 16;proc freq;tables a*b/chisq expected nocol norow nopercent;weight f;run;84理论频数理论频数理论频数理论频数 T Ti i .T..T.j j / n/ n85。