title: How to plot restricted cubic spline with hazard radio, probability of mortality, or mortality rate on y-axis?
body: <p>1) How can I change the y-axis to "odds ratio", "probability of mortality", and "mortality rate" for "fit" in the following example?</p>
  
  <p>2) How can I change the y-axis to "hazard ratio" for "fit2" in the following example?</p>
  
  <pre><code>library(Hmisc)
  library(survival)
  library(rms)
  
  data(pbc)
  d &lt;- pbc
  rm(pbc)
  d$died &lt;- ifelse(d$status == 2, 1, 0)
  d$status &lt;- ifelse(d$status != 0, 1, 0)
  
  ddist &lt;- datadist(d)
  options(datadist='ddist')
  
  fit &lt;- lrm(status ~ rcs(age, 4), data=d)
  (an &lt;- anova(fit))
  plot(Predict(fit), anova=an, pval=TRUE)
  
  fit2 &lt;- cph(Surv(time, status) ~  rcs(age, 4), data=d)
  (an2 &lt;- anova(fit2))
  plot(Predict(fit2), anova=an, pval=TRUE)
  </code></pre>
  
  <p>I am looking forward to your help!</p>
  
  <p><strong>Update 1</strong>
  Following BondedDust's answer I worte the following:</p>
  
  <pre><code># probability
  getProbability &lt;- function(x) {
       exp(x)/(1+exp(x))*100
  }
  
  fit &lt;- lrm(status ~ rcs(age, 4), data=d)
  (an &lt;- anova(fit))
  plot(Predict(fit, fun=getProbability), anova=an, pval=TRUE, ylab="Probability of death [%]")
  
  # overall probability to die
  table(d$status)
  round(table(d$status)[[2]]/sum(table(d$status))*100, digits=1) # = 44.5%
  </code></pre>
  
  <p>Since the overall probability to die is 44.5% the calculation of the predicted probability and the resulting plot seems to be correct to me as non-statistician, isn't it?</p>