小试了一下:
Modeling the cumulative incidence.
Data BMT; Set SASHelp.BMT;Run;
Data Cov;
Length Group $13.;
Group = "ALL"; Output;
Group = "AML-Low Risk"; Output;
Group = "AML-High Risk"; Output;
Run;
Ods Graphics On/ Width = 7In Height = 5In;
Proc PHReg Data = BMT Plots(Overlay) = CIF;
Ods Select CIFPlot;
Class Group(Ref = First);
Model T*Status(0) = Group / Eventcode = 1;
Baseline Covariates = Cov Out = Res CIF = _ALL_ / Seed = 1000;
Run;
Proc Lifetest Data = BMT Plots = Survival(Outside Atrisk(Maxlen = 13 Outside(0.15) ) );
Ods Select SurvivalPlot;
Time T*Status(0);
Strata Group / Order = Internal;
Run;
Ods Graphics Off;