% Analyse de la variance a un facteur

% Entree des effectifs

N=[4 7 5 8 6 6];p=length(N);n=sum(N);
% Entree des donnees

Y1=[1602 1615 1617 1624];
Y2=[1480 1482 1485 1493 1500 1507 1510];
Y3=[1548 1555 1559 1563 1575];
Y4=[1435 1438 1448 1449 1454 1458 1467 1475];
Y5=[1493 1498 1500 1502 1509 1510];
Y6=[1596 1599 1602 1604 1612 1620];

Y=[Y1 Y2 Y3 Y4 Y5 Y6]';
X=zeros(n,p);
min=0;max=n;
for i=1:p,
    X(:,i)=[zeros(min,1);ones(N(i),1);zeros(max-min-N(i),1)];
    min=min+N(i);
end

% Estimation
XX=X'*X;
m=inv(XX)*X'*Y;
Ychap=X*m;
sigma2=(norm(Y-Ychap)^2)/(n-p);
Ytild=mean(Y)*ones(n,1);

% Test de Fisher
dF=norm(Ychap-Ytild)^2;
F=dF/((p-1)*sigma2);

% Contrastes
Z=[];
for i=1:p-1,
    for j=i+1:p,
    d=m(i)-m(j);q=(1/N(i))+(1/N(j));
    Z=[Z d/sqrt(q*sigma2)];
    end,
end,