实验四方差分.docVIP

实 验 报 告 课程名称：_______________________________ 学院名称：_______________________________ 班 级：_______________________________ 姓 名：_______________________________ 学 号：_______________________________ _________________学年第 ____________学期 数 学 与 统 计 学 院 制 实验地点 课程类别 ①公共课□ ②专业课□ 实验日期 实验编组 第 组 实验所 用时间 小时 实验名称 实验四 方差分析 实验目的 1.使学生了解SAS系统 2.用SAS软件进行方差分析，并掌握其分析步骤。 3.掌握Anova过程过程和GLM过程的应用，以及结果的分析。 实验环境 SAS9.1 一、实验准备 1、SAS过程的选择： 均衡设计、拉丁方设计、正交设计等的一元、多元方差分析和重复测量的方差分析用 Anova过程；不能用Anova过程分析的数据就用GLM过程；另外，若需要对各因素间的各水平进行两两比较，则也用GLM过程。 2、结果分析的方法： 1) 单因素的方差分析： 先看该因素是否是显著的，若是则看看哪些水平是有显著性差异的，有必要的话还应该找出一个最有利的水平。 2) 双因素或多因素的方差分析： 先看交互作用是否是显著的，若是则找出一个最有利的组合； 若交互作用不显著，但存在主效应是显著的，则可以对显著的主效应分别按单因素的方法去分析。 二、实验举例及内容 （一）实验举例 例1、2002年12月进行的青贮实验，数据为7个不同品种青贮之间的可溶性有机物含量。玉米的7个不同品种variety分别为：a1高油玉米647a、a2高油玉米115a、a3高油玉米601、a4高油玉米289、a5农大80、a6玉米3138、a7特种玉米。每个品种有两个重复，实验数据一共有14个观测值。分析品种玉米青贮之间的可溶性有机物WSC的含量是否存在显著差异。 variety wsc variety wsc variety wsc variety wsc a1 20.19 a3 4.08 a5 39.38 a7 16.04 a1 21.03 a3 4.40 a5 38.95 a7 15.86 a2 36.06 a4 17.23 a6 30.04 a2 38.97 a4 16.85 a6 28.68 解： 程序如下（其中数据集为L3.anova5_1）： 主要输出结果如下： (图4.1) （图4.2） （图4.3） （图4.4） （图4.5） 结果分析： (1) 用过程univariate判断正态性 从图4.1知， Pr W 的p值为 0.1694，可判断wsc服从正态分布． (2)用过程ANOVA ，检验是否存在显著性差异 从图4.2知 ，p .0001, 　可判断存在显著性差异 (3) 用MEANS variety /　duncan，可具体的看哪些水平存在差异． 从图4.3知，(a2, a5), a6, a1, (a4,a7), a3 它们之间存在差异。 (4) 用MEANS variety /Dunnett (a1)，可分析它们与对照组‘a1’ 从图4.4知，由于检验时都出现了‘***’,说明它们与对照组‘a1’ (5) 用过程NPAR1WAY进行非参数检验．看是否存在显著性差异，以及哪个水平是最好的． 从图4.5可看出，在0.05的水平下，各水平之间存在差异，且（a5, a2）的水平是最高的。 例2 在组织培养过程中采用两种激素（KT、2,4-D）的不同浓度配比，通过观测材料的生根比例，选则合适某实验品种的配比。设计浓度及观测到的生根比例如下表： 2,4浓度（d24） 0 0.05 0.10 KT浓度（kt） 0 0.05 0.10 0 0.05 0.10 0 0.05 0.10 重复1 0.24 0.35 0.38 0.20 0.25 0.36 0.28 0.50 0.30 重复2 0.26 0.31 0.35 0.21 0.21 0.35 0.20 0.51 0.31 重复3 0.25 0.32 0.35 0.22 0.26 0.34 0.21 0.52 0.32 解： 程序如下： 主要输出结果如下： （图4.6） （图4.7）