Hi there, I'm trying to fit a Fine and Gray competing risks model and then estimate the cumulative incidence at specific time points (eg. 2 years follow-up). Currently, I'm using the proc phreg function with the baseline options, but this will output all the estimate cumulative incidence at each failure time. I've tried using the timelist option but that can only works in bayesian Cox models. My dataset is relatively big so even when running this model with one covariate, it is taking hours and hours. Ideally, I'd like to had about 8 variables in the model so currently this doesn't seem like a viable option for estimating cumulative incidence. However, producing the graph of the cumulative incidence doesn't take that long, surely there is some way to get the estimate around a certain time point rather than at all failure times. Does anyone know whether you can specifiy time points to estimate cumulative incidence outputs in the proc phreg command? ods graphics on; proc phreg data=test plots(cl overlay stratum range=0,10)=cif; class male2 (order=internal ref=first); model time_yr*outcome(0)=male2 / eventcode=1; baseline covariatesrisk out=risk2 cif=_all_; run; ods graphics off;
... View more