Hello all, I did survival analysis for my data and used log rank to test if there is significant difference between the high-risk and low-risk group. This my script and image:

k <- ggsurvplot(survfit(recsurv~T1_signature, data = MSKCC),
           data = MSKCC,
           risk.table = TRUE,
           pval = TRUE,
           conf.int = TRUE,
           ggtheme = theme_minimal(),
           risk.table.y.text.col = TRUE,
           risk.table.y.text = FALSE)

I want to add the HR in the graph as well. I did cox proportional hazards in R. This is the script:

hr <- coxph(recsurv~T1_signature, data = MSKCC)
summary(hr)

and this is the result:

Call:
coxph(formula = recsurv ~ T1_signature, data = MSKCC)

  n= 159, number of events= 6 

                    coef exp(coef) se(coef)      z Pr(>|z|)    
T1_signatureLow -2.89660   0.05521  0.87006 -3.329 0.000871 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

                exp(coef) exp(-coef) lower .95 upper .95
T1_signatureLow   0.05521      18.11   0.01003    0.3038

Concordance= 0.794  (se = 0.09 )
Likelihood ratio test= 11.17  on 1 df,   p=8e-04
Wald test            = 11.08  on 1 df,   p=9e-04
Score (logrank) test = 21.08  on 1 df,   p=4e-06

As far as I found out, for HR I need to subtract the amount of exp(coef) from 1 (1 – 0.05). In this case HR will be 0.95.

It means that the hazar ratio in my low-risk group decreased by 95%. Am I right? I want to use annotate to add HR in my plot. I am not sure should I report it as 0.05 0r 0.95?

I really appreciate any help!