Sugerencia de modelo para una regresión de Cox con covariables dependientes del tiempo

Estoy modelando el efecto del embarazo sobre el resultado de una enfermedad (vivo-muerto). Aproximadamente el 40% de las pacientes quedaron embarazadas después del momento del diagnóstico, pero en diferentes momentos. Hasta ahora, he realizado gráficos de KM que muestran un claro efecto protector del embarazo sobre la supervivencia y también un modelo regular de Cox; sin embargo, estos se han modelado utilizando solo una variable de embarazo dicotomizada y suponiendo que el efecto está presente desde el momento del diagnóstico, lo cual es claramente poco realista. ya que el tiempo medio de embarazo es de 4 años desde el diagnóstico.

¿Qué tipo de modelo absorbería el efecto de embarazos múltiples en diferentes momentos después del diagnóstico? ¿Sería correcto modelar los embarazos interactuando con el tiempo (lo que requeriría una reconstrucción de datos seria, algún software automatizado que pudiera ayudar con esto?) O existe otra estrategia de modelado preferida para estos problemas? Además, ¿cuál es la estrategia de trazado preferida para estos problemas?

Lo que necesita aquí es una covariable variable en el tiempo y no necesariamente un coeficiente variable en el tiempo . Un ejemplo conocido que podría ayudarlo con sus análisis son los datos de trasplante de corazón de Stanford .

Para presentar sus resultados, puede usar el estimador clásico de Kaplan-Meier que maneja covariables que varían en el tiempo sin problemas (recuerde, sin embargo, que este es un análisis crudo o no ajustado con todas sus limitaciones bien conocidas).

Como ejemplo, el siguiente gráfico muestra el análisis de los datos de HT de Stanford cuando se tiene en cuenta correctamente el estado del trasplante que varía en el tiempo (panel superior) y sin tenerlo en cuenta (panel inferior).

En R, esto puede abordarse con la versión de inicio / detención de un objeto Survival, p. Ej.

fit <- coxph(Surv(time1, time2, status) ~ is.pregnant + other.covariates, data=mydata)

Este documento trata esto con más detalle: http://cran.r-project.org/web/packages/survival/vignettes/timedep.pdf

Cuidado con el sesgo del tiempo inmortal en esta situación. Su grupo de embarazadas inevitablemente tendrá una mejor supervivencia que el grupo de no embarazadas, ya que no puede quedar embarazada después de su muerte (¡que yo sepa!)