Articulo de referencia

Modelo de tiempo de falla acelerado

En el ámbito estadístico del análisis de supervivencia , un modelo de tiempo de falla acelerado ( modelo AFT ) es un modelo paramétrico que ofrece una alternativa a los modelos ...

En el ámbito estadístico del análisis de supervivencia , un modelo de tiempo de falla acelerado ( modelo AFT ) es un modelo paramétrico que ofrece una alternativa a los modelos de riesgos proporcionales comúnmente utilizados . Mientras que un modelo de riesgos proporcionales supone que el efecto de una covariable es multiplicar el riesgo por una constante, un modelo AFT supone que el efecto de una covariable es acelerar o desacelerar el curso de la enfermedad por una constante. Existen sólidas evidencias científicas básicas, provenientes de experimentos con C. elegans realizados por Stroustrup et al. [ 1 ], que indican que los modelos AFT son el modelo correcto para los procesos de supervivencia biológica.

Especificación del modelo

En términos generales, el modelo de tiempo de falla acelerado se puede especificar como [ 2 ].

λ(t|θ)=θλ0(θt){\displaystyle \lambda (t|\theta )=\theta \lambda _{0}(\theta t)}

dóndeθ{\displaystyle \theta }denota el efecto conjunto de las covariables, típicamenteθ=exp([β1incógnita1++βpagincógnitapag]){\displaystyle \theta =\exp(-[\beta _{1}X_{1}+\cdots +\beta _{p}X_{p}])}(Especificar los coeficientes de regresión con un signo negativo implica que los valores altos de las covariables aumentan el tiempo de supervivencia, pero esto es simplemente una convención de signos; sin un signo negativo, aumentan el riesgo).

Esto se cumple si se toma como función de densidad de probabilidad del eventoF(t|θ)=θF0(θt){\displaystyle f(t|\theta )=\theta f_{0}(\theta t)}; entonces se deduce para la función de supervivencia queS(t|θ)=S0(θt){\displaystyle S(t|\theta )=S_{0}(\theta t)}De esto se deduce fácilmente que la vida media moderadaT{\displaystyle T}se distribuye de tal manera queTθ{\displaystyle T\theta }y la vida sin moderaciónT0{\displaystyle T_{0}}tienen la misma distribución. En consecuencia,registro(T){\displaystyle \log(T)}se puede escribir como

registro(T)=registro(θ)+registro(Tθ):=registro(θ)+ϵ{\displaystyle \log(T)=-\log(\theta )+\log(T\theta ):=-\log(\theta )+\epsilon }

donde el último término se distribuye comoregistro(T0){\displaystyle \log(T_{0})}, es decir, independientemente deθ{\displaystyle \theta }. Esto reduce el modelo de tiempo de falla acelerado a un análisis de regresión (típicamente un modelo lineal ) donderegistro(θ){\displaystyle -\log(\theta )}representa los efectos fijos yϵ{\displaystyle \epsilon }representa el ruido. Diferentes distribuciones deϵ{\displaystyle \epsilon }implican diferentes distribuciones deT0{\displaystyle T_{0}}, es decir, diferentes distribuciones de referencia del tiempo de supervivencia. Típicamente, en contextos de análisis de supervivencia, muchas de las observaciones están censuradas: solo sabemos queTi>ti{\displaystyle T_{i}>t_{i}}, noTi=ti{\displaystyle T_{i}=t_{i}}De hecho, el primer caso representa la supervivencia, mientras que el segundo caso representa un evento/muerte/censura durante el seguimiento. Estas observaciones censuradas por la derecha pueden plantear desafíos técnicos para estimar el modelo, si la distribución deT0{\displaystyle T_{0}}es inusual.

La interpretación deθ{\displaystyle \theta }En los modelos de tiempo de falla acelerado es sencillo:θ=2{\displaystyle \theta =2}significa que todo en la historia de vida relevante de un individuo sucede el doble de rápido. Por ejemplo, si el modelo se refiere al desarrollo de un tumor, significa que todas las etapas previas progresan el doble de rápido que para el individuo no expuesto, lo que implica que el tiempo esperado hasta una enfermedad clínica es 0,5 del tiempo de referencia. Sin embargo, esto no significa que la función de riesgoλ(t|θ){\displaystyle \lambda (t|\theta )}siempre es el doble de alto; ese sería el modelo de riesgos proporcionales .

Cuestiones estadísticas

A diferencia de los modelos de riesgos proporcionales, en los que el modelo semiparamétrico de riesgos proporcionales de Cox se utiliza más ampliamente que los modelos paramétricos, los modelos AFT son predominantemente totalmente paramétricos, es decir, se especifica una distribución de probabilidad pararegistro(T0){\displaystyle \log(T_{0})}(Buckley y James [ 3 ] propusieron una AFT semiparamétrica, pero su uso es relativamente poco común en la investigación aplicada; en un artículo de 1992, Wei [ 4 ] señaló que el modelo de Buckley-James carece de justificación teórica y de robustez, y revisó alternativas). Esto puede ser un problema si se requiere un cierto grado de detalle realista para modelar la distribución de la vida útil de referencia. Por lo tanto, los avances técnicos en esta dirección serían muy deseables.

Cuando se incorpora un término de fragilidad en el modelo de supervivencia, las estimaciones de los parámetros de regresión de los modelos AFT son robustas ante la omisión de covariables , a diferencia de los modelos de riesgos proporcionales. También se ven menos afectadas por la elección de la distribución de probabilidad para el término de fragilidad. [ 5 ] [ 6 ]

Los resultados de los modelos AFT son fáciles de interpretar. [ 7 ] Por ejemplo, los resultados de un ensayo clínico con la mortalidad como criterio de valoración podrían interpretarse como un cierto porcentaje de aumento en la esperanza de vida futura con el nuevo tratamiento en comparación con el grupo de control. Así, se podría informar a un paciente que se esperaría que viviera (por ejemplo) un 15 % más si tomara el nuevo tratamiento. Los cocientes de riesgo pueden ser más difíciles de explicar en términos sencillos.

Distribuciones utilizadas en los modelos AFT

La distribución log-logística proporciona el modelo AFT más utilizado . A diferencia de la distribución de Weibull , puede exhibir una función de riesgo no monótona que aumenta en los primeros tiempos y disminuye en los últimos. Es algo similar en forma a la distribución log-normal , pero tiene colas más pesadas. La función de distribución acumulativa log-logística tiene una forma cerrada simple , lo que se vuelve importante computacionalmente cuando se ajustan datos con censura . Para las observaciones censuradas se necesita la función de supervivencia, que es el complemento de la función de distribución acumulativa, es decir, se necesita poder evaluarS(t|θ)=1F(t|θ){\displaystyle S(t|\theta )=1-F(t|\theta )}.

La distribución de Weibull (incluida la distribución exponencial como caso particular) puede parametrizarse como un modelo de riesgos proporcionales o un modelo AFT, y es la única familia de distribuciones que posee esta propiedad. Por lo tanto, los resultados del ajuste de un modelo de Weibull pueden interpretarse en cualquiera de los dos marcos. Sin embargo, la aplicabilidad biológica de este modelo puede verse limitada por el hecho de que la función de riesgo es monótona, es decir, decreciente o creciente.

Cualquier distribución en un grupo multiplicativamente cerrado , como los números reales positivos , es adecuada para un modelo AFT. Otras distribuciones incluyen la log-normal , la gamma , la hipertabástica , la de Gompertz y la gaussiana inversa , aunque son menos populares que la log-logística, en parte porque sus funciones de distribución acumulativa no tienen una forma cerrada. Finalmente, la distribución gamma generalizada es una distribución de tres parámetros que incluye las distribuciones de Weibull , log-normal y gamma como casos particulares.

Referencias

  1. Stroustrup, Nicholas (16 de enero de 2016). "La escala temporal del envejecimiento de Caenorhabditis elegans" . Nature . 530 ( 7588): 103–107 . Bibcode : 2016Natur.530..103S . doi : 10.1038/nature16550 . PMC 4828198. PMID 26814965 .  
  2. Kalbfleisch y Prentice (2002). Análisis estadístico de datos de tiempo de falla (2.ª ed.) . Hoboken, NJ: Wiley Series in Probability and Statistics.
  3. Buckley, Jonathan; James, Ian (1979), "Regresión lineal con datos censurados", Biometrika , 66 (3): 429–436 , doi : 10.1093/biomet/66.3.429 , JSTOR 2335161 
  4. Wei, LJ (1992). "El modelo de tiempo de falla acelerado: una alternativa útil al modelo de regresión de Cox en el análisis de supervivencia". Statistics in Medicine . 11 ( 14–15 ): 1871–1879 . doi : 10.1002/sim.4780111409 . PMID 1480879 . 
  5. Lambert, Philippe; Collett, Dave; Kimber, Alan; Johnson, Rachel (2004), "Modelos paramétricos de tiempo de falla acelerado con efectos aleatorios y una aplicación a la supervivencia del trasplante renal" , Statistics in Medicine , 23 (20): 3177–3192 , doi : 10.1002/sim.1876 , hdl : 2268/24489 , PMID 15449337 
  6. Keiding, N.; Andersen, PK; Klein, JP (1997). "El papel de los modelos de fragilidad y los modelos de tiempo de falla acelerado en la descripción de la heterogeneidad debida a covariables omitidas". Statistics in Medicine . 16 ( 1–3 ): 215–224 . doi : 10.1002/(SICI)1097-0258(19970130)16:2 < 215::AID-SIM481 > 3.0.CO ; 2-J . PMID 9004393 . 
  7. Kay, Richard; Kinnersley, Nelson (2002), "Sobre el uso del modelo de tiempo de falla acelerado como alternativa al modelo de riesgos proporcionales en el tratamiento de datos de tiempo hasta el evento: un estudio de caso en influenza" , Drug Information Journal , 36 (3): 571–579 , doi : 10.1177/009286150203600312

Lecturas adicionales

  • Bradburn, MJ; Clark, TG; Love, SB; Altman, DG (2003), "Análisis de supervivencia Parte II: Análisis de datos multivariados: una introducción a conceptos y métodos", British Journal of Cancer , 89 (3): 431– 436, doi : 10.1038/sj.bjc.6601119 , PMC 2394368 , PMID 12888808  
  • Hougaard, Philip (1999), "Fundamentos de los datos de supervivencia", Biometrics , 55 (1): 13–22 , doi : 10.1111/j.0006-341X.1999.00013.x , PMID 11318147 
  • Collett, D. (2003), Modelado de datos de supervivencia en la investigación médica (2.ª  ed.), CRC Press, ISBN 978-1-58488-325-8
  • Cox, David Roxbee ; Oakes, D. (1984), Análisis de datos de supervivencia , CRC Press, ISBN 978-0-412-24490-2
  • Marubini, Ettore; Valsecchi, Maria Grazia (1995), Analysing Survival Data from Clinical Trials and Observational Studies , Wiley, ISBN 978-0-470-09341-2
  • Martinussen, Torben; Scheike, Thomas (2006), Modelos de regresión dinámica para datos de supervivencia, Springer, ISBN 0-387-20274-9
  • Bagdonavicius, Vilijandas; Nikulin, Mikhail (2002), Modelos de vida acelerada. Modelado y análisis estadístico, Chapman&Hall/CRC, ISBN 1-58488-186-0