proc iml ;                                                                  
n=10;k=2;                                                                       
    x={0.2   30  ,                                                            
       1.0   50  ,                                                              
       3.0   60  ,                                                              
       5.0   70  ,                                                              
       8.5   75  ,                                                              
       10.0  80  ,                                                              
       14.0  85  ,                                                             
       20.0  90  ,                                                             
       29.0  95  ,                                                             
       43.0  100 };                                                             
y=x[,1];z=j(n,1,1);                                                             
x[,1]=z;                                                                        
h=x*inv(x`*x)*x`; m=I(n)-h;m=m/(n);                                            
ylog=log(y);g=sum(ylog);                                                        
do l=.0 to .1 by .01;                                                           
 if l=0 then yl=ylog;else yl=(y##l-z)/l;                                       
 s=yl`*m*yl;                                                                                                                                
 s=log(s);                                                                     
 s=s#n/2;gg=(l-1)#g;                                                           
 s=-s+gg;print l s;                                                             
end;