proc iml;
use mult;   ßmult is the data we used in class.
read all into v;
n=41; 
a=j(n,1,1);y=j(n,3,1);b=j(n,1,1);bnorm=j(1,3,1); /**initializing matrices - needed for all matrices**/
a1=(1/n)*a; 
meanv=a1`*v; 
vcov=(1/n)*v`*(I(41) - (1/n)*a*a`)*v;
k=1:3;
count=0;
do i=1 to n;
b=v[i,k];
bnorm=b-meanv;
chisq=bnorm*inv(vcov)*bnorm`;
if chisq < 7.5 then count=count+1;
end;
print v;
print meanv;
print vcov;
print count;
run;