《概率论与数理统计.ppt》由会员分享,可在线阅读,更多相关《概率论与数理统计.ppt(157页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、1,数 理 统 计,2,第八章 假设检验,关键词:假设检验正态总体参数的假设检验 分布拟合检验 秩和检验,3,1 假设检验,统计推断的另一类重要问题是假设检验问题。它包括(1)已知总体分布的形式,但不知其参数的情况,提出参数的假设,并根据样本进行检验.(2)在总体的分布函数完全未知的情况下,提出总体服从某个已知分布的假设,并根据样本进行检验.,4,例1 设某种清漆的9个样品,其干燥时间(以小时计)分别为: 6.0 5.7 5.5 6.5 7.0 5.8 5.2 6.1 5.0根据以往经验,干燥时间的总体服从正态分布N(6.0, 0.36),现根据样本检验均值是否与以往有显著差异?,例2 一种摄
2、影药品被其制造商声称其贮藏寿命是均值180天、标准差不多于10天的正态分布。某位使用者担心标准差可能超过10天。他随机选取12个样品并测试,得到样本标准差为14天。根据样本有充分证据证明标准差大于10天吗?,例3 孟德尔遗传理论断言,当两个品种的豆杂交时,圆的和黄的、起皱的和黄的、圆的和绿的、起皱的和绿的豆的频数将以比例9:3:3:1发生。在检验这个理论时,孟德尔分别得到频数315、101、108、32、这些数据提供充分证据拒绝该理论吗?,5,例1 设某种清漆的9个样品,其干燥时间(以小时计)分别为: 6.0 5.7 5.5 6.5 7.0 5.8 5.2 6.1 5.0根据以往经验,干燥时间
3、的总体服从正态分布N(6.0, 0.36),现根据样本检验均值是否与以往有显著差异?,由于作出决策的依据是一个样本,因此,可能出现“实际上原假设成立,但根据样本作出拒绝原假设”的决策。这种错误称为“第一类错误”,实际中常常将犯第一类错误的概率控制在一定限度内,即事先给定较小的数(02)个水平,n个对象参与了试验。假定对应于因素第j个水平的组中有 个试验对象,响应变量数据为,通常假定,74,检验假设,假设等价于,75,(二)平方和分解,76,证明:,77,78,79,单因素试验方差分析表,80,81,例1 设有5种治疗荨麻疹的药,要比较它们的疗效。假设将30个病人分成5组,每组6人,令同组病人使
4、用一种药,并记录病人从使用药物开始到痊愈所需时间,得到下面的记录:(=0.05),82,这里药物是因子,共有5个水平,这是一个单因素方差分析问题,要检验的假设是“所有药物的效果都没有差别”。,83,未知参数的估计,84,85,2 双因素试验的方差分析,例 假设某药物研究者为检验a,b两种化学物质的抗癌效果,要做动物试验。通常的作法是:将一些患有某种癌的白鼠随机地分成三组。其中两组分别注射a,b两种化学物质,而第三组不作处理,作为对照。记第一组:注射a物质,第二组:注射b物质,第三组:不做处理。经过一段时间观察后,得到寿命数据。在这个药物试验中,如果白鼠的性别有可能对其寿命有显著的影响。这时应该
5、考虑将“性别”作为一个因素“双因素试验”。因素A:药物,三个水平;因素B:性别,二个水平;两个因素共有236种组合。,86,(一) 双因素等重复试验的方差分析,87,88,分别检验假设,89,90,91,92,93,双因素试验的方差分析表,94,例3 为了比较3种松树在4个不同的地区的生长情况有无差别,在每个地区对每种松树随机地选取5株,测量它们的胸径,得到的数据列表如下。,松树数据表,95,这是一批等重复的两种方式分组数据,记树种因素为A,地区因素为B,则A因素有3个水平,B因素有4个水平,总共有12个水平组合,每个组合(单元)有5个重复观测。,将树的胸径作为度量树的生长情况是否良好的数值指
6、标,我们的目标是:由以上数据来判断不同树种及不同地区对松树的生长情况是否有影响(好或坏)?,这里要考虑的影响有三种:树种的单独影响(A因素主效应),地区的单独影响(B因素主效应),以及不同树种和不同地区的结合所产生的交互影响(AB因素的交互效应)。这是一个典型的等重复双因素方差分析模型。,96,输出各单元总和及因素水平总和:,松树数据的总和表,97,双因素方差分析表,98,进一步考查A因素不同水平的均值。注意到A因素的第二水平为最大:23.55,而第三水平的均值为最小:17.65,可以认为树种2的生长情况优于树种3。能够得出这个结论,得益于观测的等重复性。,然后再来看B因素的主效应,即在扣除松
7、树种类的效应后,不同地区对树的胸径的影响。由方差分析表知, B因素的主效应不显著,即不同的地区对树的胸径没有显著影响。,最后来看AB因素的交互效应,即在扣除两种效应后,由不同树种和不同地区的结合而产生的对树的胸径的影响,这种影响可以解释为某些地区特别适合(或特别不适合)某个树种的生长。结果也不显著。,首先来看A因素主效应,即在扣除地区效应后, 松树的不同种类对树的胸径的影响。由方差分析表可以看出,A因素主效应是显著的,即松树的不同种类对树的胸径有显著影响。,99,(二) 双因素无重复试验的方差分析,100,101,分别检验假设,102,103,104,105,双因素无重复试验的方差分析表,10
8、6,例4 假定对3个小麦品种和3块试验地块进行区组设计试验,得到如下的数据:,表 小麦品种区组试验数据,107,在这个问题中我们关心的是小麦的不同品种之间在产量上的差异。由于地块不同对小麦的产量也会有影响,因此在比较试验结果时,要扣除地块的影响之后再来比较品种的差异。假定品种与地块之间无交互效应,则可对上述数据进行双因素可加效应模型的方差分析。,108,双因素无重复试验的方差分析表,109,在这个问题中我们所关心的是因素A的效应,由方差分析表知,原假设不成立,即认为小麦品种的产量之间有显著差异。,在这里,品种3的单产最高,而品种1的产量最低,因此可以断定品种3明显地优于品种1。,110,3 一
9、元线性回归分析,一、确定性关系: 当自变量给定一个值时,就确定应变量的值与之对应。 如:在自由落体中,物体下落的高度h与下落时间t之间有函数关系:,变量与变量之间的关系,111,二、相关性关系: 变量之间的关系并不确定,而是表现为具有随机性的一种“趋势”。即对自变量x的同一值,在不同的观测中,因变量Y可以取不同的值,而且取值是随机的,但对应x在一定范围的不同值,对Y进行观测时,可以观察到Y随x的变化而呈现有一定趋势的变化。,如:身高与体重,不存在这样的函数可以由身高计算出体重,但从统计意义上来说,身高者,体也重。 再如:父亲的身高与儿子的身高之间也有一定联系,通常父亲高,儿子也高。,回归分析研
10、究相关性关系的最基本,应用最广泛的方法。,112,(一)一元线性回归,113,在实际问题中,回归函数(x)一般是未知的,需要根据试验数据去估计。,114,115,一元线性回归要解决的问题:,116,(二)a,b的估计最小二乘估计,117,正规方程系数行列式,118,在误差为正态分布假定下,最小二乘估计等价于极大似然估计。,事实上,似然函数,119,120,121,例1 K.Pearson收集了大量父亲身高与儿子身高的资料。其中十对如下:,求Y关于x的线性回归方程。,122,123,(三)误差方差的估计,124,125,例2 求例1中误差方差的无偏估计。,126,(1)影响Y取值的,除了x,还有
11、其他不可忽略的因素;(2)E(Y)与x的关系不是线性关系,而是其他关系;(3)Y与x不存在关系。,(四)线性假设的显著性检验,采用最小二乘法估计参数a和b,并不需要事先知道Y与x之间一定具有相关关系,即使是平面图上一堆完全杂乱无章的散点,也可以用公式求出回归方程。因此(x)是否为x的线性函数,一要根据专业知识和实践来判断,二要根据实际观察得到的数据用假设检验方法来判断。,若原假设被拒绝,说明回归效果是显著的,否则,若接受原假设,说明Y与x不是线性关系,回归方程无意义。回归效果不显著的原因可能有以下几种:,127,128,例3 检验例1中回归效果是否显著,取=0.05。,129,(五)回归系数b
12、的置信区间,当回归效果显著时,常需要对回归系数b作区间估计。,130,(六)回归函数 函数值的点估计和置信区间,131,132,(七)Y的观察值的点预测和预测区间,133,134,135,注:在预测时, 一定要落在已有的 的数据范围内部,否则预测常常没有意义。,136,例4,在例1中F.Galton曾断言“儿子身高会受到父亲身高的影响,但身高偏离父代平均水平的父亲,其儿子身高的影响有回归到子代平均水平的趋势。”试问例1这组数据能证实这一论断吗(=0.05)?并给出x=69吋时,y的预测区间。,(1)回归到平均水平的趋势,即检验,137,138,例5 合金钢的强度y与钢材中碳的含量x有密切关系。
13、为了冶炼出符合要求强度的钢常常通过控制钢水中的碳含量来达到目的,为此需要了解y与x之间的关系。其中x:碳含量() y:钢的强度(kg/mm2)数据见下:,(1)画出散点图;(2)设(x)=a+bx,求a,b的估计;(3)求误差方差的估计,画出残差图;(4)检验回归系数b是否为零(取=0.05);(5)求回归系数b的95置信区间;(6)求在x=0.06点,回归函数的点估计和95置信区间;(7)求在x=0.06点,Y的点预测和95区间预测。,139,0.03 0.05 0.07 0.09 0.11 0.13 0.15 0.17 0.19,56 54 52 50 48 46 44 42 40 38,
14、(1)合金钢的强度y与钢材中碳的含量x的散点图,140,141,0.03 0.05 0.07 0.09 0.11 0.13 0.15 0.17 0.19 x,0,e,142,0.03 0.05 0.07 0.09 0.11 0.13 0.15 0.17 0.19,56 54 52 50 48 46 44 42 40 38,合金钢的强度y与钢材中碳的含量x的回归直线图,143,144,145,(八)可化为一元线性回归的例子,实际中常会遇到很复杂的回归问题,但在某些情况下,通过适当的变量变换,可将其化为一元线性回归来处理。下面是三种常见的可转化为一元线性回归的模型。,146,147,4 多元线性回
15、归,在实际问题中,影响Y(因变量)的因素(自变量)往往不止一个,设有,148,149,150,151,152,153,例6 某公司在各地区销售一种特殊化妆品。该公司观测了15 个城市在某月内对该化妆品的销售量Y及各地区适合使用该化妆品的人数X1和人均收入X2,得到数据如下:,表1.1.2 化妆品销售的调查数据,154,化妆品销售的调查数据(续),155,156,由回归方程可知,若固定人均收入不变,则人数每增加1千人,销售量增加0.496箱;若固定人数不变,收入每增加1元,销售量增加0.0092箱。,多元线性回归也可以像一元线性回归一样,检验模型的回归效果是否显著。所不同的是,在模型的回归效果显著的情况下,还要检验每个自变量对因变量的效应是否显著,不显著就要剔除,通常用逐步回归法可以使回归方程变得简洁、明确、显著。在此基础上可以对给定点处对应的Y进行点预测和区间预测。所有这些都可以通过SAS软件实现。,2019/10/22,课件结束!,
限制150内