单因素下的方差分析

  在方差分析中,有三个基本的假设:

  (1) 正态假设。对于因素的每个水平,其观测值都是来自正态总体的随机样本;

  (2) 方差齐次假设。各个总体的方差相同;

  (3) 独立假设。观测值之间都是独立的。

  设试验中的因素A,有r个水平A1,A2,...,An,在每个水平下进行试验得到结果xi1,xi2,...,xini,i=1,2,...,r,其被看作是来自第i个正态总体xi~N(μi,σ2),其中参数未知且每个样本都独立。从而单因素分析的数学模型可以表示为一种线性模型。

  ?????xij=μ+αi+Eij;i=1,2,...,r,j=1,2,..,niEij~N(0,σ2)且相互独立∑ri=1niαi=0

 

  其中,μ是所有总体的均值,αi=μi?μ称为第i个水平的效应,Eij是随机误差。

  1.正态性检验

  在R语言中,使用Shapiro.test(x)可以对数据x进行正态性检验,参数x是要检验的数据集,它是长度在3~5000之间的向量。

  2.方差齐次性检验

  该方法是要检验数据在不同水平下,其方差是否相等。在R语言中,使用Bartlett.test()来实现。

  方差分析的目的是,要比较因素A的r个水平下,试验结果是否有显著差异。以样本均值作为检验的标准,写出检验假设:

  H0:α1=α2=...=α1,H1:α1,α2,...,αr不全相等

  如果拒绝原假设H0,说明样本来自不同的正态总体,则由因素A的各个水平所造成均值的差异有统计意义;

  如果不能拒绝原假设H0,说明样本来自相同的正态总体,因素的不同水平之间无差异。

  案例1,某银行规定VIP客户的月均账户余额要达到100万元,并以此作为比较各个分行业绩的一项指标。这里的“分行”即为因子,账户余额是所要检验的指标,先从三个分行(对应三个水平A1、A2、A3)中,分别随机抽取7个VIP客户的账户,数据列在表(1)中。

表(1) 银行的三个分行A1、A2、A3

  (a)正态性检验

  //zheng.R  

x1=c(103,101,98,110,105,100,106)
x2=c(113,107,108,116,114,110,115)
x3=c(82,92,84,86,84,90,88)
shapiro.test(x1)
shapiro.test(x2)
shapiro.test(x3)

  效果如下:

图(1) 正态性检验的结果

  由图(1)知,P(A1)=0.948 > 0.05,不能拒绝原假设,

   P(A2)=0.4607 > 0.05,不能拒绝原假设,

   P(A3)=0.7724 > 0.05,不能拒绝原假设,

  而原假设H0是变量x服从正态分布,即A1、A2、A3都服从正态分布。

  (b)方差齐次性检测

  //qi.R  

#方差齐性检验
x=c(x1,x2,x3)
account=data.frame(x,A=factor(rep(1:3,each=7)))
bartlett.test(x~A,data=account)

效果如下:

图(2) 方差齐次性检测

  由于P=0.9341 > 0.05,不能拒绝原假设,而原假设H0是样本是“齐次的”,即三个样都是等方差的。

  (c) 单因素分析

  当数据符合正态性,和方差齐次之后,使用aov()就可以进行方差分析了。

  //fen.R  

a.aov=aov(x~A,data=account)
summary(a.aov)
plot(account$x~account$A)

如图(3)、图(4)所示:

图(3) 方差分析结果

图(4) 三个分行的箱线图

时间: 2024-10-03 13:38:45

单因素下的方差分析的相关文章

两因素重复测量方差分析,史上最详细SPSS教程!

原文地址 http://www.sohu.com/a/202657022_489312 2017-11-06 18:34 一.问题与数据 研究者想知道短期(2周)高强度锻炼是否会减少C反应蛋白(C-Reactive Protein, CRP)的浓度. 研究者招募了12名研究对象,并让研究对象参与两组试验:对照试验和干预试验.在对照试验中,研究对象照常进行日常活动:在干预试验中,研究对象每天进行45分钟的高强度锻炼,每组试验持续2周,两组试验中间间隔足够的时间. CRP的浓度在每组试验中共测量了3

SAS中用单因素ANOVA研究不同疗法对焦虑症的有效性

原文链接:http://tecdat.cn/?p=10042 本教程将介绍如何使用SAS进行单因素方差分析. 我们使用的数据可以在这里下载. 我们想研究不同疗法对焦虑症的有效性.我们收集了以下类别的75个主题的样本: 无处理(ñ1个n1 = 27). 生物反馈(ñ2n2 = 24). 认知行为治疗(n3n3 = 24). 因变量是焦虑水平.零假设是所有三个均值均相等. H0H0:μ1=μ2=μ3μ1=μ2=μ3 HAHA:至少有两个μiμi的是不同 如果甚至两种方法均存在显着差异,我们将拒绝原假

单目录下多文件 makefile编写

makefile很久就接触过了,但是一直没怎么深入的去学习和总结:在项目中我也只是看看makefile或者修改部分语句,全部自己动手写的话还真没有:知识在于沉淀,这句说的非常好,所以现在把自己理解的东西,记录下来,以便后面查阅: 这篇blog要分享的是在单目录下多文件的makefile编写,首先说明当前目录下有多少文件:fun.h   fun.c  main.c  makefile:其中*.c 文件都要依赖 *.h文件: 首先常规编译: 预处理期:gcc  -E  -o fun.i  fun.c

单因素特征选择--Univariate Feature Selection

An example showing univariate feature selection. Noisy (non informative) features are added to the iris data and univariate feature selection(单因素特征选择) is applied. For each feature, we plot the p-values for the univariate feature selection and the cor

简单的单级下拉菜单实现

效果预览 这里用纯CSS实现一个简单的单级下拉菜单,鼠标划过菜单标题时显示下拉菜单. 基本思路 在菜单处于光标之下时显示菜单,其余时候隐藏.这个过程首先得想到 :hover 伪类来实现,但是直接的利用伪类还不行,因为普通状态下菜单隐藏了,没办法触发 :hover ,不过下拉菜单隐藏了不是还有菜单标题么,菜单标题处于光标之下时显示菜单就能够显示出来了,但是鼠标移到菜单上的时候下拉菜单又会消失,这个时候想到给下拉菜单也加上 :hover 显示自身,但是万一脑洞大开,要做分体菜单怎么办,这个时候把它放

单例下并发的情况下访问私有变量的一个典型问题

一个小问题,单例下 并发的情况下访问私有变量的典型问题,以下是一个小Demo public class dssd { static void Main(string[] args) { var t = Task.Factory.StartNew(() => { Thread.Sleep(1000); var obj2 = ImportProvider.Instance; Arg ds = new Arg() { TypeSS = 2 }; obj2.Import(ds); }); var obj

第 10 章 表单元素[下]

学习要点: 1.其他元素 2.输入验证 主讲教师:李炎恢 本章主要探讨 HTML5 中表单中剩余的其他元素,然后重点了解一下表单的输入验证功能. 一.其他元素 表单元素还剩下几个元素没有讲解,包括下拉框列表 select.多行文本框 textarea.和 output 计算结果元素. 元素名称 说明 select 生成一个下拉列表进行选择 optgroup 对 select 元素进行编组 option select 元素中的项目 textarea 生成一个多行文本框 output 表示计算结果

Nginx 单IP下 配置多个server https 的问题

我们的负载均衡架构是这样的: 请求 -> 网络 -> LVS集群 -> Nginx 集群 -> APP LVS使用FULLNAT模式,每台Nginx 机器只有一个IP(内网IP),LVS也是把流量转到这个IP.如果Nginx想对多个域名使用https,比如两个域名 wandoujia.com 和 wandoulabs.com ,是可能有问题的. 看下面的配置(两个server写在不同文件中,用 include * 加载): server { listen 80; listen 44

表单(下)-EasyUI Spinner 微调器、EasyUI Numberspinner 数值微调器、EasyUI Timespinner 时间微调器、EasyUI Slider 滑块

EasyUI Spinner 微调器 扩展自 $.fn.validatebox.defaults.通过 $.fn.spinner.defaults 重写默认的 defaults. 微调器(spinner)把可编辑的文本框和两个小按钮结合起来,允许用户从某个范围的值中进行选择.与组合框(combobox)相似,微调器(spinner)允许用户输入一个值,但是它买i有下拉列表.微调器(spinner)是创建其他微调器组件(比如:数值微调器 numberspinner.时间微调器 timespinne