第9章--概率论与数理统计问题的计算机求解.pptVIP

星期日, 2013-5- 19, 07:56:11 例 9.29 用统计学工具箱生成一组G分布数据，用现成函数验证其是否为正态分布数据，显然这些数据不是正态分布的，所以假设检验结果应该是1 MATLAB求解语句： 星期日, 2013-5- 19, 07:56:11 9.4.2.3 其他分布的Kolmogorov-Smirnov检验 Kolmogorov-Smirnov检验是检验任意已知分布函数的一种有效的假设检验算法 函数调用格式： 其中，cdffun为两列的均值，第1列为自变量，第2列应该为要检验的分布函数在自变量处的值 星期日, 2013-5- 19, 07:56:11 例 9.30 生成一组G分布数据，对生成的随机数进行假设检验：该随机数满足G分布 生成G分布的数据 假设进行检验 星期日, 2013-5- 19, 07:56:11 9.5 方差分析与主成分分析 方差分析 主成分分析方法 星期日, 2013-5- 19, 07:56:11 9.5.1 方差分析 单因子方差分析 双因子方差分析 多因子方差分析 星期日, 2013-5- 19, 07:56:11 9.5.1.1 单因子方差分析 单因子方差分析就是指对一些观察来说，只有一个外界因素可能对观测的现象产生影响 求解单因子方差分析的函数调用格式 其中，X 为需要分析的数据 星期日, 2013-5- 19, 07:56:11 单因子方差分析表 星期日, 2013-5- 19, 07:56:11 接上页 星期日, 2013-5- 19, 07:56:11 例 9.31 有5种药物比较疗效，将30个病人随机地分成5组，每组使用同一种药物，并记录病人治疗时间如下表，试评价疗效 星期日, 2013-5- 19, 07:56:11 MATLAB求解语句： 星期日, 2013-5- 19, 07:56:11 9.5.1.2 双因子方差分析 如果有两种因子可能影响到某现象的统计规律，则应该引入双因子方差分析的概念 观测量 y 可以表示为一个三维数组 ，表示第 1 个因子取第 i 个水平，第 2 个因子取第 j 个水平时，组内第k个对象的观测指标。 星期日, 2013-5- 19, 07:56:11 三个假设: 为第一因子单独作用的效应 为第二因子单独作用的效应 为两个因子同时作用的效应 星期日, 2013-5- 19, 07:56:11 若 则拒绝假设 若 则拒绝假设 若 则拒绝假设 求解双因子方差分析问题的函数调用格式 3个概率的定义及意义 星期日, 2013-5- 19, 07:56:11 双因素方差表 星期日, 2013-5- 19, 07:56:11 接上页 星期日, 2013-5- 19, 07:56:11 其中 概率定义 MATLAB 求解 星期日, 2013-5- 19, 07:56:11 例 9.31 设为比较3种松树在4个不同地区的生长情况有无差别，在每个地区对每种松树随机地选择5株，测量它们的胸径，得出的数据在下文的表中给出，试对它们进行双因子方差分析 星期日, 2013-5- 19, 07:56:11 松树数据 星期日, 2013-5- 19, 07:56:11 MATLAB 求解 星期日, 2013-5- 19, 07:56:11 计算均值： 星期日, 2013-5- 19, 07:56:11 9.5.1.3 多因子方差分析 MATLAB语言的统计学工具箱还可以进行三因子甚至多因子的方差分析，可以采用manova1()函数进行多因子方差分析 星期日, 2013-5- 19, 07:56:11 9.5.2 主成分分析方法 假设某一事件发生可能受 这 N 个因素影响，而实测数据共有 M 组，这样可以假设这些数据由一个 N×M 矩阵 X 表示。记该矩阵的每一列的均值为 星期日, 2013-5- 19, 07:56:11 主成分分析方法的一般步骤 调用corr()函数，由矩阵 X 可以建立起协方差矩阵 R 星期日, 2013-5- 19, 07:56:11 由 R 矩阵可以分别得出特征向量 和对应的排序特征值 ，特征向量矩阵的每一列也都进行了相应的归一化