El dengue es una enfermedad viral transmitida al hombre por picaduras del mosquito Aedes aegypti.

El virus del dengue comprende cuatro serotipos o cepas (DEN 1-4), los cuales presentan pequeñas variaciones entre ellos. Por esta razón un individuo infectado adquiere inmunidad permanente frente a los serotipos expuestos, pero solo una protección temporal frente a los serotipos a los que no se ha expuesto (Diekmann & Heesterbeek, 2000).

El dengue tiene diversas manifestaciones que pueden conducir a estadíos graves de la enfermedad, generalmente asociados a reinfecciones que se consideran, por tanto, un factor de riesgo.

El mosquito Aedes aegypti es de hábitos típicamente domiciliarios, se cría en climas tropicales húmedos y pica con mayor frecuencia entre las 06:00 a 08:00 y 17:00 a 19:00 horas del día.

Las etapas del vector son huevo, larva, pupa y adulto (Bar, 2011).

Se estima que actualmente 3000 millones de personas están en lugares donde hay riesgo de contagiarse del virus. Además, aproximadamente ocurren 100 millones de casos por año (Bar, 2011).

Las políticas de control de la epidemia se enfocan esencialmente en la erradicación del mosquito adulto y de los hospederos de larvas. Su eficacia depende de la situación epidemiológica, reforzada en temporadas de brote y responsabilidad de las dependencias de salud de cada región. El control del mosquito se puede hacer por medios físicos, biológicos, químicos o mediante la modificación de conductas inadecuadas de la población en riesgo.

Para obtener éxito en las campañas de erradicación del mosquito por medio de insecticidas, es indispensable su uso interrumpido, de lo contrario, el mosquito puede adquirir resistencia.

Actualmente hay varias vacunas contra el dengue en fase de desarrollo, que tienen el desafío de proporcionar protección contra los 4 serotipos, al generar la producción de anticuerpos facilitadores. Es necesario que los países mejoren la vigilancia epidemiológica para contribuir al análisis económico que fundamente el costo-efectividad del uso de las vacunas en el programa de inmunización. Además, en países endémicos y no endémicos se deben considerar estrategias de aplicación adecuadas antes de la licencia (Bar, 2011).

Marco teórico

Las cadenas de Markov se deben al matemático ruso A.A.Markov, alrededor de 1905. Su intención era crear un modelo probabilístico para analizar el uso de las vocales en poemas y textos literarios. El éxito del modelo propuesto por Markov radica en que es lo suficientemente complejo como para describir ciertas características no triviales de algunos sistemas, pero al mismo tiempo sencillo para ser analizado matemáticamente (Tudor, 1994).

Definición: Una cadena de Markov es un proceso estocástico a tiempo discreto con espacio de estado discreto, y que satisface la propiedad de Markov, esto es, para cualquier , y para cualquier estado , se cumple

.

Esto quiere decir que el proceso al tiempo futuro depende únicamente del estado del proceso al tiempo , y no depende de los estados en los tiempos pasados .

En una cadena definida con posibles estados , se puede introducir la notación:

para definir la probabilidad condicional de que estando en el estado en el momento anterior inmediato pase al estado .

Si entonces el estado está comunicado con el . Si, además, entonces la comunicación es mutua.

Los se denominan probabilidades de transición en el momento y todas ellas forman la matriz de transición de dimensión :

.

Cuando las probabilidades de transición no cambian en el tiempo, la probabilidad de pasar de un estado a un estado será la misma en todo instante de observación del proceso. En este caso se tendrá una única probabilidad de transición constante en el tiempo: . Este proceso será, entonces, una cadena de Markov homogénea o estacionaria en tiempo discreto. La evolución del comportamiento al paso del tiempo de las probabilidades involucradas en la matriz de transición empleando las ecuaciones de Chapman-Kolmogorov se puede calcular multiplicando la matriz de transición por ella misma, o sea, hallando la matriz potencia, en el tiempo establecido (Taylor, 1994).

Las aplicaciones en la investigación biomédica de las cadenas de Markov han sido múltiples, tanto en la experimentación animal (Yang & Chao, 2005), como en los estudios humanos (Trajstman, 2002).

Metodología

Proponemos inicialmente, un modelo con cadenas de Markov, para un problema de transmisión de dengue, considerando un serotipo y cinco estados por los que puede transitar un individuo durante la enfermedad y dos escenarios, que favorecen o no la proliferación del vector. Para este primer modelo, se brindan tablas de resultados y gráficos que ilustran la variación en el tiempo de los principales parámetros definidos.

En un segundo momento, se seleccionan dos parámetros fundamentales en la dinámica de transmisión, control y envergadura de la epidemia, como son la probabilidad de contagio por contacto con vector infectado y la probabilidad de muerte por la enfermedad, para simular un proceso de control de estos, con vistas a estudiar su influencia en la reducción o crecimiento de la epidemia. Para este caso se procede de manera análoga a la utilizada en el modelo sin control.

Estudio sin estrategia de control

Definimos los posibles estados por los que transita una persona en la dinámica de transmisión del dengue, : susceptible, : infectado, : inmune, : muerte por la enfermedad y : muerte por causas naturales.

Se asume para el estudio que infectarse, es tener contacto con un mosquito infectado, y además se puede realizar para cualquier cepa variando las probabilidades en cuestión.

Se declara una cadena de Markov discreta homogénea con los estados antes mencionados, donde se toma un mes como unidad de tiempo. Para la cual la respectiva matriz de transición :

.

Donde:

: Probabilidad de mantenerse en el estado susceptible.

: Probabilidad de pasar de susceptible a infectado.

:

Probabilidad de muerte por causa ajena a la enfermedad.

: Probabilidad de morir de la enfermedad dado que se esté infectado.

: Probabilidad de alcanzar la inmunidad.

: Probabilidad de mantenerse en el estado infectado.

: Probabilidad de mantener la inmunidad.

El estudio se realizó para una cepa, y en dos posibles escenarios:

Escenario I: Probabilidad de contagio relativamente baja, por no existir factores favorables para el desarrollo del vector.

Escenario II: Alta probabilidad de contagio, por factores ambientales e higiénicos favorables para el desarrollo del vector transmisor.

En ambos casos los valores de las probabilidades involucradas en la dinámica se toman de forma hipotética, pero se adaptan a una situación real y fueron examinados por especialistas, ver Tabla 1.

Se trabajó con la matriz potencia para analizar la evolución de estas probabilidades al paso de tiempo. Los resultados obtenidos para el Escenario I, a partir de la matriz de transición inicial fueron:

Para v = (2119722, 500, 12505, 22063, 30), vector de prevalencia inicial, que está compuesto por la cantidad de susceptibles, infectados, recuperados, muertos por causas ajenas a la enfermedad y muertos por la enfermedad. Este expresa que, partiendo de una población inicial, las cantidades iniciales de personas, en cada uno de los compartimentos, las cuales fueron avaladas por especialistas. Valiéndonos de la matriz de transición y realizando la operación , se obtienen los siguientes resultados, ver Tabla 3:

Según los resultados alcanzados, al paso del tiempo, la probabilidad de mantenerse susceptible y la de morir por la enfermedad disminuyen; nótese que, en el primer caso, se mantiene en una vecindad de 0.99, mientras la probabilidad de contagio aumenta, también manteniéndose en un entorno de 0.00012, por lo que consideramos que, aunque el comportamiento es el deseado, cuantitativamente no lo hace de forma significativa, debido a que el escenario es no propenso al desarrollo de la epidemia.

Respecto a la evolución de las cantidades partiendo del vector de prevalencia inicial, se obtuvo que el número de susceptibles disminuye, y aumenta el número de recuperados, la cantidad de infectados también disminuye: pero en una proporción despreciable, la cantidad de muertes por la enfermedad se mantiene prácticamente constante con el paso del tiempo.

Con una metodología análoga a la utilizada en el escenario anterior, para un vector inicial de prevalencia v = (2119722, 500, 12505, 22063, 30), se obtuvo que, al paso del tiempo, la probabilidad de mantenerse susceptible disminuye, mientras aumentan las probabilidades de contagio y de morir por la enfermedad. Así mismo, el número de susceptibles disminuye y aumentan el número de recuperados y la cantidad de infectados, ambos de forma significativa, la cantidad de muertes por la enfermedad aumenta moderadamente al paso del tiempo, ver Tabla 5. Estos resultados se explican, pues en este escenario existen condiciones idóneas para el desarrollo de los vectores, por lo que se impone una estrategia de control.

Estudio con estrategia de control

Se introducen dos controles: se aplica a la probabilidad de contagiarse por contacto con el vector infectado, considera la aplicación de un tratamiento efectivo para evitar las muertes por la enfermedad logrando la inmunidad, en el Escenario II (propenso al desarrollo de una epidemia), ambos representan el porciento de eficacia lograda en la aplicación del control.

La matriz de transición para este estudio es:

.

Se consideraron los siguientes valores para la efectividad de los controles:

Estrategia I: = 0.6 (efectividad del 60 %) y = 0.4 (efectividad del 40 %).

Estrategia II: = 0.5 y = 0.7.

Estas posibilidades de estrategia fueron evaluadas y aprobadas por especialistas.

Las Tablas 6, 7 y 8 presentan una comparación de las diferentes estrategias aplicadas a las probabilidades de mantenerse susceptible, de infectarse y de morir por la enfermedad. La estrategia I resulta ser la más efectiva, pues disminuye la probabilidad de infectarse y aumenta la probabilidad de mantenerse susceptible a mayor escala, aunque la probabilidad de morirse por la enfermedad alcanza mejores resultados en la Estrategia II, en ambos casos cumple con el objetivo porque se reduce respecto al modelo sin control.

Partiendo del vector de prevalencia inicial v = (2119722, 500, 12505, 22063, 30), se estudió el comportamiento en esta población:

En las tablas 9, 10 y 11 se muestra la evolución de las cifras de susceptibles, infectados y muertos por la enfermedad en los meses 2 y 3 del estudio. Se puede verificar que los resultados de la estrategia I superan a los de la estrategia II, complementando los resultados obtenidos en las probabilidades directamente. Por tanto, el control es más efectivo al lograr mayor efectividad en el control de los contactos entre los mosquitos infectados y el ser humano y menor en la aplicación de tratamiento.

Conclusión

Se presentó un estudio basado en la matriz de transición en una cadena de Markov homogénea para la dinámica de contagio del dengue. Se observó la evolución de la probabilidad de mantenerse susceptible, de infectarse y de morir por la enfermedad en dos diferentes escenarios, uno propenso al desarrollo de la epidemia y el otro no, lo que permite que el estudio se pueda adecuar a diversas regiones. Se propusieron estrategias de control basadas en la ponderación de las probabilidades trascendentes, incorporadas en la matriz de transición.

En trabajos anteriores habíamos analizado el comportamiento de tasas y probabilidades que representan parámetros trascendentes en la dinámica de transmisión de esta enfermedad en modelos descritos por sistemas de ecuaciones diferenciales. Aunque se reportan en la bibliografía disímiles aplicaciones de las cadenas de Markov en problemas de las biociencias, es la primera vez que realizamos un estudio de la evolución de dichas probabilidades y tamaños de poblaciones, con dicha herramienta para nuestro modelo, lo que nos permitirá un análisis comparativo con los resultados obtenidos utilizando otras estrategias de modelación.

A partir de los resultados obtenidos se puede predecir, con una herramienta matemática relativamente sencilla, la evolución de un brote epidémico en un tiempo finito. Esto brindaría potencialidades de controlar la epidemia a las diferentes instancias del sistema de salud, lo que repercutiría en la salud y mejor calidad de vida de las personas. Esto representa una fortaleza de esta propuesta, si fuese aplicada. Precisamente, la limitación esencial de esta primera propuesta es que utiliza datos hipotéticos, por lo que se sugiere realizar el estudio de forma estratificada en la población, pues se asume que todos los individuos poseen la misma probabilidad de contagio.

Referencias

Fechas de Publicación

Histórico