Introducción
La dimensión fractal es una medida de la irregularidad que se puede hallar en los objetos de la naturaleza 1-11. Basada en la geometría fractal expuesta por Mandelbrot 12-14, esta dimensión es desarrollada para superar los resultados paradójicos que se obtienen al aplicar geometría euclidiana a objetos irregulares. Una de las metodologías para la determinación de la dimensión fractal es la ley de Zipf-Mandelbrot, mediante la cual se puede evaluar la complejidad de un sistema determinado a partir del análisis de la distribución de las frecuencias específicas de los eventos que lo componen cuando estas presentan un comportamiento hiperbólico 15-21. Un uso dado a esta ley se puede evidenciar en la caracterización de la complejidad del repertorio T específico en la respuesta inmune 22-24; también se utiliza para crear un diagnóstico de la monitoria fetal, diferenciando entre normalidad y enfermedad en el proceso de gestación 25.
El dengue es un virus que se transmite a través de la picadura del mosquito Aedes aegypti y, en décadas recientes, se ha convertido en una epidemia mundial, sobre todo en países tropicales y subtropicales 26. Según estimaciones de la Organización Mundial de la Salud, se producen 390 millones de infecciones al año de las cuales solo 96 millones presentan manifestaciones clínicas; 3 900 millones de personas están en riesgo de contraer dengue 27 y para 2015, solo en América, se reportan 2.35 millones de casos y 1 181 muertes por esta causa 27.
En el caso de Colombia, para el periodo 2007-2015 se registraron 669 631 casos de dengue y 814 muertes 28-30, siendo en 2014 el tercer país de América con mayor número de registros y el primero en el área andina 31; para mediados de septiembre de 2016 ya se habían reportado 36 774 casos confirmados 32. De acuerdo al Instituto Nacional de Salud, de los 37 423 casos reportados hasta la semana 36 de este mismo año, 607 correspondieron a Palmira en el Valle del Cauca, Colombia 33; teniendo en cuenta los altos valores de casos, fue evidente la necesidad de implementar sistemas de vigilancia más precisos a la hora de llevar a cabo acciones para hacer frente a esta problemática.
Rodríguez et al.34 describieron la dinámica de transmisión del dengue, tanto del denominado clásico como del hemorrágico, en el período 2001-2004 en Palmira, encontrando una disminución del 88.2% en los casos registrados de dengue clásico. Este estudio permitió observar disminuciones porcentuales en diferentes grupos generacionales de la población durante los diferentes años, así mismo estableció las distribuciones de casos de dengue clásico por grupos de edad para el mismo período de tiempo, datos que fueron analizados matemáticamente en este trabajo.
La mayoría de descripciones del comportamiento de este tipo de dinámicas se centran en metodologías de tipo epidemiológico y estadístico, aunque también se han desarrollado diversos modelos matemáticos para su comprensión. Sin embargo, desde ambas perspectivas suele analizarse el fenómeno a partir de un examen multicausal 35-40, buscando establecer relaciones cuantitativas entre el aumento o la disminución de la epidemia respecto a factores asociados como el clima 39. Esta forma de ver la dinámica de la epidemia presenta ciertas dificultades ya que el alcance de sus conclusiones se ve limitado por la cantidad de factores causales analizados; además, su aplicabilidad se restringe a los lugares y períodos temporales específicamente estudiados 41, lo que dificulta la comprensión de este fenómeno desde una perspectiva más general que permita el desarrollo de predicciones. Por ejemplo, los modelos de corredores endémicos 41 sirven para establecer los valores esperados de incidencia de una enfermedad en un tiempo y lugar determinado, pero requieren para ello información previa de 5 a 7 años.
Con el fin de resolver este problema, diferentes estudios han desarrollado varias metodologías matemáticas y físicas de predicción de epidemias de malaria y dengue desde una perspectiva acausal, es decir, sin tener en cuenta las variables que pueden influir en su variación y centrándose en las características matemáticas de esta. Entre ellas se encuentran metodologías diseñadas con base en la teoría de la probabilidad 42, la caminata al azar 43, las ecuaciones diferenciales 44,45 y las relaciones S/k de la entropía 46. Estas metodologías son útiles en la toma de decisiones en salud pública ya que evitan los inconvenientes de tener en cuenta el azar como factor predominante dentro de la dinámica y son útiles para sistemas de vigilancia y alerta temprana tanto a nivel nacional como regional.
En este contexto, el objetivo de este estudio es desarrollar una nueva metodología de caracterización desde la ley de Zipf-Mandelbrot que permita comparar la variación anual del grado de complejidad del comportamiento del dengue entre 2001 y 2004, así como entre grupos de edades para la epidemia de dengue en Palmira.
Material y métodos
Definición
La dimensión fractal estadística es la medida adimensional del grado de complejidad de un sistema que se obtiene con el inverso de la pendiente de la linealización logarítmica de la ley de Zipf, cuyo valor está entre 0 y 1 y su ecuación es:
Donde σ es el rango asignado a cada frecuencia P, P es la frecuencia de aparición, V está definida por V=1/n-1, n es el número de frecuencias medidas y F es geométricamente el punto de corte con el eje vertical de la recta obtenida en la linealización.
Procedimiento
Rodríguez et al.34 recopilaron los datos epidemiológicos del número de posibles infectados de dengue en Palmira reportados para los grupos de edad <1 año, 1 a 4, 5 a 14, 15 a 44, 45 a 59 y >60 entre los años 2001 y 2004. En principio se evaluó si el comportamiento de cada año respecto a los grupos de edad y de cada uno de los grupos de edades para los cuatro años estudiados correspondía a un comportamiento hiperbólico. Para ello las frecuencias de cada dinámica evaluada se organizaron de mayor a menor y se les asignó un rango en orden inversamente proporcional a la frecuencia, es decir, a mayor frecuencia se asigna el menor grado, que es 1, a la siguiente frecuencia se le asigna el 2, y así sucesivamente. De este modo, se construyeron gráficas donde las frecuencias se tomaron como variable dependiente y los rangos como variable independiente.
Luego, se realizó una linealización logarítmica de la hipérbole obtenida y se halló su factor de correlación R2, que indicó qué tan acertada es la aproximación lograda con la linealizacion, de tal forma que 1 correspondió a una perfecta correlación y entre más se acercó a 0 menos fue la correlación; es decir, la aproximación estuvo muy alejada del comportamiento de los datos. A continuación, se estableció la dimensión fractal estadística mediante la ecuación detallada en definiciones, que fundamentó la complejidad del sistema con el inverso de la pendiente de la linealización logarítmica. Estas medidas de complejidad fueron comparadas con el fin de establecer diferencias respecto a años y grupos de edad específicos.
Por último, se realizaron tres simulaciones en las variaciones de los seis valores que componen la dinámica de cada año con base en los valores límites hallados para las hipérboles construidas para los años 2001 a 2004 y para evidenciar diferentes formas de autoorganización fractal que puede tomar la dinámica. Además, se realizaron dos simulaciones para los rangos de edad que presentaron valores de dimensión fractal más altos y más bajos siguiendo el mismo procedimiento.
Resultados
Se encuentra que tanto las distribuciones de infectados de dengue por grupos de edad como la variación total desde el año 2001 al 2004 obedecen a un comportamiento hiperbólico. En la Figura 1 se presenta la distribución hiperbólica de infectados por grupos de edad para todos los años estudiados, en la Figura 2 la correspondiente linealización logarítmica para el año 2002 y en la Tabla 1 los cálculos necesarios para obtener la dimensión fractal estadística de dicha distribución.
p: frecuencia relativa de los valores de grupos de infectados por edad respecto a la totalidad de casos del año 2002; σ: rango asociado por la aplicación de la ley Zipf-Mandelbrot; Log: logaritmo; v: 1/n-1=0.2. Fuente: Elaboración propia
Entre los años 2001 y 2004 las dimensiones fractales encontradas estuvieron entre 0.5329 y 0.8703, con valores de R2 entre 0.7862 y 0.9518 (Tabla 2). Al observar las dimensiones fractales obtenidas para cada grupo de edad se observa que presentaron valores entre 0.4694 y 0.6689, con R2 entre 0.5583 y 0.9872 (Tabla 3).
Al realizar las tres simulaciones de las dinámicas de años, se encontraron valores de dimensión fractal entre 0.4512 y 0.6316, con R2 entre 0.9656 y 0.9881; las hipérboles halladas para cada una de estas simulaciones de las dinámicas de años se plasmaron en la Figura 3. Para la dinámica por grupos el valor de dimensión fractal más bajo fue el correspondiente al rango 15-44 años y el valor más alto fue el del rango 45-59 años; para estas dos dinámicas se desarrollaron simulaciones: la de la primera presentó un valor de dimensión fractal 0.4545 con R2 de 0.8212, mientras que la de la segunda fue de 0.6095 con R2 de 0.9242.
Discusión
Este es el primer trabajo en el que se establece una nueva metodología de caracterización de la posible epidemia de dengue en Palmira a partir de la ley de Zipf-Mandelbrot, estableciendo a su vez la complejidad del sistema para grupos de edad y para los años comprendidos entre 2001 y 2004.
En esta investigación se comprobó el comportamiento fractal estadístico de la dinámica y se calcularon dimensiones fractales que permiten comparar la variación de la complejidad de la dinámica tanto para los años como los grupos de edades evaluados sin tener en cuenta factores causales. Del mismo modo, se estableció un comportamiento finito y acotado del fenómeno, donde la determinación de los grados de complejidad resulta de utilidad en la toma de decisiones en salud pública, por ejemplo en la evaluación y seguimiento de los grupos de edad más afectados. La limitación del método dependió de que se pudieran encontrar distribuciones de frecuencias de número de infectados.
Al tratarse de una ley matemática, se tiene la gran ventaja de que la dimensión fractal estadística es independiente del tiempo y es este un primer ejemplo de su aplicación; también, al encontrarse grados de complejidad acotados se implica que la metodología es independiente de subregistros o sobreregistros, los cuales fueron sistematizados en la época del estudio original, pues la autoorganización del sistema prevé las dificultades de recolección de la información al estar trabajando con leyes matemáticas que están evaluando un fenómeno físico particular, como ya ha ocurrido en otros estudios de la física 47,48. Del mismo modo, si se implementan medidas de salud pública contra la enfermedad, el cambio en el grado de complejidad podría sugerir mejores o peores intervenciones, por lo tanto la pérdida de la complejidad permitiría evaluar la eficacia de las intervenciones, mientras que el aumento evidenciaría la falla de las mismas como se evidencia con las simulaciones realizadas (Figura 3).
Otros trabajos que tienen como objetivo el estudio del dengue o epidemias con un comportamiento similar como la malaria se centran en el análisis de factores y variables que inciden en su transmisión desde una perspectiva causal 35-41. La propagación del dengue se ha asociado a cambios en el clima -que crean condiciones ideales para la expansión y diseminación de los vectores-, al crecimiento poblacional en condiciones de infraestructura sanitaria inadecuada y a la migración poblacional 35-37,49. En este contexto se han desarrollado diferentes modelos matemáticos (biológicos-epidemiológicos) para la caracterización de estas dinámicas que, a pesar de haber contribuido a la comprensión de la dinámica de transmisión, se centran en el análisis y determinación de patrones espaciotemporales de la dinámica de extensión de la epidemia 35-37, en la estimación del número de afectados a partir de datos epidemiológicos espaciales 38 y en modelos descriptivos basados en el clima 39 desde una perspectiva causal 34-39; esto hace que sus afirmaciones queden restringidas al número de variables causales que analizan los espacios específicos o los intervalos de tiempo 39, impidiendo la generalización de los modelos o la realización de afirmaciones de carácter predictivo del fenómeno.
Con el fin de solucionar estos problemas, Rodríguez et al. 43 han desarrollado diferentes metodologías para la predicción de la dinámica de dengue y malaria independiente del análisis de variables causales, ejemplo de ello es un método para la predicción matemática de la dinámica anual de número de infectados de malaria con base en la teoría de la probabilidad, que luego fue aplicada exitosamente a la predicción del dengue.
Del mismo modo, se realizaron predicciones de las posibles trayectorias de estas epidemias a partir de ecuaciones diferenciales 44,45 y se hizo una predicción de la dinámica del dengue con base en la caminata al azar 43. Desde esta perspectiva, también se desarrolló una metodología de predicción de brotes de malaria en 820 municipios de Colombia basada en la teoría de la probabilidad y la relación S/k de la entropía, logrando un acierto de 99.86% 46. Esta fue una metodología físicomatemática acausal con la que se estableció un comportamiento matemáticamente autoorganizado de la dinámica del dengue independiente del análisis de posibles causas asociadas a partir de la ley de Zipf-Mandelbrot, lo que evidenció un comportamiento fractal estadístico independiente del período de tiempo o el rango de edad especifico evaluado, del mismo modo que ocurre en fenómenos como la distribución de palabras en los lenguajes naturales 15 o respuestas inmunes específicas 22-24.
Con la ley de Zipf-Mandelbrot se ha evaluado el grado de complejidad en el lenguaje de personas con enfermedad de Alzheimer y normales 15-18, en la monitoría cardiaca fetal 25, en el repertorio T específico contra el alérgeno Poa p9 en un paciente alérgico con presencia y ausencia de interferón α y en los clones de células Th 24.
En el presente trabajo los resultados obtenidos son independientes del tiempo o de la muestra evaluada, permitiendo establecer diferencias en la complejidad hallada para los diferentes años y rangos de edad estudiados. Además de identificar órdenes matemáticas de tipo acausal del fenómeno, este tipo de estudio puede ser aplicado en futuras investigaciones para la evaluación de otros años y otros espacios geográficos, permitiendo establecer comparaciones y facilitar la toma de decisiones de salud pública.
Esta perspectiva acausal basada en teorías físicas y matemáticas ha dado lugar al desarrollo de diagnósticos y predicciones en otras áreas de medicina como la inmunología, donde se han desarrollado predicciones de unión de péptidos y en morfometría, aplicada en células preneoplásicas y neoplásicas de cuello uterino 50 y en arterias coronarias de un modelo experimental de reestenosis. En cardiología también se han desarrollado métodos de diagnóstico de la dinámica cardiaca de aplicabilidad clínica tanto en adultos como en neonatos, mientras que en infectología se han desarrollado predicciones del número de linfocitos T CD4 en pacientes con VIH/ SIDA a partir de los valores de linfocitos y leucocitos del cuadro hemático. Estos trabajos evidencian la utilidad práctica del estudio de los fenómenos médicos a partir de teorías físicas y matemáticas, revelando órdenes matemáticas subyacentes a estos fenómenos.
Conclusiones
Se desarrolló una metodología original para la caracterización matemática de la dinámica de dengue en Palmira para el periodo 2001-2004 con base en la ley de Zipf-Mandelbrot. Por medio de esta caracterización se evidenció que la dinámica obedece a un comportamiento fractal estadístico finito y acotado tanto para la evaluación de las variaciones en los años estudiados como para los grupos de edad.
Los grados de complejidad hallados permiten comparar las variaciones de la dinámica al margen de factores causales, siendo útiles para la toma de decisiones en salud pública.