options nodate nonumber ps=55 ls=76; title 'Kuehl -- Example 7.5 -- Spectrophotometer Data'; data one; input conc day run glucose @@; datalines; 1 1 1 41.2 1 1 1 42.6 1 1 2 41.2 1 1 2 41.4 1 2 3 39.8 1 2 3 40.3 1 2 4 41.5 1 2 4 43.0 1 3 5 41.9 1 3 5 42.7 1 3 6 45.5 1 3 6 44.7 2 1 1 135.7 2 1 1 136.8 2 1 2 143.0 2 1 2 143.3 2 2 3 132.4 2 2 3 130.3 2 2 4 134.4 2 2 4 130.0 2 3 5 137.4 2 3 5 135.2 2 3 6 141.1 2 3 6 139.1 3 1 1 163.2 3 1 1 163.3 3 1 2 181.4 3 1 2 180.3 3 2 3 173.6 3 2 3 173.9 3 2 4 174.9 3 2 4 175.6 3 3 5 166.6 3 3 5 165.5 3 3 6 175.0 3 3 6 172.0 ; proc glm; class conc day run; model glucose=conc day conc*day run(day) conc*run(day); test h=conc e=conc*day; test h=conc*day e=conc*run(day); test h=run(day) e=conc*run(day); random day conc*day run(day) conc*run(day); run; proc mixed; class conc day run; model glucose=conc /solution; random day conc*day run(day) conc*run(day); run; quit;