[R] Trend test for hazard ratios

[mkleber74]

Hi Terry,

thank you for your suggestion. I tried to do it the way that you described.
I fitted my Cox model with strata and I get coefficients and P values for my
ZlnGalectin in the different strata (see summary below). So am not so
familiar with R so please excuse my perhaps stupid questions but are the P
values below like "p=2.386e-06" something similar like my P for trend (they
seem to be very low?) or if not how can I calculate it now? 


> summary(fit_mort)
Call:
coxph(formula = Surv(obstime2010, death2010) ~ ZlnGalectin3 * 
    strata(CKD_KDIGO), data = gal.df)

  n= 2573, number of events= 758 
   (743 observations deleted due to missingness)

                                                                                           
coef       exp(coef)   se(coef)      z    Pr(>|z|)
ZlnGalectin3                                                                         
0.20725   1.23029  0.07064  2.934   0.00335 **
ZlnGalectin3:strata(CKD_KDIGO)CKD_KDIGO=60,01 - 90,00  0.04165   1.04253 
0.09178  0.454   0.65000 
ZlnGalectin3:strata(CKD_KDIGO)CKD_KDIGO=90,01+           -0.06830   0.93398 
0.12113 -0.564  0.57282 
                                                                              
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

                                                                                         
exp(coef) exp(-coef) lower .95  upper .95
ZlnGalectin3                                                                            
1.230     0.8128    1.0712   1.413
ZlnGalectin3:strata(CKD_KDIGO)CKD_KDIGO=60,01 - 90,00     1.043     0.9592   
0.8709   1.248
ZlnGalectin3:strata(CKD_KDIGO)CKD_KDIGO=90,01+              0.934     1.0707   
0.7366   1.184
                                                                      
Concordance= 0.559  (se = 0.017 )
Rsquare= 0.011   (max possible= 0.978 )
Likelihood ratio test = 28.87  on 3 df,   p=2.386e-06
Wald test               = 28.64  on 3 df,   p=2.665e-06
Score (logrank) test = 28.46  on 3 df,   p=2.906e-06



--
View this message in context: http://r.789695.n4.nabble.com/Trend-test-for-hazard-ratios-tp4688729p4688837.html
Sent from the R help mailing list archive at Nabble.com.