- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Computational and Mathematical Methods in Medicine
Volume 2013 (2013), Article ID 324325, 7 pages
Multiscale Cross-Approximate Entropy Analysis as a Measure of Complexity among the Aged and Diabetic
1Department of Electrical Engineering, National Dong Hwa University, No. 1, Section 2, Da Hsueh Road, Shoufeng, Hualien 97401, Taiwan
2Research Center for Adaptive Data Analysis & Center for Dynamical Biomarkers and Translational Medicine, National Central University, Chungli 32001, Taiwan
3Department of Neurology, Buddhist Tzu Chi General Hospital and Buddhist Tzu Chi University, Hualien 97002, Taiwan
4Department of Internal Medicine, Hualien Hospital, Health Executive Yuan, Hualien 97061, Taiwan
Received 22 March 2013; Revised 27 May 2013; Accepted 1 June 2013
Academic Editor: Shengyong Chen
Copyright © 2013 Hsien-Tsai Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Complex fluctuations within physiological signals can be used to evaluate the health of the human body. This study recruited four groups of subjects: young healthy subjects (Group 1, ), healthy upper middle-aged subjects (Group 2, ), subjects with well-controlled type 2 diabetes (Group 3, ), and subjects with poorly controlled type 2 diabetes (Group 4, ). Data acquisition for each participant lasted 30 minutes. We obtained data related to consecutive time series with R-R interval (RRI) and pulse transit time (PTT). Using multiscale cross-approximate entropy (MCE), we quantified the complexity between the two series and thereby differentiated the influence of age and diabetes on the complexity of physiological signals. This study used MCE in the quantification of complexity between RRI and PTT time series. We observed changes in the influences of age and disease on the coupling effects between the heart and blood vessels in the cardiovascular system, which reduced the complexity between RRI and PTT series.
Multiple temporal and spatial scales produce complex fluctuations within the output signals of physiological systems . In recent studies on translational medicine [1–5], researchers have found that implicit information within the complex fluctuations of physiological signals can be used to evaluate health conditions.
Many recent studies [2, 3] have employed nonlinear dynamical analysis to quantify the complexity of physiological signals in the cardiovascular system. Costa et al.  were the first to propose multiscale entropy (MSE) as an approach to analyze the R-R interval (RRI) series of healthy individuals and discovered that the RRI series of young individuals were more complex than that of elderly people. Wu et al.  adopted the same method in an examination of pulse wave velocity (PWV) and found that the complexity of these series decreased with aging and/or the progression of diabetes. In addition to time and space, “coupling behavior” in the physiological system also affects the complexity of individual physiological signals, such as RRI or PWV . Drinnan et al.  indicated that pulse transit time (PTT) is influenced by RRI and other cardiovascular variables and used cross-correlation functions to quantify the phase relationship between the two time series signals in the cardiovascular system. They established that there was a strong correlation between PTT and RRI variations in healthy subjects. However, Pincus  claimed that cross-approximate entropy (Co_ApEn) is more effective than cross-correlation functions in the evaluation of complexity between the two series.
Despite the fact that Co_ApEn has been widely applied to evaluate the complexity between two time series [9–12], single-scale entropy values are not necessarily able to identify the dynamic complexity of physiological signals. Therefore, this study was an attempt to use a multiscale Co_ApEn (MCE)  to quantify the complexity between the synchronous time series of cardiac functions and the degree of atherosclerosis. We assumed that complexity would exist in RRI and PTT series of the cardiovascular system due to the mutual interaction between the heart and blood vessels. Moreover, we assumed that complexity reduces with aging and the influence of disease. We used MCE to develop an index for the quantification of complexity between the two time series capable of distinguishing between healthy individuals and those with diabetes.
2.1. Study Design
This study evaluated the influences of age and diabetes on RRI and PTT. Considering that RRI and PTT are nonlinear, cardiovascular variables, we tested the applicability of MCE in the study subjects and investigated whether this dynamic parameter could provide further information related to the clinical control of diabetes.
2.2. Subject Populations and Experiment Procedure
Between July 2009 and March 2012, four groups of subjects were recruited for this study: young healthy subjects (Group 1, age range: 18–40, ), healthy upper middle-aged subjects (Group 2, age range: 41–80, ), subjects with well-controlled type 2 diabetes (Group 3, age range: 41–80, , 6.5% glycosylated hemoglobin (HbA1c) < 8%), and subjects with poorly controlled type 2 diabetes (Group 4, age range: 41–80, , HbA1c 8%) . The other 22 subjects were excluded due to incomplete or unstable waveform data acquisition. All diabetic subjects were recruited from the Hualien Hospital Diabetic Outpatient Clinic; healthy controls were recruited from a health examination program at the same hospital. None of the healthy subjects had personal or family history of cardiovascular disease. Type 2 diabetes was diagnosed as either fasting sugar higher than 126 mg/dL or HbA1c 6.5%. All diabetic subjects had been receiving regular treatment and follow-up care in the clinic for more than two years. Regarding the use of medications, there was no significant difference in the type (i.e., antihypertensive, lipid-lowering, and hypoglycemic medications), dosage, and frequency among the well-controlled and poorly controlled diabetic subjects. This study was approved by the Institutional Review Board (IRB) of Hualien Hospital and National Dong Hwa University. All subjects refrained from caffeinated beverages and theophylline-containing medications for 8 hours prior to each hospital visit. Each subject gave informed consent, completed questionnaires on demographic data and medical history, and underwent blood sampling prior to data acquisition. Blood pressure was obtained once from the left arm of supine subjects using an automated oscillometric device (BP3AG1, Microlife, Taiwan) with a cuff of appropriate size, followed by the acquisition of waveform data from the second toe using a six-channel ECG-PWV [14, 15] as previously described.
2.3. Data Collection and Calculation of RRI and PTT Series
All subjects were permitted to rest in a supine position in a quiet, temperature-controlled room at 25 ± 1°C for 5 minutes prior to subsequent 30-minute measurements. Again, a good reproducibility of six-channel ECG-PWV system [14, 15] was used for waveform measurement from the second toe. Infrared sensors were simultaneously applied to points of reference for the acquisition of data. Electrocardiogram (ECG) measurements were obtained using the conventional method. After being processed through an analog-to-digital converter (USB-6009 DAQ, National Instruments, Austin, TX USA) at a sampling frequency of 500 Hz, the digitized signals were stored on a computer. Because of its conspicuousness, the R wave in Lead II was selected as a reference point: the time interval between the R-wave peak of the th cardiac cycle to the footpoint of the toe pulse from the left foot was defined as PTT(); the time difference between the two continues peak of ECG R wave was defined as RRI(), as shown as Figure 1.
Using ECG and photoplethysmography (PPG), we obtained the RRI series = and PTT series = from each subject. All series were retrieved from 1000 consecutive, stable ECG tracings and PPG toe pulse signals synchronous with the cardiac cycle .
Due to a trend within physiological signals [6, 16], nonzero means may be included; therefore, we used empirical mode decomposition (EMD)  to deconstruct the and series, thereby eliminating the trend from the original series. We then normalized the and series, as shown in (1). In these equations, and represent the standard deviations of series and , respectively. Complexity analysis was performed on the normalized results, and . Consider
2.4. Multiscale Cross-Approximate Entropy (MCE) Using Normalized RRI and PTT Series Together
Previous studies [1–3, 18] have employed MSE to overcome comparison difficulties at a scale factor of 1, when physiological complexity is reduced due to age or disease. However, other research  has indicated a strong relationship between variations in PTT series and RRI series; therefore, we used MCE to investigate the interactions between PTT and RRI.
2.4.1. Coarse-Grained Process and Cross-Approximate Entropy (Co_ApEn)
MSE involves the use of a scale factor , which is selected according to a 1D series of consecutive cycles. This factor enables the application of a coarse-graining process capable of deriving a new series prior to the calculation of entropy in each new individual series [1–3, 18]. Using this approach, we performed coarse-graining on the normalized 1D consecutive cycles of the series based on scale factor , thereby obtaining the series and as shown in (2). We then calculated entropy as follows:
Previous studies [19, 20] have used Co_ApEn, an improved analysis method of approximate entropy, to analyze two synchronous physiological time series, define their relationship, and calculate the complexity within that relationship [8, 21]. This method utilizes the dynamic changes between the two series to evaluate the physiological system. Similarities between changes in the two series can be used to observe the regulatory mechanisms in the physiological system. However, many studies [8, 19–21] presented their results at a scale factor of 1. To obtain a deeper understanding of the complexity of the physiological system, we utilized coarse-grained and series to calculate the Co_ApEn at each scale, using (7). We refer to this approach as multiscale cross-approximate entropy (MCE). The details of the algorithm are as follows .(1)For given , for two sets of -vectors, (2)Define the distance between the vectors , as the maximum absolute difference between their corresponding elements, as follows: (3)With the given , find the value of (where = 1 to − + 1) that is smaller than or equal to and the ratio of this number to the total number of -vectors ( − + 1). That is, let the number of satisfying the requirement , then measures the frequency of the -point pattern being similar (within a tolerance of ) to the -point pattern formed by .(4)Average the logarithm of over to obtain , as follows: (5)Increase by 1, and repeat steps 14 to obtain , .(6)Finally, take and for -point data, the estimate is where represents the chosen vector dimension, represents a tolerance range, and is the data length. To ensure efficiency and accuracy of calculation, the parameters of this study were set at = 3, = 0.15, and = 1000.
2.4.2. RRI and PTT-Based Multiscale Cross-Approximate Entropy Index (MCEI) for Small and Large Scales
The values of were obtained from a range of scale factors between 1 and 20 using the MCE data analysis method. The values of between scale factors 1 and 5 were defined as small scale; those between scale factors 6 and 20 were defined as large scale . The sum of MCE between scale factors 1 and 5 was MCEISS in (8), while the sum of MCE between scale factors 6 and 20 was MCEILS in (9). Defining and calculating these two indices of multiscale cross-approximate entropy enables the assessment and quantification of complexity in RRI and PTT between different scale factors. Consider
2.5. Multiscale Entropy Index (MEI) Using RRI or PTT Only
Sample entropy () was used to quantify the complexity of RRI or PTT series in twenty scales. The values of between scale factors 1 and 5 were defined as small scale, whereas those between scale factors 6 and 20 were defined as large scale. The sum of MSE in small scale was defined as MEISS, while the sum of MSE in large scale was MEILS .
2.6. Statistical Analysis
Average values were expressed as mean ± SD. Significant differences in anthropometric, hemodynamic, and computational parameters (i.e., RRI, PTT, MCEISS, and MCEILS) between different groups were determined using an independent sample -test. Statistical Package for the Social Science (SPSS, version 14.0 for Windows) was used for all statistical analysis. A value less than 0.05 was considered statistically significant.
3.1. Comparison of Basic Demographic and Cardiovascular Parameters in Different Groups
Table 1 presents the basic demographic parameters of Group 1 and Group 2, showing no significant difference in major demographic parameters except for age, HbA1c levels, and body height. Significant differences were observed in body mass index (BMI), waist circumference, systolic blood pressure (SBP), pulse pressure (PP), HbA1c levels, and fasting blood sugar level between Group 2 and Group 3 (Group 3 > Group 2). In addition, significant differences were also observed in HbA1c levels, triglycerides, and fasting blood sugar level between Group 3 and Group 4.
3.2. MCEILS as Parameters Indicative of Age and Diabetic Control
Figure 3 summarizes the results of the MCE analysis for the values of RRI and PTT time series over 1000 identical cardiac cycles obtained from the four groups of participants. At a scale factor of 1 (), the magnitudes of ranked as follows: Group 1/Group 3/Group 4/Group 2. The value of began dropping in all groups at a scale factor of 2 ().
Beginning at a scale factor of 3 (), the reduction in in Group 1 slowed. However, in the other groups, the values continued decreasing rapidly. Beginning at a scale factor of 5 (), the of Group 2 achieved stability with only minor fluctuations. The decline in in Group 4 remained greater than that in Group 3. When plotted against large scale factors (i.e., 6–20), the magnitudes of ranked as follows: Group 1, Group 2, Group 3, and Group 4.
MCEISS only presented a significant difference between Groups 1 and 2 (10.18 ± 0.52 versus 9.42 ± 0.70, ). The differences among Groups 2, 3, and 4 did not reach statistical significance. In comparison, MCEILS presented significant differences among all four of the groups (Group 1 versus Group 2: 28.30 ± 1.26 versus 25.96 ± 1.99, ; Group 2 versus Group 3: 25.96 ± 1.99 versus 23.14 ± 1.85, ; Group 3 versus Group 4: 23.14 ± 1.85 versus 20.13 ± 1.73, ) (Table 1).
Since Pincus and Singer’s study , Co_ApEn has generally been used to reveal similarities between two synchronous, consecutive variables within a single network. This approach has also been used to research the complexity of physiological signals [12, 19]; however, the influence of multiple temporal and spatial scales creates complexity. Thus, this study employed multiscale Co_ApEn (MCE) to evaluate the complexity between the cardiac function-related parameter, RRI, and the atherosclerosis-related parameter, PTT, in the cardiovascular systems of various subject groups.
Previous studies [1, 2, 18] have also indicated that physiological signals are generally nonlinear and exist in nonstationary states. The use of MSE to quantify complexity within the times series of a single type of physiological signal (i.e., RRI or PWV) demonstrated that the complexity of physiological signals decreases with aging  or with the influence of diabetes . In this study, although we used MSE to quantify complexity of RRI or PTT series, there were no significant differences in , , , and between well-controlled and poor-controlled diabetic subjects. Therefore, the influence of the degree of glycemic control on complexity of physiological signals might not be evaluated efficiently according to the use of MSE when analyzing single time series (i.e., RRI or PTT).
Drinnan et al.’s study  stated that cardiovascular variables such as RRI and PTT are regulated by complex physiological systems and that a strong relationship exists between variations in PTT and those in RRI. We therefore employed the Co_ApEn integrated with preprocessing coarse-graining to calculate MCEI values as well as the complexity between the synchronous time series RRI and PTT. Figure 3 shows that at small-scale factors (from 1 to 5), it is difficult to determine the influence of age, diabetes, or glycemic control based on the complexity between the time series RRI and PTT using . Similarly, MCEISS indicates only that aging reduces the complexity between the two time series. This finding is similar to that of previous studies . As the scale factor increased (from 6 to 20), began revealing significant differences between the four study groups (Figure 3). Table 1 shows that the MCEILS values of the young healthy subjects were the highest, whereas subjects with poorly controlled type 2 diabetes were the lowest. This may be due to the fact that the coupling effect between the heart and the blood vessels in the cardiovascular system varies according to age and the influence of disease [24, 25]. In other words, the complexity between the time series RRI and PTT decreases due to age and disease.
Although the MCEILS can be used to quantify the complexity of RRI and PTT and have been shown to effectively identify significant difference among study groups, limitations still exist. First, a lengthy process of data acquisition and considerable calculation and off-line processing is needed. MCE analysis involves a 30-minute measurement, as opposed to the relatively shorter duration measurement of only RRI and PTT, making the process tiring for participants. The nature of analysis postmeasurement further prevented subjects from receiving their MCEI test results immediately. Second, the medications that the diabetic patients used such as hypoglycemic, antihyperlipidemic, and antihypertensive drugs may also affect autonomic nervous activity. These effects, however, were difficult to assess. The potential effect of medications, therefore, was not considered in the statistical analysis of this study.
This study integrates cross-approximate entropy with multiple scales to analyze the complexity between two synchronous physiological signals (RRI and PTT) in the cardiovascular system. According to our results, MCEILS clearly reveals a reduction in the complexity of two physiological signals caused by aging and diabetes.
M.-T. Lo and A.-B. Liu equally contributed in this study compared with the corresponding author.
Conflict of Interests
The authors declare no conflict of interests.
The authors would like to thank the volunteers involved in this study for allowing them to collect and analyze their data. The authors are grateful for the support of Texas Instruments, Taiwan, in sponsoring the low power instrumentation amplifiers and ADC tools. The authors would also like to thank Miss Crystal J. McRae who is a native English speaker, to go over the whole paper. This research was partly supported by National Science Council under Grants NSC 100-2221-E-259-030-MY2 and NSC 101-2221-E-259-012 and National Dong Hwa University on campus interdisciplinary integration Projects no. 101T924-3 and 102T931-3. M.-T. Lo was supported by NSC (Taiwan, ROC), Grant no. 100-2221-E-008-008-MY2, joint foundation of CGH and NCU, Grant no. CNJRF-101CGH-NCU-A4, VGHUST102-G1-2-3 and NSC support for the Center for Dynamical Biomarkers and Translational Medicine, National Central University, Taiwan (NSC 101-2911-I-008-001).
- M. Costa, A. L. Goldberger, and C. K. Peng, “Multiscale entropy analysis of biological signals,” Physical Review E, vol. 71, no. 2, part 1, 2005.
- M. Costa, A. L. Goldberger, and C. K. Peng, “Multiscale entropy to distinguish physiologic and synthetic RR time series,” Computing in Cardiology, vol. 29, pp. 137–140, 2002.
- H.-T. Wu, P.-C. Hsu, C.-F. Lin et al., “Multiscale entropy analysis of pulse wave velocity for assessing atherosclerosis in the aged and diabetic,” IEEE Transactions on Biomedical Engineering, vol. 58, no. 10, pp. 2978–2981, 2011.
- R. T. Vieira, N. Brunet, S. C. Costa, S. Correia, B. G. A. Neto, and J. M. Fechine, “Combining entropy measures and cepstral analysis for pathological voices assessment,” Journal of Medical and Biological Engineering, vol. 32, no. 6, pp. 429–435, 2012.
- J. Y. Lan, M. F. Abbod, R. G. Yeh, S. Z. Fan, and J. S. Shieh, “Review: intelligent modeling and control in anesthesia,” Journal of Medical and Biological Engineering, vol. 32, no. 5, pp. 293–307, 2012.
- C.-K. Peng, M. Costa, and A. L. Goldberger, “Adaptive data analysis of complex fluctuations in physiologic time series,” Advances in Adaptive Data Analysis, vol. 1, no. 1, pp. 61–70, 2009.
- M. J. Drinnan, J. Allen, and A. Murray, “Relation between heart rate and pulse transit time during paced respiration,” Physiological Measurement, vol. 22, no. 3, pp. 425–432, 2001.
- S. M. Pincus, “Approximate entropy in cardiology,” Herzschrittmachertherapie und Elektrophysiologie, vol. 11, no. 3, pp. 139–150, 2000.
- S. M. Pincus, T. Mulligan, A. Iranmanesh, S. Gheorghiu, M. Godschalk, and J. D. Veldhuis, “Older males secrete luteinizing hormone and testosterone more irregularly, and jointly more asynchronously, than younger males,” Proceedings of the National Academy of Sciences of the United States of America, vol. 93, no. 24, pp. 14100–14105, 1996.
- F. Roelfsema, S. M. Pincus, and J. D. Veldhuis, “Patients with Cushing's disease secrete adrenocorticotropin and cortisol jointly more asynchronously than healthy subjects,” Journal of Clinical Endocrinology and Metabolism, vol. 83, no. 2, pp. 688–692, 1998.
- J. D. Veldhuis, A. Iranmanesh, T. Mulligan, and S. M. Pincus, “Disruption of the young-adult synchrony between luteinizing hormone release and oscillations in follicle-stimulating hormone, prolactin, and nocturnal penile tumescence (NPT) in healthy older men,” Journal of Clinical Endocrinology and Metabolism, vol. 84, no. 10, pp. 3498–3505, 1999.
- D.-Y. Wu, G. Cai, Y. Yuan et al., “Application of nonlinear dynamics analysis in assessing unconsciousness: a preliminary study,” Clinical Neurophysiology, vol. 122, no. 3, pp. 490–498, 2011.
- M. U. Ahmed and D. P. Mandic, “Multivariate multiscale entropy: a tool for complexity analysis of multichannel data,” Physical Review E, vol. 84, no. 6, Article ID 061918, 2011.
- A.-B. Liu, P.-C. Hsu, Z.-L. Chen, and H.-T. Wu, “Measuring pulse wave velocity using ECG and photoplethysmography,” Journal of Medical Systems, vol. 35, no. 5, pp. 771–777, 2011.
- H. T. Wu, P. C. Hsu, A. B. Liu, Z. L. Chen, R. M. Huang, C. P. Chen, et al., “Six-channel ECG-based pulse wave velocity for assessing whole-body arterial stiffness,” Blood Press, vol. 21, no. 3, pp. 167–176, 2012.
- Z. Wu, N. E. Huang, S. R. Long, and C.-K. Peng, “On the trend, detrending, and variability of nonlinear and nonstationary time series,” Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 38, pp. 14889–14894, 2007.
- N. E. Huang, Z. Shen, S. R. Long et al., “The empirical mode decomposition and the Hubert spectrum for nonlinear and non-stationary time series analysis,” Proceedings of the Royal Society A, vol. 454, no. 1971, pp. 903–995, 1998.
- M. D. Costa, C.-K. Peng, and A. L. Goldberger, “Multiscale analysis of heart rate dynamics: entropy and time irreversibility measures,” Cardiovascular Engineering, vol. 8, no. 2, pp. 88–93, 2008.
- S. Pincus and B. H. Singer, “Randomness and degrees of irregularity,” Proceedings of the National Academy of Sciences of the United States of America, vol. 93, no. 5, pp. 2083–2088, 1996.
- M. Kreuzer, H. Hentschke, B. Antkowiak, C. Schwarz, E. F. Kochs, and G. Schneider, “Cross-approximate entropy of cortical local field potentials quantifies effects of anesthesia—a pilot study in rats,” BMC Neuroscience, vol. 11, article 122, 2010.
- S. M. Pincus, “Irregularity and asynchrony in biologic network signals,” Methods in Enzymology, vol. 321, pp. 149–182, 2000.
- F. Yang, B. Hong, and Q. Tang, “Approximate entropy and its application to biosignal analysis,” in Nonlinear Biomedical Signal Processing: Dynamic Analysis and Modeling, M. Akay, Ed., vol. 2, John Wiley & Sons, Hoboken, NJ, USA, 2000.
- D. Cheng, S.-J. Tsai, C.-J. Hong, and A. C. Yang, “Reduced physiological complexity in robust elderly adults with the APOE ε4 allele,” PLoS ONE, vol. 4, no. 11, Article ID e7733, 2009.
- D. E. Vaillancourt and K. M. Newell, “Changing complexity in human behavior and physiology through aging and disease,” Neurobiology of Aging, vol. 23, no. 1, pp. 1–11, 2002.
- D. T. Kaplan, M. I. Furman, S. M. Pincus, S. M. Ryan, L. A. Lipsitz, and A. L. Goldberger, “Aging and the complexity of cardiovascular dynamics,” Biophysical Journal, vol. 59, no. 4, pp. 945–949, 1991.