数理统计方法的Matlab实现(7.0版).ppt

一、统计量的计算 例 (胃癌的鉴别) 两类Fisher判别示意图 * 数理统计的Matlab实现 * 数理统计的Matlab实现 统计量的计算 1 多元线性回归分析 3 方差分析 2 判别分析 4 主成分分析 6 聚类分析 5 因子分析 7 以向量为计算单位的统计量 sum , mean , median , var, std, max, min x=[ 909, 1091, 967,1232,1096,1164, 1086,1071,1572, 950, 808, 971, 650,789, 1120,1081,825 ,775,1224, 950, 999,1130, 914,1203,1044, 866, 1320,1336,992 ,1025,871,738,865]; mean(x ), var (x) , std(x) 注：若参数为矩阵， 默认计算各列向量的统计量，返回行向量； 若要计算各行向量的统计量，要加一个维数参数2，返回列向量。 C=[1,3,3;4,8,6;7,7,9;6,5,9] mean(C) %计算样本均值向量 mean(C,2) 计算各行向量的最值较特殊，用max(C,[ ],2) , min(C,[ ],2) 2. 计算样本的协方差矩阵和相关矩阵 cov(C ),corrcoef(C ) 3. 常用分布的分位点 nameinv (x, 参数表列) 其中函数名name的含义同前。 求分位点 。 norminv(0.025, 0,1) tinv(0.025, 10) chi2inv(0.025 ,10) finv(0.05 ,6,10) 二、方差分析 1.单因素方差分析 anova1(X) 其中，X为需要分析的数据，每一列对应于一个组（因子水平）的测试数据。 例1：某农科所为比较三种不同肥料对水稻产量的影响，进行如下试验。选一块肥沃程度较均匀的土地，将其分割成12块，随机地分成3组，所得产量如下，请评价三种不同肥料对水稻产量由无显著性差异。 78，65，80，69 A3 62，68，50，70 A2 94，91，75，85 A1 产量 水平 若p=P(F≥F值)不大于给定的显著性水平(0.05),则拒绝原假设. 相关知识点: A=[94 91 75 85; 62 68 50 70; 78 65 80 69] anova1(A) 2.双因素方差分析 anova2() 例2：某农科所为提高水稻产量，对不同水稻品种（因素A)和施肥量(因素B)进行组合试验，所得结果如下，试分析各水稻品种与各施肥量对水稻产量是否有显著影响。 72 80 88 82 B3 93 84 89 85 B2 95 90 97 87 B1 A4 A3 A2 A1 因素A 因素B AB=[87 97 90 95;85 89 84 93;82 88 80 72] anova2(AB) 2.有交互作用的方差分析 anova2(X,reps) 例3：某农科所为提高水稻产量，对不同水稻品种（因素A)和施肥量(因素B)进行组合试验,每种组合各3次，所得结果如下，试分析各水稻品种与各施肥量及它们的交互作用对水稻产量是否有显著影响。 AB2=[95 93 97 96 87 89 90 91;85 86 89 90 84 87 92 90;72 76 75 73 85 86 88 89 ] anova2(AB2,2) 88 89 85 86 75 73 72 76 B3 92 90 84 87 89 90 85 86 B2 90 91 87 89 97 96 95 93 B1 A4 A3 A2 A1 因素A 因素B 三、多元线性回归分析 例1：研究货运总量y（万吨）与工业总产值x1(亿元)、农业总产值x2(亿元)、居民非商品支出x3(亿元)的关系。数据如下： 250 275 160 275 220 240 265 210 260 160 y 3.0 42 65 10 3.2 44 70 9 2.0 36 66 8 4.0 42 78 7 1.5 45 68 6 1.2 38 72 5 3.0 42 74 4 2.0 40 65 3 2.4 40 75 2 1.0 35 70 1 x3 x2 x1 编号 相关知识点: A=load(D:\\货运总量.txt) [n,p]=size(A); X=[ones(n,1),A(:,1:p-1)] y=A(:,p); [