options ps=55 ls=76;
data one;
do year = 1870 to 1938;
input consume income price;
inc_pri=income/price;
output;
end;
cards;
1.9565 1.7669 1.9176
1.9794 1.7766 1.9059
2.0120 1.7764 1.8798
2.0449 1.7942 1.8727
2.0561 1.8156 1.8984
2.0678 1.8083 1.9137
2.0561 1.8083 1.9176
2.0428 1.8067 1.9176
2.0290 1.8166 1.9420
1.9980 1.8041 1.9547
1.9884 1.8053 1.9379
1.9835 1.8242 1.9462
1.9773 1.8395 1.9504
1.9748 1.8464 1.9504
1.9629 1.8492 1.9723
1.9396 1.8668 2.0000
1.9309 1.8783 2.0097
1.9271 1.8914 2.0146
1.9239 1.9166 2.0146
1.9414 1.9363 2.0097
1.9685 1.9548 2.0097
1.9727 1.9453 2.0097
1.9736 1.9292 2.0048
1.9499 1.9209 2.0097
1.9432 1.9510 2.0296
1.9569 1.9776 2.0399
1.9647 1.9814 2.0399
1.9710 1.9819 2.0296
1.9719 1.9828 2.0146
1.9956 2.0076 2.0245
2.0000 2.0000 2.0000
1.9904 1.9939 2.0048
1.9752 1.9933 2.0048
1.9494 1.9797 2.0000
1.9332 1.9772 1.9952
1.9139 1.9924 1.9952
1.9091 2.0117 1.9905
1.9139 2.0204 1.9813
1.8886 2.0018 1.9905
1.7945 2.0038 1.9859
1.7644 2.0099 2.0518
1.7817 2.0174 2.0474
1.7784 2.0279 2.0341
1.7945 2.0359 2.0255
1.7888 2.0216 2.0341
1.8751 1.9896 1.9445
1.7853 1.9843 1.9939
1.6075 1.9764 2.2082
1.5185 1.9965 2.2700
1.6513 2.0652 2.2430
1.6247 2.0369 2.2567
1.5391 1.9723 2.2988
1.4922 1.9797 2.3723
1.4606 2.0136 2.4105
1.4551 2.0165 2.4081
1.4425 2.0213 2.4081
1.4023 2.0206 2.4367
1.3991 2.0563 2.4284
1.3798 2.0579 2.4310
1.3782 2.0649 2.4363
1.3366 2.0582 2.4552
1.3026 2.0517 2.4838
1.2592 2.0491 2.4958
1.2365 2.0766 2.5048
1.2549 2.0890 2.5017
1.2527 2.1059 2.4958
1.2763 2.1205 2.4838
1.2906 2.1205 2.4636
1.2721 2.1182 2.4580
;

data two;
set one;
price_inc=price/income;
file 'A:\spirits.dat';
put @1 year @11 consume @21 income @31 price @41 price_inc;
run;

/*
proc plot;
plot consume*price consume*income consume*inc_pri;
run;
*/

proc reg;
model consume=income price /dw;
output out=all p=yhat r=e;

proc print;
var consume income price yhat e;
run;

symbol1 c=black l=1 i=join v=:;
symbol2 c=black v=:;

filename spirits1 'C:\PCTeX32\qmb3250\spirits1.ps';
goptions device=psepsf gsfname=spirits1 hsize=6.0in vsize=6.0in;

proc gplot;
plot e*year=1 / frame vref=0;
run;


/*
proc reg data=one;
model consume=inc_pri;
output out=sub p=yhat r=e;

proc plot;
plot e*year;
run;
*/

quit;