close all
clear all

% Définition des données pour l'ANOVA : y est la réponse, g1 et g2 sont les
% facteurs
 y = [1	7	1	9	6	2	2	6	2	8	8	2]';
% y = [2 8 3 2 2 9 2 7 2 3 2 8]';

% nom des facteurs
varnames = {'bracelet';'cadran'};
g1 = {'métal','cuir','plastique','plastique','métal','cuir','métal','cuir','plastique','plastique','métal','cuir'};
g2 = {'aiguille','digitale','aiguille','digitale','digitale','aiguille','aiguille','digitale','aiguille','digitale','digitale','aiguille'};

% ANOVA à 2 facteurs, sans interaction
[p,tbl1,stats1,terms1]=anovan(y,{g1 g2});
%[p,tbl1,stats1,terms1]=anovan(y,{g1 g2},'model','interaction','varnames',varnames);
% coefficient de détermination R2
Rsquared=1-tbl1{size(tbl1,1)-1,2}/tbl1{size(tbl1,1),2}
% coeffficients du modèle
for i=1:length(stats1.coeffs)
    disp([cell2mat(stats1.coeffnames(i)),' :  ',num2str(stats1.coeffs(i))]);
end
y_est=y+stats1.resid;
figure;
plot(y_est,y,'k*');xlabel('y_e_s_t');ylabel('y_m_e_s');
axis([1 9 1 9]);
line([1,9],[1,9],'Color','red','LineStyle','--');
title(['R2 = ',num2str(Rsquared)]);
%importances
Amp(1)=max(stats1.coeffs(2:4))-min(stats1.coeffs(2:4));
Amp(2)=max(stats1.coeffs(5:6))-min(stats1.coeffs(5:6));
Imp=Amp/sum(Amp);
disp(['fact 1 : Amp=',num2str(Amp(1)),'; Imp=',num2str(Imp(1))]);
disp(['fact 2 : Amp=',num2str(Amp(2)),'; Imp=',num2str(Imp(2))]);