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

非参数的异均值方差分析 李春红，廖娟芬 黄绍军 （广西大学 数学与信息科学学院 广西 南宁 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 0 0.222222 0.666667 0.888889 1.55556 2 0.0555556 0.416667 0.166667 0.166667 0.166667 0.0277778 表2 0 0.125 0.375 0.5 0.875 1.125 0.0694444 0.277778 0.222222 0.157407 0.148148 0.0555556 1.5 1.625 2 0.0277778 0.037037 0而为了例子的求解，只给出当时的的概率分布函数如下： 表3 0 0.0025 0.0034 0.0036 0.0081 0.0089 0.0099 0.0116 0.0135 0.0146 0.0556 0.1389 0.1389 0.1389 0.0556 0.0556 0.0556 0.0556 0.0556 0.0556 0.0178 0.0191 0.0204 0.0223 0.0225 0.0262 0.027 0.0304 0.0328 0.0278 0.0278 0.0278 0.0093 0.0278 0.0278 0.0278 0.0093 0.0093 例子：为了获得每一箱的终重为1g的虾，有三种饲料配方见文献[6]，选取第种配方下每一箱虾的尾数为，每种配方作为一实验组，每实验组有三个平行，由于精心照料每一只虾最后都存活，下面给出30天后的数据如下： 表4 饲喂实验结果 饲料配方 尾数 平均终重数据（g） 预计得到的箱重（g） 1 10 0.091 0.097 0.088 1 2 9 0