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;
![msrenee1984_1-1709124774127.png msrenee1984_1-1709124774127.png](/t5/image/serverpage/image-id/94197i482CC14485C3E327/image-size/medium?v=v2&px=400)
![msrenee1984_0-1709124730095.png msrenee1984_0-1709124730095.png](/t5/image/serverpage/image-id/94196iCF09A9970E0EC14B/image-size/medium?v=v2&px=400)