% 20061203 Jan Rohlén
%% Gauge R&R.  Data från Montogmerys bok sidan 358

X=[
37	38	37	41	41	40	41	42	41
42	41	43	42	42	42	43	42	43
30	31	31	31	31	31	29	30	28
42	43	42	43	43	43	42	42	42
28	30	29	29	30	29	31	29	29
42	42	43	45	45	45	44	46	45
25	26	27	28	28	30	29	27	27
40	40	40	43	42	42	43	43	41
25	25	25	27	29	28	26	26	26
35	34	34	35	35	34	35	34	35];


% Transponerar data för att kunna använda kommandot anova2
% Här skriver jag om matrisen för tydlighets skull. Normalt skulle man ha
% skrivit X=X' för transponeringen.
X=[
%Prodnr:
%   1      2    3      4    5      6    7      8     9    10
    37    42    30    42    28    42    25    40    25    35    %  Operatör 1 rep 1
    38    41    31    43    30    42    26    40    25    34    %  Operatör 1 rep 2
    37    43    31    42    29    43    27    40    25    34    %  Operatör 1 rep 3
    41    42    31    43    29    45    28    43    27    35    %  Operatör 2 rep 1
    41    42    31    43    30    45    28    42    29    35    %  Operatör 2 rep 2 
    40    42    31    43    29    45    30    42    28    34    %  Operatör 2 rep 3
    41    43    29    42    31    44    29    43    26    35    %  Operatör 3 rep 1
    42    42    30    42    29    46    27    43    26    34    %  Operatör 3 rep 2
    41    43    28    42    29    45    27    41    26    35    %  Operatör 3 rep 3
];


rep=3;                      % 3 test på varje produkt och operatör.

[p,table]=anova2(X, rep);     % Två faktorer, tre repetitioner

% Kvadratsummorna.
SSp=table{2,2}       % Observera krullparentesen för att få ut resultatet från table. Table är av typen Cell Array.
SSo = table{3,2}
SSop=table{4,2}
SSe=table{5,2}
SStot=table{6,2}

% Medelkvadratsummorna.
MSp=table{2,4}
MSoper=table{3,4}
MSop=table{4,4}
MSe=table{5,4}