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!