topic Re: LIFETEST CIF in SAS Software for Learning Community

LIFETEST CIF

lilies — Mon, 17 Oct 2022 01:12:20 GMT

I'm attempting to customize the CIF graph from the PROC LIFETEST statement e.g. changing the title of the graph, the thickness of the graph lines, whether the graph lines are dashed, etc. However, I'm unable to find an easy way to do so. Does anyone have any recommendations? I considered using the %CIF macro, but I also would like to display the Gray's test box on the graph, which I don't think I can do if I use the %CIF macro.

This is the code I'm using to generate the CIF graph with competing risk:

proc lifetest data=WORK.IMPORT plots=cif(test);

time Time*Status(2) / eventcode=1;

strata Variable;

run;

My question is similar to this one, which was not answered:

https://communities.sas.com/t5/Statistical-Procedures/Proc-lifetest-for-competing-risk/m-p/689519#M33245

Re: LIFETEST CIF

Ksharp — Mon, 17 Oct 2022 11:49:16 GMT

ods select none;
ods output cifPlot=cifPlot GrayTest=GrayTest;
proc lifetest data=sashelp.bmt plots=cif(test);
      time t*Status(0) / eventcode=1;
     strata Group; 
run;
ods select all;
data _null_;
 set GrayTest;
 call symputx('pvalue',put(ProbChiSq,pvalue. -l));
run;

title 'Whatever is Title';
proc sgplot data=cifPlot;
styleattrs DATALINEPATTERNS=(solid dash shortdash);
step x=Time y=CIF / group=Stratum lineattrs=( thickness=4) ;
yaxis values=(0 to 1 by 0.2);
inset "Gray Test p=&pvalue." /border;
run;

Re: LIFETEST CIF

lilies — Mon, 17 Oct 2022 14:13:46 GMT

Thank you so much @Ksharp! Is there a way to put an at risk table at the bottom of the graph like in survival plots?

Re: LIFETEST CIF

Ksharp — Tue, 18 Oct 2022 11:38:26 GMT

OK. You want this ?

ods select none;
ods output cifPlot=cifPlot GrayTest=GrayTest ;
proc lifetest data=sashelp.bmt plots=cif(test );
      time t*Status(0) / eventcode=1;
     strata Group; 
run;

ods output SurvivalPlot=SurvivalPlot;
proc lifetest data=sashelp.bmt  plots=survival(atrisk=0 to 2500 by 500);
      time t*Status(0) ;
     strata Group; 
run;
ods select all;


data _null_;
 set GrayTest;
 call symputx('pvalue',put(ProbChiSq,pvalue. -l));
run;

data want;
 merge cifPlot SurvivalPlot(keep=StratumNum Time tAtRisk AtRisk where=(tAtRisk is not missing)) ;
 by StratumNum Time;
length new_Stratum $ 80;
retain new_Stratum;
if first.StratumNum then call missing(new_Stratum);
if not missing(Stratum) then  new_Stratum=Stratum;
run;

title 'Whatever is Title';
ods graphics /ATTRPRIORITY=none ;
proc sgplot data=want;
styleattrs DATALINEPATTERNS=(solid dash shortdash);
step x=Time y=CIF / group=new_Stratum lineattrs=( thickness=4) ;
xaxistable atrisk / x=tAtRisk class=new_Stratum colorgroup=new_Stratum valueattrs=(weight=bold);
yaxis values=(0 to 1 by 0.2);
inset "Gray Test p=&pvalue." /border;
run;

Re: LIFETEST CIF

Gene24 — Fri, 04 Nov 2022 20:04:36 GMT

In a related question, I'm having difficulty editing the x and y scale and the graph title for my kaplan meier curve. This is my code here - what am I missing?

proc lifetest data=WORK.IMPORT2 atrisk plots=survival(cb=hw atrisk(outside(0.15))=0 to 1825 by 365);
time 'Time to RN'n*'Re_RN (1=RN)'n(0);
strata 'Re_ICPI_Status (2=daul; 1=single'n;
run;

title 'Kaplan meier';
proc sgplot data=WORK.IMPORT2 ;
xaxis values=(0 to 1825 by 365);
run;

Re: LIFETEST CIF

Ksharp — Sat, 05 Nov 2022 10:05:09 GMT

I didn't get any problem yet !

ods select none;
ods output SurvivalPlot=SurvivalPlot;
proc lifetest data=sashelp.bmt  plots=survival(atrisk=0 to 2500 by 500);
      time t*Status(0) ;
     strata Group; 
run;
ods select all;

title 'Whatever is Title';
ods graphics /ATTRPRIORITY=none ;
proc sgplot data=SurvivalPlot;
styleattrs DATALINEPATTERNS=(solid dash shortdash);
step x=Time y=Survival / group=Stratum lineattrs=( thickness=4) ;
yaxis values=(0 to 1 by 0.1);
xaxis values=(0 to 2500 by 200);

run;

Re: LIFETEST CIF

Ksharp — Sat, 05 Nov 2022 10:10:27 GMT

Also check Freelancer's post at
https://communities.sas.com/t5/Graphics-Programming/SGPLOT-VBAR-XAXIS-Label-Interval/m-p/842542#M23242

Re: LIFETEST CIF

Gene24 — Thu, 10 Nov 2022 18:55:37 GMT

Thank you!