Survival curves adjusted for the covariates from a Cox proportional hazard model

[Nachiketh Soodana Prakash]

I am working on a survival data on cancer patients comparing two groups of treatments. I build the model using the step by step approach of making a univariate log rank test for categorical variables

sts test variable, logrank
sts graph, by(variable)

and univariate Cox regression for continuous model

 stcox variable, nohr

When building my multivariate model from this step i wanted create KM curves for Treatment group: Chemo and Radiation

sts graph, by(treatment_group)

I wanted to adjust the survivor curves for other co-variates, which i achieved by the following

sts graph, by(treatment_group) adjustfor(varlist)

Recently i found that SAS has method of producing survival curve from a Cox Proportional Hazard Model.

You can find the article below.

http://www2.sas.com/proceedings/sugi23/Stats/p226.pdf

Is there a way to do this in STATA.

Thank you.

[Stephen Jenkins]

"adjust" is not a term I am hugely familiar with (I am an applied economist). However I note that help sts graph states:

adjustfor(varlist) adjust the estimates to zero values of varlist

Might stcurve be what you're looking for, exploiting the survivor and at() options?