WO2017144763A1 - Procedimiento para la obtención de datos útiles asociados al patrón de variabilidad de la frecuencia cardiaca - Google Patents

Abstract

La presente invención tiene por objeto un procedimiento para la descripción, representación gráfica e identificación gráfica de patrones específicos de funcionamiento de sistemas cíclicos cuasi-periódicos, como por ejemplo pero no restringido a motores de combustión alternativo, máquinas rotatorias u órganos biológicos como el corazón. También es objeto de la presente invención un procedimiento para calcular un índice de cualificación de la condición o salud cardíaca de un individuo así como su aplicación para diagnosticar y emitir juicios prognósticos sobre la funcionalidad, patología o nivel de salud de una máquina u organismo dotados con un motor u órgano de funcionamiento cíclico y para la descripción, representación gráfica compacta, e identificación gráfica de patrones específicos de funcionamiento de sistemas dinámicos, por ejemplo sistemas económicos como el mercado de valores bursátiles.

Description

DESCRIPCIÓN

PROCEDIMIENTO PARA LA OBTENCIÓN DE DATOS ÚTILES ASOCIADOS AL PATRÓN DE VARIABILIDAD DE LA FRECUENCIA CARDIACA

CAMPO DE LA INVENCIÓN

La invención que se propone consiste en un procedimiento de obtención de datos gráficos o de sus códigos multidimensionales resultantes, para caracterizar y cualificar sistemas cíclicos o cuasi-periódicos disipativos de cualquier naturaleza, ya sea natural o artificial, de una forma compacta, universal, adaptable y accesible. Dicha invención se basa, esencialmente, en los siguientes elementos:

(i) Una formulación generalizada de la variabilidad normalizada y secuencial de los intervalos temporales entre ciclos subsecuentes, utilizando una llave o clave de diversos parámetros, que determinan la expresión algorítmica particular de dicha variabilidad secuencial. Esta expresión también puede considerarse como una transformada generalizada de la serie numérica original de los intervalos temporales. Además, esta transformación es unívoca y establece una partición secuencial de la serie original, seleccionando grupos subsecuentes de N elementos en un orden preestablecido por la clave elegida. Dichos grupos se pueden considerar vectores en un espacio N-dimensional.

(ii) Una representación espacial en N-dimensiones de los vectores o puntos generados por la transformación definida por la clave elegida. En particular, para N=4, se pueden utilizar las 3 dimensiones y el color. En el caso N=5, se pueden utilizar las 3 dimensiones espaciales, el color, y una partición de la quinta dimensión que puede representarse dinámicamente, por ejemplo mediante una secuencia temporal (video) de la representación tridimensional correspondiente a cada intervalo de dicha partición, utilizando el tiempo como dimensión de representación de la quinta dimensión.

(iii) La identificación de patrones comunes o correlaciones espaciales ("clusters") entre grupos de puntos, tales como ejes o planos cristalográficos, superficies generalizadas, etc., obtenidas a través de las representaciones de un número suficiente de sistemas de la misma naturaleza.

(iv) La comparación de un sistema individual de una cierta naturaleza con dichos patrones comunes, y la determinación de la presencia de dichos patrones comunes en el sistema individual que se trata de cualificar.

Contrariamente a cualquier otro método que utilice otras transformaciones, como la transformada de Fourier o las "Wavelet Transforms" de cualquier tipo, el procedimiento propuesto es intrínsecamente adaptativo, y contempla la acumulación de conocimiento a través de la experiencia. Permite además utilizar la escala de tiempo instantánea, o cualquier otra que pueda elegirse localmente, como la escala de referencia principal del análisis; esto permite identificar patrones de variabilidad universales, independientes de otras escalas de tiempo o variables exógenas, que son inaccesibles para otras transformaciones como la de Fourier ya que éstas superponen todas las escalas temporales independientemente de su secuencia de aparición.

ANTECEDENTES DE LA INVENCIÓN

Los sistemas disipativos cíclicos se caracterizan por períodos de tiempo específicos, ya sean constantes o variables. Por lo tanto, pertenecen a una clase especial de sistemas dinámicos cuyos grados de libertad están críticamente restringidos o limitados. Al contrario de lo que ocurre con los sistemas dinámicos disipativos genéricos, que exhiben un comportamiento caótico, los atractores de sistemas disipativos cíclicos están bien

1

HOJA DE REEMPLAZO (REGLA 26) definidos y son generalmente simples. Tales sistemas son, en general, motores mecánicos o biomecánicos, que intercambian aportes o demandas de energía o trabajo con el medio ambiente, ya sea un motor de automóvil, una turbina de viento, o un corazón. Sin embargo, mientras que los periodos de tiempo constantes caracterizan ambos, relojes artificiales y motores naturales, diseñados específicamente para minimizar la variabilidad -o que no están sujetos a cambios de demandas-, los motores sujetos a demandas variables presentan también períodos de tiempo variables. Esta estrategia es la forma más económica para garantizar la capacidad de adaptación de un organismo (por ejemplo, un animal) o una máquina (por ejemplo, un coche) equipada con estos motores.

La capacidad de adaptación a los cambios en las demandas es probablemente la principal prioridad de los sistemas naturales y artificiales dotados de capacidades de movimiento. Estos sistemas son inherentemente disipativos. Por otra parte, todos los sistemas artificiales diseñados para trabajar bajo demanda constante, necesitan también acelerar desde una pausa o desacelerar hasta detenerse. El grado de adaptabilidad, o cómo un sistema responde a una demanda dada -que a su vez cuenta con su propio tiempo característico-, califica su fuerza, su robustez, o su salud, y por lo general determina su supervivencia.

En general, los sistemas móviles adaptables experimentan períodos en los que están en reposo: por ejemplo, los regímenes de ralentí de los motores de combustión interna para evitar frecuentes arranques y paradas, el sueño de los animales, etc. En estos períodos de inactividad, el sistema amplía el período de tiempo característico, exclusivamente para equilibrar la disipación interna. Sin embargo, la complejidad de los sistemas móviles más adaptables, exige que sus motores cíclicos internos tengan un número limitado de grados de libertad. Estos grados de libertad limitados, que garantizan la adaptabilidad, son casi incompatibles con periodos de tiempo constantes. La variabilidad exhibida por un sistema en modo inactivo refleja abiertamente sus características internas y describe sus compromisos, a menudo ocultos bajo un funcionamiento general. Por ejemplo, los motores de dos tiempos exhiben una fuerte variabilidad cuando están al ralentí, y la variabilidad de la frecuencia cardíaca se puede observar durante el sueño profundo o en situación de relax profundo, mucho mejor que durante el ejercicio exhaustivo en sujetos sanos.

Se ha propuesto una amplia serie de métodos para caracterizar los sistemas dinámicos disipativos: descripción de atractores extraños, exponentes de Lyapunov, análisis de entropía, leyes de potencia, análisis de Fourier, retratos de fase multidimensional, etc. Sin embargo, ninguno de estos métodos proporciona información exhaustiva o retratos unívocos para los sistemas disipativos cíclicos y adaptativos, debido a su naturaleza inherente. No hay un equivalente compacto a un código QR o un gráfico para proporcionar información completa sobre la naturaleza, la adaptabilidad, el estado de salud y las características internas de un sistema disipativo cíclico.

Un caso muy general de sistemas con elevada variabilidad y adaptabilidad es el de los organismos vivos dotados de sistema circulatorio y movilidad. El sistema circulatorio de los organismos vivos es un sistema mecánico autónomo delicadamente acoplado con el sistema respiratorio, y ambos desarrollados por la evolución como respuesta a los complejos patrones de demanda de oxígeno asociados con el movimiento. La salud circulatoria se basa en la capacidad de adaptación, lo que implica una variabilidad inherente. En el ejemplo de realización de la presente invención que se describe en lo que sigue, un gráfico N-dimensional calculado mediante el método descrito y que representa la variabilidad de la frecuencia cardíaca revela dos patrones arrítmicos universales como señas específicas de la salud: uno refleja la capacidad de adaptación cardiaca, y el otro la sintonía del ritmo cardíaco-respiratorio. Además, se identifican al menos tres patrones universales arrítmicos cuya presencia aumenta progresivamente en detrimento proporcional de los dos patrones de salud, en determinadas situaciones patológicas (infarto de miocardio, insuficiencia cardíaca, y recuperación tras muerte súbita). La presencia de las estructuras arrítmicas universales identificadas, junto

2

HOJA DE REEMPLAZO (REGLA 26) con la posición del centro de masa del gráfico de la variabilidad de la frecuencia cardíaca, proporcionan una valoración cuantitativa sin precedentes del gradiente patología-salud.

El corazón es la primera bomba volumétrica autónoma orgánica desarrollada por la naturaleza, y permitió el gran salto que supuso la movilidad para los organismos vivos. Las funciones básicas del corazón son la contractilidad y la frecuencia cardíaca. La variabilidad de la frecuencia cardíaca (heart rate variability, o HRV) representa un grado de libertad vital de la evolución de un organismo autónomo, con un sistema circulatorio (Malik, 1996) que permite la respuesta adaptativa inmediata a la demanda de oxígeno. Estas demandas pueden tener infinitos perfiles de variabilidad. Sin embargo, por principios de economía, la naturaleza responde a las demandas con la creación de un número limitado de estructuras o patrones, en lugar de dar respuestas diferentes a las infinitas soluciones posibles. Por lo tanto, cabe preguntarse: ¿en qué medida aparecen estructuras o patrones predeterminados en la HRV?, ¿podrían revelar modos generales de adaptabilidad (salud) o de fracaso (patología) algunas de esas estructuras?, ¿pueden aparecer combinadas?, ¿cuáles serían su tamaño y peso relativos generalizados respecto del resto de eventos?, ¿se trata de estructuras complementarias o antagónicas?, etc.

El sistema respiratorio, complementario al circulatorio en los organismos móviles que viven en el aire, también es una bomba volumétrica (formada por el tórax y los pulmones) con otro tiempo característico y un funcionamiento con un apreciable nivel de acoplamiento al cardíaco. La relación flexible entre la frecuencia cardíaca y la respiratoria implica una HRV específica. La frecuencia cardíaca está también sometida a otras influencias endógenas, como la digestión, la edad y el sexo, mecanismos bioquímicos, o actividad psíquica, cuyos tiempos característicos están desacoplados del control autónomo de la frecuencia cardíaca. Cuando tales influencias alcanzan o superan las demandas externas (incluyendo el ciclo circadiano), el organismo puede presentar arritmias patológicas. Sin embargo, las que más amenazan la vida son las que presentan frecuencias superiores a la frecuencia respiratoria. El ejemplo de aplicación de la presente invención que se describe aquí revela que en la especie humana pueden existir secuencias arrítmicas específicas universales (estructura interna) como producto co-evolutivo del sistema nervioso autónomo (simpático-parasimpático), y que esas secuencias pueden ser específicas del corazón sano o del corazón patológico.

La necesidad de herramientas no invasivas, precisas y concretas para el diagnóstico y el pronóstico es un factor prioritario para el avance de la Medicina. Las representaciones gráficas compactas de los sistemas fisiológicos han supuesto una gran ayuda para los médicos en cuanto a los detalles y la precisión en el diagnóstico. La obtención de imágenes de los sistemas internos y los tejidos (por ejemplo, la ecografía, la tomografía axial computarizada, o la resonancia magnética nuclear (RMN)) ha cambiado nuestras vidas. Por ello, en la cuantificación cardíaca (CQ), la ecocardiografía (echoCG) ha representado un importantísimo desarrollo. No obstante, las series de datos electrocardiográficos proporcionan información complementaria insustituible, no accesible para la echoCG en muchos aspectos, como la aparición de patrones y comportamientos patológicos temporales, ritmos ectópicos, etc. La reducción de las series temporales de un registro Holter (HR) de frecuencia cardíaca a una representación gráfica compacta (o gráfico) puede tener un valor clínico complementario al de la echoCG. Son necesarias unas consideraciones preliminares: con el fin de determinar la aparición de patrones genéricos que reflejen características normales o patológicas, y para mejorar su especificidad, el gráfico debe reducir al mínimo (i) las diferencias individuales en peso / tamaño, sexo, edad, etc., (ii) las influencias sobre la frecuencia cardíaca a largo plazo, con tiempos comparables o mayores a los de la frecuencia respiratoria, y (iii) las influencias tanto exógenas como endógenas, distintas de las de origen cardíaco o circulatorio. La representación gráfica buscada debe aportar la información completa proporcionada por la secuencia de todos los eventos cardíacos (por ejemplo, la secuencia completa latido a

3

HOJA DE REEMPLAZO (REGLA 26) latido, o los intervalos RR). Tras una revisión exhaustiva, no se han encontrado evidencias de la existencia de una representación con todas esas características buscadas, hasta la fecha.

Entre los diferentes enfoques para estudiar la HRV, incluyendo herramientas para el análisis de sistemas dinámicos y caos, análisis de entropía multiescala, análisis en el dominio de la frecuencia, y las cadenas de Markov, los mapas de Poincaré de los intervalos RR proporcionan gráficos compactos útiles para el análisis de la HRV. La alta heterogeneidad estocástica que presenta en general la HRV es indicio de la existencia de fuertes estructuras internas presumiblemente universales. Bajo esta premisa, la herramienta más adecuada para la HRV sería un análisis de las secuencias temporales normalizadas, en contra de las herramientas estocásticas globales, o las transformadas de Fourier. De hecho, la transformada de Fourier puede excluirse definitivamente como herramienta adecuada para analizar los patrones de variabilidad universales, ya que esencialmente analiza el contenido de escalas fijas de tiempo (o frecuencias) en el registro completo, de forma que mezcla los patrones de la misma naturaleza potencial para frecuencias cardíacas altas y bajas, en el mismo HR.

La representación de Poincaré (o mapas de retorno) proporciona un análisis directo de secuencia de la serie RR para identificar patrones temporales. Sin embargo, las limitaciones de la representación de los mapas de retorno de Poincaré en el plano espacial en dos dimensiones (2D) ocultan gran parte del potencial de este tipo de grafos.

Otro tipo de representación se describe, por ejemplo, en la patente US 4934374 (Ostlund et al.), donde el objeto representado es la proximidad entre vectores definidos por valores sucesivos de una señal con un retraso dado (time delay) como función del tiempo. El método descrito en este documento se basa en representar la variable tiempo en ambos ejes de un mapa bidimensional, pero cuya representación espacial sigue siendo plana (unidimensional). Asimismo, la representación descrita en dicha patente se basa en un concepto de "variabilidad enmascarada" por el módulo de un vector no universal. Por tanto, la influencia del estado físico instantáneo del paciente descartará toda forma fiable de detección de patrones universales.

Aunque esta representación permite distinguir sistemas cíclicos ordenados, cuasi-periódicos, o caóticos, sin embargo no permite distinguir secuencias complejas en el tiempo que se repiten o están específicamente asociadas a ciertas patologías. Estas características están asociadas, por ejemplo, a secuencias arrítmicas específicas que aparecen, en cualquier representación plana (como es el caso de US 493374), superpuestas e indistinguibles, impidiéndose así su identificación en las fases de evaluación o diagnóstico.

Otro ejemplo de representación plana es el registro de solicitud de patente publicado como US 2005/0171447 A1 (Esperer), en la que se reivindica el uso del mapa de retorno clásico de Poincaré en dos dimensiones. El autor de esa invención utiliza información adicional a la provista por la representación plana (en particular, la invención utiliza la densidad de local de puntos en cada zona del plano representado) para correlacionar los resultados con patologías conocidas. No obstante, esta representación oculta ineludiblemente, por superposición en el plano, muchas secuencias complejas multidimensionales que, de nuevo, la presente invención sí exhibe como característica específica.

Un ejemplo de representación multidimensional es mostrado en la solicitud de patente PCT/GB2015/050429, donde el objeto representado es un vector multidimensional consistente en los valores sucesivos de una señal cualquiera proveniente de un sistema cuasi-periódico o cíclico, separados una cantidad fija de tiempo τ entre cada dos valores sucesivos. El resultado es una órbita multidimensional más o menos compleja, cuya interpretación en términos de la identificación de determinadas secuencias complejas multidimensionales puede llegar a ser imposible, debido a la inevitable posibilidad de ocurrencia de zonas multievaluadas para una misma secuencia, lo cual impediría reconstruir dicha secuencia. Incluso variando el

4

HOJA DE REEMPLAZO (REGLA 26) valor de τ y produciendo una colección de representaciones diferentes, la identificación de secuencias complejas puede ser tan difícil que sería necesario un enorme esfuerzo computacional para compactar la información que se busca. Posteriormente, además, dicha información hay que correlacionarla con la información de patologías asociadas. En comparación, la elección de los valores de los vectores multidimensionales que la presente invención propone es la más adecuada para ofrecer la más compacta y directa representación de secuencias complejas que pueden visualizarse directamente, sin la necesidad de ningún esfuerzo computacional posterior; pueden aplicarse entonces algoritmos de identificación de forma ya conocidos (que no son objeto de la presente invención) para producir correlaciones útiles entre forma y patología de manera fiable, reproducible, compacta y rápida.

En suma, es necesario, en el presente campo técnico, desarrollar nuevos métodos de representación compacta del patrón de variabilidad de la frecuencia cardiaca, que permitan obtener una información de mayor detalle que la generada por los métodos utilizados hasta la fecha. La presente invención está orientada a resolver la citada necesidad, mediante un novedoso método de representación unívoca y compacta de sistemas cíclicos o cuasi-periódicos disipativos de cualquier naturaleza, y que resultan de especial utilidad en el análisis de patrones cardíacos. Asimismo, dicho método no depende del tiempo de obtención de datos y, mediante sus etapas de normalización y compactación, permite encontrar los patrones asociados y su evolución independientemente de la escala de tiempo y de forma multidimensional, lo que no resulta posible en otras representaciones, como las descritas en US 4934374, US 2005/0171447 o PCT/GB2015/050429. DESCRIPCIÓN BREVE DE LA INVENCIÓN

Es un primer objeto de la presente invención un procedimiento para la descripción, representación gráfica compacta, e identificación gráfica de patrones específicos de funcionamiento de sistemas cíclicos cuasi- periódicos, como por ejemplo pero no restringido a motores de combustión alternativo, máquinas rotativas, u órganos biológicos como el corazón, caracterizado por:

a) se mide y se registra un número M de intervalos de tiempo consecutivos {^_¡ }_¡=1 , correspondientes a ciclos de uno o más componentes de un complejo cardiaco "pQRSt" de un electrocardiograma, con una precisión igual o mejor que el 10% del valor medio del tiempo de ciclo, y siendo M mayor que 2;

b) se calcula la variabilidad sobre dichos M intervalos de una secuencia de vectores consecutivos de N componentes, de acuerdo con el algoritmo o transformación definida por la expresión:

,M-N

con la siguiente notación:

Los siguientes parámetros son números enteros y su elección determina la forma final de la transformación mencionada:

im N N N N s E E C c J J K K K )

donde:

- m es un índice natural que representa el orden de la variación discreta que se calcula;

5

HOJA DE REEMPLAZO (REGLA 26) - N es la dimensión o número de componentes de cada vector δ _} , siendo N > 2 ;

- N₀,N_l, yN₂ indican el número de valores que se usan para calcular el promedio local correspondiente indicado en la fórmula general del algoritmo;

- ε₀,ε_ι, ε₂ tienen valores binarios 0 ó 1 , e indican si los elementos correspondientes son, respectivamente, fijos o móviles en el cálculo de cada una de las componentes del vector δ _} ;

- ς₀, ς_ι tienen valores binarios 0 ó 1 , e indican si la media local es, respectivamente, fija o móvil;

- J₀ J₁ indican el retraso o adelanto del primer elemento que se toma en el cálculo a partir del índice y ;

- K₀, K_l, K₂ indican el retraso o adelanto del primer elemento que se toma en la serie local correspondiente para calcular el promedio local indicado.

El segundo objeto de esta invención es la representación gráfica en dos o más dimensiones de la posición del punto indicado por los valores de las componentes de cada uno de los vectores δ _} calculados según se describe anteriormente, por ejemplo, pero sin restricciones, utilizando las dos dimensiones espaciales de una pantalla, el valor numérico del color, y el tamaño de un marcador circular para representar cuatro dimensiones, o cinco si se utilizan proyecciones axonométricas o cónicas en el plano, permitiendo en este caso que el usuario pueda visualizar la extensión volumétrica de la representación tridimensional generada para la identificación de patrones gráficos específicos exhibidos por el sistema estudiado.

Un tercer objeto de la presente invención es la determinación automática de la existencia de patrones de comportamiento definidos por una función vectorial A = {α^ } ^ _n , donde los elementos a . son valores fijos concretos o funciones de una o más variables definidos por el usuario, sin restricción alguna, según el siguiente procedimiento:

Α - δ,

a. Cálculo del ángulo generalizado ^ cuyo coseno está dado por: cos(6? )

\\Α\ - \\δ\

donde los vectores δ_] son calculados según se describe anteriormente, y el símbolo significa la norma generalizada de un vector en N dimensiones, tal que:

b. Cálculo del número de eventos M' tal que el ángulo 6? es inferior a una tolerancia predeterminada s , siendo 0 < ε < 1 , preferiblemente 0 < ε < 0.1 , de tal forma que la función A debe ser explorada en su espacio de existencia para encontrar dichos eventos en que 6? . < ε ; es decir, que M' depende de la elección de los valores concretos de las variables de la función A . En particular, s A es constante,

M' es único.

c. Uso del coeficiente Μ M como un índice de cualificación del sistema según que exhiba comportamientos predeterminados definidos por el vector/4 en mayor o menor intensidad, donde dichos comportamientos pueden ser deseados o indeseados.

6

HOJA DE REEMPLAZO (REGLA 26) También se considera objeto de la presente invención un procedimiento para la determinación cuantitativa de la adaptabilidad de máquinas o sistemas cíclicos cuasi-periódicos, caracterizado por:

a. Someter a la máquina o sistema a una variación del régimen de funcionamiento, ya sea por aumento o reducción progresivos del consumo de energía, midiendo los valores consecutivos^ de los intervalos temporales correspondientes a cada ciclo.

b. El uso del coeficiente Μ Μ , calculado según el tercer objeto de esta invención mencionado, a partir de la serie {-^_¡}_¡=1 ^con la definición concreta del vector^ = t{(N + 1)/ 2- donde

N>1, y dicho vector correspondería una aceleración o deceleración homogénea del sistema si t = 1 ó t = -1 , respectivamente.

Otro objeto de esta invención es un procedimiento para calcular un índice de cualificación de la condición o salud cardíaca de un individuo, caracterizado por:

a. La medida y registro de un número de intervalos de tiempo consecutivos^}^ entre picos "R" del complejo cardíaco "pQSRt" con una precisión mejor que 10 ms, siendo mayor que dos; b. El cálculo de la serie \δ \ de vectores consecutivos de N dimensiones o componentes según se describió en el primer objeto de esta invención, calculados según las siguientes definiciones de los quince parámetros:

c. {m = 0, N = 5, N₀=5, N_t = 5, N₂=5, ε₀=\, ε_{1 =}1, ε₂=0, ς₀=1, ς_ν=\, J₀=0, ^=0,^=0,^=0,^=0},

de manera que la definición del algoritmo de cálculo según el primer objeto de esta invención queda

finalmente

En el mismo sentido, es objeto de esta invención un procedimiento para calcular un índice de cualificación de la condición o salud cardíaca de un individuo, caracterizado por:

a. El cálculo del número de eventosM' según el tercer objeto de esta invención a partir de la serie {^_¡}_¡=1 con la definición concreta del vector^ = t{(N + l)/ 2- j}_j=l ^ , donde N>1. b. El uso del índice m_sl/M , con m_sl = M' calculado en el punto anterior, directamente o en combinación con cualquier otro índice de cualificación, para determinar el nivel de salud cardíaca, por ejemplo pero no restringido al uso de m_sl / como índice directo de cualificación de la salud cardíaca.

De la misma forma, también es objeto de esta invención un procedimiento para calcular un índice de cualificación de la condición o salud cardíaca de un individuo, caracterizado por:

a. El cálculo del número de eventosM' según el tercer objeto de esta invención a partir de la serie {^_¡}_¡=1 ^con la definición concreta del vector A_N = ts> (2n j)IN}_{]=l N}, donde N puede variar entre N = 3 hasta N = 12, lo cual corresponde a una modulación sinusoidal del ritmo

7

HOJA DE REEMPLAZO (REGLA 26) cardíaco combinado con el ritmo respiratorio, donde t puede tener cualquier valor, por ejemplo, pero sin restricción, 1 o -1 .

b. El uso del índice m_S2 IM , con m_S2 = M' calculado en el punto anterior, directamente o en combinación con cualquier otro índice de cualificación, para determinar el nivel de salud cardíaca, por ejemplo pero no restringido al uso de m_S2 / como índice directo de cualificación de la salud cardíaca.

Igualmente, se considera objeto de esta invención un procedimiento para calcular un índice de cualificación de la condición o salud cardíaca de un individuo, caracterizado por:

a. El cálculo del coeficiente Μ Μ según el tercer objeto de esta invención a partir de la serie {^_¡ }_¡=1 ^ con la definición concreta del vector A_N = t j— 1,1, 0,...,θ| , donde N puede variar entre N = l hasta N = 20, correspondiente a un latido ectópico compensado, y donde t puede tener cualquier valor, por ejemplo, pero sin restricción, 1 o -1 .

b. El uso del índice m_E IM , con m_E = ' calculado en el punto anterior, directamente o en combinación con cualquier otro índice de cualificación, para determinar el nivel de salud cardíaca, por ejemplo pero no restringido al uso de m_£ / como índice directo de cualificación de la salud cardíaca.

Así mismo, es también objeto de esta invención un procedimiento para calcular un índice de cualificación de la condición o salud cardíaca de un individuo, caracterizado por:

a. El cálculo del coeficiente Μ Μ según el tercer objeto de esta invención a partir de la serie {^_¡ }_¡=1 con la definición concreta del vector A_N = t|N,— 1,...,— 1 j , donde N puede variar entre N = 2 hasta N = 20, correspondiente a una taquicardia paroxística regular, y donde t puede tener cualquier valor, por ejemplo, pero sin restricción, 1 o -1 .

b. El uso del índice m_TP I M , con m_TP = ' calculado en el punto anterior, directamente o en combinación con cualquier otro índice de cualificación, para determinar el nivel de salud cardíaca, por ejemplo pero no restringido al uso de m_TP I M como índice directo de cualificación de la salud cardíaca. Se considera también objeto de la presente invención cualquier uso para diagnosticar la funcionalidad, patología o nivel de calidad o salud de una máquina u organismo dotados con un motor u órgano de funcionamiento cíclico, que utilice cualquiera de los procedimientos descritos arriba.

Otro objeto de esta invención es cualquier uso empleado para emitir juicios prognósticos sobre la funcionalidad, patología o nivel de calidad o salud futuros de una máquina u organismo dotados con un motor u órgano de funcionamiento cíclico, que utilice cualquiera de los procedimientos descritos anteriormente.

8

HOJA DE REEMPLAZO (REGLA 26) Un último objeto de esta invención es cualquier uso para la descripción, representación gráfica compacta, e identificación gráfica de patrones específicos de funcionamiento de sistemas dinámicos, como por ejemplo pero no restringido a los sistemas económicos como el mercado de valores bursátiles, que utilice cualquiera de los procedimientos descritos anteriormente aplicados a:

a. Series secuenciales de valores {^_¡ }_¡=1 , obtenidos midiendo el valor de una cierta cantidad característica del sistema a intervalos de tiempo regulares o con una frecuencia determinada.

b. Series secuenciales de valores {^_¡ }_¡=1 , obtenidos midiendo el valor de una cierta cantidad característica del sistema siguiendo pautas temporales determinadas, no necesariamente regulares o constantes, para la adquisición o medida de los valores X_i .

DESCRIPCIÓN DE LAS FIGURAS

Figura 1 . Representación de la distancia de los vectores para un caso particular de insuficiencia cardíaca. Se muestran algunas líneas características de este caso específico y sus expresiones matemáticas. Además, se presentan ampliadas las regiones centrales de la gráfica con las zonas de interés remarcadas para apreciar los detalles.

Figura 2. (a) Gráfico en 4-dimensiones de la localización de M-3 vectores {Δ, }_{ι=1 m} __{n +1} normalizados con la media global de un HR normal (un adulto sano), (b, c) Lo mismo para dos sujetos con insuficiencia cardíaca crónica (tres ángulos de visión de cada uno). Los gráficos situados a la derecha en (b, c) muestran patrones universales cuando se proyecta con el ángulo apropiado, (d) Identificación de las líneas principales de la insuficiencia cardíaca en dos individuos diferentes con HF.

Figura 3. (a) Las diferentes regiones de actividad normal y de HF que se encuentran en este estudio: NSR (actividad normal, verde; Fantasía, cian); HF (magenta). (b1 -b3) Las diferentes regiones ocupadas por las cuatro situaciones que se encuentran en este estudio en el espacio multivariable

{log₁₀( 41), log₁₀ ( ), log₁₀(O_iV)}: NSR (n. a., verde; Fantasía, cian); MI (azul); SD (rojo); HF (magenta).

Para mayor claridad, el panel izquierdo proporciona tres situaciones (NSR n. a., Infarto de miocardio y HF), el panel central cuatro (NSR n. a., MI, SD y HF), y el panel derecho todas las situaciones estudiadas. N = 5 en todos los resultados de esta figura.

DESCRIPCIÓN DETALLADA DE LA INVENCIÓN

Un mapa de Poincaré es un grafo que consiste en la representación de un mapa de retorno o trayectoria secuencial de los valores que presenta una cierta variable en ciclos consecutivos. En concreto, un mapa de retorno bidimensional (2D) es una proyección plana donde las trayectorias complejas, con características multidimensionales (es decir, secuencias arrítmicas específicas) aparecen superpuestas e indistinguibles. Esta invención permite una generalización de dichos mapas de retorno. Entre otras posibilidades, se puede formular una variabilidad normalizada que usa un promedio móvil de orden N=5. Esta formulación se corresponde con la siguiente elección de los quince valores que definen el algoritmo base de los procedimientos descritos en esta invención:

9

HOJA DE REEMPLAZO (REGLA 26) { m = 0 , N = 5 , N₀ = 5 , N_l = 5 , N₂ = 5 , ε₀ = 1 , ε_ι = 1 , ε₂ = 0 , ς₀ = 1 , ς_ι = 1 , J₀ = 0 , J_l = 0 , Κ₀ = 0 , Κ_{ι =} 0 , Κ₂ = 0 }.

Mediante esta elección, que llamaremos NL (Normalización Local) se obtiene una expresión de la variabilidad secuencial V ^cuV^a representación

cumple las exigencias antes mencionadas para una representación ideal de la función cardiaca en términos de compacidad e integridad, y además, reduce al mínimo las influencias exógenas y las accidentales gracias a la normalización local resultante de los valores _M - El procedimiento propuesto se puede aplicar a cualquiera de los componentes del complejo cardíaco "pQRSt", aunque este ejemplo se centra en el análisis del componente principal, la serie de intervalos RR. En general, explorando el intervalo de orden desde N=2 a N=100, se observa que un medida de cuantificación de la variabilidad global, como es la norma o argumento del vector centro de masas de la distribución de puntos definidos por los vectores

S_j = _i } presenta recurrentemente un mínimo para N=5 en el caso de los

individuos sanos. Esto podría ser debido al hecho de que la frecuencia cardíaca promedio es aproximadamente cinco veces mayor que la frecuencia respiratoria en nuestra especie, hecho que da lugar a un subarmónico de orden cinco de la HRV, maximizando así la compensación global cuando N=5. Si esta hipótesis es correcta, las desviaciones de la normalidad deberían distinguirse óptimamente utilizando el orden N=5.

Afortunadamente, N=5 da lugar a la gráfica más completa, con posibilidades de representación gráfica completa, entre todos los órdenes posibles. En efecto, la expresión S_j = — lj permite

su representación en N-1 dimensiones ya que la información contenida en N-1 del total de N elementos proporciona toda la información del vector de N dimensiones. Esto es así porque el valor de cualquiera de los N elementos se puede obtener a partir de los otros N-1 ya que la suma de los N elementos es siempre nula, por definición. El vector de información secuencial en cuatro dimensiones resultante, que se puede representar gráficamente en 3D más color, da lugar a una valiosa herramienta descriptiva y comparativa desde el punto de vista clínico. Este gráfico permite la identificación de secuencias que originan patrones universales, cuya distribución y densidad proporcionan una información inmediata y completa del estado de la función cardíaca. Utilizando el procedimiento descrito, se han identificado hasta cinco secuencias universales, especialmente visibles en la representación gráfica (ver un ejemplo en la Figura 1 , con N=4, y otros ejemplos en la Figura 2, con N=5 y otra elección de parámetros), pero que no sólo están presentes cuando N = 5: todas se pueden encontrar en órdenes mayores, y de hecho se han obtenido sus expresiones generales para cualquier orden N.

Se ha verificado la universalidad de estas secuencias en una serie de bases de datos de HR, con 133 registros agrupados en cuatro situaciones básicas con características distintivas: (i) individuos en ritmo sinusal normal {normal sinus rhythm, NSR) durante su actividad normal, o en reposo en posición supina, mientras ven la película "Fantasía"; (ii) cardiopatía isquémica, concretamente infarto de miocardio (myocardial infarction, MI); (iii) insuficiencia cardíaca (heart failure, HF); y (iv) recuperación de muerte súbita (sudden death, SD). Al aplicarla a las citadas bases de datos, el porcentaje de aparición de cada una de las cinco secuencias, y la variabilidad primaria (cuya definición se verá posteriormente), proporcionan un conjunto unívoco y específico

10

HOJA DE REEMPLAZO (REGLA 26) de medidas capaces de valorar el estado del corazón en los registros incluidos en las bases de datos utilizadas y que están disponibles públicamente (ver Bases de Datos empleadas abajo).

Un registro Holter (HR) es un conjunto M de valores consecutivos que corresponden a los intervalos

RR, que puede ser expresado como x _} } ^ _m . Es fácil demostrar que el vector distancia generalizada en N dimensiones desde cualquier punto \X } formado por N intervalos RR a la línea de identidad (línea

J +k )k=0,...,N-l

de variabilidad cero, Piskorski & Guzik 2012), definida por el vector identidad {ΐ,Ι,. , .,ΐ}^ , se define como = \X , , _t - {X i _r \ . La secuencia de intervalos RR de un corazón sano "bailaría" alrededor de la línea de identidad, pero nunca llegaría a apoyarse en ella (incluso en las situaciones extremas de relax total o de ejercicio extremo, hay una variabilidad matemáticamente distinta de cero, aunque los RR sean aparentemente constantes). La expresión normalizada localmente de dicho vector distancia generalizada es precisamente la expresión de la elección de parámetros NL, es decir <5,≡ ( \ l

' N,_r D¡ = XX) la base de la normalización es la media loca ,j J / N,j X J,+,^k,,— 14| , donde

=0,...,#-1

N-l

(x)_N = ^ X _J+i■ En la Figura 1 se ha utilizado como cuarta dimensión un código simple de color que varía entre 0 y 1 (el código "Hue" del programa Mathematica® 9.0).

Otra elección de parámetros que se ha utilizado para este estudio es la siguiente:

{m = 0, N = 5 , N₀ = 5 , N_l = 5 , N₂ = 5 , ε₀ = 1 , ε_ι = 1 , ε₂ = 0 , ς₀ = 1 , ς_ι = 1 , J₀ = 0, J_l = 0 , Κ₀ = 0, Κ_{ι =} 0 , Κ₂ = 0 }, que da lugar a la siguiente expresión δ_} =

. A esta elección de parámetros la denominaremos NG (Normalización Global), y a la expresión particular de d para esta elección la denominaremos como . Esta expresión tiene las siguientes relaciones triviales con

D, = > -¾ L _Mv^∞n

IM \ I N J

La Figura 2 muestra la representación gráfica de tres HRs, analizados con el procedimiento propuesto, para N = 5: Fig. 2 (a) representa un sujeto sano, y la fig. 2 (b, c) dos pacientes con insuficiencia cardiaca crónica. Mientras que en el individuo sano da lugar a una forma densa y compacta alrededor del origen, los individuos con HF muestran líneas espaciales distintivas, que siguen secuencias preestablecidas. A la vista de esta figura se pueden extraer algunas conclusiones inmediatas. En primer lugar, las líneas expuestas son fundamentalmente rectas. En algunos casos (por ejemplo, Fig. 2 (b)), el espacio entre dos líneas está unido por un plano característico, pero hay un ángulo de visión especial (ángulo de proyección sobre el plano de vista) que sistemáticamente reduce las principales líneas rectas y planos a sólo tres en todos los casos [véase la Fig. 2 (c)].

1 1

HOJA DE REEMPLAZO (REGLA 26) Como sería esperable, las secciones de Poincaré normalizadas de la variabilidad de la base de datos NSR (aparentemente no patológicas) están relativamente centradas y homogéneamente distribuidas en torno a cero, mostrando una compensación matemática (aparentemente) con una cierta randomización, es decir, un núcleo más o menos esférico con distribución al azar, aproximadamente gaussiana. Una observación más detallada revela que todos los casos muestran una clara estructura en la cuarta dimensión (color), a lo largo de una dirección específica. Es importante destacar que esta dirección es exactamente la misma para todos los sujetos, con independencia de la forma de la distribución de puntos que representen los intervalos RR. Este hecho apunta a la existencia de un subarmónico, probablemente relacionados con el ciclo de la respiración, que refleja un automatismo del sistema simpático-parasimpático. Este efecto es relativamente raro o está casi ausente en los individuos con IM, y no se encuentra en individuos con HF y SD. La existencia de algún individuo de la base de datos NSR (alrededor del 12%) que presenta las características presentes en la HF se cuantifica y se discuten más adelante, y apoya el hecho de que la presencia del ritmo sinusal no descarta la presencia de alguna patología cardíaca.

La base de datos "Fantasía" muestra las mismas características (lyengar et al. 1996; Schimitt et al. 2007) que los individuos de la base de datos de NSR, y casi el mismo porcentaje de sujetos con características de HF. Además, los sujetos sanos mayores muestran claramente una menor variabilidad que los jóvenes.

Casi el 70% de los sujetos de la base de datos HF (Baim et al. 1986), y muchos de los de la base SD estudiados (Taddei et al. 1992) muestran las mismas líneas distintivas anteriormente mencionadas. El porcentaje de sujetos con nubes de puntos homogéneamente distribuidos es inverso proporcionalmente al compararlo con los individuos de la base de datos NSR: menos del 20%. Por otra parte, la repetición de los patrones y su extensión cumplen las características universales de la HRV. Por último, los individuos de la base de datos SD utilizada en este estudio, muestran patrones variables, totalmente distintos y visiblemente más irregulares que cualquier otro grupo. Algunos de ellos presentan las mismas líneas que los individuos con HF, pero la mayoría de ellos muestran una estructura muy compleja y aparentemente caótica.

Como una primera característica de la representación gráfica propuesta en esta invención tanto con la elección NL como con NG ( < . ο Δ_; ), es que se centra alrededor del origen por definición, de manera que puede considerarse la densidad de puntos en diferentes áreas como una firma identificativa específica de variabilidad. Por lo tanto, la norma del vector normalizado que define el centro de masa de la representación gráfica sería una medida primaria de la variabilidad para un enésimo orden subarmónico, y de su grado de compensación global (por ejemplo, en un HR de 24 horas para un ciclo circadiano). Sin embargo, nuestra experiencia demuestra que todo el registro de vectores da lugar a una compensación inadecuada, debido a la naturaleza inherente de la representación de Poincaré: puede comprobarse que el mismo intervalo RR arrítmico _y aparece como un componente común en N - 1 vectores A __k alrededor de la línea de identidad, que a menudo produce una compensación global aparente. Para evitar esto, se puede extraer un subconjunto de toda la serie tomando el índice j en saltos de N elementos (es decir, j = {l,l + N,l + 2N,...}), donde cada Χ aparece sólo una vez. La distancia desde el centro de masa de este subgrupo al origen, cuya representación gráfica es casi indistinguible de la de toda la serie excepto en sus diferentes densidades, viene dada por:

12

HOJA DE REEMPLAZO (REGLA 26)

M-N+l

donde A_{¡ k} es la componente A del vector Δ_; , y la notación del sumatorio ∑Δ_;. _¿ indica saltos de N en N en el índice i . Este coeficiente, que llamamos variabilidad primaria (primary variability, PV), se puede representar para cada individuo como una función de N. Nuestra experiencia demuestra que O^ es aproximadamente 10 a 100 veces más grande que la distancia al origen del centro de masa del conjunto original completo, hecho que amplifica enormemente la significancia de como se ha definido. También depende del número de latidos en la serie. En general, la mínima variabilidad corresponden a los individuos de NRS (Fantasía), seguidos de cerca por la variabilidades de los individuos con MI. La máxima variabilidad corresponde a SD y HF, siendo su distribución bastante similar. En valores intermedios se encuentra la variabilidad de los sujetos de la base de datos NSR con actividad normal. se muestra como una nueva medida cuantitativa con extraordinaria capacidad de diferenciación de los pacientes con distintas alteraciones de la función cardíaca, en combinación con la representación gráfica propuesta y las estructuras arrítmicas que derivan de ella.

A continuación ofrecemos una primera clasificación de los diferentes patrones que se pueden agrupar en secuencias o líneas universales. Por otra parte, vamos a demostrar que la densidad de puntos a lo largo de esas líneas proporciona una valiosa medida del estado de la función cardiaca. De hecho, cuando se aplica nuestro procedimiento a las bases de datos que hemos utilizado, que están disponibles libremente, da lugar a unos resultados con una alta especificidad descriptiva para cada grupo. Como una primera aproximación, en este trabajo nos limitamos a la identificación de las secuencias arrítmicas primarias más simples, que se pueden expresar en forma de línea, definida por el vector

parametrizada por una variable t arbitraria. Por supuesto, el vector de posición de un punto real (latido) sobre una línea arrítmica específica, puede corresponder a cualquier valor real positivo t . Esto refleja la capacidad intrínseca del procedimiento descrito: una línea arrítmica (o anomalía), que puede ser matemática y universalmente expresada, reúne todas las arritmias de la misma naturaleza, con independencia de la frecuencia cardíaca y la amplitud de la variabilidad. Se han descifrado varias secuencias primarias que se describirán a continuación y se ha calculado la densidad de puntos a lo largo de las líneas (secuencias) correspondientes. Estos resultados no excluyen la existencia de otras secuencias más complejas, con características específicas asociadas, que se determinarán en estudios futuros.

Una forma inmediata para evaluar la densidad de puntos correspondiente a una secuencia específica es cuantificar su presencia en % a lo largo del registro total de la serie de RR. Dado que la presencia de una secuencia concreta, obviamente, no es una cantidad exacta ni uniforme en todas las situaciones ni en todos los individuos, es preciso utilizar medios estadísticos para determinar la mayor o menor presencia de dicha secuencia en una condición cardíaca determinada. Una manera conveniente para representar la distribución de la presencia de una secuencia en una cierta condición es determinar el valor de F_¡ = il M_B frente y_t para la base de datos correspondiente, siendo JA el porcentaje de presencia de la secuencia

particular de un individuo, i es el rango de dicho individuo en particular, basado en su

13

HOJA DE REEMPLAZO (REGLA 26) puntuación y_t , y M_B es el número total de individuos de la base de datos, Esto vendría dado para cada secuencia, y se analizaría en combinación para cada situación.

Arritmia A1 : Un individuo sano con NSR debe exhibir una capacidad intrínseca para responder a cualquier demanda del organismo, mediante aceleraciones y desaceleraciones regulares de la frecuencia cardíaca dirigidas por el equilibrio simpático/parasimpático. Esta capacidad se debe reflejar en la apariencia de la forma más simple de HRV, que se puede expresar como una rampa lineal:

A\ = t{{N + \)l2 -j}^ _N

donde el signo positivo o negativo de t se aplica a la aceleración o deceleración de la frecuencia cardíaca respectivamente. Por ejemplo, para N = 4, con frecuencia cardíaca acelerada se formularía

^i; = t{i .5 ,0.5,-0.5,1.5}, con t > 0 ; para N = 5 y frecuencia cardíaca desacelerada

A\₅ ^~ = t{- 2,-1,0,1,2} , con t > 0 ; etc., donde los valores mayores de t indican un ascenso o descenso más brusco, o una rampa más pronunciada, de la frecuencia cardíaca, sin cambios en la estructura funcional de la variabilidad.

Esta es la forma dominante de la variabilidad cardíaca (HRV) en los individuos normales de la base de datos NSR, MIT BIH, tal como se muestra en la Tabla 1 . En efecto, esto se refleja en una forma ligeramente elipsoidal alrededor el origen en la dirección de la línea de Al₅ de cualquier gráfico en 4-dimensiones de un mapa Poincaré de quinto orden, que represente la HRV de un individuo normal con NSR, tal como hemos propuesto [véase la fig. 2 (a)].

Más importante aún: a medida que la presencia de esta forma de HRV disminuye, otras formas de arritmias patológicas, que se describirán a continuación, aumentan aproximadamente en la misma proporción. Este hallazgo apunta fundamentales a un hecho orgánico básico: esta arritmia es, en realidad, el grado de libertad básico de la HRV, y refleja la capacidad de adaptación de un organismo sano. Si alguna situación patológica deprime o limita este grado de libertad, se presentarán otras formas de HRV para compensar esta deficiencia. Es importante destacar que estas formas alternativas no son arbitrarias, especialmente en los individuos incluidos en la base de datos HF. De hecho, pueden ser consideradas como universales. En consecuencia, la forma específica de estos patrones de HRV alternativos debe estar vinculada con la afectación específica de un organismo, abriendo la puerta a nuevas formas de diagnóstico rápido no invasivo. Arritmia B1 (latido ectópico compensado): Los gráficos en 4 dimensiones de los HR pertenecientes a individuos con HF, por lo general presentan tres líneas (Figura 2 (d)) que puede ser fácilmente identificadas como (véase la figura 2 (d)):

m_{5 l} = t{- 1,1,0,0} , m_{5 2} = t{o -1,1,0}, m_{5 3} = ¿{0,0 -1,1}

En primer lugar, estas tres secuencias en realidad pertenecen a una única clase del tipo

t{.., 0,0,-1,1,0,0,...}

En segundo lugar, el valor medio de los cuatro componentes es cero para cualquier N> 2, y por lo tanto estas secuencias puede considerarse "compensadas". Esto significa que el último punto de la sección de Poincaré de enésimo orden, se asienta aproximadamente en la línea de identidad. En otras palabras, la secuencia arrítmica puede también considerarse "cerrada" o "concluida", ya que el último intervalo dura al mismo tiempo que la media local. Es destacable el hecho de que el número de intervalos con variabilidad cero, que rodean la

14

HOJA DE REEMPLAZO (REGLA 26) sucesión {- 1 , 1 } en la mayoría de los casos es mayor de dos (es decir, la línea {1 , 0, 0, - 1 } no aparece claramente). Por lo tanto, se puede concluir que hay una secuencia específica descrita por la línea siguiente:

Bl = t{..,0,0 -1,1,0,0,...} ,

que presenta una clara ubicuidad: este tipo de arritmia es definitivamente característica y claramente dominante en los individuos con HF, con una media de presencia alrededor de un orden de magnitud mayor en HF que en SD o MI, a pesar de que también es dominante en estos sujetos. Mientras que las variabilidades primarias (PV) de los individuos de bases de datos HF y SD son grandes y muy similares, lo que realmente distingue a los individuos con HF de los que se recuperan de SD, es la presencia de la arritmia B1 , mucho mayor en los primeros. Estas arritmias corresponden a latidos ectópicos aislados compensados; en este trabajo, no pretendemos identificar su origen cardíaco concreto, ya sea supraventriculares o ventriculares, ya que su naturaleza puede estar asociada a otra característica no estudiada aún.

Como puede verse, este análisis puede proporcionar una nueva base para elaborar una clasificación general basada en la naturaleza inherente, y fuertemente compensada, de estas arritmias y su capacidad de ser reducida a una sola estructura universalmente expresable. Por otra parte, se trata de arritmias relativamente raras o ausentes en los individuos de la base de datos NSR: en efecto, los latidos ectópicos pueden aparecer en los registros de los individuos normales, pero son relativamente raros. Sin embargo, este tipo de arritmias está más presente en los registros de individuos despiertos en reposo y en posición supina (registro "Fantasía" citado abajo) que los individuos normales durante la actividad normal. No obstante, el predominio de la arritmia B1 en la HF podría requerir una futura revisión del valor diagnóstico (Frolkis et al. 2003; y referencias relacionadas) de la presencia de latidos ectópicos en combinación con la capacidad de adaptarse sin problemas a las demandas de la actividad normal (adaptabilidad asociada a la presencia de arritmias A1). En la Tabla 1 se observa la fuerte correlación inversa entre la presencia de latidos ectópicos compensados (arritmias B1) y las arritmias A1 (±). Además de la Tabla 1 , la presencia opuesta de las arritmias A1 y B1 se ilustra claramente en la Figura 3.

Arritmia B2 (taquicardia regular paroxística): En la representación gráfica podemos encontrar la secuencia descrita por [Ver Figura 1 (a, b)]:

B2_{5 l} = t{ -1,-1 -1}

En muchos casos aparece una línea adicional [ver Figura 1 (c, d)]:

2¾ = í{- 1,-1 -1,-1}

En primer lugar, la aparición de estas secuencias introduce una fuente adicional de randomización, al combinarla con A1 . Esto proporciona una capacidad de adaptación adicional, que por otra parte puede estar muy limitada (aparecería una insuficiencia más intensa) al disminuir la presencia de arritmias A1 (adaptabilidad normal) en los individuos con patologías cardíacas. En segundo lugar, ninguna de estas secuencias está compensada. Esto significa que se deben considerar valores de N de orden superior para encontrar secuencias compensadas o terminadas. En este caso, tanto B2_{5 1} como B2_{5 2} son en realidad parte de la secuencia compensada:

B2 = t{4 -1,-1 -1,-1}

que en realidad pertenece a una sección de Poincaré de sexto orden (es decir, N = 6). Se observa que las secuencias compensadas se describen por líneas del tipo:

B2_N = t{(N - 2) -l,... -l} ,

15

HOJA DE REEMPLAZO (REGLA 26) por ejemplo:

£2₆ = ¿{4 -1 -1 -1 -1}

Β2_Ί = ί{5 -1 -1 -1 -1 -1}

,Β2₈ = {6 -1,-1,-1 -1,-1,-1}

B2₉ = t{l -1,-1 -1,-1 -1,-1,-1}

etc. Todas estas secuencias están casi tan presentes como la arritmia B1 en la HF, aunque la presencia de este tipo de arritmia disminuye a medida que aumenta N. Esta secuencia es difícilmente identificable si N aumenta por encima de 10. En realidad, la arritmia con mayor presencia es B2₆ , por eso llamamos a este tipo de arritmia genéricamente B2 , en lugar de B2₆ . El análisis de la presencia de esta arritmia en la HF demuestra que es tan omnipresente como B1 , con distinta dominancia de una sobre la otra dependiendo del individuo. Esta secuencia es compatible con una taquicardia paroxística regular tras una pausa, atribuible a un bloqueo auriculo-ventricular, que en algunos casos da lugar a un síndrome de Stokes-Adams. En la Tabla 1 aparece su frecuencia de presentación en cada situación. Al igual que la arritmia B1 , la arritmia B2 es característica de la HF. Es importante tener en cuenta que esta arritmia representa una pausa apreciable, seguido de una serie proporcional de frecuencia rápida para compensar la pausa.

Arritmia B3 (taquicardia paroxística regular II): Se pueden identificar varias secuencias alternativas a B2 como:

= t{- 1,4 -1 -1 -1}

^B36,2 = t{- 1,-1,4 -1,-1}

= í{- 1,-1 -1,4 -1}

B3_6A = í{- 1,-1,-1,-1,4},

o bien:

B3_{7 l} = t{- 1,5 -1,-1 -1,-1}

B3_{7 2} = í{- 1,-1,5,-1,-1,-1}

B3_{1 3} = t{- \ -1,-1,5 -1,-1}

B3_7A = t{- \ -1,-1 -1,5,-1}

B3_7A = t{- 1,-1 -1,-1 -1,5}

etc., que se pueden expresar para un índice general N como:

B3_{N ¡} = t{- 1, ... ,-1, N - 2 (en la posición i),-l, ... ,-l} .

Para un N determinado, todas las secuencias de B3 con índices / diferentes tienen la misma presencia, pero significativamente menor que B2. Esta secuencia podría asimilarse a situaciones intermedias uniendo dos secuencias consecutivas correspondientes a taquicardia paroxística regular con pausas relativas como en la arritmia B2. Llama la atención que esta arritmia es menos característica de la HF; en realidad, está tan presente en HF como en SD, con una distribución de probabilidad casi idéntica. Por otra parte, está

16

HOJA DE REEMPLAZO (REGLA 26) significativamente más presente en NSR que en los individuos con IM y en los de la base de datos NSR ("Fantasía").

En algunos casos aparece una secuencia peculiar como "sombras" de la arritmia B2 [véase, por ejemplo, Fig. 2 (d)], que puede ser identificado como las líneas:

B 3_{6 1} = t{- 4,1,1,1,1}

B3_{6 2} = t{\ -4,1,1,1}

53_{6 3} = t{l,l -4,1,1}

que puede escribirse, en general, como

B3_{N l} = t{l, ... ,1,-(N - 2)(en la posición i),l, ... ,l} ,

una secuencia complementaria en cierta forma de B3_N._2, ¡. No vamos a aportar un análisis cuantitativo de este tipo de arritmia, dada su complejidad.

Arritmia A2 (relacionada con la respiración): Es una secuencia relativamente presente, aunque mucho menos que la arritmia B1 o B2 en HF o SD, que aparece como la forma dominante de sub-arritmia en individuos normales o asintomáticos. Esta secuencia compensada se puede describir por la siguiente línea:

A2_N = t{sm(2 j)/N} _N ,

que representa una modulación sinusoidal de la frecuencia cardíaca a lo largo de un intervalo de N latidos. Este tipo de arritmia con compensación matemática, debe también reflejar una compensación fisiológica. Dado que es estadísticamente más frecuente con N = 5 que con cualquier otro orden, se puede concluir que está relacionada con la frecuencia respiratoria media en la especie humana. Si es o no más frecuente durante el sueño que en la actividad normal será objeto de futuros estudios. Su presencia en los registros se analiza en la Tabla 1 . Es relativamente dominante en NRS sobre los registros patológicos, por lo tanto es significativamente menos frecuente en individuos con HF y SD. Se puede deducir que esta arritmia es, como A1 , característicamente no patológica. En otras palabras, su presencia es compatible con un buen estado de salud cardiaca.

La combinación de las densidades de cada arritmia identificada puede constituir una valiosa firma característica de una situación específica, abriendo así el camino a nuevas exploraciones diagnósticas, cuyo significado está fuera del alcance de este ejemplo de aplicación de la invención propuesta. Un hallazgo fundamental es la relación inversa entre la presencia de ciertos tipos de arritmia, las que se pueden considerar indicativas de "salud", y las que se pueden denominar "patológicas". Concretamente, en el análisis de las bases de datos públicas que utilizadas, la arritmia A1 y la arritmia B1 son antagónicos: de hecho, la presencia relativa de una sobre la otra se invierte al pasar de una situación de normalidad a un estado patológico. En la Figura 3 se presenta un mapa universal en el que se analiza la presencia de ambas arritmias en los individuos de las bases de datos NSR y HF. Se puede observar una diferencia muy clara entre NSR y HF, basada en la presencia de A1 y B1 .

Además, el valor de la variabilidad primaria completa el conjunto de variables características para proporcionar el sello característico, distintivo, de cada situación: obsérvese en la Tabla 1 que las combinaciones de las presencias de cada arritmia y la PV configura una firma única y muy diferenciada de cada situación. La diferencia entre MI y SD es , pequeña para MI y grande para SD. La función de densidad de probabilidad multivariable para cada situación en el espacio de variables

17

HOJA DE REEMPLAZO (REGLA 26) {log₁₀( 41), log₁₀ ( ), log₁₀(O_iV)} vendría dada por la densidad de los puntos en que, en el espacio, correspondería a cada situación, como se aprecia en la Figura 3.

El valor clínico de esta representación podría ser enormemente significativo. De hecho, el efecto terapéutico de los fármacos con acción específica cardiovascular, dirigida a objetivos tales como el eje simpático/parasimpático (beta-bloqueantes adrenérgicos), quizás el eje neuro-hormonal renina-angiotensina- aldosterona (lECAs, ARA II, etc.), etc, podría dar lugar a la modificación de la presencia relativa de cada arritmia y desplazando la localización del gráfico correspondiente hacia la dirección de la región NSR. El valor pronóstico potencial de los resultados obtenidos al aplicar nuestro procedimiento a los HR es evidente. La observación de la representación gráfica en 4 dimensiones de un HR permite la fácil e inmediata identificación de alteraciones de la normalidad en personas sanas.

Futuros estudios clínicos ampliarán la profundidad de los conocimientos aprendidos con el análisis propuesto en esta invención, tales como la identificación de nuevas características y nuevos patrones arrítmicos generales, relacionados con otras patologías y situaciones, no necesariamente de origen cardíaco, por ejemplo, diabetes, hipertensión, hipotiroidismo, o incluso alteraciones psíquicas.

Tabla 1

Para concluir, hemos propuesto un procedimiento cuantitativo no sólo para proporcionar una representación universal de la variabilidad de la frecuencia cardíaca, sino también una herramienta clínica potencialmente útil para calificar la función cardiaca, los riesgos, y probablemente otros problemas de salud relacionados. En este trabajo, entre los muchos tipos de secuencia universales diferentes, probablemente existentes, hemos descrito matemáticamente algunos patrones destacados, con una estructura relativamente simple, identificando cinco tipos generales de arritmias. Las bases de datos utilizadas nos han permitido identificar las firmas características de los individuos incluidos en las bases de datos NSR, MI, HF y SD, relacionando estos tipos básicos de arritmias con las distintas situaciones de la función cardíaca. Como resultado fundamental, hemos demostrado cuantitativamente que dos de estas arritmias son característicos del estado de salud, mientras que otras tres son patológicas. Su presencia relativa en un individuo puede estar eventualmente relacionada con situaciones específicas, con la creciente evidencia clínica aportada por esta metodología que se acumulará en el futuro. Además, la evolución temporal de las estructuras arrítmicas en un paciente, visibles con la aplicación de nuestra metodología, puede proporcionar información muy valiosa acerca de su estado y/o su evolución clínica. Pensamos que se trata de una metodología de aplicación fácil, y que el potencial de sus resultados en el diagnóstico clínico no invasivo y el pronóstico podría ser vital.

18

HOJA DE REEMPLAZO (REGLA 26)

Claims

REIVINDICACIONES

1 .- Procedimiento para la obtención de datos útiles asociados al patrón de variabilidad de la frecuencia cardiaca (HRV), caracterizado por la realización de los siguientes pasos:

a) se mide y se registra un número M de intervalos de tiempo consecutivos {^, }_{¡=1 M} , correspondientes a ciclos de uno o más componentes de un complejo cardiaco "pQRSt" de un electrocardiograma, con una precisión igual o mejor que el 10% del valor medio del tiempo de ciclo, y siendo M mayor que 2;

b) se calcula la variabilidad sobre dichos M intervalos de una secuencia de vectores consecutivos } _{j m} de N componentes, de acuerdo con el algoritmo o transformación definida por la expresión:

y la siguiente notación:

^£_1 < m ) _m\

< )_£,, = ^l∑^ " ^con ) n\{m -ri)\

donde los siguientes parámetros son números enteros y su elección determina la forma final de la transformación mencionada:

{m^^^N^N^s^s^s^ g^g^J^J^K^K^K^, donde, adicionalmente:

- m es un índice natural que representa el orden de la variación discreta que se calcula;

- N es la dimensión o número de componentes de cada vector δ _} , siendo N > 2 ;

- N₀,N_l, yN₂ indican el número de valores que se usan para calcular el promedio local correspondiente indicado en la fórmula general del algoritmo;

- ε₀,ε_ι,ε₂ tienen valores binarios 0 ó 1 , e indican si los elementos correspondientes son, respectivamente, fijos o móviles en el cálculo de cada una de las componentes del vector δ _} ;

- ς_ϋ,ς tienen valores binarios 0 ó 1 , e indican si la media local es, respectivamente, fija o móvil;

- J₀ J_l indican el retraso o adelanto del primer elemento que se toma en el cálculo a partir del índice y ;

- K₀, K_l, K₂ indican el retraso o adelanto del primer elemento que se toma en la serie local correspondiente para calcular el promedio local indicado;

y donde se representa gráficamente en dos o más dimensiones la posición del punto indicado por los valores de las componentes de cada uno de los vectores <J . .

2.- Procedimiento según la reivindicación 1 , que comprende una etapa adicional de comparación de los datos útiles obtenidos con patrones de comportamiento asociados a una función vectorial ^4 = 1^ } _{^ n} , correspondiente a la parametrización de una secuencia cardíaca, donde los elementos a . son valores fijos o

19

HOJA DE REEMPLAZO (REGLA 26) funciones de una o más variables, y donde se realiza una etapa adicional de comparación de los datos útiles obtenidos con dicha función A, según los siguientes pasos:

a) se calcula el ángulo generalizado 6_¡ , cuyo coseno está dado por: cos(ft) = '„ >

\^AW\

donde el símbolo es la norma generalizada de un vector en N dimensiones, tal que:

\l

¡4=∑^{« ;}

b) se calcula el número de eventos ' tal que el ángulo^ es inferior a una tolerancia predeterminadas , siendo 0 < ε < 1 , de tal forma que la función A se explora en su espacio de existencia para encontrar dichos eventos en que #_¡ . < s ;

c) se calcula el coeficiente Μ M .

3. - Procedimiento según la reivindicación anterior, donde el cálculo de la serie \δ \ de vectores consecutivos de N dimensiones o componentes se realiza según las siguientes definiciones de parámetros:

{m = 0, N = 5, N₀ = 5, N_l = 5, N₂ = 5, ε₀ =1, ε_ι =1, ε₂ =0, ς₀ =1, ς_ι =1, J₀ =0, J_l =0, ^ =0, ^ =0, ^ =0},

de manera que la definición de la variabilidad sobre dichos M intervalos es:

4. - Procedimiento según la reivindicación anterior, donde se realizan, adicionalmente, los siguientes pasos:

a) se calcula el número de eventos a partir de la serie { _¿}.₌₁ con la definición concreta del vector^ = t{(N + \)/2-j}_{j=1 N} , donde N>1;

b) se utiliza el índice m_sl/M , con m_si = ' calculado en el punto anterior, para la determinación de la existencia de patrones de comportamiento asociados a la función vectorial A .

5. - Procedimiento según la reivindicación 3, donde se realizan, adicionalmente, los siguientes pasos: a) se calcula el número de eventos M' a partir de la serie { _¿}.₌₁ con la definición concreta del vector A_N = í{sin(2^ · j)/N}_{j=1 n}, donde N puede variar desde N = 3 hasta N = 12 , correspondiente a una modulación sinusoidal del ritmo cardíaco combinado con el ritmo respiratorio, donde t puede tener cualquier valor;

b) se utiliza el índice m_S21 M , con m_S2 = ' calculado en el punto anterior, para la determinación de la existencia de patrones de comportamiento asociados a la función vectorial A .

20

HOJA DE REEMPLAZO (REGLA 26)

6.- Procedimiento según la reivindicación 3, donde se realizan, adicionalmente, los siguientes pasos: a) se calcula el coeficiente M'/M a partir de la serie {^ }_{¡=1 M} con la definición del , donde N puede variar desde N = 1 hasta N = 20 , correspondiente a un latido

ectópico compensado, y donde i puede tener cualquier valor;

b) se utiliza el índice m_E l M , con m_E =M' calculado en el paso anterior, para la determinación de la existencia de patrones de comportamiento asociados a la función vectorial A _N.

7 - Procedimiento según la reivindicación 3, donde se realizan, adicionalmente, los siguientes pasos: a) se calcula el coeficiente M'/M a partir de la serie {^ }_{¡=1 M} con la definición concreta del vector A_N = , donde N puede variar desde N = 2 hasta N = 20 , correspondiente a una

taquicardia paroxística regular, y donde i puede tener cualquier valor;

b) se utiliza el índice m_TP I M , con m_TP =M' calculado en el punto anterior, para la determinación de la existencia de patrones de comportamiento asociados a la función vectorial A _N.

8. - Procedimiento según cualquiera de las reivindicaciones anteriores, donde el componente del complejo cardiaco "pQRSt" es el intervalo RR de un electrocardiograma.

9. - Procedimiento según cualquiera de las reivindicaciones anteriores, donde el registro de un número

M de intervalos de tiempo consecutivos

correspondientes a ciclos de un componente de un complejo cardiaco "pQRSt" se realiza con una precisión mayor que el 0.01 % del valor medio del tiempo de ciclo.

10.- Procedimiento según la reivindicación 3, donde 0 < ε < 0.1 .

1 1 .- Procedimiento según cualquiera de las reivindicaciones 6-8, donde t es 1 ó -1

21

HOJA DE REEMPLAZO (REGLA 26)