Respuestas rápidas

array:20 ["pii" => "13042666" "issn" => "02139111" "estado" => "S300" "fechaPublicacion" => "2002-12-02" "documento" => "article" "crossmark" => 0 "licencia" => "http://www.elsevier.com/open-access/userlicense/1.0/" "subdocumento" => "fla" "cita" => "Gac Sanit. 2002;16 Supl 2:59-68" "abierto" => array:3 [ "ES" => true "ES2" => true "LATM" => true ] "gratuito" => true "lecturas" => array:2 [ "total" => 6811 "formatos" => array:3 [ "EPUB" => 154 "HTML" => 5667 "PDF" => 990 ] ] "itemSiguiente" => array:16 [ "pii" => "13042667" "issn" => "02139111" "estado" => "S300" "fechaPublicacion" => "2002-12-02" "documento" => "article" "crossmark" => 0 "licencia" => "http://www.elsevier.com/open-access/userlicense/1.0/" "subdocumento" => "fla" "cita" => "Gac Sanit. 2002;16 Supl 2:70-6" "abierto" => array:3 [ "ES" => true "ES2" => true "LATM" => true ] "gratuito" => true "lecturas" => array:2 [ "total" => 19787 "formatos" => array:3 [ "EPUB" => 232 "HTML" => 18363 "PDF" => 1192 ] ] "es" => array:12 [ "idiomaDefecto" => true "titulo" => "La publicidad del tabaco, ¿incrementa el hábito de fumar entre los adolescentes? Datos desde California" "tienePdf" => "es" "tieneTextoCompleto" => "es" "tieneResumen" => array:2 [ 0 => "es" 1 => "en" ] "paginas" => array:1 [ 0 => array:2 [ "paginaInicial" => "70" "paginaFinal" => "76" ] ] "titulosAlternativos" => array:1 [ "en" => array:1 [ "titulo" => "Tobacco advertising, does it increase smoking among teenagers? Data from California" ] ] "contieneResumen" => array:2 [ "es" => true "en" => true ] "contieneTextoCompleto" => array:1 [ "es" => true ] "contienePdf" => array:1 [ "es" => true ] "resumenGrafico" => array:2 [ "original" => 0 "multimedia" => array:8 [ "identificador" => "fig1" "etiqueta" => "Figura 2" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "copyright" => "Elsevier España" "figura" => array:1 [ 0 => array:4 [ "imagen" => "138v16nSupl.2-13042667fig02.jpg" "Alto" => 443 "Ancho" => 290 "Tamanyo" => 17527 ] ] "descripcion" => array:1 [ "es" => "Anuncio de Lucky Strike de los años veinte." ] ] ] "autores" => array:1 [ 0 => array:2 [ "autoresLista" => "JP Pierce, M García, E Saltó, A Schiaffino, E Fernández" "autores" => array:5 [ 0 => array:2 [ "Iniciales" => "JP" "apellidos" => "Pierce" ] 1 => array:2 [ "Iniciales" => "M" "apellidos" => "García" ] 2 => array:2 [ "Iniciales" => "E" "apellidos" => "Saltó" ] 3 => array:2 [ "Iniciales" => "A" "apellidos" => "Schiaffino" ] 4 => array:2 [ "Iniciales" => "E" "apellidos" => "Fernández" ] ] ] ] ] "idiomaDefecto" => "es" "EPUB" => "https://multimedia.elsevier.es/PublicationsMultimediaV1/item/epub/13042667?idApp=WGSE" "url" => "/02139111/00000016000000S2/v0_201302051443/13042667/v0_201302051446/es/main.assets" ] "itemAnterior" => array:16 [ "pii" => "13042665" "issn" => "02139111" "estado" => "S300" "fechaPublicacion" => "2002-12-02" "documento" => "article" "crossmark" => 0 "licencia" => "http://www.elsevier.com/open-access/userlicense/1.0/" "subdocumento" => "fla" "cita" => "Gac Sanit. 2002;16 Supl 2:41-58" "abierto" => array:3 [ "ES" => true "ES2" => true "LATM" => true ] "gratuito" => true "lecturas" => array:2 [ "total" => 5702 "formatos" => array:3 [ "EPUB" => 160 "HTML" => 4990 "PDF" => 552 ] ] "es" => array:11 [ "idiomaDefecto" => true "titulo" => "Comportamientos adictivos no saludables: una revisión del análisis económico" "tienePdf" => "es" "tieneTextoCompleto" => "es" "tieneResumen" => array:2 [ 0 => "es" 1 => "en" ] "paginas" => array:1 [ 0 => array:2 [ "paginaInicial" => "41" "paginaFinal" => "58" ] ] "titulosAlternativos" => array:1 [ "en" => array:1 [ "titulo" => "Unhealthy addictive behaviours: a review of economic analysis" ] ] "contieneResumen" => array:2 [ "es" => true "en" => true ] "contieneTextoCompleto" => array:1 [ "es" => true ] "contienePdf" => array:1 [ "es" => true ] "autores" => array:1 [ 0 => array:2 [ "autoresLista" => "F Portillo, F Antoñanzas" "autores" => array:2 [ 0 => array:2 [ "Iniciales" => "F" "apellidos" => "Portillo" ] 1 => array:2 [ "Iniciales" => "F" "apellidos" => "Antoñanzas" ] ] ] ] ] "idiomaDefecto" => "es" "EPUB" => "https://multimedia.elsevier.es/PublicationsMultimediaV1/item/epub/13042665?idApp=WGSE" "url" => "/02139111/00000016000000S2/v0_201302051443/13042665/v0_201302051445/es/main.assets" ] "es" => array:14 [ "idiomaDefecto" => true "titulo" => "Modelos marginales y condicionales en el análisis de supervivencia multivariante" "tieneTextoCompleto" => true "paginas" => array:1 [ 0 => array:2 [ "paginaInicial" => "59" "paginaFinal" => "68" ] ] "autores" => array:1 [ 0 => array:3 [ "autoresLista" => "MA Barceló" "autores" => array:1 [ 0 => array:3 [ "Iniciales" => "MA" "apellidos" => "Barceló" "referencia" => array:1 [ 0 => array:2 [ "etiqueta" => "a" "identificador" => "affa" ] ] ] ] "afiliaciones" => array:1 [ 0 => array:3 [ "entidad" => "Grup de Recerca en Estadística, Economia Aplicada i Salut (GRECS). Departament d'Economia. Universitat de Girona. Centre de Recerca en Economia de la Salut (CRES). Universitat Pompeu Fabra. Barcelona." "etiqueta" => "a" "identificador" => "affa" ] ] ] ] "titulosAlternativos" => array:1 [ "en" => array:1 [ "titulo" => "Marginal and conditional models in multivariate survival analysis" ] ] "textoCompleto" => "Introducción En los estudios de supervivencia se está interesado por el tiempo que transcurre desde un origen (nacimiento, inicio de tratamiento, etc.) hasta que ocurre un suceso (muerte, mejora, etc.)1. En numerosas ocasiones este suceso no es único y los sujetos pueden experimentar varios sucesos (o técnicamente fallos) que pueden ser de distinto tipo o bien tratarse de recurrencias del mismo suceso de interés. Debido a esto, los individuos con más sucesos estarán «sobrerrepresentados»2 y, por otra parte, los efectos de las variables explicativas (efectos «tratamiento») podrán ser distintos entre el primer suceso y los siguientes. Esto se traducirá en una posible dependencia entre los tiempos de recurrencia, consecuencia de dos fenómenos diferentes: la existencia de heterogeneidad individual y dependencia «serial». En la primera, los tiempos de supervivencia dependen unos de otros porque los individuos comparten una o más covariables omitidas. Esta heterogeneidad individual también se denomina fragilidad, dependencia entre individuos o falso «contagio». Por lo que respecta a la dependencia «serial», que equivaldría a la autocorrelación en series temporales, la dependencia se produce en el mismo individuo, entre sus observaciones repetidas. Como consecuencia de que el individuo presenta una conducta dinámica per se, la probabilidad de que se produzca un suceso no es independiente de que se produzca uno previo. Esta dependencia también se denomina dependencia dentro de un mismo individuo o «contagio» verdadero. Con el fin de ilustrar el trabajo utilizaremos una cohorte prospectiva formada por todos los pacientes admitidos un mínimo de 24 h en la Unidad de Cuidados Intensivos (UCI) del Hospital Universitario Dr. Josep Trueta de Girona, durante el período comprendido entre el 15 de marzo y el 15 de junio de 19993. Se trata de un hospital público de nivel terciario, con 12 camas en la UCI, que cubre prácticamente toda la población de la provincia de Girona. Se pretende analizar los factores que determinan la incidencia de las infecciones nosocomiales en la UCI de dicho hospital. Las infecciones nosocomiales vienen definidas como las causadas por microorganismos adquiridos en el hospital, que afectan a los pacientes por un proceso diferente del de la infección, y que no estaban presentes en la admisión del paciente, ni tan siquiera en período de incubación. Entre los posibles factores de riesgo para la aparición de dichas infecciones distinguiremos los asociados a los pacientes (factores de riesgo intrínsecos, como su edad o sexo) y los relacionados con la UCI (factores de riesgo extrínsecos, como tratamiento antibiótico). El número total de sujetos incluidos en el estudio fue de 61 pacientes (42 varones y 19 mujeres), con un total de 40 infecciones nosocomiales. La variable dependiente (episodio de infección) consiste en el tiempo transcurrido entre el ingreso en la UCI hasta la aparición de la infección (originada por un microorganismo, bacteria u hongo), así como el tiempo entre el principio y el final de una infección. Obviamente los pacientes podrán presentar más de un episodio de infección. Sólo se analizaron infecciones nosocomiales originadas en la UCI, incluyéndose en la definición las infecciones comunitarias, las procedentes de otros hospitales y las procedentes de otras áreas hospitalarias. Para una explicación más detallada el lector interesado puede consultar en Barceló y Saez3. El objetivo principal de este artículo es presentar una revisión de la aproximación marginal y condicional en el análisis de supervivencia multivariante. Adicionalmente, y como ilustración, se presenta una aplicación a datos reales que permite valorar la utilidad de tales aproximaciones.Métodos Con el propósito de solucionar los problemas ocasionados por la existencia de más de un tiempo de supervivencia por individuo, el modelo de Cox de riesgos proporcionales fue generalizado y situado en el contexto mucho más general de los modelos de intensidad multiplicativos para procesos contadores multivariantes, a través de la aproximación de Andersen y Gill (modelo AG)4-8. En la aproximación AG cada sujeto es tratado como un proceso contador con sucesos múltiples y con «incrementos» (es decir, los tiempos entre acontecimientos sucesivos) independientes, dada la «historia» de todas las variables observables hasta el tiempo de presentación de los sucesos. Es decir, con el modelo AG el riesgo de una determinada recurrencia para un individuo sigue la hipótesis de riesgos proporcionales y no está afectado por la presentación en dicho individuo de sucesos anteriores. Dicha hipótesis, llamada de «incrementos independientes», es la hipótesis clave del modelo AG. En caso de existir alguna interrelación entre sucesos se modelizaría mediante una o más variables explicativas dependientes del tiempo9. Desde un punto de vista práctico, la aproximación AG representa a cada individuo por un conjunto de observaciones sij, tij, dij, xij, kij, j = 1, ..., ni siendo (sij, tij] el intervalo de riesgo, abierto a la izquierda y cerrado a la derecha; dij = 1 si el individuo tuvo un suceso en el instante tij y 0 en otro caso; xij es un vector de variables explicativas en el intervalo, y kij denota un estrato al que puede pertenecer el sujeto durante el intervalo. Dependiendo de la escala de tiempo, la primera observación puede o no empezar en cero. Aclararemos lo anterior mediante un sencillo ejemplo. Sean dos individuos, un varón que padeció dos episodios de infección (INF) y una mujer que no presentó ninguno (tabla 1 y fig. 1a). Los individuos se dividen en distintos registros que se corresponden con diversas bandas temporales (intervalos de observación) y cada registro es tratado como si fuesen individuos diferentes. Las bandas se construyen en función de si el individuo presenta o no un episodio de infección y de posibles cambios en las variables explicativas, no necesitando ser adyacentes. En este diseño las variables explicativas son fijas dentro de cada episodio, pero no así dentro de cada individuo (tal como suponen los modelos estándar, como el modelo de Cox de riesgos proporcionales). Nótese que es el riesgo de cada suceso concreto para un individuo determinado el que cumple la hipótesis usual de riesgos proporcionales y que el modelo supone que este riesgo no está afectado por sucesos anteriores del mismo individuo.<img src="138v16nSupl.2-13042666tab01.gif"></img><img src="138v16nSupl.2-13042666tab02.gif"></img>Figura 1. Estructura de los datos a comparar según el modelo de análisis. La aproximación AG es sencilla, pero sus hipótesis subyacentes son muy restrictivas y pueden llegar a ser irrealizables. En concreto, no es razonable suponer, tal como hemos visto, que las bandas pertenecientes a un mismo individuo no estén correlacionadas. En la práctica esta limitación puede conducir a estimadores sesgados e ineficientes. Además, puede sobrestimarse la precisión de los estimadores puesto que, habitualmente, las observaciones de un mismo individuo suelen estar correlacionadas positivamente. Como posible solución, a finales de los años ochenta empezaron a aparecer los modelos marginales y condicionales en el análisis de datos de supervivencia multivariante. Sin embargo, la aplicación de estos modelos aun hoy día no está muy extendida debido, principalmente, a una falta de sistematización en la bibliografía especializada. Es por ello por lo que consideramos especialmente importante describir estos modelos no sólo desde un punto de vista teórico10, sino también ilustrarlos sobre la base de datos utilizada anteriormente. En los modelos marginales la dependencia de las observaciones puede entenderse como un «estorbo», como algo que hay que eliminar o, en todo caso, cuyos efectos se deben controlar en el análisis, pero no como objetivo en sí. Por el contrario, los modelos condicionales consideran que la dependencia puede ser también de interés por sí misma. A grandes rasgos, los modelos marginales estiman el modelo ignorando la posible dependencia entre las observaciones, aunque corrigen posteriormente la varianza «ingenua» mediante estimadores jackknife, bootstrap o «sándwich». Los modelos condicionales estiman la dependencia (o la heterogeneidad) especificando explícitamente su distribución de probabilidad e incorporándola en el modelo. Modelos marginales La primera aproximación marginal que comentaremos es la estimación robusta de la matriz de varianzas y covarianzas. Se parte del hecho de que la aproximación AG proporciona estimadores consistentes de los parámetros asociados a las variables explicativas. Aun así, si existe dependencia, el método proporciona varianzas sesgadas. Se trataría, por tanto, de estimar de forma robusta la matriz de covarianzas, es decir, de corregirla de forma adecuada utilizando, por ejemplo, un estimador jackknife9. La idea es sencilla, se trata de ir eliminando a un individuo de la muestra cada vez y restimar el modelo resultante. Los estimadores finales se obtienen ponderando los obtenidos tras cada una de estas eliminaciones. El problema de esta aproximación marginal estriba en que se basa en la consistencia de las estimaciones. Sin embargo, las estimaciones de los parámetros serán consistentes siempre y cuando el modelo marginal esté correctamente especificado. Otro modelo marginal es el modelo de Wei, Lin y Weissfeld (WLW)11, en el que cada recurrencia del suceso (o tipo distinto de suceso) se modeliza como un estrato diferente, utilizándose los datos dentro de cada estrato de modo marginal. El resultado es que (todos) los individuos aparecen en todos los estratos. Es decir, el modelo WLW no es más que un modelo estratificado, siendo cada estrato una presentación del suceso de interés. Los estimadores finales son una media ponderada de los estimadores en cada estrato. Como vemos, se trata de un modelo marginal con respecto al conjunto de riesgo por cuanto el sujeto está en riesgo desde el inicio del estudio y, por tanto, no depende de sucesos previos. Es decir, para cada una de las k recurrencias que pueda padecer un individuo i (k = 1,2,..., K), la función de riesgo puede expresarse como:λik(t) = λi0k(t) eβ'kXik(t) t >= 0 [1] siendo λi0k(t) (k = 1,..., K) una función de riesgo subyacente no especificada, distinta para cada recurrencia, y βk = (β1k,..., βpk)' un vector de orden p ×1 de parámetros desconocidos de la regresión. En el modelo WLW los datos se organizan de forma diferente de la que veíamos en el modelo AG (tabla 1 y fig. 1b, en la que se ha supuesto que el máximo número de episodios de infección entre todos los individuos analizados es de tres). Nótese que cada intervalo de tiempo empieza en cero. Se trata de un modelo marginal con respecto al conjunto de riesgo por cuanto el sujeto está en riesgo desde el inicio del estudio y, por tanto, no depende de sucesos previos. Las limitaciones del modelo WLW pueden ser superadas, en parte, utilizando la aproximación de Prentice, Williams y Peterson (PWP)12. Se trata de un modelo marginal respecto a la estimación de los parámetros, pero condicional en relación con la construcción del conjunto de individuos en riesgo. El modelo PWP permite que el riesgo basal varíe entre diferentes recurrencias. Se trata, por tanto, de un modelo de riesgos proporcionales con estratos dependientes del tiempo, en el que la dependencia entre recurrencias se controla estratificando por el número previo de presentaciones del suceso de interés. Así, a diferencia del método WLW, en el que todos los individuos están en riesgo en todo momento, el modelo PWP incluye en el conjunto de riesgo para la recurrencia k únicamente a los individuos que han experimentado k-1 recurrencias. Prentice, Williams y Peterson12 proponen dos modelos semiparamétricos de riesgos proporcionales:1. Tiempo desde el inicio del estudio (de tiempo total). Se incluye la función de riesgo subyacente como una función del tiempo desde el inicio del estudioλik (t) = λ0k(t) eβ'kXik (t)[2]2. Tiempo desde la recurrencia inmediatamente precedente (de brecha de tiempo). La función de riesgo subyacente es una función del tiempo desde la recurrencia inmediatamente precedente al tiempo de fallo:λik (t) = λ0k (t tn(t))eβ'kXik (t)[3] La particularidad de ambos modelos reside en que se estratifica en función del número de recurrencias previas. Siguiendo con el ejemplo de antes, tanto para el modelo AG como para PWP los datos se organizarían tal como se muestra en la tabla 1 y la figura 1c. Modelos condicionales En algunos casos no sólo nos interesa el efecto marginal de las covariables sobre la función de riesgo de los tiempos de fallo de los individuos, sino también cómo los sucesos anteriores pueden influir en el riesgo de desarrollar futuros fallos. Para responder a estas cuestiones necesitamos recurrir a los modelos condicionales. En primer lugar, la hipótesis de incrementos independientes en el modelo AG puede relajarse incluyendo en el modelo variables dependientes del tiempo, como el número de recurrencias previas, a fin de capturar la posible estructura de dependencia entre las recurrencias. En concreto, se podría trabajar con el modelo autorregresivo o markoviano. Este modelo será plenamente apropiado cuando la dependencia entre observaciones sea, principalmente, «serial». Sin embargo, y aunque es un modelo muy general, supone que la dependencia entre recurrencias decae exponencialmente y que la dependencia es idéntica para todos los individuos10. Parte de las limitaciones anteriores pueden solucionarse mediante los modelos de fragilidad. Se trata de introducir una (modelos de fragilidad univariante) o más (modelos de fragilidad multivariante) variables aleatorias comunes a todos los individuos con el objetivo de permitir que los tiempos de supervivencia multivariantes, es decir, las observaciones repetidas de un mismo individuo, estén correlacionadas. Se trata, en definitiva, de que este efecto multiplicativo en la función de intensidad (o de riesgo) capture dos fuentes de variación distintas pero relacionadas13, con el propósito de modelar la dependencia entre los diferentes tiempos de supervivencia entre recurrencias. Una de las fuentes de variación surge de covariables comunes a todos los individuos, no observadas, las cuales, cuando se omiten, generan dependencia entre sucesos. La otra fuente de variación proviene de covariables individuales no observadas y, por tanto, no incluidas en el estudio, bien por circunstancias prácticas, bien porque no se asocian a ningún factor de riesgo de la supervivencia del individuo. Los efectos de esta fuente serían totalmente equivalentes a la dependencia serial citada en la introducción de este artículo. Posibles alternativas para la distribución de probabilidad de la variable aleatoria fragilidad son la distribución gamma, la distribución positiva estable y la distribución gaussiana inversa, entre otras. El modelo de fragilidad gamma es el más adecuado en nuestro caso, debido, en primer lugar, al tipo de datos de que disponemos. En concreto, bajo este tipo de distribución la fuerza de la asociación entre recurrencias es estacionaria, es decir, no depende del tiempo. Este tipo de asociación es la que más se produce en el mundo real. La ventaja más importante de la distribución gamma, sin embargo, se encuentra en la facilidad analítica que permite. En particular, el modelo AG de fragilidad14,15 puede expresarse de la siguiente manera en su caso más simple:λi(t) = ωiYi (t) α0(t) eβ'Xi (t) [4] siendo N = (Ni; i = 1, ..., n) un proceso contador multivariante que cuenta el número de sucesos que han ocurrido en cada uno de los n individuos en el tiempo t. Asociado a este proceso contador, existe un proceso de intensidad λcon componentes λi . El proceso Y con componentes Yi es un proceso no negativo, observable y predecible, es decir, su valor en cualquier momento t es conocido justo antes de t, indicando si un individuo i se observa en el conjunto a riesgo de experimentar un determinado suceso en el tiempo t-; α0 denota ahora una función de riesgo basal desconocida, y ωi (fragilidad) es un conjunto de variables aleatorias independientes e idénticamente distribuidas según una distribución de probabilidad paramétrica determinada, en este caso la distribución gamma, con media nula y varianza común. La variabilidad (medida por su varianza) de la fragilidad determina el grado de heterogeneidad. De hecho, cuando la varianza es igual a cero el modelo se reduce al modelo AG estándar. El modelo descrito es un modelo muy simple, ya que asume que la fragilidad es constante a lo largo del tiempo, es decir, la fragilidad reflejaría en este caso diferencias individuales presentes al comienzo del estudio. Asimismo, este modelo supone que la fragilidad actúa proporcionalmente sobre la función basal común a todos los individuos. Al imponer esto se asume que todos los individuos siguen un patrón uniforme. El problema es que la estimación de los parámetros del modelo utilizando directamente la máxima verosimilitud no es posible cuando la función de riesgo basal no está parametrizada, como en nuestro caso. Por otra parte, los modelos de fragilidad univariantes presentan una serie de limitaciones; entre ellas cabe destacar las siguientes: suponen que la dependencia entre observaciones es siempre positiva, no permiten la interacción entre covariables observables y la fragilidad, y más importante aún, sólo permiten la presencia de un único efecto aleatorio, siendo además éste constante a lo largo del tiempo. Con el propósito de solucionar estos problemas, Barceló y Saez16 proponen un método de especificación y estimación original basado en la verosimilitud extendida de Nelder y Lee17, así como en el método de las variables instrumentales, y estudian las ventajas analíticas de dicho método, en particular, en relación con alternativas de especificación (modelos de fragilidad multivariantes) y de estimación (algoritmo EM). De una manera muy resumida, introducen un nuevo paso en el algoritmo EM, modelizando la varianza de la fragilidad directamente. Mediante el algoritmo así modificado (algoritmo EMB), se estima no sólo el modelo de supervivencia, sino también la varianza de la fragilidad de forma simultánea. Mediante esta modificación solucionan parte de las limitaciones citadas anteriormente; en particular, permiten que la fragilidad no sea constante entre individuos. Una explicación más detallada puede encontrarse en Barceló y Saez16. Una alternativa a los modelos de fragilidad la proporcionan los modelos de Cox penalizados18. En resumen, se trata de maximizar una función de verosimilitud parcial penalizada. El término penalizador captura la variabilidad local presente en la densidad conjunta de los datos. Los modelos que se han comentado no están anidados, por lo que las usuales medidas de bondad de ajuste no pueden utilizarse. Como alternativa se suele emplear el criterio de información de Akaike (AIC). El software estadístico que permite ajustar todos los modelos vistos y que, por otra parte, se ha utilizado en la aplicación práctica, es el S-PLUS 200019. En concreto, se debe utilizar la instrucción coxph, que permite considerar diferentes tipos de truncamiento y censura. La estimación robusta de la varianza se obtiene mediante la instrucción cluster. Desafortunadamente, hay pocos programas que permitan la estimación de modelos marginales y/o condicionales en el contexto multivariante. Además del S-PLUS 2000, el lector interesado puede recurrir a otros programas como R (versión libre y gratuita de S-PLUS), SAS (versión 6 o posteriores, instrucciones proc phreg, aunque no permite la estimación robusta de la varianza) y, únicamente para el modelo AG estándar, STATA (versión 5 o posteriores, instrucciones stset y stcox).Resultados Los resultados de la estimación del modelo AG estándar aparecen en las primeras columnas de la tabla 2. En esta tabla se presentan las tasas de riesgo asociadas a toda una serie de variables explicativas y sus respectivos intervalos de confianza (IC). Por último, se ofrece una serie de medidas de bondad de ajuste del modelo. Como puede apreciarse al observar la tabla, pocas variables resultaron estar asociadas de forma significativa con el tiempo de presentación de una infección nosocomial. En concreto, ser mujer (0,258; IC del 95%, 0,067-0,997), haberse sometido a cirugía urgente (0,295; IC del 95%, 0,092-0,949), haber permanecido intubado más de 6 días mediante traqueotomía (0,182; IC del 95%, 0,041-0,814) y haber recibido tratamiento antibiótico durante más de 7 días (0,179; IC del 95%, 0,034-0,947) reducían el riesgo de tales infecciones.<img src="138v16nSupl.2-13042666tab03.gif"></img> Los resultados de la aplicación práctica del modelo AG con estimación robusta de la matriz de varianzas y covarianzas serían los que hemos visto anteriormente al hablar del modelo AG, ya que esta primera aproximación marginal sólo afecta a la matriz de varianzas y covarianzas, no a la construcción de los conjuntos a riesgo. Si nos fijamos en la tabla 2 y comparamos el modelo AG estándar y el modelo AG robusto, apreciamos una mayor precisión en los intervalos de confianza proporcionados por este último modelo, siendo este cambio muy significativo en variables como ventilación mecánica. Por el contrario, otras variables tales como cirugía urgente (0,295; IC del 95%, 0,083-1,048) y tratamiento antibiótico durante más de 7 días (0,179; IC del 95%, 0,031-1,047) llegan incluso a perder su significación estadística. En relación con los resultados del modelo PWP con estimación robusta de la varianza (tabla 3), con una importante excepción, presencia de catéter venoso central durante más de tres días (1,065; IC del 95%, 1,007-1,775), el signo de los estimadores de los parámetros coincide en los tres modelos vistos (modelo AG estándar, modelo AG robusto y modelo PWP), y no ocurre lo mismo con la magnitud ni con la significación de los parámetros. Nótese ahora que variables como la edad (1,009; IC del 95%, 0,985-1,033), la existencia previa de infecciones comunitarias (4,56; IC del 95%, 1,435-14,477), estar inestable en UCI (54,47; IC del 95%, 0,786-3.776,406), estar inestable en coma o shock (219,97, IC del 95%, 3,187-15.182,609), la presencia de ventilación mecánica de cuatro a 10 días (61,93, IC del 95%, 1,325-2.894,463), la presencia de ventilación mecánica durante más de 10 días (129,53; IC del 95%, 4,005-4.189,427) y la presencia de catéter venoso central durante más de tres días (1,07; IC del 95%, 1,007-1,775) aumentan el riesgo de padecer un episodio de infección nosocomial. Contrariamente, la presencia de sondas nasogástricas durante más de 9 días (0,13; IC del 95%, 0,023-0,748) y el tratamiento antibiótico durante más de 7 días (0,04; IC del 95%, 0,004-0,385) lo reducen. Variables tales como género y traqueotomía pierden en este modelo su significación. Resumiendo, no sólo los factores de riesgo intrínsecos (como la gravedad de la enfermedad), sino también los extrínsecos (procedimientos terapéuticos agresivos utilizados en la UCI) llevaron a un mayor riesgo de infección nosocomial durante la estancia del paciente en la UCI.<img src="138v16nSupl.2-13042666tab04.gif"></img> Los resultados del modelo AG de fragilidad gamma estimado mediante el algoritmo EMB, así como los del modelo de Cox penalizado, también se encuentran en la tabla 3. En el modelo AG de fragilidad gamma estimado mediante el algoritmo EMB, observamos que variables tales como la edad (1,009; IC del 95%, 0,990-1,028), existencia previa de infecciones comunitarias (6,923; IC del 95%, 2,578-18,594), estar inestable en UCI (49,012; IC del 95%, 1,111-2.162,135), estar inestable en coma o shock (257,040; IC del 95%, 5,997-11.017,03), la presencia de ventilación mecánica en todas sus categorías y la presencia de catéter venoso central durante más de tres días (1,082; IC del 95%, 1,011-1,822) continúan aumentando el riesgo de padecer infección nosocomial de forma estadísticamente significativa. Por el contrario, la presencia de sonda nasogástrica durante más de 9 días (0,129; IC del 95%, 0,03-0,55) y el tratamiento antibiótico en todas sus categorías lo reducen, también de forma estadísticamente significativa. Por otra parte, se observa que variables que no eran significativas en el modelo PWP, como la existencia previa de otras infecciones, traqueotomía y presencia de sonda urinaria de 5 a 12 días, pasan a serlo ahora. Estos resultados cambian un poco al observar el modelo de Cox penalizado. Mientras que variables tales como la edad (1,009; IC del 90%: 0,981-1,037), infecciones previas comunitarias (4,558; IC del 90%, 0,954-21,772), estar inestable en UCI (8.025,27; IC del 95%, 1.740-36.968,77), presencia de ventilación mecánica de cuatro a 10 días (64,406; IC del 95%, 1,35-3.065,471), presencia de ventilación mecánica durante más de 10 días (134,714; IC del 95%, 3,06-5.934,619), presencia de sonda nasogástrica durante más de 9 días (0,133; IC del 95%, 0,027-0,649) y tratamiento antibiótico durante más de 7 días (0,038; IC del 95%, 0,005-0,327) continúan manteniendo su significación, otras como traqueotomía, existencia de otras infecciones previas, presencia de sonda urinaria de 5 a 12 días y tratamiento antibiótico durante 7 días o menos la pierden. Por otra parte, el modelo de Cox penalizado parece que es incapaz de estimar los parámetros en algunas de las variables explicativas presentes en el modelo (inestable en coma o shock, ventilación mecánica inferior o igual a tres días y catéter venoso central), lo cual no ocurre en el modelo EMB. Asimismo, el modelo EMB es el que proporciona intervalos de confianza para los estimadores de los parámetros más reducidos. Adicionalmente, y a diferencia de los demás modelos, el modelo EMB permite observar los factores explicativos de la fragilidad. Como se ve, la varianza de la fragilidad depende únicamente de factores extrínsecos tales como ventilación mecánica, traqueotomía, sonda urinaria y sonda nasogástrica.Discusión Como ya se ha dicho, la problemática planteada por la existencia de datos de supervivencia multivariantes hace necesaria la generalización del modelo de Cox de riesgos proporcionales, concretándose ésta en el modelo AG. Este modelo permite trabajar con datos multivariantes; sin embargo, no contempla la posibilidad que las observaciones pertenecientes a un mismo individuo estén correlacionadas, ni la presencia de heterogeneidad individual. Con el propósito de solucionar esto, citamos los principales modelos marginales (modelo AG con estimación robusta de la varianza, modelo WLW y modelo PWP) y condicionales (modelo autorregresivo o markoviano, modelos de fragilidad y modelo de Cox penalizado). Aun cuando los modelos marginales son capaces de corregir los principales efectos de la dependencia entre observaciones (siempre suponiendo que las estimaciones de los parámetros son consistentes), son incapaces de distinguir cómo los sucesos anteriores pueden influir en el riesgo de desarrollar futuros fallos. Esto último es tratado mediante los modelos condicionales. El método WLW presenta importantes desventajas. Cook y Lawless20 señalan que el modelo WLW es válido únicamente bajo la posiblemente muy restrictiva hipótesis de censura independiente (restrictiva como consecuencia de las recurrencias). La desventaja más importante es que el método WLW supone que en el instante t todos los individuos (observados en t) constituyen el conjunto de riesgo para la correspondiente recurrencia, con independencia del número de recurrencias previas que tuviese el individuo20,21. Esto no representaría ningún problema si se tratase de diferentes sucesos de un tipo totalmente distinto, pero si se trata de recurrencias del mismo suceso podría ocurrir que algún hecho emprendido con el fin de reducir el número de tales recurrencias redujese el número de las mismas a partir de su implementación (lo mismo ocurriría en el caso de harvesting). Como consecuencia aumentaría la correlación entre los estimadores de los parámetros. Otro problema es que las estimaciones de los parámetros serán consistentes siempre y cuando la función de riesgo esté correctamente especificada. La principal desventaja del modelo PWP radica en la manera de construir los conjuntos a riesgo, ya que lleva a una «pérdida de aleatorización»22, por cuanto los individuos con más riesgo van saliendo del estudio. Esta pérdida será mayor conforme vayamos moviéndonos a estratos superiores (y cuantos más estratos se consideren). En estos casos, la consistencia de las estimaciones depende de la inclusión de todas las variables relevantes. Por otra parte, tanto el modelo AG como el modelo PWP son sensibles a los errores de especificación de la estructura de dependencia entre los tiempos de recurrencia. Aunque los métodos marginales cumplen (marginalmente) la hipótesis de riesgos proporcionales, suponen explícita (PWP) o implícitamente (WLW) que la función de riesgo basal varía de una recurrencia a otra22. Como consecuencia, son extremadamente sensibles al cumplimiento (en realidad) de la hipótesis de riesgos proporcionales9. Finalmente, y respecto a las aproximaciones condicionales, los métodos marginales no son totalmente eficientes23,24. La ineficiencia aumentará de forma proporcional al número de estratos (recurrencias) e inversamente proporcional al número de observaciones en cada estrato. En analogía a los modelos lineales es razonable suponer que la ineficiencia también aumentará conforme lo haga la correlación entre recurrencias. Es por tanto de esperar que WLW sea menos eficiente que PWP. Por lo que hace referencia a los modelos condicionales, sólo comentaremos aquí el problema en los modelos de Cox penalizados. Como argumentábamos anteriormente, la variabilidad que pretenden capturar estos modelos recoge, en realidad, dos fenómenos diferentes, fragilidad y correlación serial. Por otra parte, presentan el inconveniente de que el parámetro de suavizado debe escogerse utilizando técnicas no del todo validadas en el campo de supervivencia. Situándonos en nuestro ejemplo práctico, de entre los tres modelos marginales presentados en el artículo, el modelo PWP es posiblemente el mejor debido al tipo de datos de que disponemos. El hecho de que, tras experimentar una primera infección, aumente el riesgo de aparición de nuevas infecciones hace que el modelo PWP con estimación robusta de la varianza sea especialmente apropiado, ya que sería posible estratificar los datos para cada episodio de infección. Por lo que respecta a las estimaciones de los modelos condicionales mencionados en el artículo, el algoritmo EMB de estimación de un modelo AG de fragilidad gamma es el que presenta mejores resultados. Por un lado, es el más eficiente al presentar intervalos de confianza más reducidos, debido quizá a que el algoritmo EMB es capaz de capturar no sólo la heterogeneidad individual, sino también la dependencia serial. Por otro lado, de entre todos lo modelos estudiados, es el que proporciona un mejor ajuste en términos del criterio de Akaike. Finalmente, es el único que permite estimar los efectos explicativos de la fragilidad.Agradecimientos Agradezco al Dr. Marc Saez su gran ayuda en la elaboración de este trabajo. Agradezco, asimismo, al Dr. Santiago Pérez Hoyos y a dos evaluadores anónimos sus sugerencias. Este trabajo fue parcialmente financiado por el proyecto 07/50/1998 de la Agència d'Avaluació de Tecnologia Mèdica (AATM) del Servei Català de la Salut, Generalitat de Catalunya." "pdfFichero" => "138v16nSupl.2a13042666pdf001.pdf" "tienePdf" => true "PalabrasClave" => array:2 [ "es" => array:1 [ 0 => array:4 [ "clase" => "keyword" "titulo" => "Palabras clave" "identificador" => "xpalclavsec16355" "palabras" => array:4 [ 0 => "Análisis de supervivencia multivariante" 1 => "Enfoque marginal" 2 => "Enfoque condicional" 3 => "Modelos de fragilidad" ] ] ] "en" => array:1 [ 0 => array:4 [ "clase" => "keyword" "titulo" => "Keywords" "identificador" => "xpalclavsec16356" "palabras" => array:4 [ 0 => "Multivariate survival data" 1 => "Marginal models" 2 => "Conditional models" 3 => "Frailty models" ] ] ] ] "tieneResumen" => true "resumen" => array:2 [ "es" => array:1 [ "resumen" => "El hecho de que los individuos puedan experimentar varios sucesos (fallos) durante el período de observación (análisis de supervivencia multivariante) hace que sea necesaria la aplicación de una metodología diferente de la empleada en el análisis de supervivencia estándar univariante. Al analizar medidas repetidas, el principal problema del modelo de Cox se debe a que las observaciones no son independientes, lo cual llevará al incumplimiento de la hipótesis de riesgos proporcionales. Con el propósito de solucionar este problema, así como otros que puedan darse debido a este tipo de datos, aparece el modelo de Andersen-Gill (AG), que no es más que una generalización del modelo de Cox. Sin embargo, la hipótesis clave de este modelo (hipótesis de incrementos independientes) supone que las múltiples observaciones de un mismo individuo son independientes, condicionados a las variables explicativas. En la práctica, esta limitación puede conducir a estimadores sesgados e ineficientes, sobrestimándose normalmente la precisión de los mismos ya que, por lo común, las observaciones de un mismo individuo suelen estar positivamente correlacionadas. Con el fin de hallar una solución a esto, aparecen los modelos marginales y los modelos condicionales. La principal diferencia entre estos modelos es el trato que recibe la dependencia. Mientras que los primeros estiman el modelo ignorando la dependencia entre observaciones, corrigiéndola posteriormente mediante estimadores jackknife, bootstrap o «sándwich», los segundos la estiman especificando explícitamente la distribución de probabilidad de la misma e incorporándola en el modelo. Nuestro objetivo en este artículo consiste en ilustrar estas dos aproximaciones utilizando una base de datos de infecciones nosocomiales." ] "en" => array:1 [ "resumen" => "The fact that individuals can present more than one event (fail) in the observation period (multivariate survival analysis) demands a different methodology from that used in univariate survival analysis. The main problem of the Cox model with multivariate data is that the observations are not independent implying, among others, the violation of the proportionality hypothesis. The Andersen-Gill approximation to the Cox model (AG) overcomes in part this problem. The underlying hypotheses in the AG model, however, are very restrictive, in particular that of independent increments. Under this hypothesis the multiple observations of an individual are independent, although conditioned on the explanatory variables. In practice, this limitation could let us to biased and inefficient estimators that are usually overestimated because, in general, the observations of an individual use to be positively correlated. As a solution to this problem, marginal and conditional models turn up. The main difference between these models is the treatment that dependence receives. Whereas the first one estimates the model ignoring the dependence between recurrences and adjustes this at the end with jackknife, boostrap or «sandwich» estimates; the second one estimates the dependence specifying explicitly the probability distribution of the dependence and incorporating this into the model. Our objective in this paper is to illustrate these two approaches using nosocomial infection data." ] ] "multimedia" => array:8 [ 0 => array:8 [ "identificador" => "tbl1" "etiqueta" => "Tabla 1" "tipo" => "MULTIMEDIATABLA" "mostrarFloat" => true "mostrarDisplay" => false "copyright" => "Elsevier España" "tabla" => array:1 [ "tablatextoimagen" => array:1 [ 0 => array:1 [ "tablaImagen" => array:1 [ 0 => array:4 [ "imagenFichero" => "138v16nSupl.2-13042666tab01.gif" "imagenAlto" => 398 "imagenAncho" => 719 "imagenTamanyo" => 11875 ] ] ] ] ] "descripcion" => array:1 [ "es" => "Estructura de los datos según las diferentes aproximaciones" ] ] 1 => array:8 [ "identificador" => "tbl2" "etiqueta" => "Figura 1" "tipo" => "MULTIMEDIATABLA" "mostrarFloat" => true "mostrarDisplay" => false "copyright" => "Elsevier España" "tabla" => array:1 [ "tablatextoimagen" => array:1 [ 0 => array:1 [ "tablaImagen" => array:1 [ 0 => array:4 [ "imagenFichero" => "138v16nSupl.2-13042666tab02.gif" "imagenAlto" => 335 "imagenAncho" => 561 "imagenTamanyo" => 6798 ] ] ] ] ] "descripcion" => array:1 [ "es" => "Estructura de los datos a comparar según el modelo de análisis." ] ] 2 => array:8 [ "identificador" => "tbl3" "etiqueta" => "Tabla 2" "tipo" => "MULTIMEDIATABLA" "mostrarFloat" => true "mostrarDisplay" => false "copyright" => "Elsevier España" "tabla" => array:1 [ "tablatextoimagen" => array:1 [ 0 => array:1 [ "tablaImagen" => array:1 [ 0 => array:4 [ "imagenFichero" => "138v16nSupl.2-13042666tab03.gif" "imagenAlto" => 624 "imagenAncho" => 602 "imagenTamanyo" => 23663 ] ] ] ] ] "descripcion" => array:1 [ "es" => "Resultados de la estimación del modelo Andersen-Gill (AG) estándar y del modelo AG con estimación robusta de la matriz de variancias y covariancias" ] ] 3 => array:8 [ "identificador" => "tbl4" "etiqueta" => "Tabla 3" "tipo" => "MULTIMEDIATABLA" "mostrarFloat" => true "mostrarDisplay" => false "copyright" => "Elsevier España" "tabla" => array:1 [ "tablatextoimagen" => array:1 [ 0 => array:1 [ "tablaImagen" => array:1 [ 0 => array:4 [ "imagenFichero" => "138v16nSupl.2-13042666tab04.gif" "imagenAlto" => 810 "imagenAncho" => 599 "imagenTamanyo" => 32841 ] ] ] ] ] "descripcion" => array:1 [ "es" => "Resultados de la estimación del modelo de Prentice, Williams y Peterson (PWP) del modelo de Cox penalizado del modelo EMB. PWP y EMB con estimación robusta de la varianza" ] ] 4 => array:5 [ "identificador" => "tbl5" "tipo" => "MULTIMEDIATABLA" "mostrarFloat" => true "mostrarDisplay" => false "copyright" => "Elsevier España" ] 5 => array:5 [ "identificador" => "tbl6" "tipo" => "MULTIMEDIATABLA" "mostrarFloat" => true "mostrarDisplay" => false "copyright" => "Elsevier España" ] 6 => array:5 [ "identificador" => "tbl7" "tipo" => "MULTIMEDIATABLA" "mostrarFloat" => true "mostrarDisplay" => false "copyright" => "Elsevier España" ] 7 => array:5 [ "identificador" => "tbl8" "tipo" => "MULTIMEDIATABLA" "mostrarFloat" => true "mostrarDisplay" => false "copyright" => "Elsevier España" ] ] "bibliografia" => array:2 [ "titulo" => "Bibliografía" "seccion" => array:1 [ 0 => array:1 [ "bibliografiaReferencia" => array:24 [ 0 => array:3 [ "identificador" => "bib1" "etiqueta" => "1" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Regression models and life tables." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "Cox DR." ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "JRSS, Series B," "fecha" => "1972" "volumen" => "34" "paginaInicial" => "187" "paginaFinal" => "220" ] ] ] ] ] ] 1 => array:3 [ "identificador" => "bib2" "etiqueta" => "2" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Testing the proportional hazards assumption in the presence of unmeasured heterogeneity." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "McCall BP." ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "J Applied Econometrics" "fecha" => "1994" "volumen" => "9" "paginaInicial" => "321" "paginaFinal" => "34" ] ] ] ] ] ] 2 => array:3 [ "identificador" => "bib3" "etiqueta" => "3" "referencia" => array:1 [ 0 => array:2 [ "referenciaCompleta" => "Factors determining the incidence of nosocomial infections in an Intensive Care Unit of a public hospital in Girona, Spain. Mimeo, 2001." "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Factors determining the incidence of nosocomial infections in an Intensive Care Unit of a public hospital in Girona, Spain. Mimeo, 2001." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "Barceló MA" 1 => "Saez M." ] ] ] ] ] ] ] ] 3 => array:3 [ "identificador" => "bib4" "etiqueta" => "4" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Non-parametric inference for a family of counting processes." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "Aalen OO." ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Ann Stat" "fecha" => "1978" "volumen" => "6" "paginaInicial" => "701" "paginaFinal" => "26" ] ] ] ] ] ] 4 => array:3 [ "identificador" => "bib5" "etiqueta" => "5" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Cox's regression model for counting processes: a large sample study." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "Andersen PK" 1 => "Gill RD." ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Ann Stat" "fecha" => "1982" "volumen" => "10" "paginaInicial" => "1100" "paginaFinal" => "20" ] ] ] ] ] ] 5 => array:3 [ "identificador" => "bib6" "etiqueta" => "6" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Counting process models for life history data: a review." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "Andersen PK" 1 => "Borgan O." ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Scan J Stat" "fecha" => "1985" "volumen" => "12" "paginaInicial" => "97" "paginaFinal" => "158" ] ] ] ] ] ] 6 => array:3 [ "identificador" => "bib7" "etiqueta" => "7" "referencia" => array:1 [ 0 => array:2 [ "referenciaCompleta" => "Counting processes and survival analysis. New York: John Wiley & Sons, Inc., 1991." "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Counting processes and survival analysis. New York: John Wiley & Sons, Inc., 1991." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "Fleming TR" 1 => "Harrington DP." ] ] ] ] ] ] ] ] 7 => array:3 [ "identificador" => "bib8" "etiqueta" => "8" "referencia" => array:1 [ 0 => array:2 [ "referenciaCompleta" => "Statistical models based on counting proceses. New York: Springer-Verlag, 1993." "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Statistical models based on counting proceses. New York: Springer-Verlag, 1993." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:4 [ 0 => "Andersen PK" 1 => "Borgan O" 2 => "Gill RD" 3 => "Keiding N." ] ] ] ] ] ] ] ] 8 => array:3 [ "identificador" => "bib9" "etiqueta" => "9" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "rhDNase as an example of recurrent event analysis." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "Therneau TM" 1 => "Hamilton SA." ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:6 [ "tituloSerie" => "Stat Med" "fecha" => "1997" "volumen" => "16" "paginaInicial" => "2029" "paginaFinal" => "47" "link" => array:1 [ 0 => array:2 [ "url" => "https://www.ncbi.nlm.nih.gov/pubmed/9308130" "web" => "Medline" ] ] ] ] ] ] ] ] 9 => array:3 [ "identificador" => "bib10" "etiqueta" => "10" "referencia" => array:1 [ 0 => array:2 [ "referenciaCompleta" => "Modelos de fragilidad en el análisis multivariante de supervivencia [tesis doctoral]. Girona: Universitat de Girona, 2001." "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Modelos de fragilidad en el análisis multivariante de supervivencia [tesis doctoral]. Girona: Universitat de Girona, 2001." "idioma" => "es" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "Barceló MA." ] ] ] ] ] ] ] ] 10 => array:3 [ "identificador" => "bib11" "etiqueta" => "11" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Regression analysis of multivariate incomplete failure time data by modeling marginal distributions." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:3 [ 0 => "Wei LJ" 1 => "Lin DY" 2 => "Weissfeld L." ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "JASA" "fecha" => "1989" "volumen" => "84" "paginaInicial" => "1065" "paginaFinal" => "73" ] ] ] ] ] ] 11 => array:3 [ "identificador" => "bib12" "etiqueta" => "12" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "On the regression analysis of multivariate failure time data." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:3 [ 0 => "Prentice RL" 1 => "Williams BJ" 2 => "Peterson AL." ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Biometrika" "fecha" => "1981" "volumen" => "68" "paginaInicial" => "373" "paginaFinal" => "9" ] ] ] ] ] ] 12 => array:3 [ "identificador" => "bib13" "etiqueta" => "13" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "An additive frailty model for correlated life times." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "Petersen JH." ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:6 [ "tituloSerie" => "Biometrics" "fecha" => "1998" "volumen" => "54" "paginaInicial" => "646" "paginaFinal" => "61" "link" => array:1 [ 0 => array:2 [ "url" => "https://www.ncbi.nlm.nih.gov/pubmed/9660631" "web" => "Medline" ] ] ] ] ] ] ] ] 13 => array:3 [ "identificador" => "bib14" "etiqueta" => "14" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "A counting process approach to maximum likelihood estimation in frailty models." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:4 [ 0 => "Nielsen G" 1 => "Gill R" 2 => "Andersen P" 3 => "Sorensen T." ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Scan J Stat" "fecha" => "1992" "volumen" => "19" "paginaInicial" => "25" "paginaFinal" => "43" ] ] ] ] ] ] 14 => array:3 [ "identificador" => "bib15" "etiqueta" => "15" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "GNFCN: a computer program for gamma frailty analysis." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "Nielsen G" 1 => "Tabanera-Palacios R." ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:4 [ "titulo" => "GNFCN: a computer program for gamma frailty analysis." "fecha" => "1993" "paginaInicial" => "1" "paginaFinal" => "21" ] ] ] ] ] ] 15 => array:3 [ "identificador" => "bib16" "etiqueta" => "16" "referencia" => array:1 [ 0 => array:2 [ "referenciaCompleta" => "A modification of the EM algorithm to estimate an Andersen-Gill gamma frailty model for multivariate failure time data. Mimeo, 2001." "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "A modification of the EM algorithm to estimate an Andersen-Gill gamma frailty model for multivariate failure time data. Mimeo, 2001." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "Barceló MA" 1 => "Sáez M." ] ] ] ] ] ] ] ] 16 => array:3 [ "identificador" => "bib17" "etiqueta" => "17" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Generalized linear models for the analysis of Taguchi-type experiments." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "Nelder JA" 1 => "Lee Y." ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Applied Stochastic Models and Data Analysis" "fecha" => "1991" "volumen" => "7" "paginaInicial" => "107" "paginaFinal" => "20" ] ] ] ] ] ] 17 => array:3 [ "identificador" => "bib18" "etiqueta" => "18" "referencia" => array:1 [ 0 => array:2 [ "referenciaCompleta" => "Modeling survival data. Extending the Cox model. Statistics for biology and health. New York: Springer, 2000." "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Modeling survival data. Extending the Cox model. Statistics for biology and health. New York: Springer, 2000." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "Therneau TM" 1 => "Grambsch PM." ] ] ] ] ] ] ] ] 18 => array:3 [ "identificador" => "bib19" "etiqueta" => "19" "referencia" => array:1 [ 0 => array:2 [ "referenciaCompleta" => "Seattle: Data Analysis Products Division, MathSoft, 1999." "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Seattle: Data Analysis Products Division, MathSoft, 1999." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "S-PLUS 2000 User's Guide." ] ] ] ] ] ] ] ] 19 => array:3 [ "identificador" => "bib20" "etiqueta" => "20" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Discussion of paper by Wei and Glidden." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "Cook RJ" 1 => "Lawless JF." ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Stat Med" "fecha" => "1997" "volumen" => "16" "paginaInicial" => "841" "paginaFinal" => "3" ] ] ] ] ] ] 20 => array:3 [ "identificador" => "bib21" "etiqueta" => "21" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Discussion of paper by Wei and Glidden." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "Lipschutz KH" 1 => "Snapinn SM." ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Stat Med" "fecha" => "1997" "volumen" => "16" "paginaInicial" => "846" "paginaFinal" => "8" ] ] ] ] ] ] 21 => array:3 [ "identificador" => "bib22" "etiqueta" => "22" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Some approaches to the analysis of recurrent event data." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "Clayton D." ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Stat Med" "fecha" => "1994" "volumen" => "3" "paginaInicial" => "244" "paginaFinal" => "63" ] ] ] ] ] ] 22 => array:3 [ "identificador" => "bib23" "etiqueta" => "23" "referencia" => array:1 [ 0 => array:2 [ "referenciaCompleta" => "Frailty models for multiple event times. En: Klein YP, Goel PK, editors. Survival analysis: state of the art. Boston: Kluwer Academic Publishers, 1992; p. 415." "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Frailty models for multiple event times. En: Klein YP, Goel PK, editors. Survival analysis: state of the art. Boston: Kluwer Academic Publishers, 1992; p. 415." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "Oakes D." ] ] ] ] ] ] ] ] 23 => array:3 [ "identificador" => "bib24" "etiqueta" => "24" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:3 [ "titulo" => "Discussion of paper by Wei and Glidden." "idioma" => "en" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "Oakes D." ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:4 [ "tituloSerie" => "Stat Med" "fecha" => "1997" "volumen" => "16" "paginaInicial" => "843" ] ] ] ] ] ] ] ] ] ] ] "idiomaDefecto" => "es" "url" => "/02139111/00000016000000S2/v0_201302051443/13042666/v0_201302051445/es/main.assets" "Apartado" => array:4 [ "identificador" => "779" "tipo" => "SECCION" "es" => array:2 [ "titulo" => "Revisiones" "idiomaDefecto" => true ] "idiomaDefecto" => "es" ] "PDF" => "https://static.elsevier.es/multimedia/02139111/00000016000000S2/v0_201302051443/13042666/v0_201302051445/es/138v16nSupl.2a13042666pdf001.pdf?idApp=WGSE&text.app=https://gacetasanitaria.org/" "EPUB" => "https://multimedia.elsevier.es/PublicationsMultimediaV1/item/epub/13042666?idApp=WGSE"]

 array:20 [
  "pii" => "13042666"
  "issn" => "02139111"
  "estado" => "S300"
  "fechaPublicacion" => "2002-12-02"
  "documento" => "article"
  "crossmark" => 0
  "licencia" => "http://www.elsevier.com/open-access/userlicense/1.0/"
  "subdocumento" => "fla"
  "cita" => "Gac Sanit. 2002;16 Supl 2:59-68"
  "abierto" => array:3 [
    "ES" => true
    "ES2" => true
    "LATM" => true
  ]
  "gratuito" => true
  "lecturas" => array:2 [
    "total" => 6811
    "formatos" => array:3 [
      "EPUB" => 154
      "HTML" => 5667
      "PDF" => 990
    ]
  ]
  "itemSiguiente" => array:16 [
    "pii" => "13042667"
    "issn" => "02139111"
    "estado" => "S300"
    "fechaPublicacion" => "2002-12-02"
    "documento" => "article"
    "crossmark" => 0
    "licencia" => "http://www.elsevier.com/open-access/userlicense/1.0/"
    "subdocumento" => "fla"
    "cita" => "Gac Sanit. 2002;16 Supl 2:70-6"
    "abierto" => array:3 [
      "ES" => true
      "ES2" => true
      "LATM" => true
    ]
    "gratuito" => true
    "lecturas" => array:2 [
      "total" => 19787
      "formatos" => array:3 [
        "EPUB" => 232
        "HTML" => 18363
        "PDF" => 1192
      ]
    ]
    "es" => array:12 [
      "idiomaDefecto" => true
      "titulo" => "La publicidad del tabaco&#44; &#191;incrementa el h&#225;bito de fumar entre los adolescentes&#63; Datos desde California"
      "tienePdf" => "es"
      "tieneTextoCompleto" => "es"
      "tieneResumen" => array:2 [
        0 => "es"
        1 => "en"
      ]
      "paginas" => array:1 [
        0 => array:2 [
          "paginaInicial" => "70"
          "paginaFinal" => "76"
        ]
      ]
      "titulosAlternativos" => array:1 [
        "en" => array:1 [
          "titulo" => "Tobacco advertising&#44; does it increase smoking among teenagers&#63; Data from California"
        ]
      ]
      "contieneResumen" => array:2 [
        "es" => true
        "en" => true
      ]
      "contieneTextoCompleto" => array:1 [
        "es" => true
      ]
      "contienePdf" => array:1 [
        "es" => true
      ]
      "resumenGrafico" => array:2 [
        "original" => 0
        "multimedia" => array:8 [
          "identificador" => "fig1"
          "etiqueta" => "Figura 2"
          "tipo" => "MULTIMEDIAFIGURA"
          "mostrarFloat" => true
          "mostrarDisplay" => false
          "copyright" => "Elsevier Espa&#241;a"
          "figura" => array:1 [
            0 => array:4 [
              "imagen" => "138v16nSupl.2-13042667fig02.jpg"
              "Alto" => 443
              "Ancho" => 290
              "Tamanyo" => 17527
            ]
          ]
          "descripcion" => array:1 [
            "es" => "Anuncio de Lucky Strike de los a&#241;os veinte&#46;"
          ]
        ]
      ]
      "autores" => array:1 [
        0 => array:2 [
          "autoresLista" => "JP Pierce, M Garc&#237;a, E Salt&#243;, A Schiaffino, E Fern&#225;ndez"
          "autores" => array:5 [
            0 => array:2 [
              "Iniciales" => "JP"
              "apellidos" => "Pierce"
            ]
            1 => array:2 [
              "Iniciales" => "M"
              "apellidos" => "Garc&#237;a"
            ]
            2 => array:2 [
              "Iniciales" => "E"
              "apellidos" => "Salt&#243;"
            ]
            3 => array:2 [
              "Iniciales" => "A"
              "apellidos" => "Schiaffino"
            ]
            4 => array:2 [
              "Iniciales" => "E"
              "apellidos" => "Fern&#225;ndez"
            ]
          ]
        ]
      ]
    ]
    "idiomaDefecto" => "es"
    "EPUB" => "https://multimedia.elsevier.es/PublicationsMultimediaV1/item/epub/13042667?idApp=WGSE"
    "url" => "/02139111/00000016000000S2/v0_201302051443/13042667/v0_201302051446/es/main.assets"
  ]
  "itemAnterior" => array:16 [
    "pii" => "13042665"
    "issn" => "02139111"
    "estado" => "S300"
    "fechaPublicacion" => "2002-12-02"
    "documento" => "article"
    "crossmark" => 0
    "licencia" => "http://www.elsevier.com/open-access/userlicense/1.0/"
    "subdocumento" => "fla"
    "cita" => "Gac Sanit. 2002;16 Supl 2:41-58"
    "abierto" => array:3 [
      "ES" => true
      "ES2" => true
      "LATM" => true
    ]
    "gratuito" => true
    "lecturas" => array:2 [
      "total" => 5702
      "formatos" => array:3 [
        "EPUB" => 160
        "HTML" => 4990
        "PDF" => 552
      ]
    ]
    "es" => array:11 [
      "idiomaDefecto" => true
      "titulo" => "Comportamientos adictivos no saludables&#58; una revisi&#243;n del an&#225;lisis econ&#243;mico"
      "tienePdf" => "es"
      "tieneTextoCompleto" => "es"
      "tieneResumen" => array:2 [
        0 => "es"
        1 => "en"
      ]
      "paginas" => array:1 [
        0 => array:2 [
          "paginaInicial" => "41"
          "paginaFinal" => "58"
        ]
      ]
      "titulosAlternativos" => array:1 [
        "en" => array:1 [
          "titulo" => "Unhealthy addictive behaviours&#58; a review of economic analysis"
        ]
      ]
      "contieneResumen" => array:2 [
        "es" => true
        "en" => true
      ]
      "contieneTextoCompleto" => array:1 [
        "es" => true
      ]
      "contienePdf" => array:1 [
        "es" => true
      ]
      "autores" => array:1 [
        0 => array:2 [
          "autoresLista" => "F Portillo, F Anto&#241;anzas"
          "autores" => array:2 [
            0 => array:2 [
              "Iniciales" => "F"
              "apellidos" => "Portillo"
            ]
            1 => array:2 [
              "Iniciales" => "F"
              "apellidos" => "Anto&#241;anzas"
            ]
          ]
        ]
      ]
    ]
    "idiomaDefecto" => "es"
    "EPUB" => "https://multimedia.elsevier.es/PublicationsMultimediaV1/item/epub/13042665?idApp=WGSE"
    "url" => "/02139111/00000016000000S2/v0_201302051443/13042665/v0_201302051445/es/main.assets"
  ]
  "es" => array:14 [
    "idiomaDefecto" => true
    "titulo" => "Modelos marginales y condicionales en el an&#225;lisis de supervivencia multivariante"
    "tieneTextoCompleto" => true
    "paginas" => array:1 [
      0 => array:2 [
        "paginaInicial" => "59"
        "paginaFinal" => "68"
      ]
    ]
    "autores" => array:1 [
      0 => array:3 [
        "autoresLista" => "MA Barcel&#243;"
        "autores" => array:1 [
          0 => array:3 [
            "Iniciales" => "MA"
            "apellidos" => "Barcel&#243;"
            "referencia" => array:1 [
              0 => array:2 [
                "etiqueta" => "<span class="elsevierStyleSup">a</span>"
                "identificador" => "affa"
              ]
            ]
          ]
        ]
        "afiliaciones" => array:1 [
          0 => array:3 [
            "entidad" => "Grup de Recerca en Estadística, Economia Aplicada i Salut (GRECS). Departament d'Economia. Universitat de Girona. Centre de Recerca en Economia de la Salut (CRES). Universitat Pompeu Fabra. Barcelona."
            "etiqueta" => "<span class="elsevierStyleSup">a</span>"
            "identificador" => "affa"
          ]
        ]
      ]
    ]
    "titulosAlternativos" => array:1 [
      "en" => array:1 [
        "titulo" => "Marginal and conditional models in multivariate survival analysis"
      ]
    ]
    "textoCompleto" => "<p class="elsevierStylePara"><span class="elsevierStyleBold">Introducci&#243;n</span></p><p class="elsevierStylePara"> En los estudios de supervivencia se est&#225; interesado por el tiempo que transcurre desde un origen &#40;nacimiento&#44; inicio de tratamiento&#44; etc&#46;&#41; hasta que ocurre un suceso &#40;muerte&#44; mejora&#44; etc&#46;&#41;<span class="elsevierStyleSup">1</span>&#46; En numerosas ocasiones este suceso no es &#250;nico y los sujetos pueden experimentar varios sucesos &#40;o t&#233;cnicamente fallos&#41; que pueden ser de distinto tipo o bien tratarse de recurrencias del mismo suceso de inter&#233;s&#46; Debido a esto&#44; los individuos con m&#225;s sucesos estar&#225;n &#171;sobrerrepresentados&#187;<span class="elsevierStyleSup">2</span> y&#44; por otra parte&#44; los efectos de las variables explicativas &#40;efectos &#171;tratamiento&#187;&#41; podr&#225;n ser distintos entre el primer suceso y los siguientes&#46; Esto se traducir&#225; en una posible dependencia entre los tiempos de recurrencia&#44; consecuencia de dos fen&#243;menos diferentes&#58; la existencia de heterogeneidad individual y dependencia &#171;serial&#187;&#46; En la primera&#44; los tiempos de supervivencia dependen unos de otros porque los individuos comparten una o m&#225;s covariables omitidas&#46; Esta heterogeneidad individual tambi&#233;n se denomina fragilidad&#44; dependencia entre individuos o falso &#171;contagio&#187;&#46; Por lo que respecta a la dependencia &#171;serial&#187;&#44; que equivaldr&#237;a a la autocorrelaci&#243;n en series temporales&#44; la dependencia se produce en el mismo individuo&#44; entre sus observaciones repetidas&#46; Como consecuencia de que el individuo presenta una conducta din&#225;mica <span class="elsevierStyleItalic">per se</span>&#44; la probabilidad de que se produzca un suceso no es independiente de que se produzca uno previo&#46; Esta dependencia tambi&#233;n se denomina dependencia dentro de un mismo individuo o &#171;contagio&#187; verdadero&#46;</p><p class="elsevierStylePara"> Con el fin de ilustrar el trabajo utilizaremos una cohorte prospectiva formada por todos los pacientes admitidos un m&#237;nimo de 24 h en la Unidad de Cuidados Intensivos &#40;UCI&#41; del Hospital Universitario Dr&#46; Josep Trueta de Girona&#44; durante el per&#237;odo comprendido entre el 15 de marzo y el 15 de junio de 1999<span class="elsevierStyleSup">3</span>&#46; Se trata de un hospital p&#250;blico de nivel terciario&#44; con 12 camas en la UCI&#44; que cubre pr&#225;cticamente toda la poblaci&#243;n de la provincia de Girona&#46; Se pretende analizar los factores que determinan la incidencia de las infecciones nosocomiales en la UCI de dicho hospital&#46; Las infecciones nosocomiales vienen definidas como las causadas por microorganismos adquiridos en el hospital&#44; que afectan a los pacientes por un proceso diferente del de la infecci&#243;n&#44; y que no estaban presentes en la admisi&#243;n del paciente&#44; ni tan siquiera en per&#237;odo de incubaci&#243;n&#46; Entre los posibles factores de riesgo para la aparici&#243;n de dichas infecciones distinguiremos los asociados a los pacientes &#40;factores de riesgo intr&#237;nsecos&#44; como su edad o sexo&#41; y los relacionados con la UCI &#40;factores de riesgo extr&#237;nsecos&#44; como tratamiento antibi&#243;tico&#41;&#46; El n&#250;mero total de sujetos incluidos en el estudio fue de 61 pacientes &#40;42 varones y 19 mujeres&#41;&#44; con un total de 40 infecciones nosocomiales&#46; La variable dependiente &#40;episodio de infecci&#243;n&#41; consiste en el tiempo transcurrido entre el ingreso en la UCI hasta la aparici&#243;n de la infecci&#243;n &#40;originada por un microorganismo&#44; bacteria u hongo&#41;&#44; as&#237; como el tiempo entre el principio y el final de una infecci&#243;n&#46; Obviamente los pacientes podr&#225;n presentar m&#225;s de un episodio de infecci&#243;n&#46; S&#243;lo se analizaron infecciones nosocomiales originadas en la UCI&#44; incluy&#233;ndose en la definici&#243;n las infecciones comunitarias&#44; las procedentes de otros hospitales y las procedentes de otras &#225;reas hospitalarias&#46; Para una explicaci&#243;n m&#225;s detallada el lector interesado puede consultar en Barcel&#243; y Saez<span class="elsevierStyleSup">3</span>&#46;</p><p class="elsevierStylePara"> El objetivo principal de este art&#237;culo es presentar una revisi&#243;n de la aproximaci&#243;n marginal y condicional en el an&#225;lisis de supervivencia multivariante&#46; Adicionalmente&#44; y como ilustraci&#243;n&#44; se presenta una aplicaci&#243;n a datos reales que permite valorar la utilidad de tales aproximaciones&#46;</p><p class="elsevierStylePara"><span class="elsevierStyleBold">M&#233;todos</span></p><p class="elsevierStylePara"> Con el prop&#243;sito de solucionar los problemas ocasionados por la existencia de m&#225;s de un tiempo de supervivencia por individuo&#44; el modelo de Cox de riesgos proporcionales fue generalizado y situado en el contexto mucho m&#225;s general de los modelos de intensidad multiplicativos para procesos contadores multivariantes&#44; a trav&#233;s de la aproximaci&#243;n de Andersen y Gill &#40;modelo AG&#41;<span class="elsevierStyleSup">4-8</span>&#46; En la aproximaci&#243;n AG cada sujeto es tratado como un proceso contador con sucesos m&#250;ltiples y con &#171;incrementos&#187; &#40;es decir&#44; los tiempos entre acontecimientos sucesivos&#41; independientes&#44; dada la &#171;historia&#187; de todas las variables observables hasta el tiempo de presentaci&#243;n de los sucesos&#46; Es decir&#44; con el modelo AG el riesgo de una determinada recurrencia para un individuo sigue la hip&#243;tesis de riesgos proporcionales y no est&#225; afectado por la presentaci&#243;n en dicho individuo de sucesos anteriores&#46; Dicha hip&#243;tesis&#44; llamada de &#171;incrementos independientes&#187;&#44; es la hip&#243;tesis clave del modelo AG&#46; En caso de existir alguna interrelaci&#243;n entre sucesos se modelizar&#237;a mediante una o m&#225;s variables explicativas dependientes del tiempo<span class="elsevierStyleSup">9</span>&#46;</p><p class="elsevierStylePara"> Desde un punto de vista pr&#225;ctico&#44; la aproximaci&#243;n AG representa a cada individuo por un conjunto de observaciones s<span class="elsevierStyleInf">ij</span>&#44; t<span class="elsevierStyleInf">ij</span>&#44; d<span class="elsevierStyleInf">ij</span>&#44; x<span class="elsevierStyleInf">ij</span>&#44; k<span class="elsevierStyleInf">ij</span>&#44; j &#61; 1&#44; &#46;&#46;&#46;&#44; n<span class="elsevierStyleInf">i</span> siendo &#40;<span class="elsevierStyleItalic">s<span class="elsevierStyleInf">ij</span>&#44; t<span class="elsevierStyleInf">ij</span>&#93;</span> el intervalo de riesgo&#44; abierto a la izquierda y cerrado a la derecha&#59; d<span class="elsevierStyleInf">ij</span> &#61; 1 si el individuo tuvo un suceso en el instante t<span class="elsevierStyleInf">ij</span> y 0 en otro caso&#59; <span class="elsevierStyleItalic">x<span class="elsevierStyleInf">ij</span></span> es un vector de variables explicativas en el intervalo&#44; y <span class="elsevierStyleItalic">k<span class="elsevierStyleInf">ij</span></span> denota un estrato al que puede pertenecer el sujeto durante el intervalo&#46; Dependiendo de la escala de tiempo&#44; la primera observaci&#243;n puede o no empezar en cero&#46;</p><p class="elsevierStylePara"> Aclararemos lo anterior mediante un sencillo ejemplo&#46; Sean dos individuos&#44; un var&#243;n que padeci&#243; dos episodios de infecci&#243;n &#40;INF&#41; y una mujer que no present&#243; ninguno &#40;tabla 1 y fig&#46; 1a&#41;&#46; Los individuos se dividen en distintos registros que se corresponden con diversas bandas temporales &#40;intervalos de observaci&#243;n&#41; y cada registro es tratado como si fuesen individuos diferentes&#46; Las bandas se construyen en funci&#243;n de si el individuo presenta o no un episodio de infecci&#243;n y de posibles cambios en las variables explicativas&#44; no necesitando ser adyacentes&#46; En este dise&#241;o las variables explicativas son fijas dentro de cada episodio&#44; pero no as&#237; dentro de cada individuo &#40;tal como suponen los modelos est&#225;ndar&#44; como el modelo de Cox de riesgos proporcionales&#41;&#46; N&#243;tese que es el riesgo de cada suceso concreto para un individuo determinado el que cumple la hip&#243;tesis usual de riesgos proporcionales y que el modelo supone que este riesgo no est&#225; afectado por sucesos anteriores del mismo individuo&#46;</p><p class="elsevierStylePara"><img src="138v16nSupl.2-13042666tab01.gif"></img></p><p class="elsevierStylePara"><img src="138v16nSupl.2-13042666tab02.gif"></img></p><p class="elsevierStylePara"><span class="elsevierStyleBold">Figura 1&#46; Estructura de los datos a comparar seg&#250;n el modelo de an&#225;lisis&#46;</span></p><p class="elsevierStylePara"> La aproximaci&#243;n AG es sencilla&#44; pero sus hip&#243;tesis subyacentes son muy restrictivas y pueden llegar a ser irrealizables&#46; En concreto&#44; no es razonable suponer&#44; tal como hemos visto&#44; que las bandas pertenecientes a un mismo individuo no est&#233;n correlacionadas&#46; En la pr&#225;ctica esta limitaci&#243;n puede conducir a estimadores sesgados e ineficientes&#46; Adem&#225;s&#44; puede sobrestimarse la precisi&#243;n de los estimadores puesto que&#44; habitualmente&#44; las observaciones de un mismo individuo suelen estar correlacionadas positivamente&#46;</p><p class="elsevierStylePara"> Como posible soluci&#243;n&#44; a finales de los a&#241;os ochenta empezaron a aparecer los modelos marginales y condicionales en el an&#225;lisis de datos de supervivencia multivariante&#46; Sin embargo&#44; la aplicaci&#243;n de estos modelos aun hoy d&#237;a no est&#225; muy extendida debido&#44; principalmente&#44; a una falta de sistematizaci&#243;n en la bibliograf&#237;a especializada&#46; Es por ello por lo que consideramos especialmente importante describir estos modelos no s&#243;lo desde un punto de vista te&#243;rico<span class="elsevierStyleSup">10</span>&#44; sino tambi&#233;n ilustrarlos sobre la base de datos utilizada anteriormente&#46;</p><p class="elsevierStylePara"> En los modelos marginales la dependencia de las observaciones puede entenderse como un &#171;estorbo&#187;&#44; como algo que hay que eliminar o&#44; en todo caso&#44; cuyos efectos se deben controlar en el an&#225;lisis&#44; pero no como objetivo en s&#237;&#46; Por el contrario&#44; los modelos condicionales consideran que la dependencia puede ser tambi&#233;n de inter&#233;s por s&#237; misma&#46; A grandes rasgos&#44; los modelos marginales estiman el modelo ignorando la posible dependencia entre las observaciones&#44; aunque corrigen posteriormente la varianza &#171;ingenua&#187; mediante estimadores <span class="elsevierStyleItalic">jackknife&#44; bootstrap</span> o &#171;s&#225;ndwich&#187;&#46; Los modelos condicionales estiman la dependencia &#40;o la heterogeneidad&#41; especificando expl&#237;citamente su distribuci&#243;n de probabilidad e incorpor&#225;ndola en el modelo&#46;</p><p class="elsevierStylePara"> Modelos marginales</p><p class="elsevierStylePara"> La primera aproximaci&#243;n marginal que comentaremos es la estimaci&#243;n robusta de la matriz de varianzas y covarianzas&#46; Se parte del hecho de que la aproximaci&#243;n AG proporciona estimadores consistentes de los par&#225;metros asociados a las variables explicativas&#46; Aun as&#237;&#44; si existe dependencia&#44; el m&#233;todo proporciona varianzas sesgadas&#46; Se tratar&#237;a&#44; por tanto&#44; de estimar de forma robusta la matriz de covarianzas&#44; es decir&#44; de corregirla de forma adecuada utilizando&#44; por ejemplo&#44; un estimador <span class="elsevierStyleItalic">jackknife</span><span class="elsevierStyleSup">9</span>&#46; La idea es sencilla&#44; se trata de ir eliminando a un individuo de la muestra cada vez y restimar el modelo resultante&#46; Los estimadores finales se obtienen ponderando los obtenidos tras cada una de estas eliminaciones&#46; El problema de esta aproximaci&#243;n marginal estriba en que se basa en la consistencia de las estimaciones&#46; Sin embargo&#44; las estimaciones de los par&#225;metros ser&#225;n consistentes siempre y cuando el modelo marginal est&#233; correctamente especificado&#46;</p><p class="elsevierStylePara"> Otro modelo marginal es el modelo de Wei&#44; Lin y Weissfeld &#40;WLW&#41;<span class="elsevierStyleSup">11</span>&#44; en el que cada recurrencia del suceso &#40;o tipo distinto de suceso&#41; se modeliza como un estrato diferente&#44; utiliz&#225;ndose los datos dentro de cada estrato de modo marginal&#46; El resultado es que &#40;todos&#41; los individuos aparecen en todos los estratos&#46; Es decir&#44; el modelo WLW no es m&#225;s que un modelo estratificado&#44; siendo cada estrato una presentaci&#243;n del suceso de inter&#233;s&#46; Los estimadores finales son una media ponderada de los estimadores en cada estrato&#46; Como vemos&#44; se trata de un modelo marginal con respecto al conjunto de riesgo por cuanto el sujeto est&#225; en riesgo desde el inicio del estudio y&#44; por tanto&#44; no depende de sucesos previos&#46;</p><p class="elsevierStylePara"> Es decir&#44; para cada una de las <span class="elsevierStyleItalic">k</span> recurrencias que pueda padecer un individuo <span class="elsevierStyleItalic">i</span> &#40;k &#61; 1&#44;2&#44;&#46;&#46;&#46;&#44; K&#41;&#44; la funci&#243;n de riesgo puede expresarse como&#58;</p><p class="elsevierStylePara"><span class="elsevierStyleItalic">&#955;<span class="elsevierStyleInf">ik</span></span><span class="elsevierStyleItalic">&#40;t&#41;</span> &#61; <span class="elsevierStyleItalic">&#955;<span class="elsevierStyleInf">i</span></span><span class="elsevierStyleInf">0k</span><span class="elsevierStyleItalic">&#40;t&#41; e</span><span class="elsevierStyleItalic"><span class="elsevierStyleSup">&#946;&#39;</span></span><span class="elsevierStyleItalic"><span class="elsevierStyleInf">k</span><span class="elsevierStyleSup">X</span><span class="elsevierStyleInf">ik</span><span class="elsevierStyleSup">&#40;t&#41;</span> t</span> &#62;&#61; 0 &#91;1&#93;</p><p class="elsevierStylePara"> siendo <span class="elsevierStyleItalic">&#955;<span class="elsevierStyleInf">i0k</span></span><span class="elsevierStyleItalic">&#40;t&#41; &#40;k &#61; 1&#44;&#46;&#46;&#46;&#44; K&#41;</span> una funci&#243;n de riesgo subyacente no especificada&#44; distinta para cada recurrencia&#44; y <span class="elsevierStyleItalic">&#946;<span class="elsevierStyleInf">k</span> &#61; &#40;<span class="elsevierStyleInf">&#946;1k</span>&#44;&#46;&#46;&#46;&#44; <span class="elsevierStyleInf">  &#946;pk</span></span><span class="elsevierStyleItalic">&#41;&#39;</span> un vector de orden p &#215;1 de par&#225;metros desconocidos de la regresi&#243;n&#46;</p><p class="elsevierStylePara"> En el modelo WLW los datos se organizan de forma diferente de la que ve&#237;amos en el modelo AG &#40;tabla 1 y fig&#46; 1b&#44; en la que se ha supuesto que el m&#225;ximo n&#250;mero de episodios de infecci&#243;n entre todos los individuos analizados es de tres&#41;&#46; N&#243;tese que cada intervalo de tiempo empieza en cero&#46; Se trata de un modelo marginal con respecto al conjunto de riesgo por cuanto el sujeto est&#225; en riesgo desde el inicio del estudio y&#44; por tanto&#44; no depende de sucesos previos&#46;</p><p class="elsevierStylePara"> Las limitaciones del modelo WLW pueden ser superadas&#44; en parte&#44; utilizando la aproximaci&#243;n de Prentice&#44; Williams y Peterson &#40;PWP&#41;<span class="elsevierStyleSup">12</span>&#46; Se trata de un modelo marginal respecto a la estimaci&#243;n de los par&#225;metros&#44; pero condicional en relaci&#243;n con la construcci&#243;n del conjunto de individuos en riesgo&#46; El modelo PWP permite que el riesgo basal var&#237;e entre diferentes recurrencias&#46; Se trata&#44; por tanto&#44; de un modelo de riesgos proporcionales con estratos dependientes del tiempo&#44; en el que la dependencia entre recurrencias se controla estratificando por el n&#250;mero previo de presentaciones del suceso de inter&#233;s&#46; As&#237;&#44; a diferencia del m&#233;todo WLW&#44; en el que todos los individuos est&#225;n en riesgo en todo momento&#44; el modelo PWP incluye en el conjunto de riesgo para la recurrencia <span class="elsevierStyleItalic">k</span> &#250;nicamente a los individuos que han experimentado k-1 recurrencias&#46;</p><p class="elsevierStylePara"> Prentice&#44; Williams y Peterson<span class="elsevierStyleSup">12</span> proponen dos modelos semiparam&#233;tricos de riesgos proporcionales&#58;</p><p class="elsevierStylePara"><span class="elsevierStyleItalic">1&#46;</span> Tiempo desde el inicio del estudio &#40;de tiempo total&#41;&#46; Se incluye la funci&#243;n de riesgo subyacente como una funci&#243;n del tiempo desde el inicio del estudio</p><p class="elsevierStylePara"><span class="elsevierStyleItalic">&#955;</span><span class="elsevierStyleItalic"><span class="elsevierStyleInf">ik</span> &#40;t&#41;</span> &#61; <span class="elsevierStyleItalic">&#955;</span><span class="elsevierStyleInf">0<span class="elsevierStyleItalic">k</span></span><span class="elsevierStyleItalic">&#40;t&#41; e</span><span class="elsevierStyleItalic"><span class="elsevierStyleSup">&#946;&#39;</span></span><span class="elsevierStyleItalic"><span class="elsevierStyleInf">k</span><span class="elsevierStyleSup">X</span><span class="elsevierStyleInf">ik</span><span class="elsevierStyleSup"> &#40;t&#41;</span></span>&#91;2&#93;</p><p class="elsevierStylePara"><span class="elsevierStyleItalic">2&#46;</span> Tiempo desde la recurrencia inmediatamente precedente &#40;de brecha de tiempo&#41;&#46; La funci&#243;n de riesgo subyacente es una funci&#243;n del tiempo desde la recurrencia inmediatamente precedente al tiempo de fallo&#58;</p><p class="elsevierStylePara"><span class="elsevierStyleItalic">&#955;</span><span class="elsevierStyleItalic"><span class="elsevierStyleInf">ik</span> &#40;t&#41;</span> &#61; <span class="elsevierStyleItalic">&#955;</span><span class="elsevierStyleInf">0k</span> &#40;<span class="elsevierStyleItalic">t &#173; t<span class="elsevierStyleInf">n</span></span><span class="elsevierStyleInf">&#40;t&#41;</span>&#41;<span class="elsevierStyleItalic">e</span><span class="elsevierStyleItalic"><span class="elsevierStyleSup">&#946;&#39;</span></span><span class="elsevierStyleItalic"><span class="elsevierStyleInf">k</span><span class="elsevierStyleSup">X</span><span class="elsevierStyleInf">ik</span><span class="elsevierStyleSup"> &#40;t&#41;</span></span>&#91;3&#93;</p><p class="elsevierStylePara"> La particularidad de ambos modelos reside en que se estratifica en funci&#243;n del n&#250;mero de recurrencias previas&#46;</p><p class="elsevierStylePara"> Siguiendo con el ejemplo de antes&#44; tanto para el modelo AG como para PWP los datos se organizar&#237;an tal como se muestra en la tabla 1 y la figura 1c&#46;</p><p class="elsevierStylePara"> Modelos condicionales</p><p class="elsevierStylePara"> En algunos casos no s&#243;lo nos interesa el efecto marginal de las covariables sobre la funci&#243;n de riesgo de los tiempos de fallo de los individuos&#44; sino tambi&#233;n c&#243;mo los sucesos anteriores pueden influir en el riesgo de desarrollar futuros fallos&#46; Para responder a estas cuestiones necesitamos recurrir a los modelos condicionales&#46;</p><p class="elsevierStylePara"> En primer lugar&#44; la hip&#243;tesis de incrementos independientes en el modelo AG puede relajarse incluyendo en el modelo variables dependientes del tiempo&#44; como el n&#250;mero de recurrencias previas&#44; a fin de capturar la posible estructura de dependencia entre las recurrencias&#46; En concreto&#44; se podr&#237;a trabajar con el modelo autorregresivo o markoviano&#46; Este modelo ser&#225; plenamente apropiado cuando la dependencia entre observaciones sea&#44; principalmente&#44; &#171;serial&#187;&#46; Sin embargo&#44; y aunque es un modelo muy general&#44; supone que la dependencia entre recurrencias decae exponencialmente y que la dependencia es id&#233;ntica para todos los individuos<span class="elsevierStyleSup">10</span>&#46;</p><p class="elsevierStylePara"> Parte de las limitaciones anteriores pueden solucionarse mediante los modelos de fragilidad&#46; Se trata de introducir una &#40;modelos de fragilidad univariante&#41; o m&#225;s &#40;modelos de fragilidad multivariante&#41; variables aleatorias comunes a todos los individuos con el objetivo de permitir que los tiempos de supervivencia multivariantes&#44; es decir&#44; las observaciones repetidas de un mismo individuo&#44; est&#233;n correlacionadas&#46; Se trata&#44; en definitiva&#44; de que este efecto multiplicativo en la funci&#243;n de intensidad &#40;o de riesgo&#41; capture dos fuentes de variaci&#243;n distintas pero relacionadas<span class="elsevierStyleSup">13</span>&#44; con el prop&#243;sito de modelar la dependencia entre los diferentes tiempos de supervivencia entre recurrencias&#46; Una de las fuentes de variaci&#243;n surge de covariables comunes a todos los individuos&#44; no observadas&#44; las cuales&#44; cuando se omiten&#44; generan dependencia entre sucesos&#46; La otra fuente de variaci&#243;n proviene de covariables individuales no observadas y&#44; por tanto&#44; no incluidas en el estudio&#44; bien por circunstancias pr&#225;cticas&#44; bien porque no se asocian a ning&#250;n factor de riesgo de la supervivencia del individuo&#46; Los efectos de esta fuente ser&#237;an totalmente equivalentes a la dependencia serial citada en la introducci&#243;n de este art&#237;culo&#46;</p><p class="elsevierStylePara"> Posibles alternativas para la distribuci&#243;n de probabilidad de la variable aleatoria fragilidad son la distribuci&#243;n gamma&#44; la distribuci&#243;n positiva estable y la distribuci&#243;n gaussiana inversa&#44; entre otras&#46; El modelo de fragilidad gamma es el m&#225;s adecuado en nuestro caso&#44; debido&#44; en primer lugar&#44; al tipo de datos de que disponemos&#46; En concreto&#44; bajo este tipo de distribuci&#243;n la fuerza de la asociaci&#243;n entre recurrencias es estacionaria&#44; es decir&#44; no depende del tiempo&#46; Este tipo de asociaci&#243;n es la que m&#225;s se produce en el mundo real&#46; La ventaja m&#225;s importante de la distribuci&#243;n gamma&#44; sin embargo&#44; se encuentra en la facilidad anal&#237;tica que permite&#46;</p><p class="elsevierStylePara"> En particular&#44; el modelo AG de fragilidad<span class="elsevierStyleSup">14&#44;15</span> puede expresarse de la siguiente manera en su caso m&#225;s simple&#58;</p><p class="elsevierStylePara"><span class="elsevierStyleItalic">&#955;<span class="elsevierStyleInf">i</span></span><span class="elsevierStyleItalic">&#40;t&#41;</span> &#61; &#969;<span class="elsevierStyleItalic"><span class="elsevierStyleInf">i</span></span><span class="elsevierStyleItalic">Y<span class="elsevierStyleInf">i</span> &#40;t&#41;</span> &#945;<span class="elsevierStyleInf">0</span><span class="elsevierStyleItalic">&#40;t&#41; e</span><span class="elsevierStyleItalic"><span class="elsevierStyleSup">&#946;&#39;X</span></span><span class="elsevierStyleItalic"><span class="elsevierStyleSup">i &#40;t&#41;</span></span> &#91;4&#93;</p><p class="elsevierStylePara"> siendo <span class="elsevierStyleItalic">N &#61; &#40;N<span class="elsevierStyleInf">i</span>&#59; i &#61; 1&#44; &#46;&#46;&#46;&#44; n&#41;</span> un proceso contador multivariante que cuenta el n&#250;mero de sucesos que han ocurrido en cada uno de los <span class="elsevierStyleItalic">n</span> individuos en el tiempo <span class="elsevierStyleItalic">t</span>&#46; Asociado a este proceso contador&#44; existe un proceso de intensidad <span class="elsevierStyleItalic">&#955;</span>con componentes <span class="elsevierStyleItalic">&#955;<span class="elsevierStyleInf">i</span></span> &#46; El proceso Y con componentes Y<span class="elsevierStyleInf">i</span> es un proceso no negativo&#44; observable y predecible&#44; es decir&#44; su valor en cualquier momento <span class="elsevierStyleItalic">t</span> es conocido justo antes de <span class="elsevierStyleItalic">t</span>&#44; indicando si un individuo <span class="elsevierStyleItalic"> i</span> se observa en el conjunto a riesgo de experimentar un determinado suceso en el tiempo <span class="elsevierStyleItalic">t-</span>&#59; &#945;<span class="elsevierStyleInf">0</span> denota ahora una funci&#243;n de riesgo basal desconocida&#44; y &#969;<span class="elsevierStyleItalic"><span class="elsevierStyleInf">i</span></span> &#40;fragilidad&#41; es un conjunto de variables aleatorias independientes e id&#233;nticamente distribuidas seg&#250;n una distribuci&#243;n de probabilidad param&#233;trica determinada&#44; en este caso la distribuci&#243;n gamma&#44; con media nula y varianza com&#250;n&#46; La variabilidad &#40;medida por su varianza&#41; de la fragilidad determina el grado de heterogeneidad&#46; De hecho&#44; cuando la varianza es igual a cero el modelo se reduce al modelo AG est&#225;ndar&#46;</p><p class="elsevierStylePara"> El modelo descrito es un modelo muy simple&#44; ya que asume que la fragilidad es constante a lo largo del tiempo&#44; es decir&#44; la fragilidad reflejar&#237;a en este caso diferencias individuales presentes al comienzo del estudio&#46; Asimismo&#44; este modelo supone que la fragilidad act&#250;a proporcionalmente sobre la funci&#243;n basal com&#250;n a todos los individuos&#46; Al imponer esto se asume que todos los individuos siguen un patr&#243;n uniforme&#46;</p><p class="elsevierStylePara"> El problema es que la estimaci&#243;n de los par&#225;metros del modelo utilizando directamente la m&#225;xima verosimilitud no es posible cuando la funci&#243;n de riesgo basal no est&#225; parametrizada&#44; como en nuestro caso&#46; Por otra parte&#44; los modelos de fragilidad univariantes presentan una serie de limitaciones&#59; entre ellas cabe destacar las siguientes&#58; suponen que la dependencia entre observaciones es siempre positiva&#44; no permiten la interacci&#243;n entre covariables observables y la fragilidad&#44; y m&#225;s importante a&#250;n&#44; s&#243;lo permiten la presencia de un &#250;nico efecto aleatorio&#44; siendo adem&#225;s &#233;ste constante a lo largo del tiempo&#46; Con el prop&#243;sito de solucionar estos problemas&#44; Barcel&#243; y Saez<span class="elsevierStyleSup">16</span> proponen un m&#233;todo de especificaci&#243;n y estimaci&#243;n original basado en la verosimilitud extendida de Nelder y Lee<span class="elsevierStyleSup">17</span>&#44; as&#237; como en el m&#233;todo de las variables instrumentales&#44; y estudian las ventajas anal&#237;ticas de dicho m&#233;todo&#44; en particular&#44; en relaci&#243;n con alternativas de especificaci&#243;n &#40;modelos de fragilidad multivariantes&#41; y de estimaci&#243;n &#40;algoritmo EM&#41;&#46; De una manera muy resumida&#44; introducen un nuevo paso en el algoritmo EM&#44; modelizando la varianza de la fragilidad directamente&#46; Mediante el algoritmo as&#237; modificado &#40;algoritmo EMB&#41;&#44; se estima no s&#243;lo el modelo de supervivencia&#44; sino tambi&#233;n la varianza de la fragilidad de forma simult&#225;nea&#46; Mediante esta modificaci&#243;n solucionan parte de las limitaciones citadas anteriormente&#59; en particular&#44; permiten que la fragilidad no sea constante entre individuos&#46; Una explicaci&#243;n m&#225;s detallada puede encontrarse en Barcel&#243; y Saez<span class="elsevierStyleSup">16</span>&#46;</p><p class="elsevierStylePara"> Una alternativa a los modelos de fragilidad la proporcionan los modelos de Cox penalizados<span class="elsevierStyleSup">18</span>&#46; En resumen&#44; se trata de maximizar una funci&#243;n de verosimilitud parcial penalizada&#46; El t&#233;rmino penalizador captura la variabilidad local presente en la densidad conjunta de los datos&#46;</p><p class="elsevierStylePara"> Los modelos que se han comentado no est&#225;n anidados&#44; por lo que las usuales medidas de bondad de ajuste no pueden utilizarse&#46; Como alternativa se suele emplear el criterio de informaci&#243;n de Akaike &#40;AIC&#41;&#46;</p><p class="elsevierStylePara"> El software estad&#237;stico que permite ajustar todos los modelos vistos y que&#44; por otra parte&#44; se ha utilizado en la aplicaci&#243;n pr&#225;ctica&#44; es el S-PLUS 2000<span class="elsevierStyleSup">19</span>&#46; En concreto&#44; se debe utilizar la instrucci&#243;n <span class="elsevierStyleItalic">coxph</span>&#44; que permite considerar diferentes tipos de truncamiento y censura&#46; La estimaci&#243;n robusta de la varianza se obtiene mediante la instrucci&#243;n <span class="elsevierStyleItalic">cluster</span>&#46; Desafortunadamente&#44; hay pocos programas que permitan la estimaci&#243;n de modelos marginales y&#47;o condicionales en el contexto multivariante&#46; Adem&#225;s del S-PLUS 2000&#44; el lector interesado puede recurrir a otros programas como R &#40;versi&#243;n libre y gratuita de S-PLUS&#41;&#44; SAS &#40;versi&#243;n 6 o posteriores&#44; instrucciones <span class="elsevierStyleItalic">proc phreg</span>&#44; aunque no permite la estimaci&#243;n robusta de la varianza&#41; y&#44; &#250;nicamente para el modelo AG est&#225;ndar&#44; STATA &#40;versi&#243;n 5 o posteriores&#44; instrucciones <span class="elsevierStyleItalic">stset</span> y <span class="elsevierStyleItalic"> stcox</span>&#41;&#46;</p><p class="elsevierStylePara"><span class="elsevierStyleBold">Resultados</span></p><p class="elsevierStylePara"> Los resultados de la estimaci&#243;n del modelo AG est&#225;ndar aparecen en las primeras columnas de la tabla 2&#46; En esta tabla se presentan las tasas de riesgo asociadas a toda una serie de variables explicativas y sus respectivos intervalos de confianza &#40;IC&#41;&#46; Por &#250;ltimo&#44; se ofrece una serie de medidas de bondad de ajuste del modelo&#46; Como puede apreciarse al observar la tabla&#44; pocas variables resultaron estar asociadas de forma significativa con el tiempo de presentaci&#243;n de una infecci&#243;n nosocomial&#46; En concreto&#44; ser mujer &#40;0&#44;258&#59; IC del 95&#37;&#44; 0&#44;067-0&#44;997&#41;&#44; haberse sometido a cirug&#237;a urgente &#40;0&#44;295&#59; IC del 95&#37;&#44; 0&#44;092-0&#44;949&#41;&#44; haber permanecido intubado m&#225;s de 6 d&#237;as mediante traqueotom&#237;a &#40;0&#44;182&#59; IC del 95&#37;&#44; 0&#44;041-0&#44;814&#41; y haber recibido tratamiento antibi&#243;tico durante m&#225;s de 7 d&#237;as &#40;0&#44;179&#59; IC del 95&#37;&#44; 0&#44;034-0&#44;947&#41; reduc&#237;an el riesgo de tales infecciones&#46;</p><p class="elsevierStylePara"><img src="138v16nSupl.2-13042666tab03.gif"></img></p><p class="elsevierStylePara"> Los resultados de la aplicaci&#243;n pr&#225;ctica del modelo AG con estimaci&#243;n robusta de la matriz de varianzas y covarianzas ser&#237;an los que hemos visto anteriormente al hablar del modelo AG&#44; ya que esta primera aproximaci&#243;n marginal s&#243;lo afecta a la matriz de varianzas y covarianzas&#44; no a la construcci&#243;n de los conjuntos a riesgo&#46; Si nos fijamos en la tabla 2 y comparamos el modelo AG est&#225;ndar y el modelo AG robusto&#44; apreciamos una mayor precisi&#243;n en los intervalos de confianza proporcionados por este &#250;ltimo modelo&#44; siendo este cambio muy significativo en variables como ventilaci&#243;n mec&#225;nica&#46; Por el contrario&#44; otras variables tales como cirug&#237;a urgente &#40;0&#44;295&#59; IC del 95&#37;&#44; 0&#44;083-1&#44;048&#41; y tratamiento antibi&#243;tico durante m&#225;s de 7 d&#237;as &#40;0&#44;179&#59; IC del 95&#37;&#44; 0&#44;031-1&#44;047&#41; llegan incluso a perder su significaci&#243;n estad&#237;stica&#46;</p><p class="elsevierStylePara"> En relaci&#243;n con los resultados del modelo PWP con estimaci&#243;n robusta de la varianza &#40;tabla 3&#41;&#44; con una importante excepci&#243;n&#44; presencia de cat&#233;ter venoso central durante m&#225;s de tres d&#237;as &#40;1&#44;065&#59; IC del 95&#37;&#44; 1&#44;007-1&#44;775&#41;&#44; el signo de los estimadores de los par&#225;metros coincide en los tres modelos vistos &#40;modelo AG est&#225;ndar&#44; modelo AG robusto y modelo PWP&#41;&#44; y no ocurre lo mismo con la magnitud ni con la significaci&#243;n de los par&#225;metros&#46; N&#243;tese ahora que variables como la edad &#40;1&#44;009&#59; IC del 95&#37;&#44; 0&#44;985-1&#44;033&#41;&#44; la existencia previa de infecciones comunitarias &#40;4&#44;56&#59; IC del 95&#37;&#44; 1&#44;435-14&#44;477&#41;&#44; estar inestable en UCI &#40;54&#44;47&#59; IC del 95&#37;&#44; 0&#44;786-3&#46;776&#44;406&#41;&#44; estar inestable en coma o shock &#40;219&#44;97&#44; IC del 95&#37;&#44; 3&#44;187-15&#46;182&#44;609&#41;&#44; la presencia de ventilaci&#243;n mec&#225;nica de cuatro a 10 d&#237;as &#40;61&#44;93&#44; IC del 95&#37;&#44; 1&#44;325-2&#46;894&#44;463&#41;&#44; la presencia de ventilaci&#243;n mec&#225;nica durante m&#225;s de 10 d&#237;as &#40;129&#44;53&#59; IC del 95&#37;&#44; 4&#44;005-4&#46;189&#44;427&#41; y la presencia de cat&#233;ter venoso central durante m&#225;s de tres d&#237;as &#40;1&#44;07&#59; IC del 95&#37;&#44; 1&#44;007-1&#44;775&#41; aumentan el riesgo de padecer un episodio de infecci&#243;n nosocomial&#46; Contrariamente&#44; la presencia de sondas nasog&#225;stricas durante m&#225;s de 9 d&#237;as &#40;0&#44;13&#59; IC del 95&#37;&#44; 0&#44;023-0&#44;748&#41; y el tratamiento antibi&#243;tico durante m&#225;s de 7 d&#237;as &#40;0&#44;04&#59; IC del 95&#37;&#44; 0&#44;004-0&#44;385&#41; lo reducen&#46; Variables tales como g&#233;nero y traqueotom&#237;a pierden en este modelo su significaci&#243;n&#46; Resumiendo&#44; no s&#243;lo los factores de riesgo intr&#237;nsecos &#40;como la gravedad de la enfermedad&#41;&#44; sino tambi&#233;n los extr&#237;nsecos &#40;procedimientos terap&#233;uticos agresivos utilizados en la UCI&#41; llevaron a un mayor riesgo de infecci&#243;n nosocomial durante la estancia del paciente en la UCI&#46;</p><p class="elsevierStylePara"><img src="138v16nSupl.2-13042666tab04.gif"></img></p><p class="elsevierStylePara"> Los resultados del modelo AG de fragilidad gamma estimado mediante el algoritmo EMB&#44; as&#237; como los del modelo de Cox penalizado&#44; tambi&#233;n se encuentran en la tabla 3&#46; En el modelo AG de fragilidad gamma estimado mediante el algoritmo EMB&#44; observamos que variables tales como la edad &#40;1&#44;009&#59; IC del 95&#37;&#44; 0&#44;990-1&#44;028&#41;&#44; existencia previa de infecciones comunitarias &#40;6&#44;923&#59; IC del 95&#37;&#44; 2&#44;578-18&#44;594&#41;&#44; estar inestable en UCI &#40;49&#44;012&#59; IC del 95&#37;&#44; 1&#44;111-2&#46;162&#44;135&#41;&#44; estar inestable en coma o shock &#40;257&#44;040&#59; IC del 95&#37;&#44; 5&#44;997-11&#46;017&#44;03&#41;&#44; la presencia de ventilaci&#243;n mec&#225;nica en todas sus categor&#237;as y la presencia de cat&#233;ter venoso central durante m&#225;s de tres d&#237;as &#40;1&#44;082&#59; IC del 95&#37;&#44; 1&#44;011-1&#44;822&#41; contin&#250;an aumentando el riesgo de padecer infecci&#243;n nosocomial de forma estad&#237;sticamente significativa&#46; Por el contrario&#44; la presencia de sonda nasog&#225;strica durante m&#225;s de 9 d&#237;as &#40;0&#44;129&#59; IC del 95&#37;&#44; 0&#44;03-0&#44;55&#41; y el tratamiento antibi&#243;tico en todas sus categor&#237;as lo reducen&#44; tambi&#233;n de forma estad&#237;sticamente significativa&#46; Por otra parte&#44; se observa que variables que no eran significativas en el modelo PWP&#44; como la existencia previa de otras infecciones&#44; traqueotom&#237;a y presencia de sonda urinaria de 5 a 12 d&#237;as&#44; pasan a serlo ahora&#46; Estos resultados cambian un poco al observar el modelo de Cox penalizado&#46; Mientras que variables tales como la edad &#40;1&#44;009&#59; IC del 90&#37;&#58; 0&#44;981-1&#44;037&#41;&#44; infecciones previas comunitarias &#40;4&#44;558&#59; IC del 90&#37;&#44; 0&#44;954-21&#44;772&#41;&#44; estar inestable en UCI &#40;8&#46;025&#44;27&#59; IC del 95&#37;&#44; 1&#46;740-36&#46;968&#44;77&#41;&#44; presencia de ventilaci&#243;n mec&#225;nica de cuatro a 10 d&#237;as &#40;64&#44;406&#59; IC del 95&#37;&#44; 1&#44;35-3&#46;065&#44;471&#41;&#44; presencia de ventilaci&#243;n mec&#225;nica durante m&#225;s de 10 d&#237;as &#40;134&#44;714&#59; IC del 95&#37;&#44; 3&#44;06-5&#46;934&#44;619&#41;&#44; presencia de sonda nasog&#225;strica durante m&#225;s de 9 d&#237;as &#40;0&#44;133&#59; IC del 95&#37;&#44; 0&#44;027-0&#44;649&#41; y tratamiento antibi&#243;tico durante m&#225;s de 7 d&#237;as &#40;0&#44;038&#59; IC del 95&#37;&#44; 0&#44;005-0&#44;327&#41; contin&#250;an manteniendo su significaci&#243;n&#44; otras como traqueotom&#237;a&#44; existencia de otras infecciones previas&#44; presencia de sonda urinaria de 5 a 12 d&#237;as y tratamiento antibi&#243;tico durante 7 d&#237;as o menos la pierden&#46; Por otra parte&#44; el modelo de Cox penalizado parece que es incapaz de estimar los par&#225;metros en algunas de las variables explicativas presentes en el modelo &#40;inestable en coma o shock&#44; ventilaci&#243;n mec&#225;nica inferior o igual a tres d&#237;as y cat&#233;ter venoso central&#41;&#44; lo cual no ocurre en el modelo EMB&#46; Asimismo&#44; el modelo EMB es el que proporciona intervalos de confianza para los estimadores de los par&#225;metros m&#225;s reducidos&#46;</p><p class="elsevierStylePara"> Adicionalmente&#44; y a diferencia de los dem&#225;s modelos&#44; el modelo EMB permite observar los factores explicativos de la fragilidad&#46; Como se ve&#44; la varianza de la fragilidad depende &#250;nicamente de factores extr&#237;nsecos tales como ventilaci&#243;n mec&#225;nica&#44; traqueotom&#237;a&#44; sonda urinaria y sonda nasog&#225;strica&#46;</p><p class="elsevierStylePara"><span class="elsevierStyleBold">Discusi&#243;n</span></p><p class="elsevierStylePara"> Como ya se ha dicho&#44; la problem&#225;tica planteada por la existencia de datos de supervivencia multivariantes hace necesaria la generalizaci&#243;n del modelo de Cox de riesgos proporcionales&#44; concret&#225;ndose &#233;sta en el modelo AG&#46; Este modelo permite trabajar con datos multivariantes&#59; sin embargo&#44; no contempla la posibilidad que las observaciones pertenecientes a un mismo individuo est&#233;n correlacionadas&#44; ni la presencia de heterogeneidad individual&#46; Con el prop&#243;sito de solucionar esto&#44; citamos los principales modelos marginales &#40;modelo AG con estimaci&#243;n robusta de la varianza&#44; modelo WLW y modelo PWP&#41; y condicionales &#40;modelo autorregresivo o markoviano&#44; modelos de fragilidad y modelo de Cox penalizado&#41;&#46;</p><p class="elsevierStylePara"> Aun cuando los modelos marginales son capaces de corregir los principales efectos de la dependencia entre observaciones &#40;siempre suponiendo que las estimaciones de los par&#225;metros son consistentes&#41;&#44; son incapaces de distinguir c&#243;mo los sucesos anteriores pueden influir en el riesgo de desarrollar futuros fallos&#46; Esto &#250;ltimo es tratado mediante los modelos condicionales&#46;</p><p class="elsevierStylePara"> El m&#233;todo WLW presenta importantes desventajas&#46; Cook y Lawless<span class="elsevierStyleSup">20</span> se&#241;alan que el modelo WLW es v&#225;lido &#250;nicamente bajo la posiblemente muy restrictiva hip&#243;tesis de censura independiente &#40;restrictiva como consecuencia de las recurrencias&#41;&#46; La desventaja m&#225;s importante es que el m&#233;todo WLW supone que en el instante <span class="elsevierStyleItalic">t</span> todos los individuos &#40;observados en <span class="elsevierStyleItalic">t</span>&#41; constituyen el conjunto de riesgo para la correspondiente recurrencia&#44; con independencia del n&#250;mero de recurrencias previas que tuviese el individuo<span class="elsevierStyleSup">20&#44;21</span>&#46; Esto no representar&#237;a ning&#250;n problema si se tratase de diferentes sucesos de un tipo totalmente distinto&#44; pero si se trata de recurrencias del mismo suceso podr&#237;a ocurrir que alg&#250;n hecho emprendido con el fin de reducir el n&#250;mero de tales recurrencias redujese el n&#250;mero de las mismas a partir de su implementaci&#243;n &#40;lo mismo ocurrir&#237;a en el caso de <span class="elsevierStyleItalic"> harvesting</span>&#41;&#46; Como consecuencia aumentar&#237;a la correlaci&#243;n entre los estimadores de los par&#225;metros&#46; Otro problema es que las estimaciones de los par&#225;metros ser&#225;n consistentes siempre y cuando la funci&#243;n de riesgo est&#233; correctamente especificada&#46;</p><p class="elsevierStylePara"> La principal desventaja del modelo PWP radica en la manera de construir los conjuntos a riesgo&#44; ya que lleva a una &#171;p&#233;rdida de aleatorizaci&#243;n&#187;<span class="elsevierStyleSup">22</span>&#44; por cuanto los individuos con m&#225;s riesgo van saliendo del estudio&#46; Esta p&#233;rdida ser&#225; mayor conforme vayamos movi&#233;ndonos a estratos superiores &#40;y cuantos m&#225;s estratos se consideren&#41;&#46; En estos casos&#44; la consistencia de las estimaciones depende de la inclusi&#243;n de todas las variables relevantes&#46; Por otra parte&#44; tanto el modelo AG como el modelo PWP son sensibles a los errores de especificaci&#243;n de la estructura de dependencia entre los tiempos de recurrencia&#46;</p><p class="elsevierStylePara"> Aunque los m&#233;todos marginales cumplen &#40;marginalmente&#41; la hip&#243;tesis de riesgos proporcionales&#44; suponen expl&#237;cita &#40;PWP&#41; o impl&#237;citamente &#40;WLW&#41; que la funci&#243;n de riesgo basal var&#237;a de una recurrencia a otra<span class="elsevierStyleSup">22</span>&#46; Como consecuencia&#44; son extremadamente sensibles al cumplimiento &#40;en realidad&#41; de la hip&#243;tesis de riesgos proporcionales<span class="elsevierStyleSup">9</span>&#46; Finalmente&#44; y respecto a las aproximaciones condicionales&#44; los m&#233;todos marginales no son totalmente eficientes<span class="elsevierStyleSup">23&#44;24</span>&#46; La ineficiencia aumentar&#225; de forma proporcional al n&#250;mero de estratos &#40;recurrencias&#41; e inversamente proporcional al n&#250;mero de observaciones en cada estrato&#46; En analog&#237;a a los modelos lineales es razonable suponer que la ineficiencia tambi&#233;n aumentar&#225; conforme lo haga la correlaci&#243;n entre recurrencias&#46; Es por tanto de esperar que WLW sea menos eficiente que PWP&#46;</p><p class="elsevierStylePara"> Por lo que hace referencia a los modelos condicionales&#44; s&#243;lo comentaremos aqu&#237; el problema en los modelos de Cox penalizados&#46; Como argument&#225;bamos anteriormente&#44; la variabilidad que pretenden capturar estos modelos recoge&#44; en realidad&#44; dos fen&#243;menos diferentes&#44; fragilidad y correlaci&#243;n serial&#46; Por otra parte&#44; presentan el inconveniente de que el par&#225;metro de suavizado debe escogerse utilizando t&#233;cnicas no del todo validadas en el campo de supervivencia&#46;</p><p class="elsevierStylePara"> Situ&#225;ndonos en nuestro ejemplo pr&#225;ctico&#44; de entre los tres modelos marginales presentados en el art&#237;culo&#44; el modelo PWP es posiblemente el mejor debido al tipo de datos de que disponemos&#46; El hecho de que&#44; tras experimentar una primera infecci&#243;n&#44; aumente el riesgo de aparici&#243;n de nuevas infecciones hace que el modelo PWP con estimaci&#243;n robusta de la varianza sea especialmente apropiado&#44; ya que ser&#237;a posible estratificar los datos para cada episodio de infecci&#243;n&#46;</p><p class="elsevierStylePara"> Por lo que respecta a las estimaciones de los modelos condicionales mencionados en el art&#237;culo&#44; el algoritmo EMB de estimaci&#243;n de un modelo AG de fragilidad gamma es el que presenta mejores resultados&#46; Por un lado&#44; es el m&#225;s eficiente al presentar intervalos de confianza m&#225;s reducidos&#44; debido quiz&#225; a que el algoritmo EMB es capaz de capturar no s&#243;lo la heterogeneidad individual&#44; sino tambi&#233;n la dependencia serial&#46; Por otro lado&#44; de entre todos lo modelos estudiados&#44; es el que proporciona un mejor ajuste en t&#233;rminos del criterio de Akaike&#46; Finalmente&#44; es el &#250;nico que permite estimar los efectos explicativos de la fragilidad&#46;</p><p class="elsevierStylePara"><span class="elsevierStyleBold">Agradecimientos</span></p><p class="elsevierStylePara"> Agradezco al Dr&#46; Marc Saez su gran ayuda en la elaboraci&#243;n de este trabajo&#46; Agradezco&#44; asimismo&#44; al Dr&#46; Santiago P&#233;rez Hoyos y a dos evaluadores an&#243;nimos sus sugerencias&#46; Este trabajo fue parcialmente financiado por el proyecto 07&#47;50&#47;1998 de la Ag&#232;ncia d&#39;Avaluaci&#243; de Tecnologia M&#232;dica &#40;AATM&#41; del Servei Catal&#224; de la Salut&#44; Generalitat de Catalunya&#46;</p>"
    "pdfFichero" => "138v16nSupl.2a13042666pdf001.pdf"
    "tienePdf" => true
    "PalabrasClave" => array:2 [
      "es" => array:1 [
        0 => array:4 [
          "clase" => "keyword"
          "titulo" => "Palabras clave"
          "identificador" => "xpalclavsec16355"
          "palabras" => array:4 [
            0 => "An&#225;lisis de supervivencia multivariante"
            1 => "Enfoque marginal"
            2 => "Enfoque condicional"
            3 => "Modelos de fragilidad"
          ]
        ]
      ]
      "en" => array:1 [
        0 => array:4 [
          "clase" => "keyword"
          "titulo" => "Keywords"
          "identificador" => "xpalclavsec16356"
          "palabras" => array:4 [
            0 => "Multivariate survival data"
            1 => "Marginal models"
            2 => "Conditional models"
            3 => "Frailty models"
          ]
        ]
      ]
    ]
    "tieneResumen" => true
    "resumen" => array:2 [
      "es" => array:1 [
        "resumen" => "El hecho de que los individuos puedan experimentar varios sucesos &#40;fallos&#41; durante el per&#237;odo de observaci&#243;n &#40;an&#225;lisis de supervivencia multivariante&#41; hace que sea necesaria la aplicaci&#243;n de una metodolog&#237;a diferente de la empleada en el an&#225;lisis de supervivencia est&#225;ndar univariante&#46; Al analizar medidas repetidas&#44; el principal problema del modelo de Cox se debe a que las observaciones no son independientes&#44; lo cual llevar&#225; al incumplimiento de la hip&#243;tesis de riesgos proporcionales&#46; Con el prop&#243;sito de solucionar este problema&#44; as&#237; como otros que puedan darse debido a este tipo de datos&#44; aparece el modelo de Andersen-Gill &#40;AG&#41;&#44; que no es m&#225;s que una generalizaci&#243;n del modelo de Cox&#46; Sin embargo&#44; la hip&#243;tesis clave de este modelo &#40;hip&#243;tesis de incrementos independientes&#41; supone que las m&#250;ltiples observaciones de un mismo individuo son independientes&#44; condicionados a las variables explicativas&#46; En la pr&#225;ctica&#44; esta limitaci&#243;n puede conducir a estimadores sesgados e ineficientes&#44; sobrestim&#225;ndose normalmente la precisi&#243;n de los mismos ya que&#44; por lo com&#250;n&#44; las observaciones de un mismo individuo suelen estar positivamente correlacionadas&#46; Con el fin de hallar una soluci&#243;n a esto&#44; aparecen los modelos marginales y los modelos condicionales&#46; La principal diferencia entre estos modelos es el trato que recibe la dependencia&#46; Mientras que los primeros estiman el modelo ignorando la dependencia entre observaciones&#44; corrigi&#233;ndola posteriormente mediante estimadores jackknife&#44; bootstrap o &#171;s&#225;ndwich&#187;&#44; los segundos la estiman especificando expl&#237;citamente la distribuci&#243;n de probabilidad de la misma e incorpor&#225;ndola en el modelo&#46; Nuestro objetivo en este art&#237;culo consiste en ilustrar estas dos aproximaciones utilizando una base de datos de infecciones nosocomiales&#46;"
      ]
      "en" => array:1 [
        "resumen" => "The fact that individuals can present more than one event &#40;fail&#41; in the observation period &#40;multivariate survival analysis&#41; demands a different methodology from that used in univariate survival analysis&#46; The main problem of the Cox model with multivariate data is that the observations are not independent implying&#44; among others&#44; the violation of the proportionality hypothesis&#46; The Andersen-Gill approximation to the Cox model &#40;AG&#41; overcomes in part this problem&#46; The underlying hypotheses in the AG model&#44; however&#44; are very restrictive&#44; in particular that of independent increments&#46; Under this hypothesis the multiple observations of an individual are independent&#44; although conditioned on the explanatory variables&#46; In practice&#44; this limitation could let us to biased and inefficient estimators that are usually overestimated because&#44; in general&#44; the observations of an individual use to be positively correlated&#46; As a solution to this problem&#44; marginal and conditional models turn up&#46; The main difference between these models is the treatment that dependence receives&#46; Whereas the first one estimates the model ignoring the dependence between recurrences and adjustes this at the end with jackknife&#44; boostrap or &#171;sandwich&#187; estimates&#59; the second one estimates the dependence specifying explicitly the probability distribution of the dependence and incorporating this into the model&#46; Our objective in this paper is to illustrate these two approaches using nosocomial infection data&#46;"
      ]
    ]
    "multimedia" => array:8 [
      0 => array:8 [
        "identificador" => "tbl1"
        "etiqueta" => "Tabla 1"
        "tipo" => "MULTIMEDIATABLA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "copyright" => "Elsevier Espa&#241;a"
        "tabla" => array:1 [
          "tablatextoimagen" => array:1 [
            0 => array:1 [
              "tablaImagen" => array:1 [
                0 => array:4 [
                  "imagenFichero" => "138v16nSupl.2-13042666tab01.gif"
                  "imagenAlto" => 398
                  "imagenAncho" => 719
                  "imagenTamanyo" => 11875
                ]
              ]
            ]
          ]
        ]
        "descripcion" => array:1 [
          "es" => "Estructura de los datos seg&#250;n las diferentes aproximaciones"
        ]
      ]
      1 => array:8 [
        "identificador" => "tbl2"
        "etiqueta" => "Figura 1"
        "tipo" => "MULTIMEDIATABLA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "copyright" => "Elsevier Espa&#241;a"
        "tabla" => array:1 [
          "tablatextoimagen" => array:1 [
            0 => array:1 [
              "tablaImagen" => array:1 [
                0 => array:4 [
                  "imagenFichero" => "138v16nSupl.2-13042666tab02.gif"
                  "imagenAlto" => 335
                  "imagenAncho" => 561
                  "imagenTamanyo" => 6798
                ]
              ]
            ]
          ]
        ]
        "descripcion" => array:1 [
          "es" => "Estructura de los datos a comparar seg&#250;n el modelo de an&#225;lisis&#46;"
        ]
      ]
      2 => array:8 [
        "identificador" => "tbl3"
        "etiqueta" => "Tabla 2"
        "tipo" => "MULTIMEDIATABLA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "copyright" => "Elsevier Espa&#241;a"
        "tabla" => array:1 [
          "tablatextoimagen" => array:1 [
            0 => array:1 [
              "tablaImagen" => array:1 [
                0 => array:4 [
                  "imagenFichero" => "138v16nSupl.2-13042666tab03.gif"
                  "imagenAlto" => 624
                  "imagenAncho" => 602
                  "imagenTamanyo" => 23663
                ]
              ]
            ]
          ]
        ]
        "descripcion" => array:1 [
          "es" => "Resultados de la estimaci&#243;n del modelo Andersen-Gill &#40;AG&#41; est&#225;ndar y del modelo AG con estimaci&#243;n robusta de la matriz de variancias y covariancias"
        ]
      ]
      3 => array:8 [
        "identificador" => "tbl4"
        "etiqueta" => "Tabla 3"
        "tipo" => "MULTIMEDIATABLA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "copyright" => "Elsevier Espa&#241;a"
        "tabla" => array:1 [
          "tablatextoimagen" => array:1 [
            0 => array:1 [
              "tablaImagen" => array:1 [
                0 => array:4 [
                  "imagenFichero" => "138v16nSupl.2-13042666tab04.gif"
                  "imagenAlto" => 810
                  "imagenAncho" => 599
                  "imagenTamanyo" => 32841
                ]
              ]
            ]
          ]
        ]
        "descripcion" => array:1 [
          "es" => "Resultados de la estimaci&#243;n del modelo de Prentice&#44; Williams y Peterson &#40;PWP&#41; del modelo de Cox penalizado del modelo EMB&#46; PWP y EMB con estimaci&#243;n robusta de la varianza"
        ]
      ]
      4 => array:5 [
        "identificador" => "tbl5"
        "tipo" => "MULTIMEDIATABLA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "copyright" => "Elsevier Espa&#241;a"
      ]
      5 => array:5 [
        "identificador" => "tbl6"
        "tipo" => "MULTIMEDIATABLA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "copyright" => "Elsevier Espa&#241;a"
      ]
      6 => array:5 [
        "identificador" => "tbl7"
        "tipo" => "MULTIMEDIATABLA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "copyright" => "Elsevier Espa&#241;a"
      ]
      7 => array:5 [
        "identificador" => "tbl8"
        "tipo" => "MULTIMEDIATABLA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "copyright" => "Elsevier Espa&#241;a"
      ]
    ]
    "bibliografia" => array:2 [
      "titulo" => "Bibliograf&#237;a"
      "seccion" => array:1 [
        0 => array:1 [
          "bibliografiaReferencia" => array:24 [
            0 => array:3 [
              "identificador" => "bib1"
              "etiqueta" => "1"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Regression models and life tables&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "Cox DR&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "JRSS, Series B,"
                        "fecha" => "1972"
                        "volumen" => "34"
                        "paginaInicial" => "187"
                        "paginaFinal" => "220"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            1 => array:3 [
              "identificador" => "bib2"
              "etiqueta" => "2"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Testing the proportional hazards assumption in the presence of unmeasured heterogeneity&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "McCall BP&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "J Applied Econometrics"
                        "fecha" => "1994"
                        "volumen" => "9"
                        "paginaInicial" => "321"
                        "paginaFinal" => "34"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            2 => array:3 [
              "identificador" => "bib3"
              "etiqueta" => "3"
              "referencia" => array:1 [
                0 => array:2 [
                  "referenciaCompleta" => "Factors determining the incidence of nosocomial infections in an Intensive Care Unit of a public hospital in Girona, Spain. Mimeo, 2001."
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Factors determining the incidence of nosocomial infections in an Intensive Care Unit of a public hospital in Girona&#44; Spain&#46; Mimeo&#44; 2001&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "Barcel&#243; MA"
                            1 => "Saez M&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                ]
              ]
            ]
            3 => array:3 [
              "identificador" => "bib4"
              "etiqueta" => "4"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Non-parametric inference for a family of counting processes&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "Aalen OO&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Ann Stat"
                        "fecha" => "1978"
                        "volumen" => "6"
                        "paginaInicial" => "701"
                        "paginaFinal" => "26"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            4 => array:3 [
              "identificador" => "bib5"
              "etiqueta" => "5"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Cox&#39;s regression model for counting processes&#58; a large sample study&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "Andersen PK"
                            1 => "Gill RD&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Ann Stat"
                        "fecha" => "1982"
                        "volumen" => "10"
                        "paginaInicial" => "1100"
                        "paginaFinal" => "20"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            5 => array:3 [
              "identificador" => "bib6"
              "etiqueta" => "6"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Counting process models for life history data&#58; a review&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "Andersen PK"
                            1 => "Borgan O&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Scan J Stat"
                        "fecha" => "1985"
                        "volumen" => "12"
                        "paginaInicial" => "97"
                        "paginaFinal" => "158"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            6 => array:3 [
              "identificador" => "bib7"
              "etiqueta" => "7"
              "referencia" => array:1 [
                0 => array:2 [
                  "referenciaCompleta" => "Counting processes and survival analysis. New York: John Wiley &#38; Sons, Inc., 1991."
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Counting processes and survival analysis&#46; New York&#58; John Wiley &#38; Sons&#44; Inc&#46;&#44; 1991&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "Fleming TR"
                            1 => "Harrington DP&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                ]
              ]
            ]
            7 => array:3 [
              "identificador" => "bib8"
              "etiqueta" => "8"
              "referencia" => array:1 [
                0 => array:2 [
                  "referenciaCompleta" => "Statistical models based on counting proceses. New York: Springer-Verlag, 1993."
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Statistical models based on counting proceses&#46; New York&#58; Springer-Verlag&#44; 1993&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:4 [
                            0 => "Andersen PK"
                            1 => "Borgan O"
                            2 => "Gill RD"
                            3 => "Keiding N&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                ]
              ]
            ]
            8 => array:3 [
              "identificador" => "bib9"
              "etiqueta" => "9"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "rhDNase as an example of recurrent event analysis&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "Therneau TM"
                            1 => "Hamilton SA&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:6 [
                        "tituloSerie" => "Stat Med"
                        "fecha" => "1997"
                        "volumen" => "16"
                        "paginaInicial" => "2029"
                        "paginaFinal" => "47"
                        "link" => array:1 [
                          0 => array:2 [
                            "url" => "https://www.ncbi.nlm.nih.gov/pubmed/9308130"
                            "web" => "Medline"
                          ]
                        ]
                      ]
                    ]
                  ]
                ]
              ]
            ]
            9 => array:3 [
              "identificador" => "bib10"
              "etiqueta" => "10"
              "referencia" => array:1 [
                0 => array:2 [
                  "referenciaCompleta" => "Modelos de fragilidad en el análisis multivariante de supervivencia [tesis doctoral]. Girona: Universitat de Girona, 2001."
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Modelos de fragilidad en el an&#225;lisis multivariante de supervivencia &#91;tesis doctoral&#93;&#46; Girona&#58; Universitat de Girona&#44; 2001&#46;"
                      "idioma" => "es"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "Barcel&#243; MA&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                ]
              ]
            ]
            10 => array:3 [
              "identificador" => "bib11"
              "etiqueta" => "11"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Regression analysis of multivariate incomplete failure time data by modeling marginal distributions&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:3 [
                            0 => "Wei LJ"
                            1 => "Lin DY"
                            2 => "Weissfeld L&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "JASA"
                        "fecha" => "1989"
                        "volumen" => "84"
                        "paginaInicial" => "1065"
                        "paginaFinal" => "73"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            11 => array:3 [
              "identificador" => "bib12"
              "etiqueta" => "12"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "On the regression analysis of multivariate failure time data&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:3 [
                            0 => "Prentice RL"
                            1 => "Williams BJ"
                            2 => "Peterson AL&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Biometrika"
                        "fecha" => "1981"
                        "volumen" => "68"
                        "paginaInicial" => "373"
                        "paginaFinal" => "9"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            12 => array:3 [
              "identificador" => "bib13"
              "etiqueta" => "13"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "An additive frailty model for correlated life times&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "Petersen JH&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:6 [
                        "tituloSerie" => "Biometrics"
                        "fecha" => "1998"
                        "volumen" => "54"
                        "paginaInicial" => "646"
                        "paginaFinal" => "61"
                        "link" => array:1 [
                          0 => array:2 [
                            "url" => "https://www.ncbi.nlm.nih.gov/pubmed/9660631"
                            "web" => "Medline"
                          ]
                        ]
                      ]
                    ]
                  ]
                ]
              ]
            ]
            13 => array:3 [
              "identificador" => "bib14"
              "etiqueta" => "14"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "A counting process approach to maximum likelihood estimation in frailty models&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:4 [
                            0 => "Nielsen G"
                            1 => "Gill R"
                            2 => "Andersen P"
                            3 => "Sorensen T&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Scan J Stat"
                        "fecha" => "1992"
                        "volumen" => "19"
                        "paginaInicial" => "25"
                        "paginaFinal" => "43"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            14 => array:3 [
              "identificador" => "bib15"
              "etiqueta" => "15"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "GNFCN&#58; a computer program for gamma frailty analysis&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "Nielsen G"
                            1 => "Tabanera-Palacios R&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:4 [
                        "titulo" => "GNFCN&#58; a computer program for gamma frailty analysis&#46;"
                        "fecha" => "1993"
                        "paginaInicial" => "1"
                        "paginaFinal" => "21"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            15 => array:3 [
              "identificador" => "bib16"
              "etiqueta" => "16"
              "referencia" => array:1 [
                0 => array:2 [
                  "referenciaCompleta" => "A modification of the EM algorithm to estimate an Andersen-Gill gamma frailty model for multivariate failure time data. Mimeo, 2001."
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "A modification of the EM algorithm to estimate an Andersen-Gill gamma frailty model for multivariate failure time data&#46; Mimeo&#44; 2001&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "Barcel&#243; MA"
                            1 => "S&#225;ez M&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                ]
              ]
            ]
            16 => array:3 [
              "identificador" => "bib17"
              "etiqueta" => "17"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Generalized linear models for the analysis of Taguchi-type experiments&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "Nelder JA"
                            1 => "Lee Y&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Applied Stochastic Models and Data Analysis"
                        "fecha" => "1991"
                        "volumen" => "7"
                        "paginaInicial" => "107"
                        "paginaFinal" => "20"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            17 => array:3 [
              "identificador" => "bib18"
              "etiqueta" => "18"
              "referencia" => array:1 [
                0 => array:2 [
                  "referenciaCompleta" => "Modeling survival data. Extending the Cox model. Statistics for biology and health. New York: Springer, 2000."
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Modeling survival data&#46; Extending the Cox model&#46; Statistics for biology and health&#46; New York&#58; Springer&#44; 2000&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "Therneau TM"
                            1 => "Grambsch PM&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                ]
              ]
            ]
            18 => array:3 [
              "identificador" => "bib19"
              "etiqueta" => "19"
              "referencia" => array:1 [
                0 => array:2 [
                  "referenciaCompleta" => "Seattle: Data Analysis Products Division, MathSoft, 1999."
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Seattle&#58; Data Analysis Products Division&#44; MathSoft&#44; 1999&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "S-PLUS 2000 User&#39;s Guide&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                ]
              ]
            ]
            19 => array:3 [
              "identificador" => "bib20"
              "etiqueta" => "20"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Discussion of paper by Wei and Glidden&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "Cook RJ"
                            1 => "Lawless JF&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Stat Med"
                        "fecha" => "1997"
                        "volumen" => "16"
                        "paginaInicial" => "841"
                        "paginaFinal" => "3"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            20 => array:3 [
              "identificador" => "bib21"
              "etiqueta" => "21"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Discussion of paper by Wei and Glidden&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "Lipschutz KH"
                            1 => "Snapinn SM&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Stat Med"
                        "fecha" => "1997"
                        "volumen" => "16"
                        "paginaInicial" => "846"
                        "paginaFinal" => "8"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            21 => array:3 [
              "identificador" => "bib22"
              "etiqueta" => "22"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Some approaches to the analysis of recurrent event data&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "Clayton D&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Stat Med"
                        "fecha" => "1994"
                        "volumen" => "3"
                        "paginaInicial" => "244"
                        "paginaFinal" => "63"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            22 => array:3 [
              "identificador" => "bib23"
              "etiqueta" => "23"
              "referencia" => array:1 [
                0 => array:2 [
                  "referenciaCompleta" => "Frailty models for multiple event times. En: Klein YP, Goel PK, editors. Survival analysis: state of the art. Boston: Kluwer Academic Publishers, 1992; p. 415."
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Frailty models for multiple event times&#46; En&#58; Klein YP&#44; Goel PK&#44; editors&#46; Survival analysis&#58; state of the art&#46; Boston&#58; Kluwer Academic Publishers&#44; 1992&#59; p&#46; 415&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "Oakes D&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                ]
              ]
            ]
            23 => array:3 [
              "identificador" => "bib24"
              "etiqueta" => "24"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:3 [
                      "titulo" => "Discussion of paper by Wei and Glidden&#46;"
                      "idioma" => "en"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "Oakes D&#46;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:4 [
                        "tituloSerie" => "Stat Med"
                        "fecha" => "1997"
                        "volumen" => "16"
                        "paginaInicial" => "843"
                      ]
                    ]
                  ]
                ]
              ]
            ]
          ]
        ]
      ]
    ]
  ]
  "idiomaDefecto" => "es"
  "url" => "/02139111/00000016000000S2/v0_201302051443/13042666/v0_201302051445/es/main.assets"
  "Apartado" => array:4 [
    "identificador" => "779"
    "tipo" => "SECCION"
    "es" => array:2 [
      "titulo" => "Revisiones"
      "idiomaDefecto" => true
    ]
    "idiomaDefecto" => "es"
  ]
  "PDF" => "https://static.elsevier.es/multimedia/02139111/00000016000000S2/v0_201302051443/13042666/v0_201302051445/es/138v16nSupl.2a13042666pdf001.pdf?idApp=WGSE&text.app=https://gacetasanitaria.org/"
  "EPUB" => "https://multimedia.elsevier.es/PublicationsMultimediaV1/item/epub/13042666?idApp=WGSE"
]

Nueva respuesta rápida:

Requisitos y consideraciones

Las respuestas rápidas tienen como objetivo permitir el debate sobre observaciones científicas y de opinión sobre trabajos publicados el último año.

Deben ser enviadas desde la página web correspondiente al artículo en cuestión (no a través del sistema general de envío de artículos).

Son moderadas antes de su publicación a discreción del equipo editorial de la revista, no obstante no son peer-reviewed.

Deben ser originales y no infringir los derechos propiedad de terceros.

Los/las autores/as perderán el derecho a borrar o editar el texto, una vez sea publicado. Asimismo, la revista se reserva el derecho de borrar o editar el texto a discreción.

Algunas respuestas rápidas podrían ser seleccionadas adicionalmente para su publicación como cartas en la revista, en este caso seguirán el proceso editorial correspondiente.

Las respuestas rápidas serán publicadas en la web de la revista en 1 semana aproximadamente, después del envío del texto por parte de los/as autores/as.

La extensión máxima será de 500 palabras. Se admitirá una tabla o figura y hasta 5 referencias bibliográficas.

Se comunicará a los/las autores/ras del artículo comentado la existencia de la carta y la posibilidad de respuesta.

Al enviar esta respuesta rápida los autores/as aceptan los requisitos descritos.

Indexada en:

Síguenos:

Respuestas rápidas

Nueva respuesta rápida:

Requisitos y consideraciones

Respuesta

Conflictos de interés

Indexada en:

Síguenos:

Respuestas rápidas

Nueva respuesta rápida:

Requisitos y consideraciones

Respuesta

Conflictos de interés

Suscríbase a la newsletter