Shape of HR after an intermediate event based on time to that event, multi-state model

[albertoca]

Look also here

Hazard ratio plots with non-linear & time-varying effects in R

Good evening, I wanted to ask if someone could explain how to make in R a plot to represent a time-dependent hazard ratio with 95% CI, from a Cox model, as for example, the one I show below. I have tried hazard.ratio plot in rms or plot(cox.zph) but it is not exactly the same. Thank you very much. [image]