方差分类模型的假设检验.pdfVIP

维普资讯 第 35卷 第4期 河南师范大学学报 (自然科学版) Z．35 N o．4 2007年 11月 JournalofHenanNormalUniversity (NaturalScience) NOV．2007 文章编号 ：1O0O一2367(2007)04—0171—03 方差分类模型的假设检验 王石青，史慧娟 (华北水利水电学院 数学与信息科学学院，郑州 450008) 摘 要：主要讨论了对于广泛应用的方差分类模型，如何判断随机区组效应的存在性及是否有显著性差异．通 过对区组 内分析获得了固定效应是否存在的统计检验方法及两个随机效应是否有显著性差异的检验方法． 关键词：方差分类模型；方差分量；区组内分析；假设检验 中图分类号 ：O212．1 文献标识码 ：A 近 2O年来，线性混合效应模型在生物、医学、经济、金融、环境科学、抽样调查及工程技术领域得到愈来 愈广泛的应用 ，因为模型观测向量的协方差阵含有未知参数 ，于是关于它的参数估计、假设检验和置信区间 往往变得 比较复杂．本文考虑如下 的方差分类模型E ： 一 鄯 +U1a+U2 +￡， (1) 这里 为n×1的观测向量，X和u 分别是n× ，n×t的已知设计矩阵，为 ×1的固定效应 ，8为t× 1的随机效应 向量 ，其分布分别为N(0， ，I)，￡是n×1误差向量，其分布为N(0，2)且它们彼此独立 ， 和 是方差分量 ，且 (X)一r， (u )一m， (u )一z．于是观测 向量 的协方差阵为 cov(y)一 ；uu + ；U2U2+ I． 这样的模型，在生物、医药、经济等领域含一个个体随机效应 的区组设计、分群 (cluster)抽样 、纵 向数据 (Longitudinaldata)以及经济领域中Panal数据的分析 中，可经常遇到．文献[2]仅对含有两个方差分量的线 性混合模型，给出了固定效应的精确置信区间及随机效应的方差分量 的一致最优无偏检验 (UMPIUT)及功 效函数．对带有随机区组效应 的响应曲面模型，文献[3]获得 了区组效应是否存在的统计检验方法 ，但这种模 型也只含有一个随机效应．本文对方差分类模型，应用区组分析对固定效应 口中的元素是否存在 的假设检验 问题给出了统计检验方法 ，并且给出了两个随机效应是否有显著差异的统计检验方法． 本文将采用下列记号 ：A ，tr(A)，M(A)， (A)，A 和A 分别表示给定矩阵A的转置矩阵，迹 ，列子空 间，秩 ，广义逆和 Moore—Penrose广义逆 ，记 P 一A(AA)一A ，表示到列子空间M (A)上的正交投影阵 ， Na— l— P、，V 一 U U ． 对 固定效应的假设检验 首先给出两个引理 ： 引理 1[ 设 X～ N (，I)，A对称，则 XAX ～ 2， A幂等， (A)一 r． 引理 2c 设 X ～ N (，I)，A ，A 对称 ，X，A X和 X，A X 都服从 。分布 ，则它们相互独立 A 1A2— 0． 区组内分析是将区组效应看作 固定的，涉及到对模型(1)的 口中的元素的检验是严格 的，以下应用区组 内分析对卢中的元素是否存在进行讨论．记S：[X iu iu ]，T 一 [ ；u ；u ]，其中 是相应卢㈤ 的X 的一部分 ，卢㈤是从卢中剔除卢 之后剩余 的 一 1维参数 向量 ．于是总 的剩余平方和为 ：Q 一 VI一 收稿 日期 ：2006— 11—29 基金项 目：河南省 自然科学基金 (0611052600) 作者简介 ：王石青(1955--)，男，河北武安人 ，华北水利水电学院教授 ，主要从事数理统计方面 的研究 维普资讯 河南