MATLAB中的 anova1()函数&nbs…

Matlab 统计工具箱中单因素方差分析的命令是 anoval 。 

若各组数据个数相等，称为均衡数据。若各组数据个数不等，称非均衡数据。 

（1 ）均衡数据 

处理均衡数据的用法为： 

p=anoval(x)  

返回值p 是一个概率，当 p>a 时接受H0，x为m×r的数据矩阵，x 的每一列是一个水平的数据（这里各个水平上的样本容量nj=m）另外，还输出一个方差表和一个Box图。 

例1   为考察 5 名工人的劳动生产率是否相同，记录了每人 4 天的产量，并算出其平均值，如表3 。你能从这些数据推断出他们的生产率有无显著差别吗？ 

解   编写程序如下： 

x=[256     254     250     248     236 

    242     330     277     280     252 

    280     290     230     305     220 

    298     295     302     289     252]; 

p=anova1(x) 

求得p=0.1109>a=0.05，故接受H0，即5名工人的生产率没有显著差异。方差表对应于上面的单因素方差分析表的1~4列， F=262.2是F（4，15）分布的1-p分位数，可以验证

fcdf(2.262,4,15)=0.8891=1-p 

Box 图反映了各组数据的特征。 

注：接受H0，是将 5 名工人的生产率作为一个整体进行假设检验的结果，并不表明取其中 2 个工人的生产率作两总体的均值检验时，也一定接受均值相等的假设。实际上，读者可以用 ttest2 对本题作 的检验，看看会得到什么结果。

（2）非均衡数据  

处理非均衡数据的用法为： 

p=anova1(x,group) 

x 为向量，从第 1 组到第 r 组数据依次排列；group 为与 x 同长度的向量，标志 x 中数

据的组别（在与 x 第i 组数据相对应的位置处输入整数i（i1，2....，r））。

例 2   用 4 种工艺生产灯泡，从各种工艺制成的灯泡中各抽出了若干个测量其寿命，结果如下表，试推断这几种工艺制成的灯泡寿命是否有显著差异。 

x=[1620     1580     1460     1500 

    1670     1600     1540     1550 

    1700     1640     1620     1610 

    1750     1720     1680     1800]; 

x=[x(1:4),x(16),x(5:8),x(9:11),x(12:15)]; 

g=[ones(1,5),2*ones(1,4),3*ones(1,3),4*ones(1,4)]; 

p=anova1(x,g) 

求得 0.01

<0.05，所以几种工艺制成的灯泡寿命有显著差异。