proc mixed

data bond;
  input ingot metal $ pres @@;
  datalines;
1   n   67.0  1   i   71.9  1   c   72.2
2   n   67.5  2   i   68.8  2   c   66.4
3   n   76.0  3   i   82.6  3   c   74.5
4   n   72.7  4   i   78.1  4   c   67.3
5   n   73.1  5   i   74.2  5   c   73.2
6   n   65.8  6   i   70.8  6   c   68.7
7   n   75.6  7   i   84.9  7   c   69.0
;

proc mixed data=bond method=reml;
   class ingot metal;
   model pres=metal;
   random ingot;
   lsmeans metal / diff=control('n') adjust=dunnett;
run;

/* by simulation */
lsmeans metal / diff=control("n") cl
                adjust=simulate(report seed=4943838 cvadjust);