Dados los tiempos de supervivencia censurados por intervalos, ¿cómo realizo un modelo Cox PH censurado por intervalos R? Una búsqueda rseek muestra el paquete intcox, que ya no existe en el Rrepositorio. Estoy casi seguro de que la coxphfunción en el survivalpaquete no puede manejar datos de supervivencia censurados por intervalos.
Además, no quiero imputar los datos y luego usar la coxphfunción. Este método subestima los errores estándar de los coeficientes porque ignora la incertidumbre de la censura por intervalos.
r
survival
cox-model
interval-censoring
wcampbell
fuente
fuente
intcoxpaquete incluso si no estáCRANusando el normalinstall.packages("intcox").install.packages("intcox")sin ningún problema en particular (R-desarrollo, pero cualquier R moderna debería funcionar)Respuestas:
Como se indicó anteriormente, puede usar la función survreg. Sin embargo, una nota: este no es estrictamente un modelo Cox PH, sino modelos a escala de ubicación. Usando la transformación de registro predeterminada, este es el modelo de popa. En el caso de la distribución exponencial, los riesgos proporcionales y el modelo en popa son equivalentes, por lo que si la distribución se establece en exponencial, este es un modelo de riesgos proporcionales con una línea de base exponencial. Del mismo modo, si se utiliza un modelo de línea de base de distribución de Weibull en popa, las estimaciones de los parámetros son solo una transformación lineal de las utilizadas en el modelo de riesgos proporcionales con la distribución de línea de base de Weibull. Pero en general, survreg no se ajusta a un modelo Cox PH.
Si se desea un modelo semiparamétrico, como se encuentra implementado en intcox, una palabra de precaución: hay varios problemas con la versión actual de intcox (el algoritmo generalmente termina prematuramente de manera significativa lejos del MLE, falla por completo con observaciones sin censura, sin errores estándar) presentado automáticamente).
Una nueva alternativa que podría usar es el paquete "icenReg".
Admisión de sesgo: este es el autor de icenReg.
fuente
ic_spnecesita estimar la distribución de supervivencia de referencia (a diferencia del caso censurado correcto), que tiene tantos parámetros como tiempos únicos en sus datos. Esto crea un problema para los mini lotes; con tiempos continuos, los pasos de la línea de base no se alinearán lote a lote.Para hacer un análisis censurado por intervalos en R, debe crear un objeto Surv y luego usar survfit (). Si tiene más de una variable, el paquete intcox resuelve el problema.
fuente