data car; input weight consumption @@; datalines; 3.4 5.5 3.8 5.9 4.1 6.5 2.2 3.3 2.6 3.6 2.9 4.6 2.0 2.9 2.7 3.6 1.9 3.1 3.4 4.9 ; proc print data=car; /* Print out the data set*/ run; proc univariate data=car; /* provide the descriptive statistics of the variable(s) specified */ var consumption weight; run; proc reg data=car ; model consumption=weight; /* Model: y = a + b*x */ plot consumption*weight; /* scatter plot with fitted line */ plot consumption*weight='.' p.*weight='X' / overlay; title 'consumption vs weight'; output out=Out r=Consum_Resid; run; proc univariate normal plot data=Out; var Consum_Resid; run;