Y=[-1.70,-4.61,-5.82,-1.17,-4.23,-3.31,0.42,-2.98];
X1=[1.75,1.33,1.43,1.86,1.81,1.66,1.6,2];
X2=[0.52,0.59,0.61,0.5,0.54,0.49,0.27,0.47];
X=[(ones(1,8))',X1',X2'];

%Estimateur des moindres carres de téta
Z=lsfit(Y',X);

P=X*inv(X'*X)*X';

%Variance estimée ajustée de manière à être sans biais
V=(1/5)*Y*(Y'-P*(Y'));

K=X*Z-(1/8)*sum(Y)*ones(8,1);
R=inv(X'*X);

%Statistique correspondant au premier test
E1=(1/(2*V))*sum(K.^2);

%Statistique correspondant au deuxième test
E2=(Z(2)/sqrt(V*R(2,2)))^2;

%Quantile de la distribution d'une fisher(2,5)
f1=qf(0.95,2,5);

%Quantile de la distribution d'une fisher(1,5)
f2=qf(0.95,1,5);

if (E1<=f1) disp('La variation du taux d''hémoglobine ne dépend ni de la durée de l''opération ni du volume du sang perdu')
   else disp('La variation du taux d''hémoglobine dépend soit de la durée de l''opération soit du volume du sang perdu')
   end
   if (E2<=f2) disp('La variation du taux d''hémoglobine ne dépend pas de la durée de l''opération')
   else disp('La variation du taux d''hémoglobine dépend de la durée de l''opération')
      end