WO2019160504A1 - Système et procédé pour évaluer un risque d'événement clinique sur la base de la complexité de la fréquence cardiaque - Google Patents
Système et procédé pour évaluer un risque d'événement clinique sur la base de la complexité de la fréquence cardiaque Download PDFInfo
- Publication number
- WO2019160504A1 WO2019160504A1 PCT/SG2019/050083 SG2019050083W WO2019160504A1 WO 2019160504 A1 WO2019160504 A1 WO 2019160504A1 SG 2019050083 W SG2019050083 W SG 2019050083W WO 2019160504 A1 WO2019160504 A1 WO 2019160504A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- heart rate
- data
- subject
- hrc
- index
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 63
- 230000004044 response Effects 0.000 claims abstract description 8
- 239000002131 composite material Substances 0.000 claims description 5
- 230000000007 visual effect Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 description 11
- 238000004891 communication Methods 0.000 description 8
- 230000000747 cardiac effect Effects 0.000 description 7
- 238000001514 detection method Methods 0.000 description 7
- 230000000284 resting effect Effects 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 206010019280 Heart failures Diseases 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 5
- 230000006866 deterioration Effects 0.000 description 5
- 230000036541 health Effects 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 5
- 239000013598 vector Substances 0.000 description 5
- 208000004476 Acute Coronary Syndrome Diseases 0.000 description 4
- 0 CC(CC*(C)CC1*C1)C(C)C1[C@]2C1CC(C1)C1*2 Chemical compound CC(CC*(C)CC1*C1)C(C)C1[C@]2C1CC(C1)C1*2 0.000 description 4
- 230000002411 adverse Effects 0.000 description 4
- 238000013459 approach Methods 0.000 description 4
- 230000007423 decrease Effects 0.000 description 4
- 230000035939 shock Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 208000019622 heart disease Diseases 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000035479 physiological effects, processes and functions Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 208000020446 Cardiac disease Diseases 0.000 description 2
- 206010049418 Sudden Cardiac Death Diseases 0.000 description 2
- 210000003403 autonomic nervous system Anatomy 0.000 description 2
- 230000003906 autonomic nervous system functioning Effects 0.000 description 2
- 230000007211 cardiovascular event Effects 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000003862 health status Effects 0.000 description 2
- 230000004217 heart function Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012502 risk assessment Methods 0.000 description 2
- 206010003658 Atrial Fibrillation Diseases 0.000 description 1
- 241000282472 Canis lupus familiaris Species 0.000 description 1
- 206010007559 Cardiac failure congestive Diseases 0.000 description 1
- 208000032862 Clinical Deterioration Diseases 0.000 description 1
- 208000034657 Convalescence Diseases 0.000 description 1
- 208000003037 Diastolic Heart Failure Diseases 0.000 description 1
- 241000282326 Felis catus Species 0.000 description 1
- 208000010496 Heart Arrest Diseases 0.000 description 1
- 101000630267 Homo sapiens Probable glutamate-tRNA ligase, mitochondrial Proteins 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 102100026125 Probable glutamate-tRNA ligase, mitochondrial Human genes 0.000 description 1
- 206010039966 Senile dementia Diseases 0.000 description 1
- 208000006011 Stroke Diseases 0.000 description 1
- 206010049447 Tachyarrhythmia Diseases 0.000 description 1
- 208000001871 Tachycardia Diseases 0.000 description 1
- 208000030886 Traumatic Brain injury Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000002567 autonomic effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000005189 cardiac health Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 230000000994 depressogenic effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000002503 metabolic effect Effects 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000001734 parasympathetic effect Effects 0.000 description 1
- 230000037081 physical activity Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 208000020016 psychiatric disease Diseases 0.000 description 1
- 238000013515 script Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 210000001013 sinoatrial node Anatomy 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000002889 sympathetic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000000115 thoracic cavity Anatomy 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 230000008733 trauma Effects 0.000 description 1
- 230000009529 traumatic brain injury Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 230000002861 ventricular Effects 0.000 description 1
- 230000002618 waking effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02405—Determining heart rate variability
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7246—Details of waveform analysis using correlation, e.g. template matching or determination of similarity
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7271—Specific aspects of physiological measurement analysis
- A61B5/7275—Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
- G16H10/60—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for patient-specific data, e.g. for electronic patient records
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
Definitions
- the present disclosure generally relates to assessing clinical event risk based on heart rate complexity. More particularly, the present disclosure describes various embodiments of a system and a method for assessing the risk of clinical events based on heart rate complexity of a subject.
- HRV Heart rate variability
- ECG electrocardiogram
- HRC heart rate complexity
- perturbations of this healthy state are associated with changes in HRC and have been associated with a wide range of conditions including heart failure, senility, psychiatric disorders, traumatic brain injury, and strokes. This is postulated to be due to derangements in the balance of the autonom ic system and have been shown to correlate with severity as well as clinical outcomes in both cardiac disease and non- cardiac conditions such as trauma and stroke.
- the system comprises a physiological sensing device connectable to the subject; and an electronic device communicable with the physiological sensing device, the electronic device comprising a processor configured for performing the method.
- the method comprises: receiving, from the physiological sensing device, a physiological dataset measured from the subject, the physiological dataset comprising a set of heart rate data; generating a heart rate time series from the heart rate data; calculating a set of entropy data from the heart rate time series; calculating an HRC index from the entropy data; comparing the HRC index against predefined reference data; and generating a risk signal in response to a determination from said comparing that the HRC index is indicative of increased risks of clinical events to the subject.
- HRC index can be used as an indicator of increased risks of clinical events to the subject.
- the HRC index serves as a diagnostic modality that helps to quantify a subject’s health state.
- Clinicians can rely on the HRC index with other physiological data to look for derangements in the physiology of the subject that may predict impending significant adverse cardiovascular events. This allows the clinicians to intervene early to avert further worsening.
- Figure 1 A and Figure 1 B are schematic illustrations of a system for assessing risks of clinical events based on heart rate complexity of a subject, in accordance with embodiments of the present disclosure.
- Figure 2 is a flowchart illustration of a computerized method system for assessing risks of clinical events based on heart rate complexity of a subject, in accordance with embodiments of the present disclosure.
- Figure 3 illustrates the average profile parameters of the patients involved in the clinical study of the system and method in accordance with embodiments of the present disclosure.
- Figure 4 is a graph illustration of the results of the clinical study.
- FIG. 5 is a block diagram illustration of the technical architecture of a computer, in accordance with embodiments of the present disclosure.
- depiction of a given element or consideration or use of a particular element number in a particular figure or a reference thereto in corresponding descriptive material can encompass the same, an equivalent, or an analogous element or element number identified in another figure or descriptive material associated therewith.
- References to“an embodiment / example”,“another embodiment / example”,“some embodiments / examples”,“some other embodiments / examples”, and so on, indicate that the embodiment(s) / example(s) so described may include a particular feature, structure, characteristic, property, element, or limitation, but that not every embodiment / example necessarily includes that particular feature, structure, characteristic, property, element or limitation.
- repeated use of the phrase “in an embodiment / example” or “in another embodiment / example” does not necessarily refer to the same embodiment / example.
- the terms“a” and“an” are defined as one or more than one.
- the use of 7” in a figure or associated text is understood to mean“and/or” unless otherwise indicated.
- the term“set” is defined as a non-empty finite organization of elements that mathematically exhibits a cardinality of at least one (e.g. a set as defined herein can correspond to a unit, singlet, or single-element set, or a multiple-element set), in accordance with known mathematical definitions.
- the recitation of a particular numerical value or value range herein is understood to include or be a recitation of an approximate numerical value or value range.
- a component or a module may be, but is not limited to being, a process running on a processor, a processor, an object, an executable, a thread of execution, a program, and/or a computer.
- an application running on a controller and the controller can be a component / module.
- One or more components / modules may reside within a process and/or thread of execution.
- a component / module may be localized on one computer and/or distributed among a plurality of computers.
- FIG. 1A and Figure 1 B there is a system 100 for assessing risks of clinical events based on heart rate complexity (FIRC) of a subject 1 10.
- the subject 1 10 may be a human subject or a patient.
- the subject 1 10 may be a non-human subject or a living animal such as pet dogs and cats.
- FIRV heart rate variability
- the system 100 includes a physiological sensing device 120 connected or connectable to the subject 1 10 for measuring physiological data from the subject 1 10.
- the physiological sensing device 120 may be a heart monitoring device such as an implantable cardiac device, implantable cardioverter defibrillator (ICD), implantable loop recorder, pacemaker, or electrocardiogram (ECG) monitor.
- the physiological sensing device 120 may be a clinical device used for longer-term external electrocardiogram (ECG) or heart rate monitoring, such as an ambulatory ECG recorder or a consumer wearable (e.g. an external wearable device) with heart rate monitoring capabilities.
- the physiological sensing device 120 includes a heart rate sensor 122 for measuring heart rate data of the subject 1 10, and optionally includes a motion sensor 124 for measuring motion data of the subject 1 10.
- the physiological sensing device 120 further includes a database 126 for storing data, e.g. the heart rate data / motion data, measured by the physiological sensing device 120 from the subject 1 10.
- the system 100 further includes an electronic device 130 communicable with the physiological sensing device 120.
- the electronic device 130 includes a processor for executing program instructions, such as to perform a computerized method for assessing risks of clinical events based on HRC of the subject 1 10.
- the electronic device 130 further includes a database 132 for storing data, e.g. the heart rate data / motion data, measured by and communicated from the physiological sensing device 120.
- the electronic device 130 is a computer device 130a of the subject 1 10.
- the computer device 130a include a desktop device, mobile device, personal digital assistant (PDA), mobile phone, smartphone, tablet, phablet, and laptop computer.
- the computer device 130a may be communicable with the physiological sensing device 120 via wireless / contactless communication protocols such as Bluetooth and Wi-Fi.
- wireless / contactless communication protocols such as Bluetooth and Wi-Fi.
- the physiological sensing device 120 and the computer device 130a may share the same Wi-Fi network and data is communicable therebetween.
- the physiological sensing device 120 such as if it is an external wearable device, may also be communicable with the computer device 130a via wired communications such as USB.
- the physiological sensing device 120 stores the physiological data on the database 126 before communicating the physiological data to the computer device 130a, such as upon request by the subject 1 10, for storing on the database 132.
- the electronic device 130 is a remote server 130b and the system 100 further includes a computer device 140 of the subject 1 10.
- This computer device 140 is similar to the computer device 130a described above with reference to Figure 1A.
- the computer device 140 functions as an intermediary device between the physiological sensing device 120 and the remote server 130b, such that the physiological sensing device 120 is communicable with the remote server 130b via the computer device 140.
- the computer device 140 may be communicable with the physiological sensing device 120 via wireless / contactless communication protocols such as Bluetooth and Wi-Fi.
- the computer device 140 may be communicable with the remote server 130b across a wide area network (WAN) or the Internet.
- the computer device 140 may include a database 142 for storing data, e.g. the heart rate data / motion data, measured by and communicated from the physiological sensing device 120.
- the remote server 130b is a computer server that is located remotely away from the computer device 140.
- the remote server 130b is a physical or cloud data processing system on which a server program runs.
- the server may be implemented in hardware or software, or a combination thereof.
- the server includes computers, laptops, mini- computers, mainframe computers, any non-transient and tangible machines that can execute a machine-readable code, cloud-based servers, distributed server networks, and a network of computer systems.
- the physiological sensing device 120 communicates the physiological data to the computer device 140 for storing on the database 142.
- the computer device 140 subsequently communicates the physiological data to the remote server 130b at a later time for storing on the database 132.
- FIG 2 there is shown a computer-implemented or computerized method 200 for assessing risks of clinical events based on HRC of the subject 1 10.
- the method 200 is implemented on and performed by the electronic device 130, i.e. the computer device 130a or remote server 130b, which comprises various modules / components, including the processor for executing program instructions, for performing various steps of the method 200.
- a software application may be installed and executed on the electronic device 130 to perform the method 200.
- the software application may be implemented on Java and the electronic device 130 has suitable specification to execute the software application.
- the electronic device 130 has a 1 GFIz processor, 256 GB RAM, and 1 GB hard disk space.
- the subject 1 10 may need to setup his/her profile to include profile parameters such as gender, date of birth, age, weight, height, and health background.
- the subject profile parameters may be used to determine reference conditions for the risk assessment.
- the method 200 includes a step 202 of receiving, from the physiological sensing device 1 10, a physiological dataset measured from the subject 1 10, the physiological dataset including a set of heart rate data measured by the heart rate sensor 122.
- the physiological dataset optionally includes a set of motion data measured by the motion sensor 124.
- the method 200 further includes a step 204 of generating a heart rate time series from the set of heart rate data.
- the heart rate time series is defined as the consecutive R-R intervals on ECG records.
- the method further includes a step 206 of calculating a set of entropy data, such as multi-scale entropy (MSE) values, from the heart rate time series.
- the method 200 further includes a step 208 of calculating an FIRC index from the entropy data.
- FIRC indices are parameters or scores that can be used to assess the risks of clinical events to the subject 1 10.
- the method 200 further includes a step 210 of comparing the FIRC index against predefined reference data.
- the predefined reference data may be arbitrarily defined and/or derived from the reference conditions of the subject 1 10 and/or a population of subjects 1 10.
- the method 200 further includes a step 212 of determining, from said comparing in the step 210, if the HRC index is indicative of increased risks of clinical events to the subject 1 10. If yes, the step 212 proceeds to a step 214 of generating a risk signal in response to the determination from said comparing that the HRC index is indicative of increased risks of clinical events to the subject 1 10. Conversely, if the HRC index is not indicative, the step 212 proceeds to a step 216 of analyzing the next physiological dataset measured from the subject 1 10 and returns to the step 202.
- the method 200 may further include presenting the risk signal to alert the subject 1 10 of the increased risks of clinical events.
- the electronic device 130 is the computer device 130a having an audio and/or visual module for presenting the risk signal, such as in the form of visual graphics and sound alarms.
- the electronic device 130 is the remote server 130b configured for communicating the risk signal, such as a message, to the computer device 140 to thereby present the risk alert or message to the subject 1 10.
- the benchmark for optimal health can be defined through multi-scale entropy of a continuous heart rate time series.
- Entropy measurement of a heart rate time series can refer to sample entropy, approximate entropy, or other entropy definitions that measures the entropy data in a heart rate time series
- Multi-scale entropy refers to the entropy values at different levels of coarse-graining time scales.
- the heart rate time series x t is transformed into a new heart rate time series by coarse graining as shown in Equation 1 below for a scale m that is larger than 1 .
- m represents the level of coarse graining of x t and the coarse-grained heart rate time series is represented by y t>m ⁇
- the set of entropy data is generated from the coarse-grained heart rate time series y t,m and includes the entropy values (represented by E m ) calculated from the heart rate time series y t,m .
- the optimal healthy status is defined as the E m value from a heart rate time series y t,m having the same E m value for all the integer values of m.
- This E m value is represented by Emax and the Emax value can be calculated through numerical methods.
- Emax the E m value
- the amount of deviation in the E m values is correlated with the severity of the illness and can be a measure of the risk of future adverse clinical events.
- the amount of deviation can be transformed into the HRC index to signify the degree of deviation from a healthy state.
- the HRC index can be calculated from the set of entropy data calculated from the heart rate time series yt ,m .
- the HRC index represented by I t is calculated as the average of entropy values at scales m from 21 to 40 inclusive and weighted by a fixed value Emax, as shown in Equation 2 below. Approximate entropy may be used to calculate Emax with threshold standard deviation of 0.15 and data length of 2. The value of Emax is calculated to be 1 .8.
- the GRACE score is a well-validated cardiac risk scoring system to risk stratify patients diagnosed with acute coronary syndrome (ACS) to estimate their in-hospital and 6- month to 3-year mortality.
- ACS acute coronary syndrome
- a clinical study was performed on patients diagnosed with ACS to validate the HRC indices against the GRACE scores.
- Figure 3 illustrates a table 300 of the average profile parameters of the patients involved in the clinical study. Higher risk patients had more depressed HRC compared to lower risk patients, which was found to be most evident during waking hours and also in the first 24 hours after their presentation. Their HRC also remained largely the same while low risk patients had a decline in their HRC during their convalescence as inpatients.
- Figure 4 illustrates a graph 400 of the correlation and validation of the HRC indices against the GRACE scores.
- the R-squared or R 2 value and the p value from the graph 400 shows that the HRC indices correlates well with the GRACE scores, suggesting that the HRC indices can be used to risk stratify patients with diagnosed ACS and potentially other diseases.
- the correlation of the HRC indices with the GRACE scores was independent of left ventricular ejection fraction (LVEF).
- the HRC index I t is calculated from the entropy data and one or more profile parameters of the subject 1 10.
- the subject profile parameters may include the subject’s age derived from the initial profile setup by the subject 1 10.
- the multi-scale entropy (represent by MSE) and age (represented by Age) values are used together to calculate the HRC index l t , as shown in Equation 3 below ai and bi represent coefficients which may be arbitrarily defined or calculated through numerical methods. Incorporating the logarithmic function of the subject’s age to calculate the HRC index I t , improved the correlation from approximately 75% to approximately 80%.
- the physiological dataset includes the set of heart rate data and further includes a set of motion data measured by the motion sensor 124.
- the motion data may be used to reduce non-linearity in the heart rate data, such that the more useful non-linearity in the heart rate time series can be retrieved to improve the correlation and achieve higher prediction accuracies.
- the motion data can be used for activity context sensing to complement the heart rate data, since the heart rate data is variable in purely the time domain.
- the motion data allows for the detection of physical activity carried out by the subject 1 10.
- the motion data can be used to identify the heart rate data during certain periods, such as the resting period. Light activity such as walking would cause a linear trend in the heart rate data.
- the motion data can be used to remove the linear trend caused by the activity while retaining the linear trend intrinsic to the heart rate variability.
- the activity level (represented by m t ) is determined from the motion data and introduced to the heart rate time series x t generated from the measured set of heart rate data to thereby transform x t into a new heart rate time series (represented by yi), as shown in Equation 4 below.
- 3 2 represents a coefficient which may be arbitrarily defined or calculated through numerical methods.
- convolutional kernels may be used on the motion data instead of a simple linear correlation with the coefficient 32 as shown in Equation 4.
- the entropy data is generated from the new heart rate time series y t and the HRC index I t is calculated from the entropy data, as described above in the method 200.
- the HRC index l t is calculated from motion data together with the multi-scale entropy (represent by MSE) and age (represented by Age) values to further improve its correlation and accuracy.
- the motion data is regressed into Equation 3 which is then transformed into Equation 5 below a represents a coefficient which may be arbitrarily defined or calculated through numerical methods.
- the MSE values may change immediately after a change in the activity state of the subject 110.
- the amount of change in the MSE values can be an indicator of cardiac health and risk.
- the motion data may indicate a sudden change in physical or activity state from rest to walking.
- MSEi a first or current MSE value represented by MSEi for the rest state
- MSE2 a second or next MSE value represented by MSE2 for the walking state.
- the MSE difference MSE2 - MSE1 can be used as the indicator.
- the HRC analysis and calculation of the HRC indices can thus be modified to reflect changes in activity levels of the subject 1 10 during the awake periods and the resting / sleeping periods.
- the basal state of the human body of the subject 1 10 during the awake periods when the subject 1 10 is not engaged in extreme states of activity or inactivity is likely to best reflect the adaptability of physiological systems.
- the optimal balance of the physiological systems, or lack thereof, is most evident in the basal state during the awake periods.
- the basal state may be defined as the resting metabolic state of the body early in the morning after fasting for a minimum of 12 hours. During the resting / sleeping periods, many of the bodily functions are inactive or suppressed. Conversely, during intense physical exercise, the body’s metabolism is stretched well beyond its baseline level.
- the HRC index is significantly higher during the awake periods than the resting / sleeping periods.
- the HRC indices work optimally to indicate or predict increased risks of clinical events, such as cardiac risk, for subjects 1 10 who are awake.
- the motion data can be used to determine the physical status of the subject 1 10 and thereby identify the time periods for calculating the HRC indices.
- the set of heart rate data is divided into a plurality of subsets of heart rate data, each subset of heart rate data measured during a predefined time period. Consequently, the heart rate time series generated from the set of heart rate data may include a plurality of segments, each segment generated from one of the subsets of heart rate data. The division and segmentation would allow for calculation of the HRC indices for different time periods, as described further below.
- the heart rate time series x t generated from the measured set of heart rate data is cleaned to remove noise.
- the cleaned heart rate time series is divided into a plurality of subsets represented by X k according to the time of the day during which the respective subset of heart rate data was measured.
- the heart rate data is batched into 6-hour periods, such as 0000 to 0600 hours, 0600 to 1200 hours, 1200 to 1800 hours, and 1800 to 0000 hours, so that there are four subsets of the heart rate data per day.
- a moving window approach is used every time the FIRC index is calculated. For example, when the FIRC index is calculated at time t, the subset X k contains the heart rate data between time t hours and time t - 6 hours.
- the heart rate time series x t is transformed into a new heart rate time series y t,m by coarse graining where m represents the level of coarse graining of x t.
- the coarse graining is done mathematically through the average value of every set of m values, as shown in Equation 1 above.
- a moving window approach is used to obtain more data points in the coarse-grained data, as shown in Equation 6 below.
- the coarse-grained heart rate time series y t,m is segmented into a plurality of segments according to the subsets Xk and the segments are represented by Yk, m .
- the standard deviation of the respective heart rate time series segment Yk,m is also calculated and is represented by c3 ⁇ 4.
- the set of entropy data including the E m values is generated from the coarse-grained heart rate time series y t,m.
- the set of entropy data includes a plurality of subsets of entropy data represented by Ek,m, each subset of entropy data calculated from one of the segments of the heart rate time series Yk, m .
- sample entropy is used to calculate the entropy data or entropy values Ek,m.
- two series of vectors Uk,m and Vk, m are formed as shown in Equations 7 and 8 below.
- parameters ut,m and v t,m are denoted according to Equations 9 and 10 below so that Equations 7 and 8 can be transformed into Equations 11 and 12 below.
- each HRC index I is calculated as the average of the respective entropy values Ek,m at scales m from 21 to 40 inclusive and weighted by a fixed value E max .
- Approximate entropy may be used to calculate E max with threshold standard deviation of 0.15 and data length of 2. The value of E max is calculated to be 1 .8. [Equation 14]
- the correlation and prediction accuracy of the HRC indices may be improved by introducing other parameters such as the motion data and age of the subject 1 10. These parameters may be introduced to modify the heart rate data and/or modify the HRC indices .
- the motion data is used to exclude the sleeping / resting period from the heart rate data, since the awake period has higher correlation with cardiac risk than the sleeping / resting period.
- the motion data can also be used to reduce the non-linearity in the heart rate data.
- the motion data or activity level m t is introduced to transform the heart rate time series x t into a new heart rate time series y t , as shown in Equation 15 below.
- a3 represents a coefficient or normalization constant which may be arbitrarily defined or calculated through numerical methods.
- convolutional kernels may be used on the motion data instead of a simple linear correlation with the coefficient a ⁇ as shown in Equation 4.
- an initial HRC index I k is modified with the value of the subject’s age.
- the initial HRC index is multiplied with the inverse of the logarithmic function of the subject’s age, as shown in Equation 16 below.
- the modified HRC index is represented by I k . This modification improves the correlation of the modified HRC index I k with the GRACE score.
- the method 200 further includes calculating a composite index from the HRC index and the motion data, wherein the composite index is assistive in the determination that the HRC index is indicative of increased risks of clinical events to the subject 1 10, such as for prediction of cardiac risk.
- the composite index is represented by f(I k , m t ), where f is any function with two variables - in this case the HRC index I k and the motion data m t .
- the HRC indices are calculated for different time periods based on the respective heart rate time series segment Yk,m and respective heart rate data subset Xk. Accordingly, the HRC indices may be arranged chronologically and include a current HRC index and a previous HRC index k-i.
- the current HRC index is calculated for a first or current time period and the previous HRC index -i is calculated for a second or previous time period preceding the first time period.
- the first and second time periods may have identical durations so that the current and previous HRC indices are consistently calculated.
- the method 200 further includes determining a HRC index difference between HRC index and a sequential HRC index, such as the previous HRC index or the next HRC index.
- the HRC index difference is between the current HRC index and the previous HRC index k-i.
- the previous HRC index k-i may be associated with a basal state of the subject 1 10 and the current HRC index I k may be associated with an active state of the subject 1 10.
- the HRC index difference can be used to indicate or predict increased risks of clinical events to the subject 1 10, specifically if there is a decrease from the previous HRC index -i to the current HRC index I k .
- the decrease may be a sudden drop or a gradual decrease.
- the HRC index difference is compared to the predefined reference data, specifically a predefined difference reference.
- the HRC index difference is indicative of the increased risks if it is at least 10%, i.e. the previous HRC index -i is at least 10% more than the current HRC index Ik. Said indication would alert the subject 1 10 of possible impending adverse clinical events, such as risk of cardiac arrest.
- the risk signal is generated in response to a determination from that the HRC index difference is indicative of increased risks of clinical events to the subject 1 10.
- various embodiments herein describe the calculation of the HRC index using the method 200 wherein the HRC index can be used as an indicator of increased risks of clinical events to the subject 1 10.
- the HRC index serves as a diagnostic modality that adds a new layer of information to quantify a subject’s health state.
- the HRC index can determine the activity level and the balance of the autonomic nervous system, and allows clinicians to look for derangements in the physiology of the subject 1 10 that may predict impending significant adverse cardiovascular events.
- the HRC index may allow clinicians to detect worsening heart failure in patients before it is clinically evident by looking for derangements in the autonomic nervous system function so that there can be early intervention to avert further worsening.
- the system 100 and method 200 can be implemented with different types of physiological sensing devices 120 for various industrial applications and for assessing risks of various clinical events.
- the physiological sensing device 120 is an ICD and the method 200 can be performed for early detection of deterioration or derangement of the autonomic nervous system function that can forewarn of impending clinical events or deterioration, such as if there is a marked deterioration or change in the HRC index.
- Early detection allows for timely intervention to avert an episode of tachyarrhythmia or even a shock. Early detection may also help to reduce the incidence of inappropriate shocks by helping to detect atrial fibrillation which also markedly changes the HRC index.
- the system 100 can be used to complement the capabilities of contemporary home monitoring systems which some patients may have. Clinicians may remotely access the home monitoring systems to monitor the patients’ health status and remotely detect any deterioration in the patients’ health status. If there is any deterioration, the clinicians can be alerted quickly to the need for early intervention, thereby further reducing shocks for the patients.
- the physiological sensing device 120 may be a pacemaker and the method 200 can be performed for detection of worsening heart function.
- the physiological sensing device 120 may instead be an implantable loop recorder and the method 200 can be performed for detection of worsening heart failure in patients who do not already have an ICD, such as patients with diastolic heart failure who do not qualify for ICD implantation.
- the physiological sensing device 120 may instead be an external long-term ECG monitors and the method 200 can be performed for detection of worsening heart function and diagnosis of abnormal sinus node function.
- the calculated HRC index is advantageously used in conjunction with the physiological data measured by the physiological sensing device 120 for risk assessment of impending clinical events, particularly in cardiac diseases.
- FIG. 5 is a block diagram illustrating a technical architecture of a computer 500 in accordance with embodiments of the present disclosure.
- the computer 500 includes a processor / central processing unit (CPU) 502, memory devices 504, a database 506, a data communication module 508, and a physiological data module 510.
- CPU central processing unit
- the processor 502 executes instructions, codes, computer programs, and/or scripts which it accesses from the memory devices 504.
- the processor 502 includes suitable logic, circuitry, and/or interfaces to execute such operations or steps.
- Some non- limiting examples of the processor 502 include an application-specific integrated circuit (ASIC) processor, a reduced instruction set computing (RISC) processor, a complex instruction set computing (CISC) processor, a field-programmable gate array (FPGA), and the like. While only one processor 502 is shown, multiple processors may be present. Thus, while instructions may be discussed as executed by a processor 502, the instructions may be executed simultaneously, serially, or otherwise executed by one or multiple processors 502 (e.g. in a multi-core configuration).
- ASIC application-specific integrated circuit
- RISC reduced instruction set computing
- CISC complex instruction set computing
- FPGA field-programmable gate array
- the memory devices 504 may comprise storage devices (such as flash memory, disk drives, or memory cards), read-only memory (ROM), and random-access memory (RAM).
- the memory devices 504 store non-transitory instructions operative by the processor 502 to perform various operations or steps of the method 200 according to various embodiments of the present disclosure.
- the memory devices 504 may be referred to as computer-readable storage media and/or non-transitory computer- readable media.
- Non-transitory computer-readable media include all computer- readable media, with the sole exception being a transitory propagating signal per se.
- the database 506 is any computer-operated hardware suitable for storing and/or retrieving data. Some non-limiting examples of the database 506 are the databases 126, 132, and 142.
- the database 506 may include multiple storage units such as hard disks and/or solid-state disks in a Redundant Array of Independent Disks (RAID) configuration.
- the database 506 may include, but is not limited to, a storage area network (SAN) and/or a network attached storage (NAS) system.
- the data communication module 508 is configured for communication with other computers 500.
- the physiological data module 510 is configured to process the physiological data, such as the heart rate data and motion data, and to perform various calculations based on the physiological data to obtain the HRC index in accordance with the method 200.
- Various algorithms may be implemented in the physiological data module 510 for performing said calculations based on the physiological data and to obtain the HRC index.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Physics & Mathematics (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Cardiology (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Primary Health Care (AREA)
- Artificial Intelligence (AREA)
- Signal Processing (AREA)
- Epidemiology (AREA)
- Psychiatry (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
La présente invention concerne d'une manière générale un système (100) et un procédé (200) pour évaluer des risques d'événements cliniques sur la base de la complexité de la fréquence cardiaque d'un sujet (110). Le système (100) comprend un dispositif de détection physiologique (120) et un dispositif électronique (130). Le procédé (200) consiste à recevoir (202), à partir du dispositif de détection physiologique (120), un ensemble de données physiologiques mesuré à partir du sujet (110), l'ensemble de données physiologiques comprenant un ensemble de données de fréquence cardiaque ; à générer (204) une série chronologique de fréquence cardiaque à partir des données de fréquence cardiaque ; à calculer (206) un ensemble de données d'entropie à partir de la série chronologique de fréquence cardiaque ; à calculer (208) un indice HRC à partir des données d'entropie ; à comparer (210) l'indice HRC à des données de référence prédéfinies ; et (212) à générer un signal de risque en réponse à une détermination, à partir de ladite comparaison, du fait que l'indice HRC indique des risques accrus d'événements cliniques au sujet (110).
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
SG10201801227X | 2018-02-13 | ||
SG10201801227X | 2018-02-13 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2019160504A1 true WO2019160504A1 (fr) | 2019-08-22 |
Family
ID=67621039
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/SG2019/050083 WO2019160504A1 (fr) | 2018-02-13 | 2019-02-13 | Système et procédé pour évaluer un risque d'événement clinique sur la base de la complexité de la fréquence cardiaque |
Country Status (1)
Country | Link |
---|---|
WO (1) | WO2019160504A1 (fr) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110782989A (zh) * | 2019-09-18 | 2020-02-11 | 平安科技(深圳)有限公司 | 数据分析方法、装置、设备及计算机可读存储介质 |
CN114983371A (zh) * | 2022-05-25 | 2022-09-02 | 佳木斯大学 | 基于人工智能的心内科用心率不齐测试系统及方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012103273A2 (fr) * | 2011-01-27 | 2012-08-02 | Valencell, Inc. | Appareil et procédés permettant de surveiller des données physiologiques pendant l'interférence d'un environnement |
WO2013016290A2 (fr) * | 2011-07-22 | 2013-01-31 | President And Fellows Of Harvard College | Procédés et systèmes de detection de la dépression basés sur la complexité |
WO2015036289A1 (fr) * | 2013-09-16 | 2015-03-19 | Koninklijke Philips N.V. | Système et procédé d'estimation de la santé cardiovasculaire d'une personne |
US20150223699A1 (en) * | 2014-02-11 | 2015-08-13 | Seoul National University Bundang Hospital | System and method for assessing treatment effects on obstructive sleep apnea |
WO2016201130A1 (fr) * | 2015-06-09 | 2016-12-15 | University Of Connecticut | Procédé et appareil de surveillance de fréquence cardiaque utilisant un capteur d'électrocardiogramme |
US20170036065A1 (en) * | 2015-08-04 | 2017-02-09 | The General Hospital Corporation | System and method for assessment of cardiovascular fitness |
-
2019
- 2019-02-13 WO PCT/SG2019/050083 patent/WO2019160504A1/fr active Application Filing
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012103273A2 (fr) * | 2011-01-27 | 2012-08-02 | Valencell, Inc. | Appareil et procédés permettant de surveiller des données physiologiques pendant l'interférence d'un environnement |
WO2013016290A2 (fr) * | 2011-07-22 | 2013-01-31 | President And Fellows Of Harvard College | Procédés et systèmes de detection de la dépression basés sur la complexité |
WO2015036289A1 (fr) * | 2013-09-16 | 2015-03-19 | Koninklijke Philips N.V. | Système et procédé d'estimation de la santé cardiovasculaire d'une personne |
US20150223699A1 (en) * | 2014-02-11 | 2015-08-13 | Seoul National University Bundang Hospital | System and method for assessing treatment effects on obstructive sleep apnea |
WO2016201130A1 (fr) * | 2015-06-09 | 2016-12-15 | University Of Connecticut | Procédé et appareil de surveillance de fréquence cardiaque utilisant un capteur d'électrocardiogramme |
US20170036065A1 (en) * | 2015-08-04 | 2017-02-09 | The General Hospital Corporation | System and method for assessment of cardiovascular fitness |
Non-Patent Citations (2)
Title |
---|
LIU N. T. ET AL.: "Utility of vital signs, heart rate variability and complexity, and machine learning for identifying the need for lifesaving interventions in trauma patients", SHOCK., vol. 42, no. 2, August 2014 (2014-08-01), pages 108 - 114, XP055633468 * |
MEJADDAM A. Y. ET AL.: "Real-time heart rate entropy predicts the need for lifesaving interventions in trauma activation patients", J TRAUMA ACUTE CARE SURG., vol. 75, no. 4, October 2013 (2013-10-01), pages 607 - 612 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110782989A (zh) * | 2019-09-18 | 2020-02-11 | 平安科技(深圳)有限公司 | 数据分析方法、装置、设备及计算机可读存储介质 |
CN114983371A (zh) * | 2022-05-25 | 2022-09-02 | 佳木斯大学 | 基于人工智能的心内科用心率不齐测试系统及方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6684302B2 (ja) | 心肺機能不全に対する易罹患性(susceptibility)を判断する装置を制御する方法、コンピュータ読取可能記憶媒体 | |
Ho et al. | Predicting survival in heart failure case and control subjects by use of fully automated methods for deriving nonlinear and conventional indices of heart rate dynamics | |
US10531838B2 (en) | Apparatus, system, method and computer program for assessing the risk of an exacerbation and/or hospitalization | |
JP5723024B2 (ja) | 逐次分類器を用いる心不全検出 | |
CN105228508B (zh) | 一种测定用于分类的危险评分的系统 | |
Lake et al. | Accurate estimation of entropy in very short physiological time series: the problem of atrial fibrillation detection in implanted ventricular devices | |
JP6857612B2 (ja) | 心血管劣化の警告スコア | |
JP2020536629A (ja) | Ecgベースの心駆出率のスクリーニング | |
JP7293050B2 (ja) | 軽度認知障害判定システム | |
US11185269B1 (en) | Context-aware post traumatic stress disorder monitoring and intervention | |
US9301695B2 (en) | Stress model based on RR integral average | |
EP4193371A1 (fr) | Systèmes et procédés de génération d'évaluations de risques pour la santé | |
WO2019160504A1 (fr) | Système et procédé pour évaluer un risque d'événement clinique sur la base de la complexité de la fréquence cardiaque | |
Alves et al. | Linear and complex measures of heart rate variability during exposure to traffic noise in healthy women | |
US20150133795A1 (en) | System and method for risk stratification based on dynamic nonlinear analysis and comparison of cardiac repolarization with other physiological signals | |
TWI688371B (zh) | 心房顫動信號型態擷取及輔助診斷智能裝置 | |
US12076120B2 (en) | Systems, methods and media for estimating compensatory reserve and predicting hemodynamic decompensation using physiological data | |
US20190380661A1 (en) | Diagnostic Method And System | |
WO2023133166A2 (fr) | Système et méthode de détermination d'un état de santé cardiaque | |
US20210052217A1 (en) | Methods and systems for patient baseline estimation | |
US20230172564A1 (en) | Technologies for determining a risk of developing atrial fibrillation | |
Kane et al. | Accuracy of pulse interval timing in ambulatory blood pressure measurement | |
US20240206823A1 (en) | Method, apparatus and program for calculating disease risk based on heart rate measurement | |
US20240245345A1 (en) | Non-invasive detection and differentiation of sepsis | |
US20230326601A1 (en) | Detection of cardiac conditions from reduced lead set ecg |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 19754266 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 19754266 Country of ref document: EP Kind code of ref document: A1 |