实验二方差分析.docx
《实验二方差分析.docx》由会员分享,可在线阅读,更多相关《实验二方差分析.docx(10页珍藏版)》请在冰豆网上搜索。
实验二方差分析
实验报告
课程名称:
_______统计软件___________
学院名称:
____数学与应用数学_______
班级:
_________________
姓名:
____________________
学号:
_______________
_2015-2016_学年第___一___学期
数学与统计学院制
实验地点
三教5楼实验室
课程类别
①公共课□②专业课□
实验日期
2015.12.26
实验编组
第3组
实验所
用时间
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
三、实验步骤及程序执行过程和结果
实验分析:
考核结果
教师签名:
年月日