Introducción
Las enfermedades cardiovasculares (ECV) son un problema de salud de primer orden a nivel mundial. La Organización Mundial de la Salud ha estimado que es mayor el número de personas que muere anualmente por ECV que por cualquier otra causa y en 2015 se le atribuyó 31% de las muertes globales.1 La carga de la enfermedad y el costo relacionado con las ECV es mayor en países de ingresos bajos y medios2 y se asocia a un sinnúmero de desenlaces clínicos desfavorables para los pacientes.3
La irregularidad es una característica común a diversos objetos de la naturaleza y a algunas estructuras del cuerpo humano. Tal característica no permite su cuantificación a partir de las medidas euclidianas, lo cual permitió el surgimiento de la geometría fractal.4 Esta geometría ha permitido la comprensión de algunas dinámicas fisiológicas.5
Los fractales son objetos irregulares que pueden ser abstractos, salvajes o estadísticos6 y cuya cuantificación puede ser realizada con la dimensión fractal. Esta es una medida numérica adimensional que permite evaluar la irregularidad de un objeto.7
Los fractales estadísticos son caracterizados mediante la distinción de las frecuencias de aparición de determinadas variables; la dimensión fractal estadística establece el grado de complejidad de los fenómenos y se calcula con la Ley de Zipf-Mandelbrot, que se basa en una proporción matemática capaz de caracterizar eventos repetitivos, como la aparición de palabras en un escrito.8,9,10,11
En 1932, George Kingsley Zipf, evidenció ciertas regularidades estadísticas de los escritos y encontró que en el corpus de determinado texto existe una estructura hiperbólica al graficar la frecuencia de ocurrencia de cada palabra y un rango desde 1 hasta N, asignado a la lista de las palabras utilizadas en el texto organizadas por frecuencias decrecientes.9 El comportamiento hiperbólico encontrado, se transformó posteriormente a partir de logaritmos, y dio lugar a una linealización10,11 que permite la obtención de la dimensión fractal estadística.
En cardiología, se ha evaluado la variabilidad de la frecuencia cardíaca (VFC) en patologías como el síndrome coronario agudo o la falla cardíaca, y su relación con desenlaces como la rehospitalización y la mortalidad.12,13,14 Sin embargo, algunas investigaciones han evidenciado que la dinámica cardíaca presenta un comportamiento caótico,15,16,17 que puede analizarse desde perspectivas no lineales.18,19,20 La variabilidad de la frecuencia cardíaca se ha analizado en adultos e infantes desde diversos métodos dentro de los que figura el análisis espectral,21,22,23 siendo requeridos análisis más certeros de los sistemas cardíacos, desarrollados desde la no linealidad.
Por otro lado, fue desarrollada una metodología diagnóstica de la dinámica cardíaca a partir de la Ley de Zipf-Mandelbrot, que permite establecer la complejidad de cada dinámica y cuantificar la dimensión fractal estadística. La aplicabilidad clínica del método fue corroborada,24 diferenciando adecuadamente dinámicas cardíacas agudas de normales a partir del grado de complejidad de los sistemas cardíacos.
El objetivo de la presente investigación es corroborar la capacidad diagnóstica y utilidad clínica de un software que automatiza la metodología fundamentada en la Ley de Zipf-Mandelbrot en el marco de un estudio ciego con 80 registros Holter.24
Material y Métodos
Se realizó un estudio retrospectivo realizado en la Clínica del Country y el Grupo Insight en Bogotá, Colombia, en 2017.
Definiciones
Dimensión fractal estadística: Se realizó la linealización logarítmica para establecer la dimensión fractal estadística de Zipf-Mandelbrot, empleando el método de mínimos cuadrados y hallando el inverso multiplicativo negativo de la pendiente, correspondiente a la dimensión fractal estadística (Ecuación 1).
D, es la dimensión fractal; σ, corresponde al rango para cada frecuencia; P, es la frecuencia de aparición de la frecuencia cardiaca; V, se define por V=1/N-1 siendo N el número de frecuencias medidas y F un cofactor secundario en la linealización.
Población
Fueron tomados 80 registros Holter de mínimo 21 horas de sujetos con edad superior a los 21 años, pertenecientes a la base de datos del grupo Insight. Se dividieron en dos grupos: uno, de pacientes normales (20 registros Holter) y otro, de sujetos cursando con patologías cardíacas (60 registros Holter) a quienes se indicó una prueba de electrocardiografía ambulatoria con base en los protocolos médicos clínicos convencionales para el estudio de enfermedades cardiovasculares agudas incluyendo arritmias y cardiopatías isquémicas y estructurales.
Procedimiento
Los diagnósticos clínicos fueron enmascarados. Posteriormente, los valores de la frecuencia cardíaca mínima y máxima de los Holters durante cada hora fueron introducidos al software desarrollado, que automatiza la metodología diagnóstica de la dinámica cardíaca fundamentada en la Ley de Zipf-Mandelbrot.24 El Software determinó el número de frecuencias cardíacas que se encontraban en cada rango de 15 lat/min, efectuó el cálculo de la frecuencia de aparición de los rangos de 15 lat/min y aplicó la Ley de Zipf-Mandelbrot a estos valores.
Asimismo, se ordenó las frecuencias de aparición de los rangos de mayor a menor, relacionando cada frecuencia con un número; empezó por el uno en el caso de la mayor frecuencia y aumentó una unidad en las siguientes frecuencias de aparición organizadas de forma jerárquica. De este modo, el mayor número se asignó a la menor frecuencia.
Posteriormente, se graficaron los datos para confirmar el comportamiento hiperbólico necesario para la aplicación de la ley de Zipf-Mandelbrot, a continuación, se aplicó el método de mínimos cuadrados en estos datos; se halló la dimensión fractal estadística de las dinámicas cardíacas del estudio (Ecuación 1).
Con los hallazgos obtenidos, se realizó el análisis estadístico, para lo cual fueron desenmascarados los diagnósticos clínicos. Para efectuar la comparación entre el diagnóstico matemático obtenido por medio del software y el diagnóstico clínico convencional, se tomó este último como Gold Standard. Fue realizada una clasificación binaria, en la cual los verdaderos positivos representaban los casos categorizados matemáticamente como anormales y con diagnóstico clínico convencional de anormalidad, los falsos positivos correspondían a pacientes diagnosticados como normales por el experto clínico y con diagnóstico matemático de anormalidad, los falsos negativos eran los casos con diagnóstico matemático de normalidad y convencional de enfermedad, los verdaderos negativos eran el número de pacientes categorizados como normales mediante el diagnóstico clínico convencional y la metodología matemática incorporada en el software. El coeficiente Kappa fue calculado con la siguiente fórmula: (Ecuación 2).
Co, corresponde a las concordancias observadas, es decir, al número de casos diagnosticados igualmente por el método matemático incorporado en el software y el Gold Standard. To, representa el número total de casos y Ca, las concordancias atribuibles al azar.
Esta investigación se categoriza como de riesgo mínimo, pues se realizaron cálculos matemáticos a partir de exámenes mínimamente invasivos prescritos previamente, con base en protocolos médicos convencionales. La investigación se acoge a los lineamientos contemplados en la Resolución 8430 del Ministerio de Salud Colombiano y a los principios éticos de la Declaración de Helsinki de la Asociación Médica Mundial.
Resultados
Los valores de la dimensión fractal estadística variaron entre 0,720 y 0,913 en el caso de registros Holter normales (Tabla 1), y entre 0,454 y 0,665 en el caso de registros Holter de pacientes con enfermedad aguda (Tabla 2). A partir de los resultados de la dimensión fractal estadística, se puede observar que a medida que la dinámica cardíaca de un sujeto se va acercando a un estado de agudización, esta va perdiendo el grado de complejidad, en cambio se observa un aumento en los grados de complejidad para los sujetos con dinámica cardiaca normal.
Los valores del coeficiente de correlación R2 hallados entre la dispersión de puntos generados al graficar de manera jerárquica la frecuencia de aparición de los rangos de 15 lat/min y la línea de mayor ajuste de estos puntos oscilo entre 0,616 y 0,935 para los registros Holter normales (Tabla 1), y entre 0,667 y 0,949 en los casos de registros Holter con enfermedad aguda (Tabla 2). Los anteriores resultados muestran la posibilidad de hacer linealizaciones en comportamientos complejos, como es el caso del comportamiento de la dinámica cardiaca, en el contexto de la geometría fractal estadística. La sensibilidad y especificidad fue de 100%, mientras que el coeficiente Kappa presentó un valor de 1.
Discusión
Este es el primer trabajo en el cual es confirmada la utilidad clínica y diagnóstica de un software que automatiza una metodología diagnóstica de la dinámica cardiaca basada en la Ley de Zipf-Mandelbrot25 en el contexto de un estudio ciego con 80 registros Holter, en el cual fueron diferenciados adecuadamente los sistemas cardíacos normales de los patológicos en estado agudo mediante la organización jerárquica de las frecuencias cardíacas en los rangos y el comportamiento hiperbólico evidenciado; se obtuvo un diagnóstico matemático unívoco con base en la autoorganización matemática del fenómeno y el grado de complejidad del sistema, evidenciado a partir de la dimensión fractal estadística, al margen de consideraciones causales de cualquier orden.
El software como una herramienta de ayuda diagnóstica puede ser incorporado al ámbito clínico. Estudios han mostrado que los pacientes con alto riesgo de muerte cardíaca súbita presentan bifurcaciones o cambios espectrales abruptos y oscilaciones de frecuencia cardíaca entre los 0,1 y los 0,4 Hz en el contexto de patrones no lineales.19 La entropía aproximada26 y el análisis multifractal27 han sido empleados en la evaluación de los sistemas cardiacos. También se ha relacionado la complejidad de estos sistemas con la VFC y con dinámicas no lineales.20 Por su parte, la metodología que fue automatizada con el software presenta una gran simplicidad, pues la complejidad del sistema es lo único requerido en la distinción de dinámicas cardíacas normales de aquellas con patologías agudas.
Goldberger y Cols., establecieron que la anormalidad de un sistema se relaciona con una elevada periodicidad o aleatoriedad y que la normalidad es un estado intermedio entre los dos extremos mencionados, ofrecieron una perspectiva distinta a la planteada por la fisiología clásica;5 no obstante, no fueron propuestos parámetros matemáticos capaces de diferenciar normalidad de enfermedad, aspecto logrado en el presente estudio a partir de la dimensión fractal estadística.
La metodología diagnóstica de la dinámica cardíaca basada en la Ley de Zipf-Mandelbrot24 sigue una línea de pensamiento acausal, que ha permitido realizar diagnósticos y predicciones a partir de la incorporación de conceptos y teorías fisicomatemáticas a la Medicina. Se han desarrollado predicciones en inmunología,28 biología molecular29 y la dinámica temporal de la malaria en 820 municipios colombianos.30