How to graph HR derived from the restricted cubic spline model per procedure type

[Victor Dayan]

Hi, I am a cardiac surgeon struggling on statistics and quite new to STATA. I am evaluating the survival after aortic valve replacement with mechanical or bio prosthetic valves. It is well known that the younger patient has improved survival if a mechanical prosthesis is used, while an elder patient benefits from a bio prosthesis. In order to define the age cut-off, I plan to plot the HR and its 95%CI derived from Cox regression (y axis) and age (x-axis). Previous studies have explained this procedure as "To explore the age-dependent effect of prosthesis type on mortality, a Cox proportional-hazards model was fit to the entire weighted study population with the use of an interaction term for age (modeled as a natural spline) and prosthesis type."
I am attaching a Figure as a sample of exactly what I intend to obtain with my data.
I appreciate some help into how to proceed.
Thanks

[Marco Tomasino]

I am looking for an answer to this question as well

[Hemanshu Kumar]

You might want to take a look at

help stcox
help stcox postestimation