天天看点

SAS survival_1

ods trace on;
ods output Survivalplot=GTL_GS_SurvivalPlotData;
ods graphics / reset imagename='Survival-LifeTest';
proc lifetest data=sashelp.BMT plots=survival(atrisk=0 to 2500 by 500);
   time T * Status(0);
   strata Group / test=logrank adjust=sidak;
   run;
ods trace off;

proc format;
  value aml
  3 = 'AML-Low'
  2 = 'AML-High'
  1 = 'All';
run;
 
ods html;

proc template;
define statgraph _survivalplot;
    begingraph;
    title 'title1';
        entrytitle 'title';
        entryfootnote halign = left 'footnote' ;
    layout overlay ;
            entry 'entry';
            stepplot x = time y = survival / group = stratum lineattrs = (pattern = solid) name = "s";
            scatterplot x = time y = censored / group = stratum markerattrs = (symbol = plus) name = "c";
            innermargin;
                blockplot x = tatrisk block = atrisk / class = stratumnum display = (values label)
                                                        valueattrs = (size = 8) labelattrs = (size = 8);
            endinnermargin;
            discretelegend "s" / valueattrs = (size = 8);
    endlayout;
    endgraph;
end;
run;

proc sgrender data = GTL_GS_SurvivalPlotData template = _survivalplot;
format stratumnum aml.;
run;

ods _all_ close;      
SAS survival_1

本文来自博客园,作者:Iving,转载请注明原文链接:https://www.cnblogs.com/SAS-T/p/15369924.html