实验二方差分析.docx
- 文档编号:11543294
- 上传时间:2023-03-19
- 格式:DOCX
- 页数:12
- 大小:187.13KB
实验二方差分析.docx
《实验二方差分析.docx》由会员分享,可在线阅读,更多相关《实验二方差分析.docx(12页珍藏版)》请在冰豆网上搜索。
实验二方差分析
实验报告
课程名称:
______统计软件__________________
学院名称:
________数学与统计学院___________
班级:
_______122_______________________
姓名:
______甘天明_____________________
学号:
______1250901205___________
_2015-2016学年第_1学期
数学与统计学院制
实验地点
三教五楼实验室
课程类别
1公共课□
②专业课□
实验日期
2015-12-18
实验编组
第5组
实验所
用时间
2小时
实验名称
方差分析
实验目的
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 (2)用过程ANOVA,检验是否存在显著性差异 从图4.2知,p<.0001, 可判断存在显著性差异 (3)用MEANSvariety/ duncan,可具体的看哪些水平存在差异. 从图4.3知,(a2,a5),a6,a1,(a4,a7),a3它们之间存在差异。 (4)用MEANSvariety/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) 结果分析: (1)由anlysisofvariance中的输出结果图4.6,可得模型及d24,kt,d24*kt都是显著的。 (2)由于交互作用是显著的,所以选择最优的组合,结合图4.7可知,d24=0.1,kt=0.05的组合。 注: 若不存在交互作用,则分别选择d24和kt的最优者即可。 (二)实验题目 1.2014年下半年同一苜蓿品种接种五种生防菌进行比较实验(数据如下表)。 变量说明: deal为处理,期中CK为对照组,y代表根长。 解决的问题: 同一苜蓿品种接种不同生防菌后根长的差异。 deal y deal y deal y deal y deal y deal y XJ83073 45 NM024-2 38 XJ83060 40 NMG048-2 39 XJ83097 37 CK 39 XJ83073 41 NM024-2 33 XJ83060 38 NMG048-2 38 XJ83097 36 CK 50 XJ83073 44 NM024-2 25 XJ83060 37 NMG048-2 42 XJ83097 36 CK 40 XJ83073 36 NM024-2 35 XJ83060 32 NMG048-2 35 XJ83097 37 CK 46 XJ83073 38 NM024-2 40 XJ83060 42 NMG048-2 47 XJ83097 40 CK 47 XJ83073 39 NM024-2 41 XJ83060 40 NMG048-2 38 XJ83097 35 XJ83073 41 NM024-2 41 XJ83060 42 NMG048-2 39 XJ83097 32 XJ83073 36 NM024-2 35 XJ83060 36 NMG048-2 36 XJ83097 37 XJ83073 38 NM024-2 46 XJ83060 39 NMG048-2 31 XJ83097 38 XJ83073 45 NM024-2 38 XJ83060 40 NMG048-2 40 XJ83097 36 XJ83073 37 NM024-2 45 XJ83060 42 NMG048-2 38 XJ83097 33 XJ83073 40 NM024-2 40 XJ83060 39 NMG048-2 37 XJ83097 42 2.为了研究某生长素不同浓度对豌豆生长的影响,从而选择推荐使用浓度而进行了此试验(数据见下表)。 m1、m2、m3、m4、m5分别代表不同的浓度,即5个处理(treat)。 Weight为豌豆植株干重,单位为g/盆。 每个处理有4个重复,共20个观测值。 选择一种对豌豆生长有利的生长素浓度。 处理(treat) 豌豆植株干(单位: g/盆) 重复1重复2重复3重复4 m1 m2 m3 m4 m5 10 9 10 8 12 12 13 14 14 15 13 16 6 7 6 8 5 4 3 3 三、实验步骤及程序执行过程和结果 1.程序: datagc; inputdeal$y@@; cards; XJ8307345NM024-238XJ8306040NMG048-239XJ8309737CK39 XJ8307341NM024-233XJ8306038NMG048-238XJ8309736CK50 XJ8307344NM024-225XJ8306037NMG048-242XJ8309736CK40 XJ8307336NM024-235XJ8306032NMG048-235XJ8309737CK46 XJ8307338NM024-240XJ8306042NMG048-247XJ8309740CK47 XJ8307339NM024-241XJ8306040NMG048-238XJ8309735 XJ8307341NM024-241XJ8306042NMG048-239XJ8309732 XJ8307336NM024-235XJ8306036NMG048-236XJ8309737 XJ8307338NM024-246XJ8306039NMG048-231XJ8309738 XJ8307345NM024-238XJ8306040NMG048-240XJ8309736 XJ830007337NM024-245XJ8306042NMG048-238XJ8309733 XJ8307340NM024-240XJ8306039NMG048-237XJ8309742 ; run; procanovadata=gc; classdeal; modely=deal; meansdeal/duncanalpha=0.05; meansdeal/dunnett('deal'); run; quit; 2.程序: PROCNPAR1WAYdata=li_2WILCOXON; classdeal; vary; run; quit; datali_2; dotreat=1to5; doj=1to4; inputy@@; output; end; end; cards; 109108 12121314 14151316 6768 5433 run; procanova; classtreat; modely=treat; meanstreat; run; quit; 实验分析: 考核结果 教师签名: 年月日
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 方差分析