Options ls=78 nodate;

FILENAME instore 'beefpork.dat';

DATA one;
INFILE instore DLM='09'x DSD MISSOVER;
INPUT year pbe cbe  ppo cpo pfo dinc cfo rdinc rfp;

PROC REG DATA=one;
  MODEL pbe = cbe cfo rfp year / P R INFLUENCE;
  output out=resid1 R=resid;
  paint student.<-2 or student.>+2 / symbol='*';
  plot R.*P. student.*P. / hplots=2 vplots=2;
  plot pbe*cbe STUDENT.*cbe
       pbe*cfo STUDENT.*cfo
       pbe*rfp STUDENT.*rfp
       pbe*year STUDENT.*year/hplots=2 vplots=2;
  TITLE 'Best model Residuals plot';


Proc plot data=resid1 vpct=50 hpct=50;
  plot resid*pbe resid*cbe resid*cfo resid*rfp resid*year;
  Title 'seperate plot';

Proc reg corr data=one;
 MODEL pbe = cbe cfo rfp year / VIF PARTIAL;
 PLOT pbe*cbe pbe*cfo pbe*rfp pbe*year;

proc UNIVARIATE DATA=resid1 NORMAL PLOT;
  VAR resid;
  TITLE 'NORMAL PLOT OF RESIDUAL FROM THE BEEF PRICE REDUCED MODEL';

Proc chart data=resid1;
  var resid;

Submit to see output