平台型网站如何推广,网站首页图片轮转,网络培训心得体会5篇,网页制作用什么工具文章目录 前言方差分析方差分析的原理什么是方差分析误差分解 单因子方差分析数学模型效应检验 练习 前言
本片开始介绍有关方差分析的知识。 方差分析
方差分析的基本原理是在20世纪20年代由英国统计学家Ronald A.Fisher在进行实验设计时为解释实验数据而首先引入的。方差分… 文章目录 前言方差分析方差分析的原理什么是方差分析误差分解 单因子方差分析数学模型效应检验 练习 前言
本片开始介绍有关方差分析的知识。 方差分析
方差分析的基本原理是在20世纪20年代由英国统计学家Ronald A.Fisher在进行实验设计时为解释实验数据而首先引入的。方差分析是一种统计方法。目前方差分析方法广泛应用于分析心理学、生物学、工程和医药的实验数据。本章首先介绍方差分析的基本原理然后介绍只涉及一个因变量的单因子独立样本和双因子独立样本的方差分析方法
方差分析的原理
什么是方差分析
方差分析(analysis of varianceANOVA) 是分析类别自变量对数值因变量影响的一种统计方法 。自变量对因变量的影响也称为自变量效应effect。由于影响效应的大小表现在因变量的误差里有多少是由自变量造成的因此方差分析是通过对数据误差的分析来检验这种效应是否显著。 研究分类型自变量对数值型因变量的影响 一个或多个分类型自变量 两个或多个 (k 个) 处理水平或分类 一个数值型因变量 有单因子方差分析和双因子方差分析 单因子方差分析涉及一个分类的自变量 双因子方差分析涉及两个分类的自变量
例题 数据example8_1.RData为分析小麦品种对产量的影响一家研究机构挑选了3个小麦品种品种1、品种2、品种3然后选择条件和面积相同的30个地块每个品种在10个地块上试种实验获得的产量数据如下表所示。 分析“小麦品种” 对“产量”的影响 在上表中“小麦品种”是类别变量称为实验的因子factor品种1、品种2、品种3是因子的3个不同取值称为处理treatment或水平level 这里的“地块”就是接受处理的对象或实体称为实验单元experiment unit 产量则是因变量每个地块上获得的产量是样本观察值。分析小麦品种对产量影响的统计方法就是方差分析。 如果只分析品种一个因子对产量的影响则称为单因子方差分析(one-way analysis of variance) 如果两个因子对产量的单独影响但不考虑它们对产量的交互效应(interaction)则称为只考虑主效应(main effect)的双因子方差分析或称为无重复双因子分析(two-factor without replication) 如果除了考虑两个因子对产量的单独影响外还考虑二者对产量的交互效应则称为考虑交互效应的双因子方差分析或称为可重复双因子分析(two-factor with replication)
误差分解
怎样分析小麦品种对产量是否有显著影响呢由于品种对产量的影响效应体现在产量数值的误差里因此分析时首先应从数据误差的分析入手。方差分析的原理就是通过对数据误差的分析来判断类别自变量小麦品种对数值因变量产量的影响效应是否显著。 怎样分析数据的误差呢从上表可以看出每个品种每种处理各有10个实验数据这些数据实际上是从每个品种的产量总体中抽出来的样本量为10的随机样本共获得3个样本的30个实验数据.可以看出这30个产量数据是不同的。 因此有如下定义 总误差(total error) 反映全部观测数据的误差 所抽取的全部30个地块的产量之间的误差
总误差可能是由于不同处理小麦的不同品种造成的也可能是由于其他随机因素造成的。 处理误差(treatment error)—组间误差(between-group error)处理误差来自不同的处理之间 由于不同处理造成的误差它反映了处理(品种)对观测数据(产量)的影响因此称为处理效应(treatment effect) 随机误差(random error)—组内误差(within-group error) 由于随机因子造成的误差也简称为误差(error)
由其他随机因素对观测数据造成的误差称为随机误差random error也简称为误差error。本例中随机误差反映了除品种外其他随机因素对产量的影响。 由于随机误差主要存在于每种处理的内部因此有时也称为组内误差within-group error)。 随机误差(random error)—组内误差(within-group error) 由于随机因子造成的误差也简称为误差(error)
数据的误差用平方和(sum of squares)表示记为SS 总平方和(sum of squares for total)记为SST 反映全部数据总误差大小的平方和 例如全部30个地块的产量之间的误差平方和就是总平方和它反映了全部产量的总离散程度。 处理平方和(treatment sum of squares)记为SSA 反映处理误差大小的平方和。(例如不同品种之间产量的误差平方和就是处理平方和) 也称为组间平方和(between-group sum of squares) 误差平方和(sum of squares of error)记为SSE 反映随机误差大小的平方和称为误差平方和 也称为组内平方和(within-group sum of squares)
全部数据的总误差平方和被分解成两部分一部分是处理平方和SSA另一部分是误差平方和SSE。很显然这三个误差平方和的关系为SSTSSASSE。 数据误差的来源及其分解过程可用下面的图来表示 方差分析的基本原理就是要分析数据的总误差中处理误差是否显著存在。如果处理(小麦品种)对观测数据(产量)没有显著影响意味着处理误差不显著。这时每种处理所对应的总体均值(i)应该相等。 如果存在处理误差每种处理所对应的总体均值(i)至少有一对不相等 就本例而言在只考虑品种一个因子的情况下方差分析也就是要检验下面的假设
单因子方差分析
数学模型
设因子A有I种处理比如品种有(“品种1”、“品种2”、“品种3”)3种处理单因子方差分析可用下面的线性模型来表示 效应检验
设因子A有I种处理单因子方差分析要检验的假设为 H0 ai 0 (i1,2,…,I) 没有处理效应 H1 ai 至少有一个不等于0 有处理效应 为获得上述检验的统计量首先需要计算处理平方和SSA、误差平方和SSE然后将各平方和除以相应的自由度df以消除观测数据多少对平方和大小的影响其结果称为均方mean square也称为方差。 最后将处理均方MSA除以误差均方MSE即得到用于检验处理效应的统计量F。这一计算过程可以用方差分析表的形式来表示。下表显示了单因子方差分析表的一般形式 有了统计量就可以根据P值做出决策若Pα拒绝H0αi不全为0表示处理效应显著因子对观察值有显著影响。
例题 数据example8_2. RData沿用例7-1。检验小麦品种对产量的影响是否显著α0.05 解设小麦品种对产量的影响效应分别为a1品种1、a2品种2、a3品种3 为使用R做方差分析需要将短格式数据example8_1转为长格式数据并存为 example8_2。R代码和结果如下所示
# 将表8-1的短格式数据转为长格式数据
load(C:/example/ch8/example8_1.RData)
example8_1-cbind(example8_1,idfactor(1:10))
library(reshape)
example8_2-melt(example8_1,id.varsc(id), variable_name品种)
example8_2-rename(example8_2,c(id地块,value产量))
save(example8_2,fileC:/example/ch8/example8_2.RData)
example8_2melt()函数用法把数据转变成长格式数据 melt(data,id.vars,measure.vars,variable.name,value.name) id.vars:标识变量依旧在列上位置保持不变的变量 measure.vars度量变量需要放进同一列的变量 variable.name为新列变量取名 value.name对应值所在的变量名
为便于理解方差分析的结果可以先对样本数据做一些描述性分析。首先绘制3个品种产量的箱线图并计算各样本的均值和标准差观察不同品种产量之间的差异。R代码和结果如下所示
#绘制3个品种数据产量的箱线图
load(C:/example/ch8/example8_2.RData)
attach(example8_2)
boxplot(产量~品种,dataexample8_2,colgold,main,ylab产量, xlab品种)#计算描述统计量
my_summary-function(x){with(x,data.frame(均值mean(产量),标准差sd(产量),nlength(产量)))}
library(plyr)
ddply(example8_2,.(品种),my_summary)#ddply()函数用于对data.frame进行分组统计 从箱线图中可以看出品种2的产量明显低于品种1和品种3而品种1和品种3的产量差异似乎不大。从3个品种产量的均值也可见其差异。但要想知道不同品种之间的产量差异是否显著还需要做方差分析。R代码和结果如下所示 #方差分析表 attach(example8_2)model_1w-aov(产量~品种)summary(model_1w)#方差分析模型的参数估计
model_1w$coefficients注函数aov formuladataNULL用于拟合方差分析模型。参数 formula用于指定模型形式data为数据框。单因子方差分析模型的R表达式为y~A其中y为因变量A为因子。
方差分析表显示了品种效应和随机效应的平方和Sum Sq、自由度Df、均方Mean Sq、检验统计量值F value、检验的P值PrF。由于P0.0001580.05,拒绝原假设表示αii123至少有一个不等于0这意味着品种对产量的影响效应显著。 练习
1、一家牛奶公司有4台机器装填牛奶每桶的容量为4L。下面是从4台机器中抽取的装填量样本数据。 检验机器对装填量是否有显著差异 解提出假设 H0: 机器对装填量的影响不显著 H1机器对装填量的影响显著 #绘制出箱线图
load(C:/example/ch8/exercise8_1.RData)
exercise8_1
attach(exercise8_1)
boxplot(填装量~机器,dataexercise8_1,colgold,main,ylab填装量, xlab机器)#方差分析表
attach(example8_2)
model_1w-aov(填装量~机器)
summary(model_1w)方差分析表显示了机器效应和随机效应的平方和Sum Sq、自由度Df、均方Mean Sq、检验统计量值F value、检验的P值PrF。由于P0.01350.05,拒绝原假设机器对填装量的影响效应显著。 #方差分析模型的参数估计
model_1w$coefficients
