1、p=anova1（X）含义：比较样本 mn的矩阵X中两列或多列数据的均值。其中，每一列表示一个具有m个相互独立测量的独立样本。返回：它返回X中所有样本取自同一总体（或者取自均值相等的不同总体）的零假设成立的概率p。解释：若p值接近0（接近程度有解释这自己设定），则认为零假设可疑并认为至少有一个样本均值与其它样本均值存在显著差异。Matlab程序：Score=75 62 71 58 73;81 85 68 92 90;73 79 60 75 81;P=anova1（Score）输出结果：方差分析表和箱形图ANOVA TableSource SS df MS F ProbFColumns 604.
2、9333 2 302.4667 4.2561 0.040088Error 852.8 12 71.0667 Total 1457.7333 14 由于p值小于0.05，拒绝零假设，认为3种教学方法存在显著差异。例2（双因素方差分析）为了考察4种不同燃料与3种不同型号的推进器对火箭射程（单位：海里）的影响，做了12次试验，得数据如表2所示。表2 燃料-推进器-射程数据表推进器1推进器2推进器3燃料158.256.265.3燃料249.154.151.6燃料360.170.939.2燃料475.848.7在Matlab中利用函数 anova2函数进行双因素方差分析。p=anova2（X,reps）
3、比较样本X中两列或两列以上和两行或两行以上数据的均值。不同列的数据代表因素A的变化，不同行的数据代表因素B的变化。若在每个行-列匹配点上有一个以上的观测量，则参数reps指示每个单元中观测量的个数。当 reps=1（默认值）时，anova2将两个p值返回到向量p中。 H0A：因素A的所有样本（X中的所有列样本）取自相同的总体； H0B：因素B的所有样本（X中的所有行样本）取自相同的总体。 当reps1时，anova2还返回第三个p值： H0AB：因素A与因素B没有交互效应。如果任意一个p值接近于0，则认为相关的零假设不成立。disp1=58.2 56.2 65.3;49.1 54.1 51.6
4、;60.1 70.9 39.2;75.8 58.2 48.7;p=anova2（disp1,1）方差分析表Columns 157.59 3 52.53 0.43059 0.73875Rows 223.8467 2 111.9233 0.91743 0.44912Error 731.98 6 12 1.9967 Total 1113.4167 11 由于燃料和推进器对应的p值均大于0.05，所以可以接受零假设H0A和H0B，认为燃料和推进器对火箭的射程没有显著影响。例3（双因素方差分析）设火箭的射程在其它条件基本相同时与燃料种类及推进器型号有关。现在考虑4种不同的燃料及3种不同型号的推进器，对于
5、每种搭配个发射了火箭两次，得数据见表3。问各自变量和自变量的交互效应是否对火箭的射程有显著影响？表3 燃料-推进器-射程数据表52.641.260.842.850.548.458.373.240.771.551.041.4disp2=58.2 52.6 49.1 42.8 60.1 58.3 75.8 71.5;56.2 41.2 54.1 50.5 70.9 73.2 58.2 51.0;65.3 60.8 51.6 48.4 39.2 40.7 48.7 41.4;p=anova2（disp2,2）Columns 370.9808 2 185.4904 9.3939 0.003506Row
6、s 261.675 3 87.225 4.4174 0.025969Interaction 1768.6925 6 294.7821 14.9288 6.1511e-005Error 236.95 12 19.7458 Total 2638.2983 23 显著。方差分析上机练习为研究广告的效果，考察4种广告方式：当地报纸（paper）、当地广播（radio）、店内销售员（people）和店内展示（display）的效果。共设有144个销售点，每种广告随机抽取36个销售点记录销售额，分布在6个地区的144个销售点的销售情况生成的数据集ADS见下表。数据集ADS中有3个变量：AD表示广告的类型、
7、AREA表示地区、SALES表示销售额（单位：千元）。请完成以下练习：（1） 概括下列数据：用箱形图、条形图直观地呈现四种广告方式下销售量的分布情况；计算四种广告方式下销售量的均值、方差、标准差、最大和最小值；（2） 进行单因素方差分析：检验四种广告方式下销售量数据是否服从正态分布，方差是否相等；检验四种广告方式下的销售量是否有显著差异（）；若四种广告方式下的销售量有显著差异，指出哪些类型的广告效果有显著的不同？（3） 在设计广告效果的试验时，虽然地区差异对销售量的影响并不是我们感兴趣的，但希望排除这一因素的影响。数据集ADS记录了各个销售点所在的地区AREA。试用双因素方差分析方法分析销售数
8、据，并指出广告方式和地区对销售量是否有显著影响（）？广告方式（AD）与地区（AREA）之间有无交互效应？表 ADS数据集中的数据广告方式（变量:AD）销售额（单位：千元）（变量SALES）地区1地区2地区3地区4地区5地区6当地报纸（paper）75 57 7668 75 8377 75 7266 66 7676 81 6370 86 6294 54 7088 56 8687 65 6584 77 7879 62 7580 62 70当地广播（radio）69 51 10054 78 7990 77 6083 74 6933 79 7368 75 65100 61 6870 53 7368 6
9、3 8379 66 6576 73 7481 57 65店内销售员（people）63 67 8558 82 7880 87 6287 70 7770 75 4068 61 5564 40 6776 70 7751 61 7542 71 6564 50 6278 37 83店内展示（display）52 61 6141 44 8676 57 5275 75 6333 69 6052 61 4361 66 4169 43 5165 58 5060 52 5544 45 5852 45 60参考答案（1）箱形图：boxplot（ads） 结果：有异常值。（其它：略）（2）正态性检验 Paper：
10、Hist（X1,6）频数直方图分布的正态性检验： normplot（X1） 均服从正态分布。单因素方差分析Columns 5866.0833 3 1955.3611 13.4831 8.8495e-008Error 20303.2222 140 145.023 Total 26169.3056 143 P=8.8495e-008Columns 1444.2222 5 288.8444 1.9582 0.089763Rows 5866.0833 3 1955.3611 13.2559 1.5637e-007Interaction 1158 15 77.2 0.52336 0.92341Error
11、 17701 120 147.5083 从以上分析结果可知：0.05P1=0.0897630.1，地区对检验水平有一定影响，但不显著。P2=1.5637e-0070.010.1，地区和广告方式对销售量无交互效应。二 Matlab作回归分析回归分析的相关数学理论可以参见概率论与数理统计教程，下面仅以示例说明如何利用matlab处理回归分析。1.一元线性回归分析【例1】为了了解百货商店销售额x与流通费率（反映商业活动的一个质量指标，指每元商品流转额所分摊的流通费用）y之间的关系，收集了九个商店的有关数据，见下表1.试建立流通费率y与销售额x的回归方程。表1 销售额与流通费率数据样本点销售额x（万元）流通费率y11.57.024.54.837.53.6410.53.1513.52.7616.52.5719.52.4822.52.3925.52.2【分析】：首先绘制散点图以直观地选择拟合曲线，这项工作可结合相关专业领域的知识和经验进行，有时可能需要多种尝试。选定目标函数后进行线性化变换，针对变换后的线性目标函数进行回归建模与评价，然后还原为非线性回归方程。【Matlab数据处理】：【Step1】：绘制散点图以直观地选择拟合曲线x=1.5 4.5 7.5 10.5 13.5 16.5 19.5 22.5 25.5;y=7.0 4.8 3.6 3.1 2.7
