理学方差分析建模.pptx

严谨笃实尚善立人统计与应用数学学院单因素方差分析1单因素方差分析模型设指标变量为Y,因素为A，它可以取几个不同的水平，记为A1,…,Aa,在水平Ai上作若干次试验，试验结果为yi1,…,yini,需要考虑因素A对指标的影响是否显著.将A在某个水平下的试验结果看成一个随机变量(总体)，因此考虑A取不同水平时的指标有无显著差异，就是检验几个总体的均值是否相等.假定在水平Ai下的总体服从正态分布，他们可以有不同的均值?i,但有相同的方差?2，需要检验假设：不全相等(3.1)

严谨笃实尚善立人统计与应用数学学院将每次试验结果yij分解为其中?表示总的均值，?i=?i-?是水平Ai对指标的效应,满足?ij是相互独立的随机误差，服从正态分布N(0,?2).建立方差分析模型如下：(3.2)

严谨笃实尚善立人统计与应用数学学院2.单因素方差分析的Matlab实现[p,c,s]=anova1(X,group)输入：X是一个向量，从第一个总体的样本到第r个总体的样本依次排列，group是与X有相同长度的向量，表示X中的元素是如何分组的.group中某元素等于i,表示X中这个位置的数据来自第i个总体.因此group中分量必须取正整数，从1直到r.[p,c,s]=anova1(X)%比较X各列数据的均值是否相等输出p是零假设成立时概率，对给定的?，若p?，则有显著差异；c是方差分析表，s用于多重比较的输入.输入X各列的元素相同，即各总体的样本大小相等，称为均衡数据的方差分析，不均衡时用下面的命令

严谨笃实尚善立人统计与应用数学学院例3.1.某水产研究所为了比较四种不同配合饲料对鱼的饲喂效果，选取了条件基本相同的鱼20尾，随机分成四组，投喂不同饲料，经一个月试验以后，各组鱼的增重结果列于下表。表3-1饲喂不同饲料的鱼的增（单位：10g）饲料鱼的增重（xij）A131.927.931.828.435.9A224.825.726.827.926.2A322.123.627.324.925.8A427.030.829.024.528.5四种不同饲料对鱼的增重效果是否显著？

严谨笃实尚善立人统计与应用数学学院解：这是单因素均衡数据的方差分析，Matlab程序如下：A=[31.9 27.9 31.8 28.4 35.924.8 25.7 26.8 27.9 26.222.1 23.6 27.3 24.9 25.827.0 30.8 29.0 24.5 28.5];%原始数据输入B=A;%将矩阵转置,Matlab中要求各列为不同水平[p,c,s]=anova1(B)运行后得到一表一图，表是方差分析表（重要）；图是各列数据的盒子图，离盒子图中心线较远的对应于较大的F值，较小的概率p.

严谨笃实尚善立人统计与应用数学学院Source方差来源SS平方和df自由度MS均方差F统计量P值Columns(因素A组间)SSAr-1SS/(r-1)7.140.0029Error误差（组内）SSEn-rSS/(n-r)Total总和SSTn-1表中所列出的各项意义如下：因为p=0.00290.01，故不同饲料对鱼的增重效果极为显著.

严谨笃实尚善立人统计与应用数学学院四种不同饲料对鱼的增重效果极为显著，那么哪一种最好呢？请看下图此时，第一个图对应第一种饲料且离盒子图中心线较远，效果最突出。如果从原始数据中去掉第一种饲料的试验数据，得到的结果为各种饲料之间对鱼的增重效果不显著.

严谨笃实尚善立人统计与应用数学学院p=anova1(B(:,2:4))

严谨笃实尚善立人统计与应用数学学院例3.2为比较同一类型的三种不同食谱的营养效果，将19支幼鼠随机分为三组，各采用三种食谱喂养.12周后测得体重，三种食谱营养效果是否有显著差异？食谱体重增加量甲164190203205206214228257乙185197201231丙187212215220248265281解：这是单因素非均衡数据的方差分析A=[164190203205206214228257185197201231187212215220248265281];group=[ones(1,8),2*ones(1,4),3*ones(1,7)];p=anova1(A,group)表3-2饲喂不同食谱鼠的增重

严谨笃实尚善立人统计与应用数学学院方差分析表均值盒子图由于概率p=0.1863比较大，故认为三种食料没有显著差异.

严谨笃实尚善立人统计与应用