利用SAS软件对方差不同质试验数据的分析.docxVIP

青海大学学报 (自然科学版)Journal of Qinghai UniversityVol12 青海大学学报 (自然科学版) Journal of Qinghai University Vol122 No14 Aug12004 第 22 卷 第 4 期 2004 年 8 月 利用 SAS 软件对方差不同质试验数据的分析 乔有明 ( 青海大学农牧学院 ,青海 西宁 810003) 摘要 :文中利用了一组试验数据说明如何利用 SAS 软件解决试验数据分析中经常碰到的问 题 ,如正态性检验 、同质性检验 、方差异质数据的分析等 。并对应用不同统计模型的结果进行 了比较 。 关键词 :试验数据分析 ;方差异质 ; SAS 软件 中图分类号 :O21211 文献标识码 :B 文章编号 :1006 - 8996 (2004) 04 - 0073 - 03 标准的方差分析对原始数据有一些基本假定 ,其中之一就是假设方差是同质的 。如果数据不符合 这些假设 ,也未经过适当的数据转换 ,则方差分析所得到的结果就值得怀疑 ,甚至可以说是错误的 。所 以方差同质性的假设就成为这些假设中最为重要的一个 ,而这个假设也常常是最容易被忽略或最不容 易满足的 。Peddada 和 Smith 分析了方差异质情况下数据的线性模型和固定效应1 ,2 。Smith 和 Hickling 推导出了一个适于小样本且方差不同质时的单向分组资料方差分析公式3 。这些研究成果有的侧重于 数理统计方面 ,有的还未得到进一步的验证 ,所以在教科书或文章中被引用的情况极少 。事实上统计分 析软件 SAS 已提供了方差不同质的数据的分析方法 ,遗憾的是这一方法很少引起人们的注意 。本文通 过一个实例来说明 : ①利用 SAS 软件检验数据及误差的正态性 ; ②利用 SAS 软件对方差不等数据进行 分析 ; ③比较不同的统计模型对同一数据的处理结果 。 试验数据 本文所采用的数据 (见表 1) 引自 Steel 和 Torrie1980 年的著作4 。试验测定了红三叶草接种 6 种不 同的根瘤菌后植株体中的氮含量 ,以比较不同根瘤菌的接种效果 。这是一个单项分组资料 ,每个处理含 1 氮量重复测定 5 次 。每个处理代表一个总体 ,分别有其平均数μi 和方差σ2 。对这样一个资料作 i 方差分析 ,至少认为在理论上它满足两个基本 假定 : ①试验误差是独立的 ,并作正态分布 。即 可以假设研究者在进行此试验时是将 6 个处理 随机地分配给 30 株 ,且各株之间互不影响 ; ② 试验误差方差是同质的 ,因为测试的 30 株最初 来自同一个总体 ,所以处理之间的误差方差理 论上没 有 差 异 。但 实 际 存 在 的 一 种 可 能 情 况 是 ,由于处理的作用 ,各个处理的平均数和方差 都受到影响而产生差异 。 表 1 不同根瘤菌接种红三叶草后植株体中的氮含量 含氮量测定值 (mg) 菌种 1914 1717 1710 2017 1413 1713 3216 2418 1914 2110 1414 1914 2710 2719 911 2015 1118 1911 3211 2512 1518 1818 1116 1619 3310 2413 2017 1816 1412 2018 1 2 3 4 5 6 2 SAS 程序及其运行结果 2 . 1 方差分析程序及其结果 对表 1 的数据进行方差分析的 SAS (版本 :8e) 程序如下 ,结果见表 2 、表 3 。 options nodate ; data ;infile’clover . dat’;/ 3 表 1 数据文件 3 /input trtmt N ;proc anova ;/ 3 方差分析程序 3 data ; infile’clover . dat’;/ 3 表 1 数据文件 3 / input trtmt N ; proc anova ;/ 3 方差分析程序 3 / class trtmt ; model N = trtmt ; means trtmt/ 1sd ; run ; 表 2 表 1 资料的方差分析 变异来源 自由度 平方和 均方 F 值 概率 ( Pr F) 8471046667 2821928000 11291974667 1691409333 111788667 14137 010001 模 型 误 差 总变异 5 24 29 表 3 6 种接种处理植株含氮量的差异显著性 这一结果表示各处理之间的植株含氮量有极 显著的差异 ,需进一步比较各处理平均数的差异 。 这里采用 LSD 法 ,其中 LSD0105 = 414818 。 2 . 2 误差的正态性检