[R-es] Interpretación de coeficientes en un cox proportional hazards con variable strata

José Luis Cañadas Reche canadasreche en gmail.com
Mie Feb 4 13:54:04 CET 2015


Buenas.

Abajo pongo la  salida de un modelo de cox , dónde he estratificado por 
una variable de país (Countryb)  y por otra (Q6). Además hay interacción 
entre la variable mobilityPDurG2 (es una variable 0,1, y 0 es la 
categoría de referencia)  país.
  La categoría de referencia para país es "united kingdom".

Mi duda surge si quiero calcular el hazard ratio para los que tienen un 
1 en mobilityPDurG2 en los diferentes países respecto a los que tienen 0.

En principio sería con los exp(coef), de forma que por ejemplo para el 
Reino unido sería 0.78592 y para  Bélgica 0.65489, pero no estoy seguro 
y no he encontrado por ahí como se interpreta cuando tengo interacción 
entre la variable de estrato y una covariable.

Gracias..

summary(mod.cph.strat3)
Call:
coxph(formula = S1.nuevo ~ mobilityPDurG2 * strata(Countryb) +
     strata(Q6), data = datos[filtro, ], method = "breslow")

   n= 5885, number of events= 4397

coef exp(coef) se(coef)      z Pr(>|z|)
mobilityPDurG2 -0.24090   0.78592  0.06930 -3.476 0.000509 ***
mobilityPDurG2:strata(Countryb)Countryb=belgium -0.42328   0.65489  
0.21738 -1.947 0.051514 .
mobilityPDurG2:strata(Countryb)Countryb=france -0.13352   0.87501  
0.13945 -0.957 0.338326
mobilityPDurG2:strata(Countryb)Countryb=germany -0.01344   0.98665  
0.13199 -0.102 0.918871
mobilityPDurG2:strata(Countryb)Countryb=italy 0.19129   1.21081  
0.11428  1.674 0.094149 .
mobilityPDurG2:strata(Countryb)Countryb=netherlands -0.10238   0.90269  
0.13224 -0.774 0.438819
mobilityPDurG2:strata(Countryb)Countryb=poland 0.17324   1.18916  
0.25840  0.670 0.502579
mobilityPDurG2:strata(Countryb)Countryb=spain -0.49122   0.61188  
0.14938 -3.288 0.001008 **
mobilityPDurG2:strata(Countryb)Countryb=sweden -0.28491   0.75208  
0.12522 -2.275 0.022884 *
mobilityPDurG2:strata(Countryb)Countryb=switzerland -0.14876   0.86178  
0.16448 -0.904 0.365771
---



Más información sobre la lista de distribución R-help-es