非参数异均值方差分析1.docVIP

作者简介：李春红（1971- ），女，广西南宁，广西大学，副教授，主要从事应用统计的研究. 非参数的异均值方差分析 李春红，廖娟芬 黄绍军 （广西大学 数学与信息科学学院 广西 南宁 530004） 摘要 在方差齐性且分布未知的条件下，检验单因素等重复实验中的多个母体是否具有某种不同的均值，文章给出了样本中心化和秩统计量两种方法对其进行检验. 关键词 非参数； 异均值； 方差分析 中图分类号：O212.7 文献标识码：A 随着科学技术的发展，对于根据技术和实验得到的结果的正确性和可信性都需要验证.在文献[1] 提出了一种新的统计量来进行单因素等重复实验的非参数方差分析，文献[2] 提出了在正态分布下的异方差分析法，本文给出了单因素等重复实验的非参数异均值方差分析的两种方法. 1模型的引进 考虑单因素异方差分析问题，因素为，假定有个水平，目的是比较个水平的差异是否存在某种关系: (1) 令 (2) 是由仿真数字化设计等方法求得，要检验其是否满足式(1)，相当于检验假设 于是可得异均值方差分析的统计模型为 2模型的假设检验 2.1样本中心化 设抽取的因素的个水平的样本为.令 其中. 对进行排序，以表示在在合样本中的秩. 记.根据文献[1]引进统计量： 在原假设成立时及由引理得：，而当原假设不成立时，倾向于取较大的值，由此可得检验的否定域为： 其中是适当选定的常数. 常数的确定要根据的分布，而在原假设成立时的分布可参看文献,从而可实现本问题中的假设检验. 引进记号： R A C 引理1[1] C’(R- A)’(R- A)C.. 引理2[1] R为随机阵，当确定，成立时，R有种取法，且每种取法的概率均等，为，则有个值（可能有些相等，看成是有差别），且取每个值的概率为. 2.2秩统计量 先对进行排序，以表示在中的秩， 于是引进统计量： 原假设成立及由引理知：，则当原假设不成立时，倾向于取较大的值，因此检验否定域可取为： 是适当选定的常数. 为确定常数，要计算在原假设成立时的分布，原则上可以按照引理算出又由以下定理保证，??以比较方便的由计算机算出在原假设成立时的分布. 引进记号： R* A* C 则有 定理1 C’(R*- A*)’(R*- A*)C. 证明：假定确定的值，就有 . 定理2 R*为随机阵，当确定，成立时， R*有种取法，且每种取法的概率均等，为，则有个值（可能有些相等，看成是有差别），且取每个值的概率为. 证明：与引理2的证明类似，参考文献[1]. 通过计算机编程很容易算出的分布函数，从而实现本问题中的假设检验. 3模拟计算部分的分布函数及例子 下面给出用Matlab编程计算当的概率分布函数如下： 表1 00.2222220.6666670.8888891.5555620.05555560.4166670.1666670.1666670.1666670.0277778表2 00.1250.3750.50.8751.1250.06944440.2777780.2222220.1574070.1481480.0555556 1.51.62520.02777780.0370370 而为了例子的求解，只给出当时的的概率分布函数如下： 表3 00.00250.00340.00360.00810.00890.00990.01160.01350.01460.05560.13890.13890.13890.05560.05560.05560.05560.05560.0556 0.01780.01910.02040.02230.02250.02620.0270.03040.03280.02780.02780.02780.00930.02780.02780.02780.00930.0093 例子：为了获得每一箱的终重为1g的虾，有三种饲料配方见文献[6]，选取第种配方下每一箱虾的尾数为，每种配方作为一实验组，每实验组有三个平行，由于精心照料每一只虾最后都存活，下面给出30天后的数据如下： 表4 饲喂实验结果 饲料配方尾数平均终重数据（g）预计得到的箱重（g）1