Code: 
options linesize=76;
options pagesize=53;
data gpadat;
 infile 'a:\anal2.dat';
 input gpa inc;
 label inc='Income'
       gpa='GPA';
title 'Income/GPA Analysis';
proc plot;
 plot inc*gpa='*';
proc reg;
 model inc=gpa;
 output out=new1 r=resid;
proc univariate normal;
 var resid;
run;


Output: 

                   Plot of INC*GPA.  Symbol used is '*'.

         |
         |
   35000 +
         |                                                              *
         |
         |                                               *    *
         |                                                    *    *
         |                                                    *    *    *
   30000 +
         |                                               *    *    *
         |                                               *    *    *
         |                                          *    *    *
         |                                     *
         |                                *    *    *    *         *
 I 25000 +                           *                   *
 n       |                      *         *    *
 c       |                                *
 o       |                      *    *         *
 m       |                           *              *
 e       |                           *         *    *
   20000 +                 *
         |            *              *
         |            *
         |
         |            *    *    *
         |       *
   15000 +       *    *    *
         |  *
         |  *
         |       *
         |       *
         |  *
   10000 +  *
         |
         ---+----+----+----+----+----+----+----+----+----+----+----+----+--
           1.6  1.8  2.0  2.2  2.4  2.6  2.8  3.0  3.2  3.4  3.6  3.8  4.0

                                         GPA



                              Income/GPA Analysis                              


Model: MODEL1  
Dependent Variable: INC        Income                                  

                            Analysis of Variance

                            Sum of         Mean
   Source          DF      Squares       Square      F Value       Prob>F

   Model            1 2039676985.8 2039676985.8      406.776       0.0001
   Error           62 310883236.89 5014245.7564
   C Total         63 2350560222.7

       Root MSE    2239.25116     R-square       0.8677
       Dep Mean   22757.41708     Adj R-sq       0.8656
       C.V.           9.83965

                            Parameter Estimates

                     Parameter      Standard    T for H0:               
    Variable  DF      Estimate         Error   Parameter=0    Prob > |T|

    INTERCEP   1    271.540275  1149.4898110         0.236        0.8140
    GPA        1   7994.978418  396.40525736        20.169        0.0001

                  Variable
    Variable  DF     Label

    INTERCEP   1  Intercept                               
    GPA        1                                          


                            Univariate Procedure

Variable=RESID         Residual

                                  Moments

                  N                64  Sum Wgts         64
                  Mean              0  Sum               0
                  Std Dev    2221.408  Variance    4934655
                  Skewness    -0.0318  Kurtosis   -0.43214
                  USS        3.1088E8  CSS        3.1088E8
                  CV                .  Std Mean    277.676
                  T:Mean=0          0  Pr>|T|       1.0000
                  Num ^= 0         64  Num > 0          34
                  M(Sign)           2  Pr>=|M|      0.7080
                  Sgn Rank          7  Pr>=|S|      0.9631
                  W:Normal   0.980895  Pr W         0.6960


                              Quantiles(Def=5)

                   100% Max  4840.703       99%  4840.703
                    75% Q3    1637.36       95%  3381.108
                    50% Med  131.0351       90%   2738.16
                    25% Q1   -1642.27       10%  -2740.63
                     0% Min  -4911.31        5%  -3583.21
                                             1%  -4911.31
                   Range     9752.011                    
                   Q3-Q1     3279.633                    
                   Mode      -4911.31                    


                                  Extremes

                     Lowest    Obs     Highest    Obs
                   -4911.31(      42) 3097.725(      29)
                    -4543.8(      61) 3381.108(      52)
                   -4075.93(      43) 4159.418(      28)
                   -3583.21(      40) 4663.046(      22)
                   -3063.48(      21) 4840.703(      48)