anova10.m

统计学习的Matlab程序:包括各种假设检验,模型识别和线性回归等.

function [h,P,Q,R,SA,Se,F,F1]=anova10(X,alpha)
%ANOVA10 单因素试验方差分析
%参考文献：《应用数理统计》（叶慈南编）p179.单因素试验方差分析
%举例：p207习题四第2题的数据X=
%[138,144,135,14,143;163,148,152,146,157;155,144,159,147,153];

%By：Ji Lin
%Email: [e-mail]
%Blog: http://linji526.spaces.live.com
%Date: 2007，2008/08/18

[r,c]=size(X);%r为单因素的不同水平的个数，c为样本的个数
n=r*c;
P=sum(sum(X.'))^2/n;
Q=sum(sum(X.').^2/c);
R=sum(sum(X.^2));
SA=Q-P;
Se=R-Q;
MSA=SA/(r-1);
MSe=Se/(n-r);
F=MSA/MSe;
F1=finv(1-alpha,r-1,n-r);
if F<F1
    h='无显著差异';
else
    h='有显著差异';
end