clear all
aa=load('../Données/hurruc.data');
t=0.5*(aa(:,1)+aa(:,2));
h=aa(:,8);
xb=mean(h);
s=std(h);
nb=4;
N=[normq(0.25),normq(0.5),normq(0.75)]*s+xb;
obs=[length(h(h<=N(1))),length(h(h>N(1)&h<=N(2))),...
    length(h(h>N(2)&h<=N(3))),length(h(h>N(3)))];
att=[4,4,4,4];
U=sum((obs-att).^2/att)
chisqp(U,1)