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

[Francisco Viciana]

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
+--------------------------------------------------------------