data crdreg1;
input n wt;
nc = n;
cards;
10  20.5
10  21.6
10  22.9
20  36.4
20  38.5
20  38.0
30  43.1
30  44.4
30  44.2
40  45.0
40  45.5
40  46.1
50  46.4
50  47.3
50  48.4
;
proc glm data=crdreg1;
classes n;
model wt = n;
lsmeans n/stderr pdiff;
estimate ' n = 10' intercept 1 n 1 0 0 0 0;
estimate ' n = 20' intercept 1 n 0 1 0 0 0;
estimate ' n = 30' intercept 1 n 0 0 1 0 0;
estimate ' n = 40' intercept 1 n 0 0 0 1 0;
estimate ' n = 50' intercept 1 n 0 0 0 0 1;
contrast 'Cubic+Quartic' n -1 2 0 -2 1,
                         n 1 -4 6 -4 1;
run;
quit;

proc glm data=crdreg1;
model wt = n n*n n*n*n n*n*n*n/solution;
estimate 'n, linear ' n 1;
estimate 'n, quadratic ' n*n 1;
estimate 'n, cubic ' n*n*n 1;
estimate 'n, quartic ' n*n*n*n 1;
contrast 'Regressions' n  1,
                       n*n 1,
                       n*n*n 1,
                       n*n*n*n 1;
contrast 'N3 and N4' n*n*n 1,
                     n*n*n*n 1;
estimate 'N = 10' intercept 1 n 10 n*n 100 n*n*n 1000 n*n*n*n 10000;
estimate 'N = 20' intercept 1 n 20 n*n 400 n*n*n 8000 n*n*n*n 160000;
estimate 'N = 30' intercept 1 n 30 n*n 900 n*n*n 27000 n*n*n*n 810000;
estimate 'N = 40' intercept 1 n 40 n*n 1600 n*n*n 64000 n*n*n*n 2560000;
estimate 'N = 50' intercept 1 n 50 n*n 2500 n*n*n 125000 n*n*n*n 6250000;
run;
quit;

proc glm data=crdreg1;
model wt = n n*n/solution;
run;
quit;

proc glm data=crdreg1;
classes nc;
model wt = n n*n nc;
run;
quit;