Para la proyección de brotes epidemiológicos se
involucra la aplicación de métodos estocásticos para encontrar
patrones, describir frecuencias de enfermedades que son de
interés del sector salud y para el monitoreo de alertas epidémicas
en la población. En el presente trabajo se propone el uso de un
modelo computacional basado en un Autómata Celular estocástico
global para la simulación del modelo epidemiológico SEIR
(Susceptibles, Expuestos, Infectados, Recuperados) con el que se
pueda observar la morfología de eventos epidemiológicos y poder
simular estrategias de vacunación (Inmunización). Para validar la
morfología se contrasta el modelo computacional con un modelo
matemático clásico.
Presentacion: Modelado para estudio de brotes epidémicos usando un Autómata Celular Estocástico Global.
1. Modelado para estudio de brotes epid´micos usando un
e
Aut´mata celular Estoc´stico Global.
o a
´
Autores: Hector Cuesta-Arvizu, Angel Bravo-Salgado
Armin R. Mikler y Adri´n Trueba-Espinosa
a
Centro Universitario UAEM Texcoco
Center for Computational Epidemiology and Response Analysis - University of North Texas
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 1 / 19
3. Contenido
Introducci´n.
o
Modelo Epidemiol´gico SEIR.
o
Modelado Matem´tico para representar el SEIR.
a
Aut´mata celular.
o
Modelo Epidemiol´gico Estocastico Global.
o
Simulador de Brotes Epidemiol´gicos.
o
Estrategias de Vacunaci´n.
o
Conclusiones.
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 3 / 19
4. Introducci´n
o
Motivaci´n
o
Los Epidemi´logos y los sistemas de salud p´blica utilizan modelos
o u
para estudiar la propagaci´n de enfermedades infecciosas durante un
o
brote epid´mico.
e
Dichos Modelos incluyen: Modelos Matem´ticos, Estad´
a ısticos y
Computacionales.
Simulando dichos modelos es una forma en la que se puede observar
diferentes evoluciones en diferentes escenarios que en otro caso no se
podr´ ya que puede ser muy costoso, no ´tico o simplemente no
ıa e
existen los medios para reproducir el escenario.
Otra ventaja de la simulaci´n es el tener la capacidad de crear
o
estrategias de intervenci´n que permitan manipular el contexto de la
o
epidemia.
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 4 / 19
5. Modelado de Enfermedades Infecciosas
Modelado de Enfermedades Infecciosas
SEIR Susceptibles-Expuestos-Infectados-Recuperados
El Modelo Epidemiol´gico SEIR mantiene el ciclo de vida de una
o
enfermedad infecciosa a trav´s de cuatro estados: Susceptible (S),
e
Expuesto (E), Infectado (I), Recuperado o removido(R). Cada uno de
esos estados representa el n´mero de individuos en dicho grupo.
u
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 5 / 19
6. Modelado de Enfermedades Infecciosas
Modelado de Enfermedades Infecciosas
Modelado Matem´tico para SEIR
a
El Modelo SEIR se puede representar con el siguiente sistema de
Ecuaciones Diferenciales:
dS
= −β ∗ S ∗ I
dt
dE
=β∗S ∗I −σ∗E
dt
dI
=σ∗E −γ∗I
dt
dR
=γ∗I
dt
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 6 / 19
7. Modelado de Enfermedades Infecciosas
Modelado de Enfermedades Infecciosas
Modelado Matem´tico para SEIR
a
Se realiz´ la simulaci´n en el sistema de ecuaciones diferenciales
o o
en Lenguaje R:
En la grafica se observa la curva de un brote epid´mico.
e
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 7 / 19
8. Aut´mata celular
o
Aut´mata celular
o
¿Que es un Aut´mata celular?
o
Modelo Discreto estudiado en teor´ de la computaci´n y
ıa o
matem´ticas. Para problemas no lineales.
a
Facts:
Consiste en un n´mero infinito de c´lulas en un grid regular donde cada
u e
c´lula tiene un numero de estados finitos.
e
El grid puede constar de cualquier n´mero finito de dimensiones.
u
Cada c´lula representa a un individuo de la poblaci´n.
e o
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 8 / 19
9. Aut´mata celular
o
Aut´mata celular
o
Vecindarios
El vecindario es una seleccion de c´lulas relativas a cierta c´lula
e e
especifica cuya posici´n en el Grid no cambia.
o
Cada c´lula tiene el mismo set de reglas para actualizar su estado
e
basado en los valores de su vecindario.
Cada vez que las reglas son aplicadas a todos las c´lulas del grid una
e
nueva generaci´n es producida.
o
Vecindarios Locales y Globales, Vecindarios de Von Neumann y Moore.
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 9 / 19
10. Modelo
Idealizaci´n para el Estudio
o
Se asume para este modelo de contagio (SEIR) lo siguiente:
Se supone una poblaci´n cerrada (que no cambia a trav´s del tiempo).
o e
Una mezcla homog´nea de contactos entre individuos de la poblaci´n.
e o
Cada individuo tiene en promedio el mismo numero de contactos.
No se consideran las variables demograficas o distancias geograficas
como factores para este modelo.
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 10 / 19
11. Modelo
El Modelo de Contacto Estoc´stico Global
a
El objetivo de este modelo es describir la din´mica de una enfermedad
a
infecciosa en una poblaci´n cerrada.
o
Es un modelo Global de Interacci´n Humano-Humano
o
Su proposito es el simular la din´mica de contacto entre individuos de
a
la poblaci´n. Facilitando el an´lisis de la propagaci´n de cierta
o a o
enfermedad.
El Aut´mata celular es representado en un grafo de cayley que muestra la
o
interacci´n entre c´lulas.
o e
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 11 / 19
12. Modelo
Interacci´n de Contacto Global
o
Contactos por generaci´n (Time Step):
o
CR∗N
C= 2
Total de contactos en el evento:
Ctot = Σtπ CR∗N donde te = (1, 2, 3, ..., n)
t=1 2
C = N´mero de interacciones por generaci´n.
u o
CR = Promedio de Contacto.
N = N´mero de individuos en la poblaci´n.
u o
tπ = N´mero de generaciones.
u
π = Restricci´n de Termino, cuando E + I = 0
o
Ctot = N´mero total de interacciones en el evento.
u
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 12 / 19
13. Software de Simulacion
Simulador de Brotes Epidemicos
Opciones tecnol´gicas para el Simulador
o
La mayor contribuci´n del presente trabajo es el desarrollo de un
o
software para simular brotes epid´micos incorporando un modelo de
e
aut´mata celular estoc´stico global.
o a
Opciones Tecnologicas:
C# .NET (como lenguaje de programaci´n)
o
WindowsForms y MonoDesktop (para crear interfaces graficas
multiplataforma)
Background Worker (Algoritmo de paralelizaci´n, sincronizando un
o
pool de hilos)
M´dulos:
o
M´dulo de especificaci´n.
o o
M´dulo de simulaci´n.
o o
M´dulo de visualizaci´n.
o o
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 13 / 19
14. Software de Simulacion
Simulador de Brotes Epid´micos
e
M´dulos de Especiaci´n y Simulaci´n
o o o
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 14 / 19
15. Software de Simulacion
Simulador de Brotes Epid´micos
e
M´dulo de Visualizaci´n
o o
En la Figura A podemos observar la curva de una epidemia tipo SEIR.
(A)
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 15 / 19
16. Estrategias de Intervenci´n
o
Estrategias de Vacunaci´n
o
Modelo de Vacunaci´n
o
Figura B.- Modelo de Vacunaci´n para SEIR
o
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 16 / 19
17. Estrategias de Intervenci´n
o
Estrategias de Vacunaci´n
o
Tipos de Estrategias de Vacunaci´n
o
Vacunaci´n programada.
o
Vacunaci´n disparada por picos en la poblaci´n infectada.
o o
Figura C.- Grafica de la estrategia de vacunaci´n programada.
o
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 17 / 19
18. Conclusiones
Conclusiones y Trabajo Futuro
Conclusiones y Trabajo Futuro
La Simulaci´n ayuda a entender la propagaci´n de una enfermedad
o o
infecciosa.
Se puede observar diferentes salidas en escenarios donde se aplican
estrategias de intervenci´n.
o
Trabajo Futuro:
Usar diferentes tipos de modelos de contacto.
Extender el modelo a SEIRS (Donde la poblaci´n Recuperada pierde
o
despu´s de cierto periodo su resistencia y regresa a ser Susceptible)
e
Integrar Estacionalidad.
Integrar Aspectos Demograficos y Geograficos.
Hector Cuesta-Arvizu (UAEM-UNT) Modelado para estudio de brotes epid´micos
e November 25, 2011 18 / 19