程序如下：
data IleumSucraseActivity;                                                                                                                                                                                                                                                      
  length treatment $20;
  input treatment$ block IleumKm IleumVmax;                                                                                 
  Non-Antibiotics 5 46.18 28.06
  Non-Antibiotics 3 55.61 24.64
  Non-Antibiotics 4 58.71 27.35
  Non-Antibiotics 3 47.32 22.15
  Non-Antibiotics 4 43.61 12.53
  Non-Antibiotics 8 52.00 25.94
  Non-Antibiotics 11 54.65 21.34
  Non-Antibiotics 5 64.54 25.73
  Non-Antibiotics 10 51.27 24.36
  Non-Antibiotics 5 57.61 23.04
  Non-Antibiotics 5 48.37 32.42
  Non-Antibiotics 9 56.04 40.59
  Non-Antibiotics 6 53.89 11.33
  Antibiotics 8 39.02 27.74
  Antibiotics 4 55.56 16.92
  Antibiotics 7 56.97 25.41
  Antibiotics 9 58.29 31.49
  Antibiotics 4 53.18 28.04
  Antibiotics 11 41.89 43.48
  Antibiotics 5 40.35 38.12
  Antibiotics 3 54.84 40.96
  Antibiotics 10 47.66 32.51
  Antibiotics 6 61.14 35.47
  Antibiotics 5 52.42 38.02
  ;
  run;
Proc print data=Activity;
run;
proc glimmix data=Activity nobound;
  class treatment block;
  model IleumKm = treatment;
  random block/subject=treatment*block type=cs residual;
  random _residual_/subject=treatment*block group=treatment;
  lsmeans treatment/tdiff lines;
  covtest homogeneity;
run;


结果：
ERROR: QUANEW Optimization cannot be completed.