%newsas(piecewrk);
title 'Piecework Operation (nwk 9.7)';
data piecewrk;
   Label Y = 'Productivity'
         X = 'Employee age' X2='Age**2' X3='Age**3'
         D='Age-40'  D2='D**2' D3='D**3';
   Input X Y;
   X2  =  X**2;
   X3  =  X**3;
   D   =  X - 40;     /* deviation from mean=40 */
   D2  =  D**2;
   D3  =  D**3;
datalines;
20 97
20 93
25 99
25 105
30 109
30 106
35 109
35 111
40 100
40 105
45 97
45 101
50 105
50 103
55 105
55 109
60 112
60 110
;
proc print;
Proc Reg data=Piecewrk;
Raw:     Model Y = X X2 X3 /ss1 ss2;
Dev:     Model Y = D D2 D3 /ss1 ss2 ;
         plot y  * d = '*'
              p. * d = '#' / overlay;
run; 
/*  9.7(c,d)
Proc GLM data=Piecewrk;
         Model Y = X X2 X3 /ss1 ss2;
         Estimate 'Age 53' INTERCEPT 1  X 53  X2 2809  X3 148877 ;
         Estimate 'Age 58' INTERCEPT 1  X 58  X2 3364  X3 195112 ;
         Estimate 'Age 62' INTERCEPT 1  X 62  X2 3844  X3 238328 ;
*/