goptions reset=all; Title " "; Footnote " "; 
axis1 label=(font=swiss h=1.5 angle=90); 
*******************************************;


TITLE 'CONSTRUCTION REVIEW';
TITLE2 'WORKERS IN THOUSANDS';
DATA CONSTRT;
   INPUT ELECTRIC DATE :DATE7. MASONRY CARPENT
      ROOFERS PAPER HEATING;
elec1 = lag(electric); roof1=lag(roofers);
   FORMAT DATE MONYY5.;  mnth=month(date);
CARDS;
241.9 1JAN77  196.6 77.9 77.0 84.5 318.4
239.3 1FEB77  216.3 78.7 85.7 83.0 318.1
245.5 1MAR77  235.1 82.4 94.7 88.7 326.5
249.6 1APR77  252.2 87.9 102.8 96.4 337.8
255.7 1MAY77  261.5 92.7 107.8 103.3 348.4
262.8 1JUN77  278.0 99.5 111.0 111.9 359.7
268.1 1JUL77  280.8 103.1 115.6 116.8 370.5
276.1 1AUG77  290.5 106.5 116.4 120.6 374.6
277.1 1SEP77  287.4 103.9 117.7 115.5 374.7
278.6 1OCT77  288.7 105.8 122.4 114.3 376.7
276.2 1NOV77  286.2 102.8 123.5 112.1 377.1
274.6 1DEC77  276.7 101.3 112.5 104.3 376.2
262.1 1JAN78  234.3 91.4 93.1 89.6 361.0
256.2 1FEB78  237.4 84.8 94.2 87.8 350.9
259.9 1MAR78  270.6 93.9 109.7 94.7 362.3
263.6 1APR78  292.7 98.6 126.1 106.5 375.7
278.0 1MAY78  295.5 103.8 128.7 113.5 384.0
289.5 1JUN78  308.3 113.0 135.5 126.6 399.2
296.6 1JUL78  306.9 115.6 140.3 133.3 413.0
305.1 1AUG78  315.3 115.1 141.3 135.4 423.7
307.8 1SEP78  310.6 110.5 139.2 129.1 423.4
308.2 1OCT78  316.3 110.4 142.0 130.5 425.7
310.1 1NOV78  310.5 107.3 136.7 127.4 423.1
307.6 1DEC78  304.4 105.2 127.7 118.7 420.6
298.2 1JAN79  263.7 97.7 103.4 101.6 396.9
292.1 1FEB79  259.8 96.3 98.4 100.4 388.5
294.5 1MAR79  303.5 102.0 124.1 104.0 395.0
296.9 1APR79  307.2 98.9 128.9 112.8 398.3
303.3 1MAY79  316.8 101.3 133.5 124.1 407.3
312.1 1JUN79  325.4 103.2 139.2 134.1 416.8
322.8 1JUL79  328.6 107.4 144.0 143.3 425.3
327.7 1AUG79  335.2 107.2 145.9 145.1 430.9
327.0 1SEP79  327.9 106.2 145.6 137.9 426.9
330.4 1OCT79  334.5 105.5 146.3 136.4 426.7
329.7 1NOV79  328.7 100.4 144.9 129.5 424.2
329.5 1DEC79  321.5 101.5 141.1 121.5 421.9
320.3 1JAN80  293.8 91.6 119.8 106.8 409.8
315.7 1FEB80  285.8 89.7 114.6 104.7 398.8
312.6 1MAR80  292.0 88.6 119.0 104.1 394.6
306.5 1APR80  299.3 86.4 133.9 114.8 388.3
308.6 1MAY80  301.7 88.6 135.7 122.9 390.3
316.3 1JUN80  307.9 80.0 138.1 128.7 396.7
319.5 1JUL80  310.7 93.0 140.1 131.8 395.9
326.4 1AUG80  314.9 96.9 142.9 130.1 405.1
330.8 1SEP80  312.7 95.8 142.9 126.1 407.9
329.3 1OCT80  318.5 96.3 144.7 126.8 410.5
330.6 1NOV80  307.7 94.2 144.2 119.2 406.9
327.2 1DEC80  296.2 93.0 132.6 110.0 399.6
316.2 1JAN81  259.2 84.8 107.4 93.2 388.2
310.1 1FEB81  258.8 80.9 102.0 90.7 379.8
308.5 1MAR81  271.5 84.0 122 94.2 378.8
311.1 1APR81  281.0 92.0 120.4 105.1 382.2
313.6 1MAY81  283.7 95.4 122.1 112.2 379.4
318.3 1JUN81  289.3 95.7 129.2 120.8 389.0
321.3 1JUL81  291.1 95.7 132.6 126.2 393.7
327.4 1AUG81  295.9 96.7 135.7 126.1 397.1
326.7 1SEP81  292.7 92.4 135.4 122.3 393.2
326.4 1OCT81  282.6 88.9 136.4 118.5 393.7
322.5 1NOV81  275.5 84.9 135.9 113.7 387.1
318.6 1DEC81  260.2 80.5 125.2 104.2 379.7
301.9 1JAN82  214.3 70.9 95.8 89.2 363.2
296.1 1FEB82  224.8 71.3 96.0 86.9 354.4
298.3 1MAR82  228.7 73.4 101.6 90.2 350.6
297.7 1APR82  244.7 76.3 115.0 100.4 351.5
303.5 1MAY82  260.4 83.3 120.2 112.2 358.5
305.0 1JUN82  262.2 85.9 125.3 117.2 362.1
307.6 1JUL82  270.4 85.8 129.4 123.1 367.2
;

proc print;
proc glm; class mnth; model electric = elec1 roof1 mnth/solution;
proc glm; class mnth; model roofers  = elec1 roof1 mnth/solution;

data a;
  do t=1 to 30;  t1=t+.02;
     if t < 5 then do; e=0; r=0; end;
     if t=5 then do; e =  0; r=30; end;
     if t>5 then do; e = .8707*e1 + .2268*r1;
                     r = .0064*e1 + .8738*r1;
     end;
electric = 301 + e; roofers=124+r;
output; r1=r;e1=e; end;
proc print; var e r;
proc gplot; plot e*t  r*t1/overlay legend;
  symbol1 v=dot i=needle;
run;