options nodate nonumber ps=55 ls=80; title 'Kuehl -- Example 3.1 Grain/Density p. 85'; data grain; input density grain @@; datalines; 10 12.2 10 11.4 10 12.4 20 16.0 20 15.5 20 16.5 30 18.6 30 20.2 30 18.2 40 17.6 40 19.3 40 17.1 50 18.0 50 16.4 50 16.6 ; proc print; run; proc means mean std; class density; var grain; run; proc glm; class density; model grain=density; contrast 'Linear ' density -2 -1 0 1 2; contrast 'Quadratic ' density 2 -1 -2 -1 2; contrast 'Cubic ' density -1 2 0 -2 1; contrast 'Quartic ' density 1 -4 6 -4 1; run; proc glm; class density; model grain=density; estimate 'Linear ' density -2 -1 0 1 2; estimate 'Quadratic ' density 2 -1 -2 -1 2; estimate 'Cubic ' density -1 2 0 -2 1; estimate 'Quartic ' density 1 -4 6 -4 1; run; quit;