I'm working on some stcox analyses where I want to test different thresholds of exposure and I'd like to show the survival graphs all on one plot.
My syntax looks like this for the stcox, where threshold30 is a binary variable:
svyset [pweight=weight]
stset deathdate [pw=weight], id(idauniq) failure(Entry_death==1) origin(interviewdate) scale(365.25)
stcox threshold30 $covariates
stcurve, survival at1(threshold30=0) at2(threshold30=1) xlabel(0(2)14) ylabel(0.7(0.05)1) name(Survival_Cox30)
stcox threshold40 $covariates
stcurve, survival at1(threshold40=0) at2(threshold40=1) xlabel(0(2)14) ylabel(0.7(0.05)1) name(Survival_Cox40)
stcox threshold50 $covariates
stcurve, survival at1(threshold50=0) at2(threshold50=1) xlabel(0(2)14) ylabel(0.7(0.05)1) name(Survival_Cox50)
I'd be grateful for advice on how to combine these stcurve graphs onto one graph.
Best wishes,
Daisy
My syntax looks like this for the stcox, where threshold30 is a binary variable:
svyset [pweight=weight]
stset deathdate [pw=weight], id(idauniq) failure(Entry_death==1) origin(interviewdate) scale(365.25)
stcox threshold30 $covariates
stcurve, survival at1(threshold30=0) at2(threshold30=1) xlabel(0(2)14) ylabel(0.7(0.05)1) name(Survival_Cox30)
stcox threshold40 $covariates
stcurve, survival at1(threshold40=0) at2(threshold40=1) xlabel(0(2)14) ylabel(0.7(0.05)1) name(Survival_Cox40)
stcox threshold50 $covariates
stcurve, survival at1(threshold50=0) at2(threshold50=1) xlabel(0(2)14) ylabel(0.7(0.05)1) name(Survival_Cox50)
I'd be grateful for advice on how to combine these stcurve graphs onto one graph.
Best wishes,
Daisy
