[R-es] Interpretación de coeficientes en un cox proportional hazards con variable strata
Lo que me resulta estraño del modelo que planteas es la combinación de la "interacion" con la "estratificación": S1.nuevo ~ mobilityPDurG2 * strata(Countryb) + strata(Q6), Esta claro que "coxph" se lo ha tragado, pero ni idea de los que puede significar. En principio, yo creo que la estratificación se usan con varaibles "confusores", con efecto causal conocido, pero que no tiene interés directo en el análisis (por ejemplo el sexo en el análisis de las desigualdades sociales en supervivencia) , ya que para cada estrato se estima una función empÃrica basal de "hazars rates", sobre las que hay que aplicar las "hazard ratios" derivadas de los coeficientes del modelo. Por lo tanto si estas interesado en disciriminal el efecto de la movilidad en cada uno de los paises, creo que o debes de estratificar por "Countryb". Para eso el modelo que yo plantearia, seria este: S1.nuevo ~ mobilityPDurG2 * Countryb + strata(Q6) # puede que con -1, para facilitar la interpretación del efecto paÃs A partir de este modelo, para cada paÃs tendremos dos coeficientes, una el efecto del "Countryb" y otro la interacción con "mobility". En principio, y obviando la cuestión de la categorÃa paÃs de referencia, exponenciando el coeficiente de interacción "exp(mobilityPDurG2:Countryb(pais.uno))" se obtendrÃa la "hazar ratio" (Riesgo Relativo) de los individuos con mobilityPDurG2=1 en el "paÃs.uno" El 04/02/15 a las 13:54, José Luis Cañadas Reche escribió:
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
---
_______________________________________________ R-help-es mailing list R-help-es en r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es
+-------------------------------------------------------------- | Francisco J. Viciana Fernández | Coordinador del Registro de Población | Servicio de EstadÃsticas Demográficas y Sociales | Instituto de EstadÃstica y CartografÃa de AndalucÃa | Leonardo Da Vinci, nº 21. Isla de La Cartuja. | 41071 SEVILLA. | franciscoj.viciana en juntadeandalucia.es +--------------------------------------------------------------