用STATA做基于RCS的COX回归预测HR值置信区间图要怎么做呢,类似于这种,

用STATA做基于RCS的COX回归预测HR值置信区间图要怎么做呢,类似于这种,谢谢!