[数学]第1章方差分析及MATLAB实现2012.ppt

第1章 方差分析 例. 某水产研究所为了比较四种不同配合饲料对鱼的饲喂效果，选取了条件基本相同的鱼20尾，随机分成四组，投喂不同饲料，经一个月试验以后，各组鱼的增重结果列于下表。 表 饲喂不同饲料的鱼的增（单位：10g） 28.5 24.5 29.0 30.8 27.0 A4 25.8 24.9 27.3 23.6 22.1 A3 26.2 27.9 26.8 25.7 24.8 A2 35.9 28.4 31.8 27.9 31.9 A1 鱼的增重（xij） 饲料 四种不同饲料对鱼的增重效果是否显著 ？ 解：这是单因素均衡数据的方差分析，Matlab程序如下： A=[31.9 27.9 31.8 28.4 35.9 24.8 25.7 26.8 27.9 26.2 22.1 23.6 27.3 24.9 25.8 27.0 30.8 29.0 24.5 28.5]; %原始数据输入 B=A; % 将矩阵转置,Matlab中要求各列为不同水平 p=anova1(B) 运行后得到一表一图，表是方差分析表（重要）；图是各列数据的盒子图，离盒子图中心线较远的对应于较大的F值，较小的概率p. n-1 SST Total 总和 SS/(n-r) n-r SSE Error误差（组内） 0.0029 7.14 SS/(r-1) r-1 SSA Columns (因素A组间) P值 F统计量 MS均方差 df 自由度 SS 平方和 Source 方差来源 表中所列出的各项意义如下： 因为p=0.00290.01，故不同饲料对鱼的增重效果极为显著 .如果没有给出概率。 四种不同饲料对鱼的增重效果极为显著 ，那么哪一种最好呢？请看下图 此时，第一个图对应第一种饲料且离盒子图中心线较远，效果最突出。如果从原始数据中去掉第一种饲料的试验数据，得到的结果为各种饲料之间对鱼的增重效果不显著 . p=anova1(B(:,2:4)) 例.为比较同一类型的三种不同食谱的营养效果，将19支幼鼠随机分为三组，各采用三种食谱喂养. 12周后测得体重，三种食谱营养效果是否有显著差异？ 187 212 215 220 248 265 281 丙 185 197 201 231 乙 164 190 203 205 206 214 228 257 甲 体重增加量 食谱 解：这是单因素非均衡数据的方差分析 A=[164 190 203 205 206 214 228 257 185 197 201 231 187 212 215 220 248 265 281]; group=[ones(1,8),2*ones(1,4),3*ones(1,7)]; p=anova1(A, group) 方差分析表: 均值盒子图 由于概率p=0.1863比较大，故认为三种食料没有显著差异. (3) 多重比较的MATLAB实现 为了便于解决实际问题，我们给出多重比较的MATLAB命令。 c=multcompare(s) 其中输入s，由[p,c,s]=anova1(B);得到 输出C共有6列，每一行给出均值差的置信区间 例. 四个实验室试制同一型号纸张，为了比较光滑度每个实验室测量了8张纸，进行方差分析 42.8 41.2 37.8 37.2 35.4 34.8 34.8 34 A4 45 44 43 43 40 39 35 34 A3 45.8 43.3 42.9 41.8 41.4 39.7 39.3 39.2 A2 58 47.7 46 45.5 44.5 43.8 41.5 38.7 A1 纸张光滑度 实验室 解： A=[38.7,41.5,43.8,44.5,45.5,46,47.7,58 39.2,39.3,39.7,41.4,41.8,42.9,43.3,45.8 34,35,39,40,43,43,44,45 34,34.8,34.8,35.4,37.2,37.8,41.2,42.8]; %输入数据 B=A’; % MATLAB只对各列进行分析 [p,c,s]=anova1(b); % 方差分析 c=multcompare(s) % 多重比较 从方差分析表可知：四个实验室生产有差异，那么如何比较？软件输出c如下所示：1,2列表示比较的实验室号码，3，5 列分别为置信区间左右端点 ，第4 列是均值差的统计量观测值. 1.0000 2.0000 -1.4753 4.0375 9.5503 1.0000 3.0000 -0.1753 5.3375