1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75
| clear all
% Lecture des données et création des vecteurs
[club,v1,v2,v3]=textread('totbut.txt','%s %d %d %d');
% Création des vecteurs contenant les moyennes des but
% de chaque pays
x2007=[mean(v1(v3==1)),mean(v1(v3==2)),mean(v1(v3==3)),mean(v1(v3==4)),mean(v1(v3==5))]
x2008=[mean(v2(v3==1)),mean(v2(v3==2)),mean(v2(v3==3)),mean(v2(v3==4)),mean(v1(v3==5))]
% pareilavec les écarts types
s2007=[std(v1(v3==1)),std(v1(v3==2)),std(v1(v3==3)),std(v1(v3==4)),std(v1(v3==5))]
s2008=[std(v2(v3==1)),std(v2(v3==2)),std(v2(v3==3)),std(v2(v3==4)),std(v1(v3==5))]
% Vecteur contenant le nombre but par pays
n=[sum(v3==1),sum(v3==2),sum(v3==3),sum(v3==4),sum(v3==5);
% Moyennes globales pour 2007 et 2008
xb2007=sum(x2007.*n)/sum(n);
xb2008=sum(x2008.*n)/sum(n);
% Test ANOVA pour 2007
I=length(n); % Nombre de groupes
DFG2007=I-1;
DFE2007=sum(n)-I;
DFT2007=sum(n)-1;
SSG2007=sum(n.*((x2007-xb2007).^2));
SSE2007=sum((n-1).*(s2007.^2));
SST2007=sum((x2007-xb2007).^2);
MSG2007=SSG2007/(length(n)-1);
MSE2007=SSE2007/(sum(n)-length(n));
MST2007=SST2007/DFT2007;
F2007=MSG2007/MSE2007; % Statistique
% Comparaison à la variable aléatoire de Fisher
F2007=fp(F2007,length(n)-1,sum(n)-length(n))
% Tableau ANOVA
fprintf(1,'Source\t | DF\t | SS\t | MS\n')
fprintf(1,'Groupe\t | %g\t | %g\t | %g\n',DFG2007,SSG2007,MSG2007)
fprintf(1,'Erreur\t | %g\t | %g\t | %g\n',DFE2007,SSE2007,MSE2007)
fprintf(1,'Total\t | %g\t | %g\t | %g\n',DFT2007,SST2007,MST2007)
% On compare à un taux d'erreur acceptable de 5%
if F2007>0.95
display('Avec moins de 5% de chance de se tromper, en 2007 le nombre de but dépend du pays')
else
display('Avec moins de 5% de chance de se tromper, en en 2007 le nombre de but ne dépend pas du pays')
end
% Test ANOVA pour 2008
I=length(n);
DFG2008=I-1;
DFE2008=sum(n)-I;
DFT2008=sum(n)-1;
SSG2008=sum(n.*((x2008-xb2008).^2));
SSE2008=sum((n-1).*(s2008.^2));
SST2008=sum((x2008-xb2008).^2);
MSG2008=SSG2008/(length(n)-1);
MSE2008=SSE2008/(sum(n)-length(n));
MST2008=SST2008/DFT2008;
F2008=MSG2008/MSE2008;
F2008=fp(F2008,length(n)-1,sum(n)-length(n))
% Tableau ANOVA
fprintf(1,'Source\t | DF\t | SS\t | MS\n')
fprintf(1,'Groupe\t | %g\t | %g\t | %g\n',DFG2008,SSG2008,MSG2008)
fprintf(1,'Erreur\t | %g\t | %g\t | %g\n',DFE2008,SSE2008,MSE2008)
fprintf(1,'Total\t | %g\t | %g\t | %g\n',DFT2008,SST2008,MST2008)
if F2008>0.95
display('Avec moins de 5% de chance de se tromper, en en 2008 le nombre de but dépend du pays')
else
display('Avec moins de 5% de chance de se tromper, en en 2008 le nombre de but ne dépend pas du pays')
end |
Partager