PLOTS=BOXPLOT option in PROC GLIMMIX fails to create plots for random intercept models


If you specify a random intercepts model using the SUBJECT= effect syntax in the RANDOM statement and request PLOTS=BOXPLOT, the GLIMMIX procedure does not generate box plots in the ODS Graphics results. An example of commonly used syntax is:

    proc glimmix data=test plots=boxplot;
     class trt id;
     model bp = trt / dist=normal solution;
     random intecept / solution subject=id;
  run;

There are no log messages indicating any problems related to the ODS Graphics.

To circumvent the problem, do not use the SUBJECT= effect syntax. An equivalent analysis can be obtained, but in a less efficient manner requiring more memory and time to process, with the following code:

  proc glimmix data=test plots=boxplot;
     class trt id;
     model bp = trt / dist=normal solution;
     random id / solution;
  run;