matlab方差分析可编辑文档.docx

4.9方差分析

4.9.1单因素方差分析

单因素方差分析是比较两组或多组数据的均值，它返回原假设一一均值相等的概率.

函数anoval

格式p=anoval(X)%X的各列为彼此独立的样本观察值，其元素个数相同，p为各列均值相等的概率值，若P值接近于0,则原假设受到怀疑，说明至少有一列均值与其余列均值有明显不同。

p=anoval(X,group) %X和group为向量且group要与X对应

p=anova1(X,group,displayopt,) %displayopt=on∕off表示显示与隐藏方差分析

表图和盒图

[p,table]=anoval(???)%table为方差分析表

[p,table,stats]=anoval(???) %stats为分析结果的构造

说明anoval函数产生两个图：标准的方差分析表图和盒图。

方差分析表中有6歹∣J:第1列(SOUrCe)显示：X中数据可变性的来源；第2列(SS)显示：用于每一列的平方和；第3列(df)显示：与每一种可变性来源有关的自由度；第4列(MS)显示：是SS/df的比值；第5歹IJ(F)显示：F统计量数值，它是MS的比率；第6列显示：从F累积分布中得到的概率，当F增加时，P值减少。

例4-84设有3台机器，用来生产规格相同的铝合金薄板。取样测量薄板的厚度，精确至%。厘米。得结果如下：

TOC\o1-5\h\z机器1:0.236 0.238 0.248 0.245 0.243

机器2：0.257 0.253 0.255 0.254 0.261

机器3：0.258 0.264 0.259 0.267 0.262

检验各台机器所生产的薄板的厚度有无显著的差异？

解：

?X=∣0.2360.2380.2480.2450.243;0.2570.2530.2550.2540.261;-

0.2580.2640.2590.2670.262];

?P=anoval(X)

结果为：

P=

l.3431e-005

还有两个图，即图4-22和图4-23o

ANOVATableCOIUJmS 0 00105Error 0 00019Total 0 0012S图422-ProbF T]

ANOVATable

COIUJmS 0 00105

Error 0 00019

Total 0 0012S

图422

-ProbF T]

134305e-005

图4-23

例4-85建筑横梁强度的研究：3000磅力量作用在一英寸的横梁上来测量横梁的挠度,

钢筋横梁的测试强度是：82 86 79 83 84858687；其余两种更贵的合金横梁强度

测试为合金1：748278 75 76 77；合金2： 79 7977788279]。

检验这些合金强度有无明显差异？

解：

?strength=[8286798384858687748278757677797977788279];

alloy={,st,,st,st,st7st,,st?,st,,,st,al1,al1?,al1,al1?,al1,alΓ,???

,al27al2,,,al27al2,,,al27al2,);

?[ρ.(able,s(ats]=anova1(s(rength,alloy,on)

结果为

1.5264e-004table=,Source ,

1.5264e-004

table=

,Source ,SS

Groups1184.8000]

Error [1()2.0000]

Total 1286.80001

stats=

gnames:{3xlcell)n:∣866]

source:,anovaΓmeans:[847779]df:17

s:2.4495

￡27191l?H

MS[92.4000][6.00()01

[]

F[15.4000]

[I

11

,ProbF[1.5264e-004]

[]

[]

ANOVATableSourceSSdfMSFProbFGroups 184 8 2 92 4154 0.0002Error 102 17 6Tosl 286 8 19

ANOVATableSourceSSdfMSFProbF

Groups 184 8 2 92 4154 0.0002

Error 102 17 6

Tosl 286 8 19

图4-24

图425

说明P

说明P值显示，3种合金是明显不同