DM 'LOG;CLEAR;OUT;CLEAR;';
options nodate nonumber ls=76 ps=66;

data in; input gatorid gatorlen gatorwgt @@;
     gatorsq = gatorlen**2;
cards;
1   94 130     2  74  51     3  82  80     4  58  28     5  86  80
6   94 110     7  63  33     8  86  90     9  69  36    10  72  38    
11  85  84    12  86  83    13  88  70    14  72  61    15  74  54 
16  61  44    17  90 106    18  89  84    19  68  39    20  76  42 
21  78  57    22  90 102   
;
/* Use proc surveyreg to estimate the average alligator weight */

proc surveyreg data=in ; 
    model gatorwgt = gatorlen gatorsq / clparm solution df=19; 
 
/* To estimate a mean substitute 1 for intercept, mean(x) for gatorwgt */
   estimate 'Average gator weight' intercept 1 gatorlen 90 gatorsq 8100;
run;