Residual_data=[
99.6673    30.2942    2.18    1    2.35    1    0    18.1
100.278    29.9931    2.51    1    1.28    1    0    18.4
100.237    29.9858    1.73    1    1.19    1    0    18.3
100.39     29.6955    -0.46    1    2.65    1    0    18.1
100.101    29.9061    3.16    1    8.55    1    0    18
99.757     28.3195    -1.34    1    0.79    1    0    15.3
100.121    29.1754     0.1    1    1.55    1    0    17
100.653    27.7479    -1.63    1    2.04    1    0    16.1
99.708     27.8292    -2.24    1    -1.44    1    0    14.6
99.9325    26.4496    -1.28    1    -2.48    1    0    13.6
];
[m,n]=size(Residual_data);
 i=1:m
    a(i)=Residual_data(i,3);
    b(i)=Residual_data(i,5);
    av=sum(a(i))/m;
    bv=sum(b(i))/m;
    a2(i)=a(i)-av;
    b2(i)=b(i)-bv;
    all1=sum(a2(i).^2);
    bll1=sum(b2(i).^2);
  
for i=1:m
       temp1=Residual_data;  
        Residual_data(i,:)=[];
        temp2=Residual_data;
[k,da]=size(temp2);
 t=1:k
   c(t)=temp2(t,3);
   d(t)=temp2(t,5);
    av2=sum(c(t))./(m-1);
    bv2=sum(d(t))./(m-1);
    a22(t)=c(t)-av2;
    b22(t)=d(t)-bv2;
    all2(t)=sum(a22(t).^2);
    bll2(t)=sum(b22(t).^2);
end
Residual_data=temp1;
return
   x1(t)=((2*m-5).*all1)./((2*m-3).*all2);
   y1(t)=((2*m-5).*bll1)./((2*m-3).*bll2);
   Residual_data=temp1; 
   
   Ftest=[x1(t);y1(t)]'

 想利用循环算出每出去一个后的f检验值,总是不对,我是新手,只有上来问问好心人! 谢谢!