多元正态分布及其抽样分布.pptx

《多元正态分布及其抽样分布.pptx》由会员分享，可在线阅读，更多相关《多元正态分布及其抽样分布.pptx（87页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第二章第二章第二章第二章多元正态分布及其抽样分布多元正态分布及其抽样分布多元正态分布及其抽样分布多元正态分布及其抽样分布内 容第一节 多元正态分布的定义第二节 多元正态的性质第三节 多元正态参数的极大似然估计第四节 多元正态的样本分布第一节第一节 多元正态分布的定义多元正态分布的定义一、标准多元正态分布一、标准多元正态分布 则设随机向量其分量独立同分布于密度函数为其中的均值为 协方差矩阵为 二、一般的正态分布二、一般的正态分布 设随机向量 ,若其的密度函数为其中 的均值为 协方差为称 服从均值为E(X)，协方差为 的正态分布。三、一般的三、一般的p p维正态和维正态和p p维标准正态的关系维标

2、准正态的关系 设 ，其中 是一个 阶非退化矩阵，服从 维标准正态分布，则 服从p维正态分布，且均值向量为 x x的协方差矩阵为的协方差矩阵为其密度函数为 若 ，则1存在，是非退化 元正态分布；若 ，则 不存在，是退化 元正态分布，不存在密度函数。值得注意值得注意 设随机向量 ，是常数向量，是一个 的常数矩阵，则 服从正态分布，记为 ，其中 例：设随机向量 ，则 的分布是退化的三元正态分布。第二节第二节 多元正态分布的性质多元正态分布的性质二、x x是一个服从p维正态分布，当且仅当它的任何线性函数 服从一元正态分布 。一、多元正态分布的特征函数一、多元正态分布的特征函数 三、X X服从 维正态分

3、布，则 ，其中 为 常数矩阵，为 维的常数向量，则 四、设 ，则 的任何子向量也服从多元正态分布，其均值为 的相应子向量，协方差为 的相应子矩阵。五、设 ,，相互独立，且，则对任意 个常数 ，有 六、，则 分布。七、将 作如下的分块：子 向量相互独立，当且仅当 。证：必要性 八、设 ，其中 是 阶矩阵，是 阶矩阵，则 与 相互独立，当且仅当 。九、设 ，其中 是 阶矩阵，是 阶矩阵，则 与 相互独立，当且仅当 。同上可证。十、将 作如下的分块：则 与 相互独立，与相互独立 。证：则给定 时 的条件分布为 ，其中 十一、将 作如下的分块：为 给定的条件下 数学期望。十二、偏相关系数十二、偏相关系

4、数 矩阵 称为条件协方差矩阵，它的元素用表示。是当 给定的条件下，与 （）的偏相关系数，定义为 它度量了在值 给定的条件下，与（）相关性的强弱。例 设XN6(,)，其协方差矩阵为，计算偏相关系数。求x7给定的条件下，x1，x6的偏协方差矩阵3 3 实例分析及实例分析及SAS/CORRSAS/CORR 例1 今对31人进行人体测试，考察的7个指标是：x1：年龄 x2：体重 x3：肺活量 x4：1.5英里跑所需时间 x5：休息时的脉搏 x6：跑步时的脉搏 x7：跑步时记录的最大的脉搏 对这些指标进行一些相关分析。SAS的程序的程序data a;input x1-x7;cards;44 89.47

5、44.609 11.37 62 178 18240 75.07 45.313 10.07 62 185 18538 89.02 49.874 9.22 55 178 18047 48 61.24 47.920 11.50 52 170 17652 82.78 47.467 10.50 53 170 172;proc corr nosimpl cov;var x1;with x7;partial x3;run;proc corr nosimpl cov;分析相关系数nosimpl是要求不打印描述性统计量。var x1;指定分析相关系数的变量。with x7;with指定变量与var指定的变量之间

6、的相关系数。partial x3;当指定的变量给定时，计算偏相关系数。x1x2x3x4x5x6x7x11.00000-0.23354-0.304590.18875-0.16410-0.33787-0.43292P值值0.20610.09570.30920.37770.06300.0150 x2-0.233541.00000-0.162750.143510.043970.181520.24938P值值0.20610.38170.44120.81430.32840.1761x3-0.30459-0.162751.00000-.86219-0.39936-0.39797-0.23674P值值0.09

7、570.3817.00010.02600.02660.1997x40.188750.14351-0.862191.000000.450380.313650.22610P值值0.30920.4412.00010.01100.08580.2213x5-0.164100.04397-0.399360.450381.000000.352460.30512P值值0.37770.81430.02600.01100.05180.0951x6-0.337870.18152-0.397970.313650.352461.000000.92975P值值0.06300.32840.02660.08580.0518.

8、0001x7-0.432920.24938-0.236740.226100.305120.929751.00000P值值0.01500.17610.19970.22130.0951|r|underH0:PartialRho=0 x1x7-0.545730.0018第三节 极大似然估计及其性质则总体的密度函数为X1，X2，Xn是从总体中抽取的一个简单随机样本，满足X1，X2，Xn相互独立，且同正态分布称X X为样本数据矩阵。一、样本的联合密度函数一、样本的联合密度函数为样本联合密度函数。所以，似然函数还可以表示为：二、和 的极大似然估计的极大似然估计 所谓和的极大似然估计，是寻找 和 满足条件

9、令 可以证明和的极大似然估计为三、相关系数的极大似然估计 （一）极大似然估计的不变性质 设 是的极大似然估计是 ，而且变换f()是一一对应的，则f()的极大似然估计就是（二）简单相关系数的极大似然估计 其中Sij是样本协方差矩阵S中相应位置上的元素（三）偏相关系数的极大似然估计则偏相关系数的极大似然估计其中，。（四）复相关系数的极大似然估计将x和S作如下的分块 的线性函数为 定义定义(复相关系数复相关系数)一个变量y与一组变量X1,X2,XK的负相关系数是以y为被解释变量，X1,X2,XK为自变量的回归方程的可决系数。为自变量的回归方程的可决系数。为了研究四川经济增长的影响因素，欲建立四川省经

10、济增长模型。主要经济指标采用国内生产总值增长率（Y），投资指标资本形成总额增长率（X1），人口指标用自然增长率（X2），就业指标失业率（X3）和消费指标居民消费水平增长率（X4）。分析指标之间的关系。data a;input y x1-x4;cards;数据行;proc corr nosimpl noprob cov;run;proc iml;sigma22=76.58605619 2.59407381-3.45807619 49.03157071,2.59407381 5.14447619 -0.78252381 4.24046429,-3.45807619-0.78252381 3.637

11、47619-2.32063571,49.03157071 4.24046429-2.32063571 53.90793143;sigma12=57.79053524 4.91975476-2.98844524 52.41117214;fcorr=sigma12*inv(sigma22)*t(sigma12)/54.8989690;print fcorr;proc reg;model y=x1-x4;run;Analysis of Variance Sum of MeanSource DF Squares Square F Value Pr FModel 4 1089.28592 272.321

12、48 501.20 FWilksLambda0.545616206.874330.0004PillaisTrace0.454383806.874330.0004Hotelling-LawleyTrace0.832790156.874330.0004RoysGreatestRoot0.832790156.874330.0004直接检验两个总体的均值向量是否相等。DependentVariable:x1（对（对X1进行的检验）进行的检验）SumofSourceDFSquaresMeanSquareFValuePrFModel10.874667910.8746679116.900.0002Error

13、361.863008400.05175023CorrectedTotal372.73767632X1在类间有显著性差异。DependentVariable:x2（对（对X2进行的检验）进行的检验）SumofSourceDFSquaresMeanSquareFValuePrFModel10.083120770.083120771.950.1710Error361.533700280.04260279CorrectedTotal371.61682105X2在类间没有显著性差异。DependentVariable:x3（对（对X3进行的检验）进行的检验）SumofSourceDFSquaresMea

14、nSquareFValuePrFModel116.4695844316.4695844321.45FModel10.001126940.001126940.030.8643Error361.369780950.03804947 CorrectedTotal371.37090789X4在类间没有显著性差异。第四节第四节 抽样分布抽样分布 一、维希特（一、维希特（WishartWishart）1 1、定义、定义随机矩阵的分布 矩阵中的每一个元素均为随机变量，则矩阵X的分布是其列向量拉长，组成一个长向量 特别当 是 阶对称阵，则 的分布为的下三角部分组成的长向量 在一元正态随机变量中，我们曾经讨论了

15、 分布，在多元正态随机变量也有类似的样本分布。维希特分布（Wishart）相当于一元统计中的 分布。定义定义 维希特（维希特（WishartWishart）分布的统计量）分布的统计量 设 个随机向量 独立同分布于，则随机矩阵服从自由度为的非中心维斯特分布，记为。二、维斯特(Wishart)分布有如下的性质：（1）若A1和A2独立，其分布分别和，则的分布为，即维斯特分布有可加性。（2），C为mp阶的矩阵，则的分布为分布。三、三、抽样分布抽样分布 定理1：设X1,X2,Xn是来自多元正态总体Np(,)的简单随机样本，有则有证明：独立故且相互独立。独立当 ，时，由卡方分布的定义可知可见维希特分布是由

16、卡方分布在多元下的推广。服从自由度为 的卡方分布。定理定理2 2 设 独立同正态分布，则统计量 证：由于样本均值 相互独立的标准正态分布的平方和为自由度为 的卡方分布。在一元正态的情形下，我们有样本的统计量当总体的方差未知时，我们必须用样本的方差来代替总体的方差，则那么在多元正态的情形下，是否有相同的问题呢？回答时肯定的。定义：称T2服从参数为P和n的非中心霍特林（Hotelling)分布，当。当时，服从自由度为n的中心霍特林分布，记为。定理定理：定定理理：设 是来自多元正态总体 的简单随机样本，有 定理：设 是来自多元正态总体 的简单随机样本，设 是来自多元正态总体 的简单随机样本，（1）Wilks分布 定义：设和，且相互独立，和，则称服从Wilks分布，记。可以证明，当和时，Wilks分布可以用分布近似。四、基于四、基于维斯特维斯特(Wishart)分布的统计量分布的统计量在一元方差分析中，常常遇到基于独立的分布随机变量比值的统计量。在多元统计分析中，起到相同作用的是统计量和分布。2、统计量和分布设k个总体，它们服从。分别抽出如下的样本：W=E+B当当K个总体的均值相等时个总体的均值相等时,服从服从Wilks分布。分布。