options ls=80;
data toluca;
input x y;
label x=Lot Size y=Work Hours;
cards;
   80  399
   30  121
   50  221
   90  376
   70  361
   60  224
  120  546
   80  352
  100  353
   50  157
   40  160
   70  252
   90  389
   20  113
  110  435
  100  420
   30  212
   50  268
   90  377
  110  421
   30  273
   90  468
   40  244
   80  342
   70  323
;

proc reg data=toluca;
  model y=x/cli clm;
output out=howdy p=pred u95=upper95i l95 =lower95i u95m=u95mean
l95m=l95mean;

Test1: test x=0;

proc print data=howdy;


proc plot data=howdy;
  plot  upper95i*x='U'
	pred*x='P'
	lower95i*x='L'
	y*x='*'/overlay;

proc plot data=howdy;
  plot  u95mean*x='U'
        pred*x='P'
        l95mean*x='L'
        y*x='*'/overlay;