data cubic;
call streaminit(8675309);
  do nrep=1 to 3;
    do x=1 to 10;
       response = 3 + 1*x - .5*x**2 + .05*x**3 +
                                      rand('normal',0,1.5);
       output;
    end;
  end;
run;

ods rtf bodytitle file="C:\Users\baileraj\Desktop\fig6p1.rtf" ;
ods trace on;
* lists names of output objects, tables, graphs produced;
ods graphics on;
proc reg data=cubic;
  model response = x;
run;
ods graphics off;
ods trace off;
ods rtf close;

ods rtf bodytitle file="C:\Users\baileraj\Desktop\fig6p1.rtf";
ods graphics on;
proc reg data=cubic plots(only)=DFBETAS;
  model response = x;
run;
ods graphics off;
ods rtf close;

ods rtf bodytitle file="C:\Users\baileraj\Desktop\fig6p1.rtf";
ods graphics on;
proc reg data=cubic;
ods select ParameterEstimates FitPlot;
   model response = x;
run;
ods graphics off;
ods rtf close;

ods listing sge=on;  * start graphics editor;
ods graphics on;
proc reg data=cubic;
  model response = x;
run;
ods graphics off;