[R-es] transformación a formato largo en modelos multietapa

Hola,

Estoy intentando realizar una estimación de riesgos de evento mediante un
modelo de cox multietapa. El trabajo es un modelo básico de cancer: estado
inicial --> evento intermedio --> evento final (normalmente la muerte).

Para ello estoy usando el paquete mstate de Hein Putter. Así, he definido
mi matriz de transición y mis covariables de interés, pero al intentar
poner los datos en formato largo (una fila = un estado a riesgo de una
transición, en lugar de una fila un paciente) usando el comando msprep, me
escupe un error , a saber :

Error in hlpsrt[, 1] - hlpsrt[, 2] :
  argumento no-numérico para operador binario

Mi código para esta parte es:
formato_largo<-msprep(time = c(*NA*, "tiempo_evento_intermedio",
"tiempo_evento_final"),
              status = c(*NA*, "evento_intermedio", "evento_final"),
              data = mydata,
              trans = tmat,
              keep = covs)

Por más que modifico las variables no consigo salir del error. Si alguien
tiene una sugerencia o me puede ayudar, estaré enormemente agradecido.

Un saludo y gracias de antemano

Jose

Jose Miguel Torregrosa Díaz
Hématoloie Oncologique et Thérapie Cellulaire
CHU de Poitiers, France




************************* AVISO LEGAL *************************
Este mensaje electrónico está dirigido exclusivamente a sus
destinatarios, pudiendo contener documentos anexos de carácter privado
y confidencial.
Si por error, ha recibido este mensaje y no se encuentra entre los
destinatarios, por favor, no use, informe, distribuya, imprima o copie
su contenido por ningún medio. Le ruego lo comunique al remitente y
borre completamente el mensaje y sus anexos.

	[[alternative HTML version deleted]]