clear all                                             % On vide la mémoire
aa=load('../Données/hurruc.data');                    % On lit les données
h=aa(:,8);                                            % le nombre d'ouragans est dans la 8e colonne
xb=mean(h);                                           % moyenne expérimentale
s=std(h);                                             % écart-type expérimental
nb=4;                                                 % Nombre de catégories équiprobables
N=[normq(0.25),normq(0.5),normq(0.75)]*s+xb;          % Définition des catégories
obs=[length(h(h<=N(1))),length(h(h>N(1)&h<=N(2))),... % On compte combien sont dans chaque catégorie
    length(h(h>N(2)&h<=N(3))),length(h(h>N(3)))];     % cela définit les "obs"
att=[4,4,4,4];                                        % On définit les attendus
U=sum((obs-att).^2/att)                               % Calcul de la statistique U
chisqp(U,1)                                           % Et on compare à la table de la v.a. chi^2