今天分享一个SAS画图小技巧。画图时，会出现某个分组没有数据的情况，这时对于缺失的数据会输出一个灰色的图例（这通常不是我们想要的），如何将这个图例去掉呢？

以sashelp.class数据集为例，我们以性别为分组，X轴为height，Y轴为weight画一个散点图。

#随机置空几行数据#;
data class;
 set sashelp.class;
 if _n_ in (1,5,7) then call missing(sex,height,weight);
run;

#画散点图#;
proc sgplot data=class;
  styleattrs datacontrastcolors=(green blue);
  scatter x=height y=weight/ group=sex markerattrs=(symbol=starfilled);
  xaxis values=(55 to 75 by 5);
  yaxis values=(40 to 160 by 20);
run;

筛掉缺失组别的数据

第一种方法，非常简单，只要将分组缺失的数据where排除掉就可以了。

proc sgplot data=class;
  where sex>'';
  styleattrs datacontrastcolors=(green blue);
  scatter x=height y=weight/ group=sex markerattrs=(symbol=starfilled);
  xaxis values=(55 to 75 by 5);
  yaxis values=(40 to 160 by 20);
run;

可能有一些其他的情况，比如，我们画的图比较复杂，一张图里需要好几个xxplot语句，这些缺失分组的数据，在xxplot1中没有用到，但在xxplot2中又是需要画出来的数据，这种情况下，就不能采用这种筛选数据的方式了。

plot语句选项控制

SAS有相应的选项来控制不输出分组缺失的数据。在GTL中，使用includemissinggroup=false；在SGPLOT中，使用nomissinggroup。

proc sgplot data=class;
  styleattrs datacontrastcolors=(green blue);
  scatter x=height y=weight/ group=sex markerattrs=(symbol=starfilled) nomissinggroup;
  xaxis values=(55 to 75 by 5);
  yaxis values=(40 to 160 by 20);
run;

proc template;
define statgraph scatterplot;
 begingraph /backgroundcolor = white designwidth = 800px designheight = 440px;;
  layout lattice / columns=1; 
   layout overlay / walldisplay=all
        yaxisopts = (linearopts = (viewmin=40 viewmax=160
                     tickvaluesequence = (start = 40 end = 160 increment = 20))
                     display = (tickvalues ticks label line))

        xaxisopts = (linearopts = (viewmin=55 viewmax=75
                     tickvaluesequence = (start = 55 end = 75 increment = 5) tickvaluefitpolicy=none) 
                     display = (tickvalues ticks label line));
    scatterplot x=height y=weight/ name="scatter"  group=sex markerattrs=(symbol=starfilled) includemissinggroup=false;
        discretelegend "scatter"/ title="性别" location = outside;
      endlayout;
    endlayout;
  endgraph;
end;
run;

proc sgrender data=class template=scatterplot;
run;

Legend选项控制

在legend语句后加上exclude= (" ")，GTL同理。

proc sgplot data=class;
  styleattrs datacontrastcolors=(green blue);
  scatter x=height y=weight/ group=sex markerattrs=(symbol=starfilled);
  xaxis values=(55 to 75 by 5);
  yaxis values=(40 to 160 by 20);
  discretelegend/exclude= (" ");
run;

使用legenditem自定义图例

SGPLOT代码示例

proc sgplot data=class;
  styleattrs datacontrastcolors=(green blue);
  scatter x=height y=weight/ group=sex markerattrs=(symbol=starfilled);
  xaxis values=(55 to 75 by 5);
  yaxis values=(40 to 160 by 20);
  legenditem type=marker name="m_marker"/ markerattrs=(color=blue  symbol=starfilled) label="男";
  legenditem type=marker name="f_marker"/ markerattrs=(color=green symbol=starfilled) label="女";
  discretelegend "m_marker" "f_marker"/title="性别";
run;

GTL代码示例

data attrmap;
 length id $9 value markercolor markersymbol $100;
 id='sex';value='男';markercolor='blue';markersymbol='starfilled';output;
 id='sex';value='女';markercolor='green';markersymbol='starfilled';output;
run;

proc template;
define statgraph scatterplot;
  begingraph /backgroundcolor = white designwidth = 800px designheight = 440px;
    legenditem type=marker name="m_marker"/ markerattrs=(color=blue symbol=starfilled) label="男";
    legenditem type=marker name="f_marker"/ markerattrs=(color=green symbol=starfilled) label="女";
    layout lattice / columns=1; 
      layout overlay / walldisplay=all
        yaxisopts = (linearopts = (viewmin=40 viewmax=160
                     tickvaluesequence = (start = 40 end = 160 increment = 20))
                     display = (tickvalues ticks label line))

        xaxisopts = (linearopts = (viewmin=55 viewmax=75
                     tickvaluesequence = (start = 55 end = 75 increment = 5) tickvaluefitpolicy=none) 
                     display = (tickvalues ticks label line));

        scatterplot x=height y=weight/ name="scatter"  group=sex;
        discretelegend "m_marker" "f_marker"/title="性别" location = outside;
      endlayout;
    endlayout;
  endgraph;
end;
run;

proc sgrender data=class template=scatterplot dattrmap=attrmap;
  dattrvar sex='sex';
run;

legenditem可以完全按照我们的需求定制图例，适用于一些特殊情况，除了去掉缺失分组的图例，正常数据的图例也可以不让它显示，或者对于某些组别，我们知道只有x种可能取值，但是有时候数据中只有x-2，x-1种，这时SAS其他语句输出的图例就是图中出现的组别才显示相应的图例，如果需要将这x种的图例都强制画出来，也可以使用legenditem语句。

最终输出的效果如图：

大家可以根据自己的需求灵活选用以上方法，有其他方法也可以留言或者私信哦！