options nocenter nodate;

data one;
  input spacing density rep y;
  cards;
1     1     1      0.241
1     1     2      0.233
1     1     3      0.157
1     1     4      0.378
1     2     1      0.304
1     2     2      0.204
1     2     3      0.191
1     2     4      0.302
1     3     1      0.207
1     3     2      0.235
1     3     3      0.233
1     3     4      0.212
1     4     1      0.307
1     4     2      0.268
1     4     3      0.172
1     4     4      0.237
2     1     1      0.256
2     1     2      0.189
2     1     3      0.153
2     1     4      0.195
2     2     1      0.226
2     2     2      0.249
2     2     3      0.191
2     2     4      0.182
2     3     1      0.347
2     3     2      0.309
2     3     3      0.343
2     3     4      0.237
2     4     1      0.288
2     4     2      0.180
2     4     3      0.252
2     4     4      0.314
;

proc glm data=one;
  class rep spacing density;
  model y=rep spacing rep*spacing
          density spacing*density;
  random rep rep*spacing;
  test h=spacing e=rep*spacing;
run;

proc mixed method=type3 data=one;
  class rep spacing density;
  model y=spacing density spacing*density / ddfm=satterthwaite;
  random rep rep*spacing;
  lsmeans spacing density;
  lsmeans spacing*density / slice=spacing;
  lsmeans spacing*density / slice=density;
run;