y = mu + e

data test;
  input y @@;
  datalines;
59.9 91.8 82.4 59.1 82.9 70.6 77.2 69.5 78.5 85.6
run;

proc sql;
  select mean(y) as mean,  css(y)/count(y) as MLE, css(y)/(count(y)-1) as REMLE
  from test
;quit;

proc mixed data=test;
  model y= / s;
run;

proc mixed data=test method=ml;
  model y= / s;
run;

proc genmod data=test;
  model y=;
run;

/* REML: default */
proc glimmix data=test method=rspl;
  model y= / s;
run;

/* ML */
proc glimmix data=test method=mspl;
  model y= / s;
run;

proc nlmixed data=test;
  parms mu=100 s2e=100;
  model y~normal(mu,s2e);
run;