[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