Tengo una pregunta metodológica, y no tengo ningún conjunto de datos.
Supongamos que quiero ajustar un modelo de riesgos proporcionales (Cox) para los datos de supervivencia. Tengo múltiples observaciones para cada individuo (datos en formato largo). El interés particular reside en un predictor continuo (como los niveles de lípidos en sangre) y estoy examinando la asociación/efecto sobre el riesgo de infarto de miocardio.
Soy un novato en esto, pero utilizo la regresión de Cox dependiente del tiempo (sin agrupación ya que sólo se analiza un evento [los eventos repetidos no son de interés]). Me parece que este es el método estándar.
Ahora, los paquetes JMbayes, JM, joinR, lcmm pueden ajustar modelos conjuntos ( http://www.r-bloggers.com/joint-models-for-longitudinal-and-survival-data/ ) que parecen ser una fusión entre los modelos de efectos mixtos y la regresión de Cox.
Esto parece una buena idea, para combinar con métodos robustos... Un par de ventajas se reportan para los modelos conjuntos, de los cuales la "precisión en la trayectoria de cada paciente" se menciona repetidamente. Sin embargo, he buscado en pubmed, google scholar y google las publicaciones que utilizan este enfoque y no he encontrado mucho.
¿Debería seguir con la "habitual" regresión de Cox dependiente del tiempo (proceso de recuento)? ¿Ventajas? ¿Inconvenientes?