clear all;
clf;

y = [59 61 50 54 50 58 46 44 74 70 81 85 69 67 79 81]';
g1 = [160 160 160 160 160 160 160 160 180 180 180 180 180 180 180 180]; 
g2 = [20 20 20 20 40 40 40 40 20 20 20 20 40 40 40 40]; 
g3 = {'c1'; 'c1'; 'c2'; 'c2'; 'c1'; 'c1'; 'c2'; 'c2';'c1'; 'c1'; 'c2'; 'c2';'c1'; 'c1'; 'c2'; 'c2'};

%[p, table, stats]  = anovan(y, {g1 g2 g3}, 'varnames',{'Temperature' 'Concentration' 'Catalyst'}, 'model', 'linear');
%[p, table, stats]  = anovan(y, {g1 g2 g3}, 'varnames',{'Temperature' 'Concentration' 'Catalyst'}, 'model', 'interaction');
[p, table, stats]  = anovan(y, {g1 g2 g3}, 'varnames',{'Temperature' 'Concentration' 'Catalyst'}, 'model', 'full');
%[p, table, stats]  = anovan(y, {g1 g2 g3}, 'varnames',{'Temperature' 'Concentration' 'Catalyst'}, 'model', [0 1 0; 1 0 1]);

[c, m, h, gnames] = multcompare(stats, 'alpha', 0.05, 'ctype', 'lsd', 'dimension', [2]);
%[c, m, h, gnames] = multcompare(stats, 'alpha', 0.05, 'ctype', 'lsd','dimension', [1,3]);
%[c, m, h, gnames] = multcompare(stats, 'alpha', 0.05, 'ctype', 'lsd', 'dimension', [1,2]);

disp(sprintf('\n'));
disp([gnames, num2cell(m)]);
for i = 1 : size(c,1)
    disp(sprintf('%s\t%s\t%f\t%f\t%f', char(gnames(c(i,1),1)), char(gnames(c(i,2),1)), c(i, 3), c(i, 4), c(i, 5)));
end;