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

Thanks for the contribution.
There is also an example code by Chris Jackson in this tutorial (example 2.1):