So brilliant and simple. Thank you! ... View more

Thank you @Rick_SAS! I ended up studying the box= dataset, formatting my data accordingly, and then using proc sgpanel as I am more familiar with manipulating the visuals in sgpanel. I think that the results turned out fantastic!   proc sgpanel data=fig1 noautolegend sganno=anno; panelby dpt / rows=2 columns=1 novarname uniscale=column; vbox _value_ / group=sys nomean lineattrs=(color=black) medianattrs=(color=black) whiskerattrs=(color=black) capshape=none; refline 0 / lineattrs=(color=crimson); rowaxis label=" " offsetmax=.1; keylegend / title=" " position=bottom noborder; run;    ... View more

Normally, I would use a TEXT plot for this sort of thing, as it works well in most cases. However, with a box plot, you do not have access to the numbers used to plot each box to correctly place the text around the box. I would consider using an "inside" AXISTABLE to put the significance text above or below the box plot. ... View more

To expand on @PaigeMiller 's response with your example code: The following sort and BY statement will run ALL the combinations of depth, year and version in sort order. If you want all of the same versions together in output then use a sort order in both By statements with the Version first, or the year though that may not be critical.   Proc sort data=SOC; by depth year version; run; proc glimmix data=SOC plots=studentpanel; by depth year version; nloptions maxiter=500; class year block system depth plot subplot version; model deltaSOC_Mg_ha = system / ddfm=kr; random block; random _residual_ / subject=plot group=system type=arh(1); run; Also note that the values of the By variables can be addressed in Title statements and even used in some ODS destinations to control some output aspects.   If you don't want all the combinations then you still sort and use the BY statement but modify the WHERE clause to use IN with value lists. For example (assuming some likely values for example). where depth in ("0_15" "15_30") and year in (2008 2009 2010) and version in ("new" "proposed"); Or if you want to exclude some of a long list of values:   variable NOT IN (<list>)     ... View more

@yabwon, this is fantastic! Thank you for the upate.  ... View more

For #1, use the NOVARNAME option on the PANELBY statement to drop the left-side labels (e.g. depth=). For your custom row values, create a user-defined format and set it on the SGPANEL procedure using the FORMAT statement: FORMAT resp myfmt.;   For #2, set ROWHEADERPOS=LEFT on the PANELBY statement. You might also want to set NOHEADERBORDER and HEADERBACKCOLOR to make the labels appear more like axis labels.   Hope this helps! ... View more

That is correct. ... View more