options ls=75 nodate;

data one; 
   input x y;
   label y="height(cm)";
   if x=1 then numpass=0;
   else if x=2 then numpass=25;
   else if x=3 then numpass=75;
   else if x=4 then numpass=200;
   else if x=5 then numpass=500;
   lnumpass=log(1+numpass);
   cnumpass=numpass;  
   yid=_n_;
/* a duplicate version of numpass (for use in a class statement) */
   cards;
            1         20.7
            2         12.9
            3         11.8
            4          7.6
            5          7.8
            1         15.9
            2         13.4
            3         12.6
            4          9.5
            5          9.0
            1         17.8
            2         12.7
            3         11.4
            4          9.9
            5          8.5
            1         17.6
            2          9.0
            3         12.1
            4          9.0
            5          6.7
;
run;
proc glm;
   class numpass;
   model y=numpass;
   output out=anova p=ymean;
run;   
proc reg;
   model y=numpass;
   output out=slr p=yhat;
run;  
data both;
   merge anova slr; 
   by yid;
   lof=(yhat-ymean);
   lof2=lof*lof;
run;
proc sort;  by numpass; run;
proc print;
   sum lof lof2 ymean yhat;
run;