function OrthExpAnalysis3
% 正交試驗的方差分析Variance Analysis of Orthogonal experiment
%
%   Author: HUANG Huajiang
%   Copyright 2003 UNILAB Research Center, 
%   East China University of Science and Technology, Shanghai, PRC
%   $Revision: 1.0 $  $Date: 2002/08/18 $
%
%   [Ref] 胡上序和陳德釗編著，實驗數(shù)據(jù)的分析和處理，杭州：浙江大學出版社，
%   1996，p.273
%
%   DATA: 正交試驗方案及結果---最后一列為試驗指標（收率）,其余是以水平表示的正交試驗方案
%   NF: 因素數(shù)（列數(shù)）
%   NL: 水平數(shù)
%   r:  每個水平有r次試驗

clear all
clc

DATA = ...              % 最后一列為試驗指標（收率）,其余是以水平表示的正交試驗方案
   [1  1  1  1  31;
    1  2  2  2  54;
    1  3  3  3  38;
    2  1  2  3  53;
    2  2  3  1  49;
    2  3  1  2  42;
    3  1  3  2  57;
    3  2  1  3  62;
    3  3  2  1  64];
[m,n] = size(DATA);
NF = n-1;              % 因素數(shù)（列數(shù)）
NL = max(DATA(:,1))    % 水平數(shù)
r = 0;
for i=1:m;
    if DATA(i,1)==1
        r = r+1;        % 每個水平有r次試驗
    end
end

A = DATA(:,1:NF);
Y = DATA(:,NF+1);
% 計算每個因素（列）相同水平的指標總和K
for i=1:NL              % i:水平
    for j=1:NF          % j:因素（列號）
        K(i,j) = sum(Y(find(A(:,j)==i)))
    end
end

xa = sum(Y)/length(Y)
xT = sum(Y)
P = 1/(NL*r)*xT^2
j=1:NF;
Q(j) = sum(K(:,j).^2)/r
S = Q - P               % 各因素的平方和
ST = sum(S)             % 總平方和
Sav = S/(NL-1)          % 各因素的均方
% 顯著性檢驗
F = Sav(1:NF-1)/Sav(end)