/* Example 4.3.1.------- testing mu=mu0 */
proc iml;
sugar=  {30 90 -10 10 30 60  0 40};
systol= {-8  7  -2  0 -2  0 -2  1};
diastol={-1  6   4  2  5  3  4  2};
n=8; f=8-1; p=3;
x=sugar`||systol`||diastol`;
xb=J(n,1,1)*J(1,n,1/n)*x;
s=(x-xb)`*(x-xb)/f;
Tsq=n*J(1,n,1/n)*x*inv(s)*x`*J(n,1,1/n);
F0=(f-p+1)/(f*p)*Tsq;
Fa=finv(0.95,3,5);
pv=1-probf(F0,3,5);
print F0 FA pv;
if F0> Fa then print "F0 > Fa: reject";
         else print "F0 < Fa: accept"
quit;
          F0          Fa           PV
   18.824668    5.4094513      0.00373
        F0 > Fa: reject