Este es un conjunto de datos simulado que pretende responder a una pregunta sobre si los sujetos con aumento de value
entre trip
(de v0
a v2
, 3 visitas) tendrán un tiempo de supervivencia significativamente más corto/largo ( time
) . El conjunto de datos simulado se muestra a continuación (puede reproducirse en R):
subj = rep(1:50, each=3)
trip = rep(c("v0", "v1", "v2"), 50)
value = round(abs(rnorm(n=150, mean=6, sd=2)), 3)
time = round(rexp(n=150, rate=0.3), 3)
censor = rep(rbinom(n=50, size=1, prob=0.5), each=3)
type = c(rep("type 1", 12), rep("type 2", 6), rep("type 3", 5), rep("type 4", 16), rep("type 5", 11))
dat = data.frame(subj, trip, value, time, censor, type)
Hay un total de 50 temas ( subj
), cada una con 3 visitas ( trip
), y para cada visita, hay una medida value
tomada. Cada sujeto tuvo un tiempo de supervivencia ( time
), el estado de censura ( censor
0=censor, 1=muerte), y el tipo de enfermedad ( type
). Le rogamos que nos indique algún método estadístico bien establecido para este tipo de análisis. Lo ideal es que dichos métodos puedan ampliarse para incorporar covariables adicionales a nivel de sujeto, con la implementación del paquete R disponible. Probablemente el modelo PH de Cox con efecto aleatorio para el tipo de enfermedad ( type
) puede ser una forma de considerar? Gracias por las sugerencias.