Thanks ballardw, That is very helpful. I was able to get onto the right track I think, though I need to alter the min max values to get a better-looking figure (in terms of the color transitions). My SGPLOT code differed slightly, mainly in terms of where the colorresponse= and rattrid= options are located. *Create custom RATTRMAP; data myrattrmap; retain id "myID"; length min $ 6 max $ 6 color altcolor colormodel1 colormodel2 colormodel3 $ 15; input min $ max $ color $ altcolor $ colormodel1 $ colormodel2 $ colormodel3 $ ; datalines; -0.99 -0.29 . . blue lightblue verylightblue -0.289 0.02 white white . . . 0.02 0.50 . . verylightred lightred red ; run; PROC SGPLOT Data=work.mycorr_sexF_2 rattrmap=myrattrmap; Heatmapparm x=x y=y colorresponse=value / rattrid=myID outline nomissingcolor; *coloresponse statement tells SAS that the third variable is numeric value (correlation stat); Text x=x y=y text=value; RUN; Doing this, I was able to generate the following: It does not look great. I think I can change the min and max values and also determine if I need more colormodel statements? I will read up on how to define these colormodels at the following: https://documentation.sas.com/doc/en/pgmsascdc/9.4_3.5/grstatproc/n1fxl54aizx1tjn1nfder5bsqqco.htm & https://documentation.sas.com/doc/en/pgmsascdc/9.4_3.5/grstatproc/p0edl20cvxxmm9n1i9ht3n21eict.htm
... View more