Abstract: Implementación de un algoritmo para la detección de brotes aplicado a un sistema microbiológio

Introduccion

Un algoritmo estadístico para el cálculo del valor umbral debe tener en cuenta las fluctuaciones en los datos debidas a ciclos estacionales, tendencia y brotes pasados, que varían de un microorganismo a otro. Además, el sistem a debe ser suficientemente robusto para ser aplicado a una gran variedad de microorganismos. Farrington et al. (1996) desarrollaron un algoritmo que tiene en cuenta estos factores, calculando el valor umbral mediante un modelo de regresión log-lin eal ajustado por sobredispersión, estacionalidad, tendencia y brotes pasados.

El objetivo de este trabajo es la adaptación a este sistema de vigilancia del modelo propuesto por Farrington et al. (1996). El paquete estadístico Stata 5.0 fue utilizado para el desarrollo del algoritmo.

Asumimos que los datos siguen una distribución de Poisson de media m _i excepto cuando existe sobredispersión; en ese caso la media de la distribución de Poisson se supone distribuída de acuerdo a una Gamma con media m _i y varianza m _i(f -1), f parámetro de sobredispersión, lo que resulta en una Binomial-Negativa con media m _i y varianza f m _i. La componente sistemática del modelo es: ln m _i =a +b t_i +g n_i, donde t_i es la variable tiempo medida en semanas y n_i es el número de hospitales en los que se produjeron los casos declarados. Asumimos también que existe correlación entre los datos de un mismo año, y que los datos de años distintos son independientes.

Descripción del algoritmo

Para cada uno de los microorganismos notificados se siguen los siguientes pasos:

Los datos históricos para estimar el modelo de regresión para la semana t₀ son los correspondientes a las semanas t₀-3 a t₀+3 de los cinco años anteriores, resultando un total de 35 datos:

Se estima un modelo inicial utilizando poisson y se contrasta la sobredispersión de los datos.

Si hay sobredispersión, se estima un modelo de regresión Binomial-Negativa haciendo uso del ado ml con el métodolf para maximizar la función de verosimilitud y tomando como valores iniciales los obtenidos con el modelo de Poisson. Este nuevo modelo proporciona estimaciones iniciales y permite definir una función de ponderación que asigna pesos bajos a las observaciones atípicas. A continuación se reajusta el modelo teniendo en cuenta los pesos def inidos, y de este modo se corrige el efecto de posibles brotes pasados. En este modelo la correlación existente entre datos pertenecientes a un mismo año se corrigió en el cálculo del valor umbral.

En otro caso, se reestima el modelo de Poisson con poisml utilizando las opciones robust, que corrige el efecto de brotes pasados, y cluster, que tiene en cuenta la correlación existe nte entre los datos del mismo año.

En ambos casos, tanto en el modelo de Poisson como en el de la Binomial-Negativa, se contrasta la significación de la tendencia, de tal forma que si esta no es significativa se elimina del modelo y se repite la estimación sin esta variable.

En cada caso el ajuste del modelo nos proporciona una predicción por intervalo con una confianza del 97.5% para la semana t₀. El valor umbral será el límite superior del intervalo de predicción.

El resultado final del algoritmo es una tabla con la siguienta estructura:

        ------------------------+-----------------------------------
                      Organismo |      Casos      Umbral       Brote
        ------------------------+-----------------------------------
                       A.caviae |          5        2.92       Brote
                  Aeromonas sp. |          0        0.02           -
                         C.coli |          4        8.78           -
                    C.difficile |          1        3.20           -
                      C.laridis |          0        2.18           -
              Campylobacter sp. |         16       24.16           -
               Y.enterocolitica |          8        8.53           -
        Y.enterocolitica ser.03 |          1        3.45           -
        ------------------------+-----------------------------------

Referencias