/* Example 4.4.2---Paired T^2-Test */
proc iml;
x={2600 2350, 2850 2350, 2800 2250, 2750 2200, 2400 2250, 2450 2300};
y={1450 1300, 1750 1450, 1700 1350, 1650 1300, 1350 1250, 1400 1300};
n=6; p=2; f=n-1;
d=x-y;
db=J(1,n,1/n)*d;
ddb=J(n,1,1)*db;
sd=(d-ddb)`*(d-ddb)/f;
F0=n*(f-p+1)/(f*p)*db*inv(sd)*db`;
Fa=finv(1-0.05,p,f-p+1);
print F0 Fa ;
if F0>Fa then print "F0 > Fa:  reject";
         else print "F0 < Fa:  accept";
Ta=sqrt(fa*10/4);
Tfa4=tinv(1-0.05/4, f);
print Ta Tfa4;
e1=Tfa4*sqrt(sd[1,1]/n);
e2=Tfa4*sqrt(sd[2,2]/n);
l1=db[1,1]-e1;
u1=db[1,1]+e1;
l2=db[1,2]-e2;
u2=db[1,2]+e2;
print "delt1=("l1 ","u1"),",
      "delt2=("l2 ","u2"),";
quit;
           F0             FA
      2453.84      6.9442719
          F0 > Fa: reject
           TA           TFA4
    4.1666149      3.1633814
   delt1=( 1043.0584 , 1140.2749 ),
   delt2=( 872.50647 , 1044.1602 ),