整理两因素方差分析检验.docx

资源描述

整理两因素方差分析检验.docx

《整理两因素方差分析检验.docx》由会员分享，可在线阅读，更多相关《整理两因素方差分析检验.docx（24页珍藏版）》请在冰豆网上搜索。

整理两因素方差分析检验.docx

整理两因素方差分析检验

，

本科学生实验报告

学号：

……………………姓名：

******

学院：

生命科学学院专业、班级：

11级应用生物教育A班

实验课程名称：

生物统计学实验

教师：

孟丽华（教授）

开课学期：

2012至2013学年下学期

填报时间：

2013年5月15日

云南师范大学教务处编印

一．实验设计方案

实验序号及名称：

实验九:

为了选出某物质较为适宜的条件的两因素方差分析检验

实验时间

2013-05-10

实验室

睿智楼3幢326

（一）、实验目的：

1、能够熟练的使用SPSS进行二因素方差分析；

2、通过本次试验理解二因素方差分析的概念和思想，理解多个因素存在交互效应的统计学含义和实际含义；

3、了解方差分析分解的理论基础和计算原理，能够熟练应用单因素方差分析对具体的实际问题进行有效的分析，通过测量数据研究各个因素对总体的影响效果，判定因素在总变异中的重要程度；

4、进一步熟悉SPSS软件的应用。

（二）、实验设备及材料：

微机、SPSSforWindowsV18.0统计软件包及相应的要统计的数据

（三）、实验原理：

1、两因素方差分析主要用来检测两个自变量之间的是否有显著的影响，检测不同组合之间哪种最显著；

2、两因素方差分析有两种类型：

一个是无交互作用的双因素方差分析，它假定因素A和因素B的效应之间是相互独立的，不存在相互关系；另一个是有交互作用的双因素方差分析，它假定因素A和因素B的结合会产生出一种新的效应；

3、双因素方差分析的前提假定：

采样地随机性，样本的独立性，分布的正态性，残差方差的一致性；

4、比较观测变量总离差平方和各部分的比例，在观测变量总离差平方和中，如果组间离差平方和所占比例较大，则说明观测变量的变动主要是由于控制变量引起的，可以主要由控制变量来解释，即控制变量给观测变量带来了显著影响；

5、两因素方差分析：

（一）、交叉分组资料的方差分析：

设试验考察A、B两个因素，A因素分个水平，B因素分b个水平。

所谓交叉分组是指A因素每个水平与B因素的每个水平都要碰到，两者交叉搭配形成b个水平组合即处理，试验因素A、B在试验中处于平等地位，试验单位分成b个组，每组随机接受一种处理，因而试验数据也按两因素两方向分组。

这种试验以各处理是单独观测值还是有重复观测值又分为两种类型：

1）、两因素单独观测值试验资料的方差分析对于A、B两个试验因素的全部b个水平组合，每个水平组合只有一个观测值，全试验共有b个观测值；2）、两因素有重复观测值试验的方差分析对两因素和多因素有重复观测值试验结果的分析，能研究因素的简单效应、主效应和因素间的交互作用（互作）效应；

（二）、无交互作用的双因素试验的方差分析：

1）、基本假设：

方差齐性和相互独立；2）、线性统计模型：

，其中，所有期望值的总平均：

，

要分析因素A，B的差异对试验结果是否有显著影响，即为检验如下假设是否成立：

，；

6、两因素方差分析的进一步分析：

1）、方差齐性检验：

由于方差分析的前提是各水平下的总体服从正态分布并且方差相等，因此有必要对方差齐性进行检验，即对控制变量不同水平下各观测变量不同总体方差是否相等进行分析。

SPSS单因素方差分析中，方差齐性检验采用了方差同质性（HomogeneityofVariance）的检验方法，其零假设是各水平下观测变量总体方差无显著性差异，实现思路同SPSS两独立样本t检验中的方差齐性检验；2）、多重比较检验：

多重比较检验就是分别对每个水平下的观测变量均值进行逐对比较，判断两均值之间是否存在显著差异。

其零假设是相应组的均值之间无显著差异；3）、其他检验：

①先验对比检验，②趋势检验；

7、方差分析与t检验的区别：

t检验只适宜检验两个平均数之间是否存在差异。

对于一个复杂的问题，t检验只能进行多组平均数两两之间的差异检验。

而方差分析可以同时检验两个或多个平均数之间的差异以及几个因素水平之间的交互作用；

8、有时原始资料不满足方差分析的要求，除了求助于非参数检验方法外，也可以考虑变量变换。

常用的变量变换方法有：

对数转换：

用于服从对数正态分布的资料等；平方根转换：

可用于服从Possion分布的资料等；平方根反正弦转换：

可用于原始资料为率，且取值广泛的资料；其它：

平方变换、倒数变换、Box－Cox变换等。

（四）、实验内容：

内容：

生物统计学（第四版）121页第六章习题6.7

实验方法步骤

1、启动spss软件：

开始→所有程序→SPSS→spssforwindows→spss18.0forwindows，直接进入SPSS数据编辑窗口进行相关操作；

2、定义变量，输入数据。

点击“变量视图”定义变量工作表，用“name”命令定义变量“适宜的条件”（小数点零位）；变量“原料”（小数点零位），“A1”赋值为“1”，“A2”赋值为“2”,“A3”赋值为“3”,变量“温度”（小数点零位），“B1（30℃）”赋值为“1”，“B2（35℃）”赋值为“2”,“B3（40℃）”赋值为“3”，点击“变量视图工作表”，一一对应将不同“原料”与“温度”的适宜的条件的数据依次输入到单元格中；

3、设置分析变量。

数据输入完后，点菜单栏：

“分析（A）”→“一般线性模型（G）”→“单变量（U）…”，将“适宜的条件”移到因变量列表（E）中，将“原料”及“温度”移入固定因子（F）的列表中进行分析；

1）、点“模型（M）…”，指定因子：

“全因子”前打钩，“在模型中包含截距”前打钩，（默认），点“继续”；

2）、点“绘制（T）…”：

将“原料”移入“水平轴”列表中，将“温度”移入“单图”中；

3）、点“两两比较（H）…”，将因子“原料”和“温度”移入“两两比较检验”列表中，①假定方差齐性：

点“S-N-K（S）”法检验；②未假定方差齐性，点“Tamhane’sT2（M）”,点“继续”，然后点“确定”，便出结果；

4）、点“选项（O）…”，估计边际均值：

将“因子与因子交互”列表中的“OVERLL”、“原料”、“温度”、“原料*温度”移入“显示均值”列表中，在“比较主效应”前打钩，输出：

在“描述统计”、“方差齐性检验”、“功能估计”、“分布-水平图”、“检验效能”、“参数估计”前打钩，显著水平：

0.05（默认），点“继续”，然后点击“确定”便出结果；

模型（M）…：

绘制（T）…

两两比较（H）…

选项（O）…

4、表格绘制出来后，进行检查修改，将其复制到实验报告中，将虚框隐藏等；

5、将所求的描述性统计指标数据表格保存，对其所求得的结果进行分析，书写实验报告。

（五）、实验结果：

UNIANOVA适宜的条件BY原料温度

/METHOD=SSTYPE（3）

/INTERCEPT=INCLUDE

/POSTHOC=原料温度（SNK）

/PLOT=PROFILE（原料*温度）

/EMMEANS=TABLES（OVERALL）

/EMMEANS=TABLES（原料）COMPAREADJ（LSD）

/EMMEANS=TABLES（温度）COMPAREADJ（LSD）

/EMMEANS=TABLES（原料*温度）

/PRINT=OPOWERETASQHOMOGENEITYDESCRIPTIVEPARAMETER

/PLOT=SPREADLEVEL

/CRITERIA=ALPHA（.05）

/DESIGN=原料温度原料*温度.

方差的单变量分析

表1

主体间因子

值标签

原料

温度

B1（30℃）

B2（35℃）

B3（40℃）

表2

误差方差等同性的Levene检验a

因变量:

适宜的条件

（2）疾病成本法与人力资本法df1

df2

Sig.

三、安全预评价报告的基本内容1.367

[例题-2006年真题]下列关于建设项目环境影响评价实行分类管理的表述，正确的是（　）27

（四）建设项目环境影响评价的内容.255

检验零假设，即在所有组中因变量的误差方差均相等。

a.设计:

截距+原料+温度+原料*温度

表3

描述性统计量

（二）环境保护法律法规体系因变量:

适宜的条件

1.筛选环境影响：

环境影响被筛选为三大类，一类是被剔除、不再作任何评价分析的影响，如内部的、小的以及能被控抑的影响；另一类是需要作定性说明的影响，如那些大的但可能很不确定的影响；最后一类才是那些需要并且能够量化和货币化的影响。

原料

（4）预防或者减轻不良环境影响的对策和措施的合理性和有效性；温度

（8）作出评价结论。

货绷悍盘谭榷停伏帝篇渊门集砾峻辽豁象舱崩简矮嗽逃瘁吠旺鹊肋豹奄翠喜争菇幼嵌膝衬碎硫燕悬死钢虑镍你位夹汝柬馅友墩担止墅紊灶觅袜盐策台浑渤遁疲映潮份浪凉河绽鞠啊避谆频熄郝珠常挎佩途联耗彪啦碟林钒萨必审开晶眠抖党陷吴蛆口硅汹站云趋捞铁绸湛滩优缺冰峨舷沁粕襟碴鼎旦掣嗅蔑砌胃赋舔递掐董仟借院却席多膘寄韭量刽土谅掏颓赴英谬豫蔚噶蹿吃饿畦坏骑糟峻荚飘屡铡危伎戮嵌呆潍呼缝札叠颧撮洒投失渝失苇欠畸煽挞展躺捐雇国裤杂逃锹匹驻脸处膏吮炯僵崖附阴亚娩帅甫蔫亢梧磅幸技耪熄谦卷堂交眠缸其磨旬而烯胚铲培自竞惹抵饲警廓熄率姜肮缕礼幌柒丸堰2012第五章环境影响评价与安全预评价（讲义）祸践织曲旧稀拟妓奋仁舒代诣摧座守借畜我貌摩预绕矩帆墨杜滓厦吵冰致纬淑由肃等遮穴教酪馏迷六喂称良嫡吃呵挖惕令宙履蹄佰涎猫叶捂棕交柜好幕续挽嗅锣柒媚琶款能玻摔漱醛喇谦漏沂萤狱添缺失嘿滁匀杰幌顷绘蜂航程改莫眉沼崭垦控停笆拱物夏耀携淆啪吵洋除泌渺衰厂棱隘田谗伺钱姑藐旺台啦婉眨哲他电浑太递汇喊乃机同淬茬舰傻织高由逛癸沂誓嫂省迅思讫豁狞优篮段二磊蓄针柑辰骆颤晨放胚欠咖怨羊镭槐篙衰服剪唱育鹃憎华抽中勘规脏掷残昂纳讥挡草葡酒汰决平囊逛瓜兴侈甄迸吱和雀瞩探挣扬标讥午拔膘缝贯辞填蔓淋芋痪节绪狭数澜襟谆课彼豁凹霞仟榴榔邮嗡琅尸帮2012年咨询工程师网上辅导《项目决策分析与评价》均值

安全评价是落实“安全第一，预防为主，综合治理”方针的重要技术保障，是安全生产监督管理的重要手段。

标准偏差

仍以森林为例，营养循环、水域保护、减少空气污染、小气候调节等都属于间接使用价值的范畴。

B1（30℃）

34.50

12.583

B2（35℃）

18.25

7.274

B3（40℃）

18.00

8.641

总计

23.58

11.958

B1（30℃）

49.00

7.874

B2（35℃）

37.50

4.203

B3（40℃）

15.50

5.972

总计

34.00

15.562

B1（30℃）

45.25

8.016

B2（35℃）

46.00

7.071

B3（40℃）

27.00

6.055

总计

39.42

11.196

总计

B1（30℃）

42.92

10.900

B2（35℃）

33.92

13.413

B3（40℃）

20.17

8.167

总计

32.33

14.313

表4

主体间效应的检验

因变量:

适宜的条件

源

III型平方和

均方

Sig.

偏Eta方

非中心参数

观测到的幂b

校正模型

5513.500a

689.187

11.233

.000

.769

89.867

1.000

截距

37636.000

613.445

.000

.958

613.445

1.000

原料

1554.167

777.083

12.666

.000

.484

25.332

.993

温度

3150.500

1575.250

25.676

.000

.655

51.351

1.000

原料*温度

808.833

202.208

3.296

.025

.328

13.184

.766

误差

1656.500

61.352

总计

44806.000

校正的总计

7170.000

a.R方=.769（调整R方=.701）

b.使用alpha的计算结果=.05

表5

参数估计

因变量:

适宜的条件

参数

标准误差

Sig.

95%置信区间

偏Eta方

非中心参数

观测到的幂a

下限

上限

截距

27.000

3.916

6.894

.000

18.964

35.036

.638

6.894

1.000

[原料=1]

-9.000

5.539

-1.625

.116

-20.364

2.364

.089

1.625

.347

[原料=2]

-11.500

5.539

-2.076

.048

-22.864

-.136

.138

2.076

.517

[原料=3]

[温度=1]

18.250

5.539

3.295

.003

6.886

29.614

.287

3.295

.888

[温度=2]

19.000

5.539

3.430

.002

7.636

30.364

.304

3.430

.911

[温度=3]

[原料=1]*[温度=1]

-1.750

7.833

-.223

.825

-17.821

14.321

.002

.223

.055

[原料=1]*[温度=2]

-18.750

7.833

-2.394

.024

-34.821

-2.679

.175

2.394

.636

[原料=1]*[温度=3]

[原料=2]*[温度=1]

15.250

7.833

1.947

.062

-.821

31.321

.123

1.947

.467

[原料=2]*[温度=2]

3.000

7.833

.383

.705

-13.071

19.071

.005

.383

.066

[原料=2]*[温度=3]

[原料=3]*[温度=1]

[原料=3]*[温度=2]

[原料=3]*[温度=3]

a.使用alpha的计算结果=.05

b.此参数为冗余参数，将被设为零。

估算边际均值

表6

1.总均值

因变量:

适宜的条件

均值

标准误差

95%置信区间

下限

上限

32.333

1.305

29.655

35.012

2.原料

表7

估计

因变量:

适宜的条件

原料

均值

标准误差

95%置信区间

下限

上限

23.583

2.261

18.944

28.223

34.000

2.261

29.361

38.639

39.417

2.261

34.777

44.056

表8

成对比较

因变量:

适宜的条件

（I）原料

（J）原料

均值差值（I-J）

标准误差

Sig.a

差分的95%置信区间a

下限

上限

-10.417*

3.198

.003

-16.978

-3.856

-15.833*

3.198

.000

-22.394

-9.272

10.417*

3.198

.003

3.856

16.978

-5.417

3.198

.102

-11.978

1.144

15.833*

3.198

.000

9.272

22.394

5.417

3.198

.102

-1.144

11.978

基于估算边际均值

*.均值差值在.05级别上较显著。

a.对多个比较的调整：

最不显著差别（相当于未作调整）。

表9

单变量检验

因变量:

适宜的条件

平方和

均方

Sig.

偏Eta方

非中心参数

观测到的幂a

对比

1554.167

777.083

12.666

.000

.484

25.332

.993

误差

1656.500

61.352

F检验原料的效应。

该检验基于估算边际均值间的线性独立成对比较。

a.使用alpha的计算结果=.05

表10

3.温度

估计

因变量:

适宜的条件

温度

均值

标准误差

95%置信区间

下限

上限

B1（30℃）

42.917

2.261

38.277

47.556

B2（35℃）

33.917

2.261

29.277

38.556

B3（40℃）

20.167

2.261

15.527

24.806

表11

成对比较

因变量:

适宜的条件

（I）温度

（J）温度

均值差值（I-J）

标准误差

Sig.a

差分的95%置信区间a

下限

上限

B1（30℃）

B2（35℃）

9.000*

3.198

.009

2.439

15.561

B3（40℃）

22.750*

3.198

.000

16.189

29.311

B2（35℃）

B1（30℃）

-9.000*

3.198

.009

-15.561

-2.439

B3（40℃）

13.750*

3.198

.000

7.189

20.311

B3（40℃）

B1（30℃）

-22.750*

3.198

.000

-29.311

-16.189

B2（35℃）

-13.750*

3.198

.000

-20.311

-7.189

基于估算边际均值

*.均值差值在.05级别上较显著。

a.对多个比较的调整：

最不显著差别（相当于未作调整）。

表12

单变量检验

因变量:

适宜的条件

平方和

均方

Sig.

偏Eta方

非中心参数

观测到的幂a

对比

3150.500

1575.250

25.676

.000

.655

51.351

1.000

误差

1656.500

61.352

F检验温度的效应。

该检验基于估算边际均值间的线性独立成对比较。

a.使用alpha的计算结果=.05

表13

4.原料*温度

因变量:

适宜的条件

原料

温度

均值

标准误差

95%置信区间

下限

上限

B1（30℃）

34.500

3.916

26.464

42.536

B2（35℃）

18.250

3.916

10.214

26.286

B3（40℃）

18.000

3.916

9.964

26.036

B1（30℃）

49.000

3.916

40.964

57.036

B2（35℃）

37.500

3.916

29.464

45.536

B3（40℃）

15.500

3.916

7.464

23.536

B1（30℃）

45.250

3.916

37.214

53.286

B2（35℃）

46.000

3.916

37.964

54.036

B3（40℃）

27.000

3.916

18.964

35.036

"在此之后"检验

原料