Code below is the example from SAS GUIDE. I am working on a competing risk survival analysis.
But this example only has one covariate affecting the event. But my model has more than one categorical variables. How can I create dataset ("risk" in my example) in the "baseline covariates =" syntax to plot CIF plot according to the categorical variable "class5"?
All the variables are categorical except Age. Thanks very much!
/*My code*/
proc phreg data=pat plots(overlay) = cif;
class class5(order = internal ref=first) bmi race drink smoke cvd;
model time*status(0) = class5 age bmi race drink smoke cvd / risklimits = WALD eventcode=1;
baseline covariates = risk / rowid = class5;
run;
/*Example from SAS guide*/
proc format;
value DiseaseGroup 1='ALL'
2='AML-Low Risk'
3='AML-High Risk';
data Bmt;
input Disease T Status @@;
label T='Disease-Free Survival in Days';
format Disease DiseaseGroup.;
datalines;
1 2081 0 1 1602 0 1 1496 0 1 1462 0 1 1433 0
1 1377 0 1 1330 0 1 996 0 1 226 0 1 1199 0
1 1111 0 1 530 0 1 1182 0 1 1167 0 1 418 2
1 383 1 1 276 2 1 104 1 1 609 1 1 172 2
1 487 2 1 662 1 1 194 2 1 230 1 1 526 2
1 122 2 1 129 1 1 74 1 1 122 1 1 86 2
1 466 2 1 192 1 1 109 1 1 55 1 1 1 2
1 107 2 1 110 1 1 332 2 2 2569 0 2 2506 0
2 2409 0 2 2218 0 2 1857 0 2 1829 0 2 1562 0
2 1470 0 2 1363 0 2 1030 0 2 860 0 2 1258 0
2 2246 0 2 1870 0 2 1799 0 2 1709 0 2 1674 0
2 1568 0 2 1527 0 2 1324 0 2 957 0 2 932 0
2 847 0 2 848 0 2 1850 0 2 1843 0 2 1535 0
2 1447 0 2 1384 0 2 414 2 2 2204 2 2 1063 2
2 481 2 2 105 2 2 641 2 2 390 2 2 288 2
2 421 1 2 79 2 2 748 1 2 486 1 2 48 2
2 272 1 2 1074 2 2 381 1 2 10 2 2 53 2
2 80 2 2 35 2 2 248 1 2 704 2 2 211 1
2 219 1 2 606 1 3 2640 0 3 2430 0 3 2252 0
3 2140 0 3 2133 0 3 1238 0 3 1631 0 3 2024 0
3 1345 0 3 1136 0 3 845 0 3 422 1 3 162 2
3 84 1 3 100 1 3 2 2 3 47 1 3 242 1
3 456 1 3 268 1 3 318 2 3 32 1 3 467 1
3 47 1 3 390 1 3 183 2 3 105 2 3 115 1
3 164 2 3 93 1 3 120 1 3 80 2 3 677 2
3 64 1 3 168 2 3 74 2 3 16 2 3 157 1
3 625 1 3 48 1 3 273 1 3 63 2 3 76 1
3 113 1 3 363 2
;
data Risk;
Disease=1; output;
Disease=2; output;
Disease=3; output;
format Disease DiseaseGroup.;
run;
ods graphics on;
proc phreg data=Bmt plots(overlay=stratum)=cif;
class Disease (order=internal ref=first);
model T*Status(0)=Disease / eventcode=1;
Hazardratio 'Pairwise' Disease / diff=pairwise;
baseline covariates=Risk out=out1 cif=_all_ / seed=191;
run;
ods graphics on;
title 'Cause-specific Analysis';
proc phreg data=Bmt plots(overlay)=cif;
class Disease (order=internal ref=first);
model T*Status(0)=Disease / eventcode(cox)=1;
baseline covariates=Risk out=out2 cif=_all_;
run;
