proc freq data=r.fac2;
tables facility*year / norow nocol nopct
sparse out=FreqOut(where=(percent^=.));
run;
proc print data=FreqOut;
var facility year Count;
run;
/* Create heat maps of the tabular results */
ods graphics / width=400px;
title "shoes Facility by Year ";
title2 "sneakers (N=5952)";
/* basic heat map */
proc sgplot data=FreqOut;
heatmap x=year y=facility / freq=Count
discretex discretey
colormodel=(green redTwoColorRamp)outline;
run;