《分布检验和拟合优度检验》由会员分享，可在线阅读，更多相关《分布检验和拟合优度检验（43页珍藏版）》请在金锄头文库上搜索。
1、第八章 分布检验和拟合优度2检验分布检验和拟合优度检验第八章 分布检验和拟合优度2检验Kolmogorov-Smirnov单样本检验及一些正态性检验1235Kolmogorov-Smirnov两样本分布检验Pearson2 拟合优度检验分布检验和拟合优度分布检验和拟合优度检验检验第一节 KS单样本分布检验v一、适用范围vKolmogorov-Smirnov检验常译为柯尔莫哥洛夫-斯米尔诺夫检验，简写为K-S检验，亦称D检验法，也是一种拟合优度检验法。K-S单样本检验主要用来检验一组样本数据的实际分布是否与某一指定的理论分布相符合。v二、基本原理和方法v1、基本原理：v这种检验主要是将理论分布下
2、的累计频数分布与观察到的累计频数分布相比较，找出它们间最大的差异点，并参照抽样分布，定出这样大的差异是否处于偶然。分布检验和拟合优度分布检验和拟合优度检验检验2、方法v用Fn(x)表示样本量为n的随机样本观察值的累计分布函数，且Fn(x)=i/n（i是等于或小于x的所有观察结果的数目，i=1，2，n）。F(x)表示理论分布的累计概率分布函数。K-S单样本检验通过样本的累计分布函数Fn(x)和理论分布函数F(x)的比较来做拟合优度检验。检验统计量是F(x)与Fn(x)间的最大偏差Dn：v若对每一个x值来说，Fn(x)与F(x)都十分接近，则表明实际样本的分布函数与理论分布函数的拟合程度很高。分布
3、检验和拟合优度分布检验和拟合优度检验检验三、检验步骤v1.建立假设组：vH0：Fn(x)=F(x)vH1：Fn(x)F(x)v2.计算样本累计频率与理论分布累计概率的绝对差，令最大的绝对差为Dn；v3.用样本容量n和显著水平a在附表11中查出临界值Dna；v4.通过Dn与Dna的比较做出判断，若DnDna，则认为拟合是满意的。分布检验和拟合优度分布检验和拟合优度检验检验四、实例v例8.1：正态拟合。v某织布厂工人执行的生产定额（织机每小时生产织物的米物）情况如表8-1，试检验这些样本数据能否作正态拟合？v表8-1工人执行生产定额情况分组表按定额执行情况分组工人数3.753.754.252020
4、4.254.254.753723724.754.755.254984985.255.255.751031035.755.756.257 710001000分布检验和拟合优度分布检验和拟合优度检验检验例8.1 正态拟合v解：首先，由于做正态拟合的均值、标准差未知，因此，先计算样本均值和标准差，再做正态拟合。通过对样本资料的计算得：=4.85；s=0.352,分别作为和 的估计值，建立假设：vH0：样本数据服从均值为4.85，标准差为0.352的正态分布vH1：样本数据不服从均值为4.85，标准差为0.352的正态分布v计算资料列如表8-2：分布检验和拟合优度分布检验和拟合优度检验检验表8-2v表
5、8-2正态拟合计算表X的组限标 准 化标准正态概率累计概率 (理论概率) 累计工人数实 际 累计 频 率(2)-(4)的绝 对值甲乙（1）（2）（3）（4）（5）不足4.25-1.700.0450.0450.0450.04520200.0200.0200.0250.0254.25-4.75-1.70-0.280.3450.3900.3903923920.3920.3920.0020.0024.75-5.25-0.281.140.4830.4830.8730.8738908900.8900.8900.0170.0175.25-5.751.142.560.1220.1220.9950.995993
6、9930.9930.9930.0020.0025.75-6.252.56-+0.0050.0051.0001.000100010001.0001.0000.0000.000合 计1.0001.000-分布检验和拟合优度分布检验和拟合优度检验检验例8.1v根据表8-2中第(5)列数据，取最大绝对差数1ooo=0.025作为检验统计量。若取a=0.05，n=1000，从临界值表中查得：v。v因为1ooo0.043，故认为样本数据所提供的信息无法拒绝H0，即接受H0，认为可做正态分布的拟合。v检验法是一种精确分布的方法，不受观察次数多少的限制。这个方法可应用于分组或不分组的情形。检验量Dn也可用于检
7、验随机样本是否抽自某特定的总体的问题。分布检验和拟合优度分布检验和拟合优度检验检验第二节 K-S双样本分布检验v一、适用范围vK-S双样本检验主要用来检验两个独立样本是否来自同一总体（或两样本的总体分布是否相同）。其单尾检验主要用来检验某一样本的总体值是否随机地大于（或小于）另一样本的总体值。v二、理论依据和方法v1、理论依据：v与K-S单样本检验相似，K-S双样本检验是通过两个样本的累计频数分布是否相当接近来判断Ho是否为真。如果两个样本间的累计概率分布的离差很大，这就意味着两样本来自不同的总体，就应拒绝Ho。分布检验和拟合优度分布检验和拟合优度检验检验2、方法v如果令S1（x）表示第一个样
8、本观察值的累计概率分布函数，S2（x）表示另一个样本观察值的累计概率分布函数，那么K-S双样本的单尾检验统计量为：vK-S双样本的双尾检验统计量为：分布检验和拟合优度分布检验和拟合优度检验检验三、检验步骤v1、双尾检验假设：H0：S1(x)=S2(x)vH1：S1(x)S2(x)v单尾检验假设：vH0：S1(x)=S2(x)或H0：S1(x)=S2(x)vH1：S1(x)S2(x)H1：S1(x)S2(x)v2、把两组样本分别排成累计频数分布（对两个分布用相同的间隔或分类，并利用尽可能多的间隔。v3、计算检验统计量D值，如是单尾检验，应按H1的方向计算D值。分布检验和拟合优度分布检验和拟合优度
9、检验检验三、检验步骤v4、显著性检验：v小样本情况下，及n1=n2=n，n30，用附表12。对于单尾检验和双尾检验，该表列出了不同显著性水平下的临界值。v大样本情况下，n1不一定等于n2，但都小于40的双尾检验，可用附表12续表中的公式算出D的临界值。v当n1和n2都较大，但又是单尾检验时，用算式分布检验和拟合优度分布检验和拟合优度检验检验四、实例v例8.2（小样本）检验两矿的金属含量率是否相同。在甲、乙两矿坑中各抽取10个矿石样本，矿石中含有某种金属含量率（%）的资料如表8-3所示：v表7-3v解：这是一个双样本的K-S检验，根据题意，建立双侧检验假设组：v甲矿3.1 1.2 2.9 3.0
10、 0.6 2.8 1.6 1.7 3.2 1.7乙矿3.8 2.1 3.2 7.2 2.3 3.5 3.0 4.6 3.1 3.2分布检验和拟合优度分布检验和拟合优度检验检验四、实例v1、列等距分组表，计算各组次数f甲、f乙，累计次数F甲、F乙，累计频率、及其差额。计算结果列如表8-4所示。v表7-4例7.2的计算表v金属含量率（%）次数累计次数累计频率F甲/10- F乙/10f甲f乙F甲F乙F甲/10F乙/100.00.910101/1001/101.01.940505/1005/102.02.922727/102/105/103.03.93610810/108/102/104.04.901
11、10910/109/101/105.05.90010910/109/101/106.06.90010910/109/101/107.07.901101010/1010/100分布检验和拟合优度分布检验和拟合优度检验检验四、实例v2、确定检验统计量：v本例D=5/10v3、检验与判断。由于n1=n2=10，属小样本，查附表12得临界值0.05=7/10，因为D=5/10C，故在5%的显著性水平下拒绝Ho，即高分组的学生智力显著高于低分组的学生。分布检验和拟合优度分布检验和拟合优度检验检验第三节 卡方(2)拟合优度检验v一、什么是卡方(2)拟合优度检验v人们通常关心随机变量的概率分布，如：“随机变
12、量服从参数为n=10和p=2的二项分布”，这样的命题假设可以用“拟合优度检验”来检验。即设计一个检验来比较从假设的分布中抽取的样本，看所假设的分布函数与样本数据是否“拟合”。v所以，拟合检验就是检验抽取样本的总体分布与某种特定分布的符合程度，也就是检验观察值与理论数之间的紧密程度。以2分布为依据的这种检验，称为2拟合优度检验分布检验和拟合优度分布检验和拟合优度检验检验v英国统计学家Pearson(皮尔逊)于1900年首先提出了卡方统计量。v1 1、数据：由随机变量X X的N N个观测组成。这N N个观测可划分为k k类，即把X的样本空间S划分成k个互不相交的部分S1，S2，Sk，且Si与Sj相
13、互独立。即vSiSj=,(ij)，v记Oi为类i中的观测数，i=1,2,k.则第三节 卡方(2)拟合优度检验分布检验和拟合优度分布检验和拟合优度检验检验v2 2、假设条件 1.1.样本是随机的 2. 2.度量尺度至少是名义的3 3、检验统计量 在零假设为真的条件下，令X X的一个随机观测落入类i i的概率为p pi。定义E Ei为H H0为真时观测值落入类i i的期望观测数，即E Ei= =piN，i=1,2,k.给出如下卡方检验统计量：第三节 卡方(2)拟合优度检验分布检验和拟合优度分布检验和拟合优度检验检验4 4、零分布：由于的精确分布难以求得，所以我们用自由度为k-1的卡方分布来近似。5
14、、假设组：H0：pi=p（i=1,2,k.）H1：pip（对某个i.）若(自由度为k-1的卡方分布的1-a分位数），则拒绝H0，p-值近似等于p(X2(c-1)Q)，这个概率可由附表10获得。第三节 卡方(2)拟合优度检验分布检验和拟合优度分布检验和拟合优度检验检验三、2检验的具体步骤v1 1. .数据分组. .根据样本观测值的范围划分为组；v2.2.求落在各组的频数i和频率Yi/n。v3.3.求理论概率i。当0成立时,出现在(bi-1,bi)内的概率iv4.计算检验统计量2。v5. .求出拒绝域.根据给定的显著性水平和自由度k-r-1查2分布表（附表10），可得临界值C，统计量2的拒绝域为2
15、C。v6. .作出判断.若2C则拒绝0，否则接受0。分布检验和拟合优度分布检验和拟合优度检验检验实例 单样本拟合检验v一、检验某固定比率的假设v例8.4：据标准规定，某批工业产品中不良品的比例为10%，则可检验如下假设：0：P=0.1；v1：P0.1。v为此，我们在产品批中抽出100个作为样本，发现不合格品数（Y1）为16，则合格品数Y2=100-16=84。v当0成立时，不合格品的期望数应为nP1=10个，相应地，合格品的期望数n（1-P1）=90。则：vk=，自由度为k-1=1，显著水平a0.05，查表10得临界值为3.841。由于n2a2，所以拒绝o假设。分布检验和拟合优度分布检验和拟合
16、优度检验检验二、检验某固定比率的假设v例8.5：检验随机变量在(,)区间是否为均匀分布。假设如下：v0：在（0,1）区间为均匀分布（假设分10类，pi=1/10）；v1：在（0,1）区间不是均匀分布(pip1/10)；v从未知总体中抽取50个样本。为了检验，我们可以将（0,1）区间分为10等份，即00.1，0.10.2，0.91.0。如果o为真，那么任何观察值落入类i的概率为1/10，任何小区间的期望观测数为(1/10)50=5。实例 单样本拟合检验分布检验和拟合优度分布检验和拟合优度检验检验例8.5v实际50个样本落入类i的观察如下：v区间 0- -0.10.1- -0.2 0.2-0.3
17、0.3- -0.40.4- -0.5 0.5-0.6v观测数645674v区间 0.6- -0.7 0.7-0.80.8-0.90.9- -1.0v观测数6534v检验统计量v查2分布表，自由度为，显著水平a=0.05时，查得a2=16.92，因2=2.80a2，所以接受o假设，即观察值取自均匀分布。分布检验和拟合优度分布检验和拟合优度检验检验三、检验多面体无偏性的假设v例8.6：以六面体的骰子为例。如果将一颗骰子抛掷120次，其结果如表8-4-4所示：v表8-48-4v根据题意，检验假设如下：vo：这颗骰子是无偏的(pi=1/6)；v1：这颗骰子是有偏的(pi1/6)；v如果零假设为真，各点
18、出现的期望次数vnP1=1/6120=20点数123456合计观测次数132816103221120实例 单样本拟合检验分布检验和拟合优度分布检验和拟合优度检验检验三、检验多面体无偏性的假设v例8.6：v检验统计量为：v查2分布表，自由度k-1=6-1=5，取a=0.05，查得a2=11.07因20.052，故应拒绝o假设，认为这颗骰子是有偏的。分布检验和拟合优度分布检验和拟合优度检验检验实例 总体分布拟合检验v一、正态分布拟合v例8.7：一家钟表厂把检验钟表的精确度作为质量控制的一部分。该厂将700只手表效准后使之走24小时，然后记下每只表走快或走慢的秒数( (数据见表5- -3)。这些数据
19、是否提供了充分的证据，说明观察值并非来自正态总体。v解：假设：0:样本数据来自正态总体分布；v1:样本数据并非来自正态总体分布。v表8- -3中，K=11=11，实际观察频数Oi已知，预期频数Ei则尚需确定。分布检验和拟合优度分布检验和拟合优度检验检验表53 700只手表时间误差的频数分布2424小时内走快或走慢的秒数表的数目O Oi0 09.999.993838101019.9919.995151202029.9929.996262303039.9939.997474404049.9949.998383505059.9959.999191606069.9969.998181707079.99
20、79.997272808089.9989.996161909099.9999.995252100100109.99109.993535合计700700分布检验和拟合优度分布检验和拟合优度检验检验一、正态分布拟合v1.1.预期频数Ei的计算v根据概率分布原理，我们可以通过求正态分布曲线下的面积来确定理论预期频数。为了计算正态分布曲线下的面积，利用公式Zo=(Xo-u)u)/将Xo标准化，求标准正态表上相应的面积（即频率）。v因为零假设中并没指定总体分布的均值(u)(u)和标准差()()。所以只有将样本均值=54.71和标准差S=27.61分别作为u和的估计值。v如在区间10-19.99内的预期频
21、数，可按如下步骤计算：分布检验和拟合优度分布检验和拟合优度检验检验v1. .预期频数Ei的计算v(1)(1)分别对x=10和x=20标准化：Z=(=(10-54.71)/27.61=-1.62和Z=(=(20-54.71) )/27.61=-1.26。其余类推。v(2)(2)查标准正态分布表，介于0和-1.62之间的面积（概率）为0.4474，介于0和-1.26之间的面积为0.3962，所以介于-1.62和-1.26之间的面积等于0.4474-0.3962=0.0512。其余类推。v(3)(3)于是落在10与20之间的预期频数为0.0512700=35.84。其余类推。一、正态分布拟合分布检验
22、和拟合优度分布检验和拟合优度检验检验v2.2.约束条件r r的确定v预期频数之和必等于700，即等于样本容量，这就构成一个约束；v又由于我们必须通过样本来估计u和，所以对数据还须增加两个约束。于是=3，自由度k-=11-3=8。v注意：如果u和在零假设中已被指定，那就不必再用样本数据来估计，这时=1。v检验统计量X2= (Oi-i) )2/i=20.3558，2的临界值为15.507。由于2的计算值大于临界值，所以否定零假设，样本数据并非来自正态分布。一、正态分布拟合分布检验和拟合优度分布检验和拟合优度检验检验小预期频数v在应用卡方检验时，有可能遇到预期频数很小的情形，这时将随机分布取作卡方的
23、近似分布并不完全正确。v对于什么样的预期频数才算小预期频数，学者们的意见并不一致。较保守的学者一般要求预期频数至少应大于等于5。而科库兰ochran( (1952，1954) )主张预期频数小于1的就算小预期频数，很多学者同意这种意见。本教材也采用了这一观点。v对小预期频数的处理一般采用的科克兰的法则。如将相邻类目的频数合并（前提是不破坏其分类意义），以达到所要求的最小频数。合并后的类数应相应地减小。分布检验和拟合优度分布检验和拟合优度检验检验二、二项分布拟合v例8.88.8:一个市场分析员想研究食品店的顾客对待信用卡付款方式的态度。研究员从100家超级市场各抽选了25名经常性顾客作为随机样本
24、，并对其中每一个进行访问以确定此人是否喜欢除信用卡付款方式以外的别的某种付款方式。调查结果列于表8-5：v解：这个分析员应先提出如下假设：v0：在这些容量为25的样本中，喜欢另外某种付款方式的顾客数服从二项分布；v1：不服从二项分布。（取a=0.05）分布检验和拟合优度分布检验和拟合优度检验检验表5-5 例5.5中的抽样结果喜欢另外某一种付款方式的顾客数商店数0 04 41 15 52 28 83 310104 414145 515156 612127 716168 810109 96 61010或更多0 0合计100100二、二项分布拟合分布检验和拟合优度分布检验和拟合优度检验检验v首先求得
25、P的估计值如下：P=4(0)+5(1)+8(2)+6(9)/2500=0.20。v由于二项分布的参数P没有指定，必须通过样本数据对它作出估计，因此要损失一个自由度。v1.1.预期频数iv通过计算函数f(x)=C25x(0.2)x(0.8)25-x(其中x为某一特定商店中喜欢另外某种付款方式的顾客数,x=0，1，2，25)或查的二项分布表，可以得到所需的相对预期频数。二、二项分布拟合分布检验和拟合优度分布检验和拟合优度检验检验v2.2.自由度的确定v合并后的类目数10，但由于预期频数之和必须与观察频数之和一致，这个自由度应减去，又由于必须通过样本数据来估计，自由度再减。于是真正的自由度应为10-
26、2=8。v检验统计量=28.1=28.1v显著性水平a=0.05和自由度8相对应的2的临界值为15.507。因为2的计算值大于的临界值，所以否定零假设，从而得出这些数据并非来自二项分布总体的结论（P0.05）。二、二项分布拟合分布检验和拟合优度分布检验和拟合优度检验检验三、泊松分布拟合v例8.98.9:旅馆管理人员对90天内房间预定和注销的格局进行了研究，其观察结果如表5-7。他想了解“每日注销的房间数是否服从泊松分布”？(取a=0.05)v解：由于泊松分布的参数并未给出，我们只有利用表中数据对它进行估计：v =0(9)+1(17)+8(2)/90=2.6 =0(9)+1(17)+8(2)/9
27、0=2.6v泊松分布的函数表达式为 (x=0,1,2,)分布检验和拟合优度分布检验和拟合优度检验检验表5-7 该旅游区旅馆注销的房间数注销的房间数观察到这一注销数的天数0 01 12 23 34 45 56 67 78 89 9或更多9 917172525151511117 72 22 22 20 0合计9090三、泊松分布拟合分布检验和拟合优度分布检验和拟合优度检验检验v利用表8-88-8的数据可算出：X2=(Oi-i) )2/i=6.674v本例中, ,自由度为6，因为经过小预期频数修正后的类目数为8，约束数为2：预期频数之和等于观察频数之和，样本数据来估计。v查2表，在显著水平a=0.0
28、5下，我们不能否定“数据来自泊松分布”这一零假设P0.10。三、泊松分布拟合分布检验和拟合优度分布检验和拟合优度检验检验五、KS检验与X拟合检验的比较vK-S检验与X检验均属拟合优度检验，但X拟合检验常用于对定类尺度测量的数据，而K-S检验还可用于对定序尺度测量的数据进行拟合检验。v另外对于X拟合检验来说，一般是要求一个大样本；在将所有的观察值分为n个组进行检验时，原则上还要求每组出现的观察值数目不能少于5，观察值数目少于5的组一般要加以合并，这样也就使样本信息受到了一些损失；而KS检验则无上述要求。第三，对于特别小的样本数目，X拟合检验不能应用，而KS检验不受限制。因此KS检验比X拟合检验显得更为有用。如果两种检验的样本容量相等时，KS检验对于Ho假设为不真可提供一个更高的拒绝率，因而与X拟合检验相比较具有更大的检出力。分布检验和拟合优度分布检验和拟合优度检验检验
