**180**wrote:

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!

**10**• written 13 months ago by Rahil •

**180**