Just going to make a plug for my survival graphing macro that I've shared on here:  https://communities.sas.com/t5/SAS-Communities-Library/Kaplan-Meier-Survival-Plotting-Macro-NEWSURV/ta-p/479747   Maybe you'll find it helpful.   It has all sorts of features for survival graphs.  I think you could make your graph with the following macro call: %newsurv(data=file, time=os_time_km, cens=vitalstatusos_km, cen_vl=1, summary=0, display=, risklist=0 to 22.5 by 2.5, risklocation=inside, ylabel=Overall Survival Probability (%), yincrement=25, classdesc=At Risk, parheader=, riskdivider=0, height=550px, linesize=2.5, censormarkers=0)     ; ... View more

I did this by making a format for the p-values and applying it.  I added an ENTRYFOOTNOTE for a footnote clarifying the symbols, and you can change the text to be what you'd like.   data work.Almtarfi_2; input Cult Biomass_BRC Biomass_RB yield_BRC Yield_RB Height_BRC Heigh_RB; CARDS; 1930 62.125 75.255 2563.16 2196.06 123.075 142.8 1933 60.62 77.91 2396.49 1965.79 106.233 124.325 1935 51.275 77.2 2317.14 2490.56 114.75 137.25 1940 63.31 73.025 2738.19 2888.38 142.85 169.5 ; run; proc format; value pvaldots 0-<0.0001 = '***' 0.0001-<0.001 = '**' 0.001-<0.05 = '*' 0.05-high =' ' other = ' '; run; * for P <0.05 and and ** for P <0.001 , *** for P <0.0001, and nothing when P > 0.05; /*Code below by Chris Hemedinger, found on the SAS blogs*/ ods trace off; /* Prepare the correlations coeff matrix: Pearson's r method */ %macro prepCorrData(in=,out=); /* Run corr matrix for input data, all numeric vars */ ods select none; ods noresults; ods output PearsonCorr=blah; proc corr data=&in. pearson outp=work._tmpCorr vardef=df ; run; ods results; ods select all; proc sql noprint; select variable into :vlist separated by '|' from blah; quit; /* prep data for heat map */ data &out.; keep x y r p; format p pvaldots.; set blah (rename=(variable=_name_)); *set work._tmpCorr(where=(_TYPE_="CORR")); array v{*} %do i = 1 %to %sysfunc(countw(%superq(vlist),|)); %scan(%superq(vlist),&i,|) %end;; array pval{*} %do i = 1 %to %sysfunc(countw(%superq(vlist),|)); p%scan(%superq(vlist),&i,|) %end;; x = _NAME_; do i = dim(v) to 1 by -1; y = vname(v(i)); r = v(i); p = pval(i); /* creates a lower triangular matrix */ if (i<_n_) then do; r=.;p=.; end; output; end; run; proc datasets lib=work nolist nowarn; delete _tmpcorr; quit; %mend; /* Create a heat map implementation of a correlation matrix */ ods path work.mystore(update) sashelp.tmplmst(read); proc template; define statgraph corrHeatmap; dynamic _Title; begingraph; entrytitle _Title; entryfootnote halign=left '*: P <0.05;**: P <0.001;***: P <0.0001;Blank: P > 0.05' / textattrs=(size=10pt); rangeattrmap name='map'; /* select a series of colors that represent a "diverging" */ /* range of values: stronger on the ends, weaker in middle */ /* Get ideas from http://colorbrewer.org */ range -1 - 1 / rangecolormodel=(cxD8B365 cxF5F5F5 cx5AB4AC); endrangeattrmap; rangeattrvar var=r attrvar=r attrmap='map'; layout overlay / xaxisopts=(display=(line ticks tickvalues)) yaxisopts=(display=(line ticks tickvalues)); heatmapparm x = x y = y colorresponse = r / xbinaxis=false ybinaxis=false name = "heatmap" display=all; textplot x=x y=y text=p / textattrs=(size=14pt); continuouslegend "heatmap" / orient = vertical location = outside title="Pearson Correlation"; endlayout; endgraph; end; run; /* Build the graphs */ ods graphics /height=600 width=800 imagemap; %prepCorrData(in=work.Almtarfi_2,out=work.Almtarfi_r); proc sgrender data=work.Almtarfi_r template=corrHeatmap; dynamic _title="Correlation matrix for all variables"; run; ... View more

OK I ran what you had and got the same errors.  I fixed it by changing the order before the proc corr.  Put the ODS OUTPUT after the ODS NORESULTS statement. ods select none; ods noresults; ods output PearsonCorr=blah; ... View more

Can you copy and paste the exact code you're running now, just so I can check if you're running the same thing that I posted? ... View more

It looks like the actual PROC CORR isn't running. Try deleting the comments around the procedure to make sure the procedure isn't accidentally being commented out. ... View more

What happens when you run the PROC CORR with the ODS OUTPUT statement that I put before it? And then the PROC SQL code that I had listed? ... View more

Here is my solution.  I used ODS OUTPUT to get the p-values and TEXTPLOT to graph them.   data work.Almtarfi_2; input Cult Biomass_BRC Biomass_RB yield_BRC Yield_RB Height_BRC Heigh_RB; CARDS; 1930 62.125 75.255 2563.16 2196.06 123.075 142.8 1933 60.62 77.91 2396.49 1965.79 106.233 124.325 1935 51.275 77.2 2317.14 2490.56 114.75 137.25 1940 63.31 73.025 2738.19 2888.38 142.85 169.5 /*Code below by Chris Hemedinger, found on the SAS blogs*/ ods trace off; /* Prepare the correlations coeff matrix: Pearson's r method */ %macro prepCorrData(in=,out=); /* Run corr matrix for input data, all numeric vars */ ods select none; ods noresults; ods output PearsonCorr=blah; proc corr data=&in. pearson outp=work._tmpCorr vardef=df ; run; ods results; ods select all; proc sql noprint; select variable into :vlist separated by '|' from blah; quit; /* prep data for heat map */ data &out.; keep x y r p; set blah (rename=(variable=_name_)); *set work._tmpCorr(where=(_TYPE_="CORR")); array v{*} %do i = 1 %to %sysfunc(countw(%superq(vlist),|)); %scan(%superq(vlist),&i,|) %end;; array pval{*} %do i = 1 %to %sysfunc(countw(%superq(vlist),|)); p%scan(%superq(vlist),&i,|) %end;; x = _NAME_; do i = dim(v) to 1 by -1; y = vname(v(i)); r = v(i); p = pval(i); /* creates a lower triangular matrix */ if (i<_n_) then do; r=.;p=.; end; output; end; run; proc datasets lib=work nolist nowarn; delete _tmpcorr; quit; %mend; options mprint; /* Create a heat map implementation of a correlation matrix */ ods path work.mystore(update) sashelp.tmplmst(read); proc template; define statgraph corrHeatmap; dynamic _Title; begingraph; entrytitle _Title; rangeattrmap name='map'; /* select a series of colors that represent a "diverging" */ /* range of values: stronger on the ends, weaker in middle */ /* Get ideas from http://colorbrewer.org */ range -1 - 1 / rangecolormodel=(cxD8B365 cxF5F5F5 cx5AB4AC); endrangeattrmap; rangeattrvar var=r attrvar=r attrmap='map'; layout overlay / xaxisopts=(display=(line ticks tickvalues)) yaxisopts=(display=(line ticks tickvalues)); heatmapparm x = x y = y colorresponse = r / xbinaxis=false ybinaxis=false name = "heatmap" display=all; textplot x=x y=y text=p; continuouslegend "heatmap" / orient = vertical location = outside title="Pearson Correlation"; endlayout; endgraph; end; run; /* Build the graphs */ ods graphics /height=600 width=800 imagemap; %prepCorrData(in=work.Almtarfi_2,out=work.Almtarfi_r); proc sgrender data=work.Almtarfi_r template=corrHeatmap; dynamic _title="Correlation matrix for all variables"; run;     ... View more

Hello @ballardw, The 'A0'x is an unbreakable space. This is similar to a space in that it doesn't show anything (it's a blank space), but SAS considers it a character. In situations where a regular space would be ignored by SAS procedures an unbreakable space will be treated like any other character ('A', 'W', anything). In this case, when the label is set to be a blank (e.g. label=' ') I believe SAS ignores it and prints the x variable's label. But with an unbreakable space it leaves it alone so you can add a blank label and still allocate the graph space for it. ... View more