I am attempting to fit a Cox proportional hazard model to my data. I think I have the formula correct but am having trouble understanding the output. I have tried looked through the documentation and it is hard for me to understand. Any help would be greatly appreciated, thank you!
coxfit1 <- coxph(Surv(days, status)~GENE1, data=dataset1) summary(coxfit1)
Where "days" is days until an event occurred (or last known followup if no event), "status" is an event (recurrence), GENE1 is expression data of a gene that I am testing if it has an effect on recurrence.
Call: coxph(formula = Surv(days, status) ~ GENE1, data = dataset1) n= 34, number of events= 22 coef exp(coef) se(coef) z Pr(>|z|) GENE1 0.6370 1.8908 0.2362 2.697 0.00699 ** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 exp(coef) exp(-coef) lower .95 upper .95 GENE1 1.891 0.5289 1.19 3.004 Concordance= 0.618 (se = 0.068 ) Rsquare= 0.166 (max possible= 0.98 ) Likelihood ratio test= 6.17 on 1 df, p=0.01298 Wald test = 7.27 on 1 df, p=0.006993 Score (logrank) test = 7.81 on 1 df, p=0.005198
Now, this is one that is obviously significant, but what do the different parts of this output mean? Where is the hazard ratio??? And which of this information is appropriate for reporting?