多元统计分析第三章假设检验及方差分析

文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持. 1 第3章 多元正态总体的假设检验与方差分析 从本章开始，我们开始转入多元统计方法和统计模型的学习。统计学分析处理的对象是带有随机性的数据。按照随机排列、重复、局部控制、正交等原则设计一个试验，通过试验结果形成样本信息（通常以数据的形式），再根据样本进行统计推断，是自然科学和工程技术领域常用的一种研究方法。由于试验指标常为多个数量指标，故常设试验结果所形成的总体为多元正态总体，这是本章理论方法研究的出发点。 所谓统计推断就是根据从总体中观测到的部分数据对总体中我们感兴趣的未知部分作出推测，这种推测必然伴有某种程度的不确定性，需要用概率来表明其可靠程度。统计推断的任务是“观察现象，提取信息，建立模型，作出推断”。 统计推断有参数估计和假设检验两大类问题，其统计推断目的不同。参数估计问题回答诸如“未知参数的值有多大?”之类的问题,而假设检验回答诸如“未知参数的值是0吗?”之类的问题。本章主要讨论多元正态总体的假设检验方法及其实际应用，我们将对一元正态总体情形作一简单回顾，然后将介绍单个总体均值的推断， 两个总体均值的比较推断，多个总体均值的比较检验和协方差阵的推断等。 3.1一元正态总体情形的回顾

一、 假设检验 在假设检验问题中通常有两个统计假设（简称假设）,一个作为原假设（或称零假设），另一个作为备择假设（或称对立假设），分别记为0H和1H。 1、显著性检验 为便于表述，假定考虑假设检验问题：设1X，2X，…,nX 来自总体),(2N的样本，我们要检验假设

0100:,:HH （3.1） 原假设0H与备择假设1H应相互排斥，两者有且只有一个正确。备择假设的意思是，一旦否定原假设0H，我们就选择已准备的假设1H。

当2已知时，用统计量nXz 文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.

2 在原假设0H成立下，统计量z服从正态分布z)1,0(~N，通过查表，查得)1,0(N的上分位点2z。 对于检验问题（，我们制定这样一个检验规则（简称检验）： 当2zz时，拒绝0H； 当2zz时，接受0H。 （3.2） 我们称2z为临界值，是)1,0(N的上分位点，不同的临界值代表不同的检验。称拒绝原假设0H的统计量z的范围为拒绝域，称接受0H的统计量z的范围为接受域，因此给出一个检验，就是给出一个拒绝域。 2、两类错误 由于样本具有随机性，因此在根据样本进行判断时，有可能犯两种类型的错误。一类错误是，原假设0H本来正确，但按检验规则却作出了拒绝0H的判断，这类错误称为第一类错误（弃真错误），其发生的概率2zzP称为犯第一类错误的概率；另一类错误时，原假设0

H

本来不正确，但按检验规则却作出了接收0H的判断，这类错误称为第二类错误（存伪错误），其发生的概率称为犯第二类错误的概率，记为。 同时控制这两类错误是困难的，当时在样本容量n固定的条件下，要使和同时减小，通常是不可能的。在假设检验的应用中,由奈曼(NEYMAN)与皮尔逊(PEARSON)提出了一个原则，即在控制犯第一类错误的概率条件下，尽量使犯第二类错误的概率小,这种检验问题, 称为显著性检验问题。根据这一原则，原假设受到保护，不至于被轻易拒绝，一旦检验结果拒绝了原假设，则表明拒绝的理由是充分的，如果接受了原假设，则只是表明拒绝的理由还不充分，未必意味着原假设就是正确的。所以，在实际问题中，为了通过样本观测值对某一猜测取得强有力的支持，通称我们把这一猜测的否定作为原假设，而把猜测本身作为备择假设。 3、关于检验的p值 下面，我们再介绍进行检验的另一种方式——p值，我们就以(，对于样本，我们通过统

计量，计算出nxz00，是一确定值，这里的x是样本观测值的均值，再由统计量z服

从正态分布z)1,0(~N,计算}{0zzP为检验的p值。 由于2zz等价于p=}{0zzP2zzP，所以检验规则可以表述为： 文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持. 3 当p时，拒绝0H； 当p时，接受0H。接受0H。 （3.3） 上述p值的检验规则与（，p值越小，拒绝原假设的理由就充分。通常SAS等软件的计算机输出一般只给出p值，由你自己给定的值来判断检验结果 二、单一变量假设检验的回顾 1、 单个正态总体均值的检验 考虑假设检验问题：设1X，2X，…,nX 来自总体),(2N的样本，我们要检验假设 （1） 总体方差2已知 构造统计量

在原假设0H成立下， z服从正态分布z)1,0(~N，可得这样一个检验规则： 当2zz时，拒绝0H； 当2zz时，接受0H。 （2） 总体方差2未知 构造统计量

在原假设0H成立下，t服从自由度为1n的t分布t)1(~nt可得这样一个检验规则：

当)1(2ntt时，拒绝0H； 当)1(2ntt时，接受0H。 （ 2、 两个正态总体均值的比较检验 考虑假设检验问题 211210:,:HH （

设121,,,nXXX是取自总体),(211N的容量为1n的样本，221,,,nYYY是取自),(222N的容量为2n的样本，给定显著性水平。

（1） 两个总体方差21和22已知 文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持. 4 构造检验统计量222121nnYXz （

在原假设0H成立下， z服从正态分布z)1,0(~N，检验规则为： 当2zz时，拒绝0H； 当2zz时，接受0H。 （2） 两个总体方差21和22都未知，但21=22=2 用样本方差s代替，构造检验统计量 在原假设0H成立下，t服从正态分布t)2(~21nnt，检验规则为： 当)2(212nntt时，拒绝0H； 当)2(212nntt时，接受0H。 3、多个正态总体均值的比较检验（方差分析） 设k个正态总体分别为),(21N，),(22N，…, ),(2kN从k个总体取in

个独立样本如下: 考虑假设检验问题 假设0H成立条件下,构造检验统计量为:

这里kiiiXXnSSA12)(称为组间平方和；kiiijnjXXSSEi12)(1)(称为组内平方和； kiijnjXXSSTi12)(1)(称为总平方和。其中iXnjijiXn1)(1,XkinjijXn11)(1

knnnn21 给定检验水平，查F分布表，使FFP，可确定出临界值F，

再利用样本值计算出F值，若FF，则拒绝0H，否则不能拒绝0H。 附注：多元假设检验与SAS过程

本章的主要内容是多元假设检验和方差分析，其中的计算一般都很复杂，可用国际上著名的专业软件——SAS软件计算。SAS中有GLM，ANOVA和NESTED等过程可用方差分析。其中GLM过程最常用。 SAS的GLM过程采用了一般线性模型： 文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持. 5 在方差分析问题中，变量 mxx...1是示性变量，即只取0或1的变量。GLM过程对每一因子的每一水平，通过CLASS语句产生1个示性变量，也称分类变量。 GLM过程主要有四个语句：PROC GLM，CLASS，MODEL和LSMEANS语句。 PROC GLM语句 用以调用GLM过程，有许多选项，一般形式是： Proc glm [data=数据集名称] [outstat=输出的统计量] [order=formatted|freq|data|internal]； CLASS语句 说明哪些变量是分类变量。方差分析中的因素都是分类变量，如： Class V1 V2 V3； 此语句指示计算机把因子V1，V2 ，V3作为分类变量，可以是字符型变量或数字型变量。如果是字符型变量，长度限于10个字符以内。 MODEL语句 语句中等号前是响应变量，如： Model Y=A； 单因子ANOVA Model Y=A B C； 主效应模型 Model Y=A B A*B； 含交互效应的因子模型 Model Y1 Y2=A B； 多因子方差模型MANOVA LSMEANS语句 用以求待估参数的最小二乘估计。 Lsmeans A B A*B； MANOVA语句 用以说明是做多元方差分析。 3.2 均值等于常数向量的检验

在经济生产、管理决策中的很多实际问题，通常要选取多个指标进行考察，根据历史数据，将p项指标的历史平均水平记作0，考虑新的p项指标平均值是否与历史数据记载的平均值有明显差异？若有差异，进一步分析差异主要在哪些指标上，先看下面的实例： 例3.1测量20名健康女性排汗量1x、钠含量2x、钾含量3x得表3.1。问健康女性1x、2x、

3x的均值是不是4、50、10？

表3-1 20名健康女性排汗量1x、钠含量2x、钾含量3x数据 排汗量1x 钠含量2x 钾含量3x 3.7 48.5 9.3 5.7 65.1 8.0 3.8 47.2 10.9 3.2 53.2 12.0 3.1 55.5 9.7 4.6 36.1 7.9 2.4 24.8 14.0 7.2 33.1 7.6 6.7 47.4 8.5 5.4 54.1 11.3 3.9 36.9 12.7 4.5 58.8 12.3 3.5 27.8 9.8 4.5 40.2 8.4