/*--Fig 4.6.2 Heat map--*/
data GTL_GS_HeatmapParm;
  drop pi val;
  pi=constant('Pi');
  do x=1 to 20;
    do y=1 to 10;
      value=sin(x*pi/10) + cos(y*pi/10);
      val=ranuni(2);
      group=ifc(val < 0.5, 'M', 'F');
      output;
    end;
  end;
run;
proc print data=GTL_GS_HeatmapParm; run;

ods html;

proc template;
define statgraph _heatmapparm;
    begingraph;
        layout overlay / xaxisopts = (display = (ticks tickvalues)) yaxisopts = (display = (ticks tickvalues));
            heatmapparm x = x y = y colorresponse = value / name = 'a';
            continuouslegend 'a';
        endlayout;
    endgraph;
end;
run;

proc sgrender data = GTL_GS_HeatmapParm template = _heatmapparm;
run;