The dataset is presented below,

insheet using "https://ift.tt/aXA9elG", comma
foreach i of varlist familyhx smokinghx sex cancerstage school {
encode `i', gen(`i'2)
drop `i'
rename `i'2 `i'
}
ssc inst center,all replace
center co2 il6 crp lengthofstay

I want to fit a survival model with the interaction term between the quadratic time variable (canerstage) and a continuous variable (co2) as follows,

xtmelogit remission c.c_il6 c.c_crp c.c_lengthofstay i.sex i.sex#c.cancerstage c.cancerstage c.c_co2 c.cancerstage#c.c_co2 c.cancerstage#c.cancerstage#c.c_co2 c.cancerstage#c.cancerstage#i.sex|| did:, intpoints(10) or

Now, I want to graph the estimated hazard curve showing the effect of co2 in the manner (Mean-1*SD,Mean, Mean+1*SD) while holding all other variables at their mean in Stata.

Thanks for your code or any hint!