Journal of Healthcare Engineering
Volume 2017 (2017), Article ID 7406896, 10 pages
https://doi.org/10.1155/2017/7406896

Research Article

A Real-Time Analysis Method for Pulse Rate Variability Based on Improved Basic Scale Entropy

Yongxin Chou,¹ Ruilei Zhang,¹ Yufeng Feng,² Mingli Lu,¹ Zhenli Lu,1,3 and Benlian Xu¹

¹School of Electrical and Automatic Engineering, Changshu Institute of Technology, Changshu 215500, China
²Changshu No. 1 People’s Hospital, Changshu, China
³State Key Laboratory of Robotics, Shenyang Institute of Automation, Chinese Academy of Sciences, Shenyang 110014, China

Correspondence should be addressed to Yongxin Chou; lutchouyx@163.com

Received 17 November 2016; Revised 18 February 2017; Accepted 7 March 2017; Published 9 May 2017

Academic Editor: Valentina Camomilla

Copyright © 2017 Yongxin Chou 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.

Abstract.

Base scale entropy analysis (BSEA) is a nonlinear method to analyze heart rate variability (HRV) signal. However, the time consumption of BSEA is too long, and it is unknown whether the BSEA is suitable for analyzing pulse rate variability (PRV) signal. Therefore, we proposed a method named sliding window iterative base scale entropy analysis (SWIBSEA) by combining BSEA and sliding window iterative theory. The blood pressure signals of healthy young and old subjects are chosen from the authoritative international database MIT/PhysioNet/Fantasia to generate PRV signals as the experimental data. Then, the BSEA and the SWIBSEA are used to analyze the experimental data; the results show that the SWIBSEA reduces the time consumption and the buffer cache space while it gets the same entropy as BSEA. Meanwhile, the changes of base scale entropy (BSE) for healthy young and old subjects are the same as that of HRV signal. Therefore, the SWIBSEA can be used for deriving some information from long-term and short-term PRV signals in real time, which has the potential for dynamic PRV signal analysis in some portable and wearable medical devices.

1. Introduction

Electrocardiogram (ECG) signal has been used for many diseases to assist in diagnosis in a clinic. The subtle changes of heart beat periods are called heart rate variability (HRV). The continuous heart rate or continuous RR wave intervals extracted from ECG signal are denoted as heart rate variability (HRV) signal [1]. An increasing number of studies have shown that HRV is a useful quantitative indicator for assessing the balance between the cardiac sympathetic nervous system and the parasympathetic nervous system and can be engaged in the diagnosis and prevention of some cardiovascular diseases such as sudden cardiac death and arrhythmia [2–4]. Pulse signal or continuous blood pressure signal generated by the systolic and diastolic of heart contains abundant physiological and pathological information of the cardiovascular system [5, 6]. The subtle change of vessel pulse periods is denoted as pulse rate variability (PRV). The continuous pulse rate or continuous PP wave intervals extracted from pulse signal or continuous blood pressure signal are defined as PRV signal [7]. Because a heartbeat produces a vessel pulse, many studies show that PRV is a substitute for HRV to present the physiological and pathological changes of the cardiovascular system when the subjects are sleeping or testing, as well as in some nonstationary states [8–10]. In addition, due to the wide distribution of human vessels, the acquisition of a pulse signal is easier than that of an ECG signal. Therefore, the pulse signal is employed in many wearable and portable medical devices such as smart watches, wristbands, and smart glasses but not ECG signal [11, 12], and PRV signal has more practical values than HRV signal.

Because PRV signal has similar characteristics with HRV signal, the analysis methods of HRV signal are often employed to analyze PRV signal. These methods are divided into time domain methods, frequency domain methods, time-frequency domain methods, and nonlinear methods [13]. HRV signal and PRV signal generated by heartbeat are neither stochastic nor periodic; they are the results of many independent factors and have nonlinear properties. Thus, the nonlinear methods are more useful for analyzing HRV signal and PRV signal, and there are many nonlinear methods such as recurrence quantification analysis, detrended fluctuation analysis, the Lyapunov exponent, and information entropy analysis [14, 15]. Among them, the information entropy analysis is an effective tool to present the complexity of the nonlinear signal. The sample entropy (SampEn), the approximate entropy (ApEn), the sign series entropy analysis (SSEA), the base scale entropy analysis (BSEA), and so on are been used for analyzing HRV signal [13, 16–18]. However, because of the long time consumption of these methods, they are not suitable for the PRV signal in real time. The BSEA, proposed by Li and Ning, can effectively detect the complexity dissimilarity of short-term HRV signal (about 5 minutes) in different physiological or pathological states [17], while it is unknown whether the BSEA is suitable for analyzing pulse rate variability (PRV) signal, so far. In addition, the 5 minutes of HRV signal analysis is too long for some acute cardiovascular disease (ACVD), and its time consumption still needs to be improved.

Therefore, this study proposed an improved basic scale entropy on the basis of BSEA with the theory of sliding window iterative; we denote it as sliding window iterative basic scale entropy analysis (SWIBSEA). The BSEA and SWIBSEA are engaged in analyzing the measured PRV signals, and by the results of the experiments, the accuracy and time consumption are compared between BSEA and SWIBSEA. In addition, the structure of this paper is as follows: in Section 2, the theories of BSEA and SWIBSEA are presented and then the experimental data are introduced. The results are shown in Section 3. Then, the results are discussed in Section 4. The conclusion is given in the last section.

2. Methods and Materials

2.1. Basic Scale Entropy Analysis

The process of BSEA is as follows [16]: (1) a series of vectors are constructed from PRV signal, and for each vector, we compute their basic scale (BS). (2) The vectors are symbolized and classified according to BS, each of these categories is a heart or pulse beat mode. (3) Computing the probability of each beat mode, and getting entropy of their probabilities, the entropy is denoted as BSE.

For a PRV signal with the length of , the consecutive data points are used to construct a vector:

Thus, we will get vectors which are denoted as temporal sequence vectors (TSVs). is the length of TSV; the larger the value of , the more complex of the beat mode that TSV expresses. For each TSV, the BS is defined by the root mean square (RMS) of the difference for two adjacent data points:

where is the BS of the ith TSV.

Then, the BS is multiplied by a constant α; the result is as the standard for the vector symbolization. The are symbolized, and the results are named symbol sequence vectors (SSVs) and denoted as . The symbolization process is as follows:

where is the mean of the ith TSV. is the (i + j)th data points of or is the (j + 1)th datum in the ith TSV. The symbols 0, 1, 2, and 3 are the labels of different scopes for PRV amplitude and are employed for probability calculation; their values are of no practical significance. is the (j + 1)th datum in the ith SSV. is used to control the value of BS and to adjust the division range of PRV amplitude, the way to choose the value of is as [19].

After getting the SSVs, we compute the probability of each vector. There are 4 symbols, 0, 1, 2, and 3, to express the vector, so we can get 4m kinds of different SSVs, denoted by . Each SSV is a heart or pulse beat mode. Then, we compute the probability of each beat mode in SSVs:

where , and # is the number of . The beat state with probability 0 is denoted as “disabled mode.”

Therefore, we define BSE as

The BSE can be used to describe the change of heartbeat mode. Obviously, . When there is only one pulse mode, . When there are 4m pulse modes, and each mode has equal probability, is the maximum. The larger the entropy value, the more complicated the heartbeat mode, whereas the smaller the entropy value, the simpler the heartbeat mode.

2.2. Sliding Window Iterative Basic Scale Entropy Analysis

We improve the BSEA with the theory of sliding window iterative, and define the improved method as sliding window iterative base scale entropy analysis. The process is shown in Figure 1.

Figure 1: The process of SWIBSEA.

In Figure 1, a data buffer with the length of is set to store PRV data points which are extracted from dynamic pulse signal. The process of SWIBSEA contains data updating and iterative. The PRV datum will be updated by the sliding window theory, and the BSE will be calculated with the iterative theory.

During data updating, we set 1 byte in buffer to store the latest PRV datum and denote it as . Then, we delete the oldest PRV datum , and the data in higher addresses move to lower addresses, . If we image the PRV signal as fixed, the buffer is like sliding forward in the data. This is the process of sliding window.

During iterative, the data in and is only changed in data updating, so it is only needed to compute the entropy of SSVs corresponding to and , and then we will get the new BSE of all PRV data in buffer.

2.2.1. Construct the Temporal Sequence Vectors

According to the theory of BSEA, we need to reconstruct the PRV data with the length of to TSVs. In the process of updating data, only the data in and are changed; thus, for speeding up the calculation, we only construct the TSVs corresponding to them.

where is the TSV of is the TSV of .

2.2.2. The Symbolization of Temporal Sequence Vectors and

For computation, by (3),

where is the mean of . is the base scale of , by (2):

Similarly, is symbolized and denoted as , .

2.2.3. Encode the Symbol Sequence Vector

During the data updating, only the numbers of and are changed; thus we only need to update their numbers. Moreover, for speeding up the calculation, we encode and to generate their storage addresses. If we denote the addresses as and (as is shown in Figure 1), then

and

2.2.4. Entropy Calculation with Iterative Theory

We update the PRV data in buffer by sliding window theory; when getting a new PRV datum, the oldest PRV datum should be deleted. Thus, when computing entropy, we should just subtract the entropy of and add the entropy of from the entropy before data updating. Meanwhile, the intermediate variables are also need update by iterative after each data updating.

As shown in Figure 1, we denote the numbers of and as and , the probabilities of and as and , and the entropies of and as and before and after updating PRV data, respectively. After data updating, should subtract 1 and should add 1. Then, . The initial value of . In the iterative process of computing entropy, the antilogarithm of logarithm must be over 0; according to the changes of , there are four kinds of iterative methods to computing BSE as follows:

(1) , which means the beat modes that and expressed all exist. Thus, during computing BSE, the antilogarithm of logarithm will be over 0. By (5),

where .

Because only the modes of and are changed, thus, . With (12) and (13),

Simplify (14):

Note that when the beat mode of is the same as that of before and after data updating, the number of the beat mode that expressed should subtract 1: , and the number of the beat mode that expressed should add 1: . Then, (14) is

That means the entropy is not changed when data updating. Because the beat modes of and are the same, the total number of beat modes remains unchanged.

(2) , which means the beat mode disappeared after data updating, and . By (14),

(3) , which means the beat mode of appeared the first time after data updating, and . By (14),

(4) , which means the beat mode of disappeared, the beat mode of appeared the first time after updating data, and the total number of beat modes is unchanged. Then, . By (14),

By (15), (16), (17), (18), and (19), we will obtain the BSE of PRV signal based on sliding window and iterative theory.

2.3. Experimental Data

The experimental data we used are from the international authority of the database: PhysioNet/Fantasia [20]. In this database, there are 40 health subjects which have the same proportion of men and women, 20 of them are the elderly (65–85 years old, data name: f2o01m~f2o20m), the remaining subjects are the young (21–34 years old, data name: f2y01m~f2y20m). The ECG signal, continuous blood pressure signal, and respiration signal are recorded when the subjects under rest and watching the Fantasia movies to keep awake. However, only half of the subjects’ blood pressure signals are recorded (data name: f2y01m-f2y10m, f2o01m-f2o10m). The data sampling frequency is 250 Hz, and the duration is 66 minutes.

The experimental data we used are the continuous blood pressure signals. Compared with the ECG signal, the blood pressure signal is uncalibrated. Therefore, the dynamic difference threshold method is used to calibrate the P waves [21], and the accuracy of the calibration is determined manually. Then, by making a first-order difference for the locations of the calibrated P waves, we will obtain a set of continuous PP intervals, which is PRV signal.

In reality, there are some singularities in PRV signal that have a bad effect on the signal processing results. A pretreatment method is employed to delete the singularities from these short-term PRV signals. The steps are as follows [22]:

Step 1. For the first datum of PRV signal, if

then is a singularity and thus deleted. In (20), are the mean and the standard deviation of a short-term PRV signal, respectively.

Step 2. For the ith datum of PRV signal, if

then is singular and deleted, where is the datum before .

3. Experimental Results

For simulating the process of PRV analysis in microcontroller system, the length of buffer is set to , which is the length of sliding window. In this study, the sliding window theory is used to data updating. When obtaining a new PRV datum, the data in buffer are analyzed with BSEA and SWIBSEA, and the performance of the two methods are compared.

The experimental data are continuous blood pressure signals that are used to generate PRV signals; the results of a young subject and an old subject are selected randomly and shown in Figure 2. With the individual differences in heart rate, the length of two PRV signals are different. For 66 minutes of data, the young one is 2793 points corresponding to the mean pulse rate is 56.1 beat per minute (bpm) and the old one is 4849 points corresponding to the mean pulse rate is 73.8 bpm.

(a)

(b)

Figure 2: The PRV signal of a young subject (a) and an old subject (b).

In this study, the performance of SSEA and SWISSEA are quantitatively evaluated by mean square error (MSE) and program running time.

The MSE is defined as

where is the length of entropy, is the BSE extracted by BSEA, and is the BSE extracted by SWIBSEA.

For the program running time, we test the programs of SSEA and SWISSEA in MATLAB 2016a on a PC with i7-6700HQ CPU (2.60 GHz, 16 GB buffered RAM).

3.1. Comparison of BSEA and SWIBSEA

The BSEA and the SWIBSEA are used to compute the entropy of the PRV signals in Figure 2. The results are shown in Figure 3.

(a)

(b)

Figure 3: The comparison of BSEA and SWIBSEA for a young subject (a) and an old subject (b), when , , and .

When the total number of PRV data points in buffer is less than window width, , for BSEA, the BSE = 0 based on its theory; for SWIBSEA, this is the initial process, and we compute the entropy by iterative, and the value increases with . When , the BSEA updates PRV data by sliding window and calculates BSE, while the SWIBSEA computes BSE by sliding window iteration; their entropies are shown in Figure 2. From the figure, it can be seen that their values are the same and their MSE = 0. The BSE of the young subject is 3.784 ± 0.050 (mean ± std), and the old subject is 4.056 ± 0.053, except for the entropy of the initial stage.

The time consumption of BSEA and SWIBSEA are 0.132 s and 4.769 s for the young subject and 0.192 s and 8.438 s for the old subject, respectively. For the two 66 minutes PRV signals, the time BSEA cost are 36 times and 44 times for SWIBSEA, respectively. Although the blood pressure signals of the old and the young subjects have the same length, but because of the individual differences, their pulse rates are different, and the PRV signal lengths are different. Thus, the time consumptions of the young and the old are difference. While compared with BSEA, SWIBSEA saves a lot of running time in computing BSE and keeps its values unchanged.

3.2. Comparison of BSEA and SWIBSEA under the Different Lengths of SSV

According to the process of BSEA and SWIBSEA, the length of the SSV has great influence on the running time of their program. The longer the length is, the more the heartbeat modes are represented by the SSVs. Therefore, keeping the width of sliding window and α unchanged (here, they are assigned to 300 and 0.5, respectively), we increase the value of from 2 to 10, and the time consumption of these two methods for a young subject and an old subject are as shown in Figures 4 and 5.

(a)

(b)

Figure 4: The time consumption of SWIBSEA (b) and BSEA (a) under the different length of SSVs for a young subject, when and .

(a)

(b)

Figure 5: The time consumption of SWIBSEA (b) and BSEA (a) under the different length of SSV for an old subject, when and .

Figure 4 shows the time consumption of the two method for the 66 minutes PRV signal of a young subject; Figure 4(a) is comparison of two method, and Figure 4(b) is the time consumption of SWIBSEA. When increases from 2 to 10, the time they cost are all increased. The time consumption of SWIBSEA are from 0.182 s to 0.218 s, and that of BSEA are 42, 45, 45, 46, 50, 51, 62, 134, and 426 times for SWIBSEA. The growth rate of BSEA is 85.317 s, while that of SWIBSEA is only 0.036 s.

Figure 5 is the time consumption of an old subject. When increases from 2 to 10, the time consumption of SWIBSEA are from 0.115 s to 0.137 s, and that of BSEA are 40, 42, 42, 44, 46, 49, 60, 126, and 388 times for SWIBSEA. The growth rate of BSEA is 48.524 s, while that of SWIBSE is only 0.022 s. Although the time consumption of the two methods are all increased, but the increase of SWIBSEA is far less than that of BSEA, SWIBSEA will save much more running time of a program.

The variation of will cause some changes of the SSVs in the window, which inevitably causes the change of the BSE value. As is shown in Figure 6, when increases from 2 to 10, the entropies of young subjects and old subjects are increased. There are significant differences between BSEA and SWIBSEA (, two-sample t-test), and the increase of does not affect the difference between the young and the elderly.

Figure 6: The BSE under different ’s, the results are shown with mean ± std, and . “” is the BSE of old subjects; “o” is the BSE of young subjects.

3.3. Comparison of BSEA and SWIBSEA under Different

The width of sliding window is corresponding to the length of buffer in a microcontroller system, and the range of is varied from several minutes PRV data points (short-term PRV signal) to several hours data points (long-term PRV signal), or even to 24 hours PRV data points. The short-term PRV signal is used to derive some changes in the autonomic nervous system within a short time, and the long-term PRV signal is used to reflect the long time and slow changes of the autonomic nervous system. They have potential to apply in a clinic. However, with the increase of , the time consumption will increase inevitably. The results are shown in Figures 7 and 8. The values of and are chosen randomly; here, , . Then, the increases from 100 to 1000 data points with the interval of 100. For the PRV signal of the young subject (in Figure 7), the time SWIBSEA used are from 0.1301 s to 0.1093 s, respectively. The time BSEA used are 15, 30, 42, 56, 67, 80, 93, 104, 112, and 121 times for SWIBSEA. For the PRV signal of the old subject (in Figure 8), the time SWIBSEA used are from 0.199 s to 0.185 s, respectively. The time BSEA used are 17, 33, 49, 63, 79, 93, 110, 124, 138, and 153 times for SWIBSEA. The increase of BSEA from to are 25.188 s, but that of SWIBSEA decreases and the decrease is only 0.014 s. Because the SWIBSEA computes entropy with iterative and only updates the related variables of and , the BSEA needs to update all the variables in sliding window.

(a)

(b)

Figure 7: The time consumption of SWIBSEA (a) and BSEA (b) under the different for a young subject, when .

(a)

(b)

Figure 8: The time consumption of SWIBSEA (a) and BSEA (b) under the different for an old subject, when , .

The change of will have influence on BSE, as is shown in Figure 9. When increases from 100 to 1000 data points, the difference of BSE between the young subject and the old subject are more and more larger.

Figure 9: The BSE under different . The results are shown with mean ± std. and . “” is the BSE of an old subject; “o” is the BSE of a young subject.

4. Discussion

BSEA, as a nonlinear method, has been employed to HRV signal analysis. Li and Ning [17] used BSEA to analyze the short-term HRV (500 data points) extracted from ECG signals of the PhysioNet/fantasia. When , the results of BSEA show that the entropy increases with aging. The entropies are tested by two-sample t-test and . If the BSE of healthy young subjects represent the best physiological state of the human body, the BSE of healthy old subjects deviated from that of the old subjects, which indicates that normal aging can lead to some function degradation of the body’s control system. We use SWIBSEA to analyze the HRV signals in [17] and get the same results as [17]. Similar to the HRV signal, the PRV signal is also generated from heartbeat. When , , and , the BSE of PRV signals are shown in Figure 10, and the entropies in the figure are the mean of BSE which are derived by SWIBSEA. The two-sample t-test of and <0.05. We do the linear fitting of the entropies by aging, as shown on the solid line, the entropies increase with aging, the same as the result of short-term HRV signals analysis.

Figure 10: The SSE results of young subjects and old subjects. “o” is the young subjects, and “” is the old subjects. Solid line is the linear fitting result of BSE changed with age. Chain dotted line is the 95% confidential region.

BSE has been effectively used for HRV signals analysis. In this study, on the basis of BSEA and with the theory of sliding window iterative, we proposed SWIBSEA for improving the computing efficiency of BSEA. Different from BSEA, it is not necessary to process all the data in buffer after date updating with sliding window for SWIBSEA. For example, when the buffer cache is 1024 bytes, . For BSEA, it needs more than 8169 bytes of memory space to store intermediate variables, that is, bytes for storing PRV data, bytes for storing bytes for storing SSVs, 4m = 1024 bytes for storing the number of different SSVs, and 1 byte for storing BSE. Meanwhile, it also needs times shifting operations to update data; times RMS operations to compute BS; 1020 times loops and times comparison operations to construct SSVs; times comparison operations to updating ; 1024 times multiplications for computing ; 1024 times multiplications, 1024 times logarithms, and 1023 times additions for getting BSE. However, For SWIBSEA, it only needs 2061 bytes, that is, 1024 bytes for buffering PRV data, 2 bytes for storing and , bytes for storing and , 4m = 1024 bytes for storing , and 1 byte for storing BSE. Meanwhile, it needs 1023 times shifting operations to update PRV signal, and 2 times RMS operations to compute and , times comparison operations to compute and , at most 9 times multiplications, 3 times logarithms, and 6 times additions to compute BSE (sometimes, it does not need to update the BSE). Since the storing addresses of are generated by encoding, no loops and comparison operations are required when update . Compared with BSEA, the SWIBSEA saves 6108 bytes buffer space. In addition, the sliding window iterative theory is used by the SWIBSEA, thus, its time consumption are reduced, and it can be engaged in PRV signal analysis in real time.

Moreover, the proposed method significantly reduces the time consumption and buffer cache, and can be used for both short-term and long-term PRV signals (in Section 3.3) by adjusting the width of the sliding window. The length of PRV signal only determines the initialization time; after the initialization, the BSE can be calculated by iterative, and the time consumption between the short-term PRV signal and the long-term PRV signal has not significantly increased.

5. Conclusion

In this study, the sliding window iterative theory is used to improve the BSEA, and the SWIBSEA is proposed and employed to analyze the data of healthy young and old subjects from MIT/PhysioNet/Fantasia database. The results show that compared with BSEA, the SWIBSEA reduces the computing time and saves the buffer cache while keeping the BSE unchanged. Meanwhile, by adjusting the width of sliding window, the SWIBSEA can analyze the long-term and short-term PRV signals in real time. The experimental results show that the BSE increases with aging, and normal aging leads to some functions degradation of the control system. Therefore, the SWIBSEA could be employed in some wearable and portable devices for analyzing dynamic PRV signal in real time.

Conflicts of Interest

The authors indicated no potential conflicts of interest.

Acknowledgments

This work was supported by the State Key Laboratory of Robotics of China (Grant 2014-O08), the Natural Science Fundamental Research Program of Higher Education Colleges in Jiangsu Province of China (Grant 14KJB510001), the Special Research Foundation of Young Teachers (Grant XZ1613), the Suzhou Municipal Science and Technology Plan Project (Grants SYG201548, SYG201504), and the National Natural Science Foundation of China (Grant 61273312, Grant 61673075, and Grant 61473283).

References

R. E. Kleiger, P. K. Stein, and J. T. Bigger, “Heart rate variability: measurement and clinical utility,” Annals of Noninvasive Electrocardiology, vol. 10, no. 1, pp. 88–101, 2005.

J. L. Hamilton and L. B. Alloy, “Atypical reactivity of heart rate variability to stress and depression across development: systematic review of the literature and directions for future research,” Clinical Psychology Review, vol. 50, pp. 67–79, 2016.

A. K. Schuster, J. E. Fisher, J. F. Thayer, D. Mauss, and M. N. Jarczok, “Decreased heart rate variability correlates to increased cardiovascular risk,” International Journal of Cardiology, vol. 203, pp. 728–730, 2015.

E. Ebrahimzadeh, M. Pooyan, and A. Bijar, “A novel approach to predict sudden cardiac death (SCD) using nonlinear and time-frequency analyses from HRV signals,” PLoS One, vol. 9, no. 2, article e81896, 2014.

E. Zahedi, V. Sohani, M. A. Ali, K. Chellappan, and G. K. Beng, “Experimental feasibility study of estimation of the normalized central blood pressure waveform from radial photoplethysmogram,” Journal of Healthcare Engineering, vol. 6, no. 1, pp. 121–144, 2015.

S. Hu, V. Azorinperis, and J. Zheng, “Opto-physiological modeling applied to photoplethysmographic cardiovascular assessment,” Journal of Healthcare Engineering, vol. 4, no. 4, pp. 505–528, 2013.

Y. X. Chou, A. Zhang, P. Wang, and J. Gu, “Pulse rate variability estimation method based on sliding window iterative DFT and Hilbert transform,” Journal of Medical and Biological Engineering, vol. 34, no. 4, pp. 347–355, 2014.

E. Yu, D. He, Y. Su, L. Zheng, Z. Yin, and L. Xu, “Feasibility analysis for pulse rate variability to replace heart rate variability of the healthy subjects,” in 2013 IEEE International Conference on Proceedings of the Robotics and Biomimetics (ROBIO), pp. 1065–1070, Shenzhen, December 2013.

A. Schäfer and J. Vagedes, “How accurate is pulse rate variability as an estimate of heart rate variability?: a review on studies comparing photoplethysmographic technology with an electrocardiogram,” International Journal of Cardiology, vol. 166, no. 1, pp. 15–29, 2013.

C. C. Chuang, J. J. Ye, W. C. Lin, K. T. Lee, and Y. T. Tai, “Photoplethysmography variability as an alternative approach to obtain heart rate variability information in chronic pain patient,” Journal of Clinical Monitoring and Computing, vol. 29, no. 6, pp. 801–806, 2015.

M. B. Mashhadi, E. Asadi, M. Eskandari, S. Kiani, and F. Marvasti, “Heart rate tracking using wrist-type photoplethysmographic (PPG) signals during physical exercise with simultaneous accelerometry,” IEEE Signal Processing Letters, vol. 23, no. 2, pp. 227–231, 2015.

J. Guyomard and R. Stortelder, Heart Rate Measurement through PPG: Heartbeat Measurement in a Wireless Headset, The Delft University of Technology, Delft, Netherlands, 2015.

U. R. Acharya, O. Faust, V. Sree et al., “Linear and nonlinear analysis of normal and CAD-affected heart rate signals,” Computer Methods and Programs in Biomedicine, vol. 113, no. 1, pp. 55–68, 2014.

M. P. Tarvainen, J. P. Niskanen, J. A. Lipponen, P. O. Ranta-Aho, and P. A. Karjalainen, “Kubios HRV-heart rate variability analysis software,” Computer Methods and Programs in Biomedicine, vol. 113, no. 1, pp. 210–220, 2014.

D. Cornforth, H. Jelinek, and M. Tarvainen, “A comparison of nonlinear measures for the detection of cardiac autonomic neuropathy from heart rate variability,” Entropy, vol. 17, no. 3, pp. 1425–1440, 2015.

J. Li and X. Ning, “Dynamical complexity detection in short-term physiological series using base-scale entropy,” Physical Review E Statistical Nonlinear & Soft Matter Physics, vol. 73, no. 5, p. 052902, 2006.

J. Li and X. B. Ning, “The base-scale entropy analysis of short-term heart rate variability signal,” Chinese Science Bulletin, vol. 50, no. 12, pp. 1269–1273, 2005.

C. H. Bian, Q. L. Ma, J. Si et al., “Sign series entropy analysis of short-term heart rate variability,” Chinese Science Bulletin, vol. 54, no. 24, pp. 4610–4615, 2009.

A. L. Song, “Optimum parameters setting in symbolic dynamics of heart rate variability analysis,” Acta Physica Sinica, vol. 60, no. 2, pp. 307–351, 2011.

“PhysioNet/fantasia database,” December 2016, http://www.physionet.org/physiobank/database/fantasia/.

A. H. Zhang, P. Wang, and Y. X. Chou, “Peak detection of pulse signal based on dynamic difference threshold,” Journal of Jilin University, vol. 44, no. 3, pp. 847–853, 2014.

Y. X. Chou and A. H. Zhang, “Age-related alterations in the sign series entropy of short-term pulse rate variability,” in 11th ICIC Conference on Advanced Intelligent Computing Theories and Applications, pp. 723–729, Fuzhou, China, 2015.