Next Article in Journal
Hierarchy of Relaxation Times and Residual Entropy: A Nonequilibrium Approach
Next Article in Special Issue
Complexity Changes in Brain Activity in Healthy Ageing: A Permutation Lempel-Ziv Complexity Study of Magnetoencephalograms
Previous Article in Journal
Mixture and Exponential Arcs on Generalized Statistical Manifold
Previous Article in Special Issue
Complexity of Simple, Switched and Skipped Chaotic Maps in Finite Precision
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Permutation Entropy and Permutation Min-Entropy in Multiple Emotional States Analysis of RRI Time Series

1
School of Control Science and Engineering, Shandong University, Jinan, 250061, China
2
Centro de Investigaciones Ópticas (CONICET La Plata—CIC), C.C. 3, 1897 Gonnet, Argentina
3
Departamento de Ciencias Básicas, Facultad de Ingeniería, Universidad Nacional de La Plata (UNLP), 1900 La Plata, Argentina
4
School of Instrument Science and Engineering, Southeast University, Nanjing 210018, China
*
Authors to whom correspondence should be addressed.
Entropy 2018, 20(3), 148; https://doi.org/10.3390/e20030148
Submission received: 14 January 2018 / Revised: 14 February 2018 / Accepted: 15 February 2018 / Published: 26 February 2018
(This article belongs to the Special Issue Permutation Entropy & Its Interdisciplinary Applications)

Abstract

:
This study’s aim was to apply permutation entropy (PE) and permutation min-entropy (PME) over an RR interval time series to quantify the changes in cardiac activity among multiple emotional states. Electrocardiogram (ECG) signals were recorded under six emotional states (neutral, happiness, sadness, anger, fear, and disgust) in 60 healthy subjects at a rate of 1000 Hz. For each emotional state, ECGs were recorded for 5 min and the RR interval time series was extracted from these ECGs. The obtained results confirm that PE and PME increase significantly during the emotional states of happiness, sadness, anger, and disgust. Both symbolic quantifiers also increase but not in a significant way for the emotional state of fear. Moreover, it is found that PME is more sensitive than PE for discriminating non-neutral from neutral emotional states.
PACS Codes:
87.85.Ng; 05.45.Tp; 87.19.Hh; 87.19.ug; 87.19.uj

1. Introduction

Emotion recognition by using physiological signals, including electroencephalogram (EEG) [1], photoplethysmography (PPG) [2], skin temperature [3], and electrocardiogram (ECG) [4], has attracted increasing attention. Among them, ECG signals and their related RR interval (RRI) time series are mostly used.
Time-domain and frequency-domain indices were conventionally used for emotional state recognition. Guo et al. [5] applied time-domain indices, mean, coefficient of variation of RR intervals, standard deviation of the RR intervals, and standard deviation of the successive differences of the RR intervals together with a support vector machine classifier, reaching 50.3% of correct rate to discriminate negative and positive emotional states. With frequency-domain indices—low-frequency power (LF), high-frequency power (HF) and low-frequency power to high-frequency power ratio (LF/HF)—54.3% of correct rate has been obtained. Kim et al. [6] used spectrum, amplitude, mean, maximum, and standard deviation together with a support vector machine classifier to achieve an accuracy of 61.8% for classifying the emotional states of sadness, stress, surprise, and anger. Jang et al. [7] used mean, HF, and LF of RRI time series and other physiological signals to classify boredom, pain, and surprise and reached 84.7%, the best classification rate. Ma et al. [8] used mean, median, standard deviation, maxima, minima, maximal interval, and spectrum average of heart rate variability (HRV) and obtained an average of 88.88% recognition rate for classifying joy and sadness.
The use of nonlinear dynamic methods can enhance the accuracy of emotion recognition. Indeed, Valenza et al. [9] evaluated deterministic chaos, recurrence plots, and detrended fluctuation analysis from HRV, respiration activity, and electrodermal response series and attained an accuracy of 90%. Features from the instantaneous spectrum and bispectrum, as well as the dominant Lyapunov exponent, were extracted from RRI time series and an overall accuracy of 79.29% was achieved in recognizing four emotional states (happiness, sadness, anger, and relaxation) with 79.15% on the valence axis and 83.55% on the arousal axis [10]. Goshvarpour et al. [11] used Poincare plot indices, an efficient nonlinear tool to extract features from RRI time series, and obtained the best classification rate of 97.45% to distinguish the neutral, happiness, relaxation, sadness, and fear emotional states.
Among the nonlinear approaches, permutation entropy (PE) plays an important role because of its robust, fast, and easy computation. This entropic measure, introduced by Bandt and Pompe [12], is based on the comparison between the neighboring values of a time series. PE quantifies the diversity of the ordinal pattern distribution. It is maximal in the case of a uniform distribution (white noise). When some ordinal patterns appear more frequently than others, then the PE decreases. It has been applied in different fields [13,14,15,16]. Within the physiological field, PE has been frequently used to analyze EEG signals for different purposes [17,18,19]. PE has also been used to characterize the RRI time series. Actually, PE of HRV was confirmed to separate workload levels between different stages of anesthesia [20]. Aziz and Arif [21] found that the PE of the RRI time series of normal sinus rhythm subjects becomes greater than that of the RRI time series from congestive heart failure subjects. Graff et al. [22] reported that PE is significantly different between groups with cardio depressive or negative results for the head-up tilt test.
However, little attention has been paid so far to PE of the RRI time series for emotion recognition. To the best of our knowledge, the research developed by García-Martínez et al. [23] is the only one in which this symbolic entropic measure has been considered for automatic recognition of emotions. More precisely, EEG records have been analyzed for discerning between emotions of calm and stress. Results obtained confirm that PE is useful for this purpose. In addition, permutation min-entropy (PME) is an improved symbolic alternative for identifying the existence of hidden temporal correlations in time series and it allows a better discrimination of time series with similar temporal correlations [24]. It is defined as the Rényi PE when the order q goes to infinity. It can be shown that, in this limit, the min-entropy is a function of the highest probability only. Since it is estimated by implementing the Bandt and Pompe scheme, all the advantages associated with this symbolic method are naturally inherited [24]. PME improves PE in this sense. It has not been previously implemented for characterizing different emotional states. Thus, this study aimed to analyze how the PE and PME of the RRI time series change during the happiness, sadness, anger, fear, and disgust emotional states when comparing with the neutral emotional state. The effect of using different embedding delays to estimate PE and PME has also been considered in this research.

2. Methods

2.1. Subjects

Sixty subjects with a mean age of 23 years, a height of 169 centimeters and a weight of 61 kilograms were enrolled in this study (30 females and 30 males). All subjects were healthy graduate students with no history of psychotropic medication or abuse of alcohol, confirmed by a medical examination at the University Hospital. All subjects signed the consent forms in the study. Before and after the signal recording, heart rate (HR), systolic blood pressure (SBP), and diastolic blood pressure (DBP) were measured using the OMRON HEM-7051 electronic sphygmomanometer device. The details of the subjects are summarized in Table 1.

2.2. Experimental Procedure and Data Collection

ECG signals were collected under neutral and five non-neutral emotional states of happiness, sadness, anger, fear, and disgust. Video clips (each about ten minutes) of ‘World Heritage In China’ (a documentary), ‘Top Funny Comedian’ (a situation comedy), ‘Nuan Chun’ (a touching movie), ‘Never Talk to Strangers’ (a domestic violence television series), ‘The Conjuring’ (a horror movie), and ‘The Fly’ were used to induce the emotional states of neutral, happiness, sadness, anger, fear, and disgust, respectively. Figure 1 shows still images from the six emotion-stimulating video clips.
Subjects were tested in a quiet room and while sitting in a comfortable chair. Ten minutes were given to the subjects to calm themselves. Before playing the emotion-stimulating videos, subjects with ECG electrodes watched a 1-minute video to adapt to the experimental conditions. Then, emotion-stimulating videos were played in randomized order to the subjects. A gap of 10 min between videos was used to make sure the subjects could recover from the previous emotional state. ECG signals were recorded by a RM6240B signal acquirement system with standard leads at a rate of 1000 Hz. For each emotional state, ECG signals were acquired from the fourth to eighth minute (for 5 min) of each video clip.
After signal collection, a sym8 wavelet filter was used to remove the high-frequency interference and baseline drift in the ECG signal [25]. R peaks were detected by the difference threshold algorithm [26]. Additionally, we checked each ECG signal and manually removed the error R peaks before generating the RRI time series. The RRI time series has a length of about 300 to 400 intervals for each emotional state of each subject.

2.3. Permutation Entropy and Permutation Min-Entropy

A time series {x(t), t = 1, 2, …, N } can be partitioned into subsets of length D and lag τ as {X(t) = [x(t), x(t + τ),…, x(t + (D − 1)τ)], t = 1, 2,…, N − (D − 1)τ)}, where D is the embedding dimension and τ is the embedding delay. X(t) is rearranged in an ascending order and an ordinal pattern π is created. For D different numbers, there will be D! possible ordinal patterns. For example, when D = 3, there are 3! = 6 ordinal patterns. They are denoted as (supposing τ = 1):
x ( t ) < x ( t + 1 ) < x ( t + 2 ) { 123 } x ( t ) < x ( t + 2 ) < x ( t + 1 ) { 132 } x ( t + 1 ) < x ( t ) < x ( t + 2 ) { 213 } x ( t + 1 ) < x ( t + 2 ) < x ( t ) { 231 } x ( t + 2 ) < x ( t ) < x ( t + 1 ) { 312 } x ( t + 2 ) < x ( t + 1 ) < x ( t ) { 321 }
The above definitions apply to time series that have a continuous distribution so that equal values are very rare. Otherwise, equalities may cause a deviation from true PE [27]. In order to avoid this drawback, a very small amount of noise can be added to break ties.
The relative frequency of each π (# means number) is defined as:
p ( π ) = # { X ( t )   has   type   π   } N ( D 1 ) τ
The permutation entropy is defined as:
PE = p ( π ) ln   p ( π )
It is clear that the bound of PE is [0, lnD!]. It is also worth remarking here that the condition N >> D! should be satisfied in order to achieve a reliable estimation of the relative frequency for the different ordinal patterns. As an illustrative example, an artificial time series x(t) ={3, 5, 2, 1, 4, 8, 5, 6} can be symbolized with ordinal patterns as follows:
[ 3   5   2 ] { 312 } [ 5   2   1 ] { 321 } [ 2   1   4 ] { 213 } [ 1   4   8 ] { 123 } [ 4   8   5 ] { 132 } [ 8   5   6 ] { 231 }
The PE of this toy time series is calculated as:
PE = 1 / 6 ln ( 1 / 6 ) 1 / 6 ln ( 1 / 6 ) 1 / 6 ln ( 1 / 6 ) 1 / 6 ln ( 1 / 6 ) 1 / 6 ln ( 1 / 6 ) 1 / 6 ln ( 1 / 6 ) 1.7918
For the purpose of identifying the presence of hidden correlational structures in complex time series, Zunino et al. [24] defined a new entropy, the permutation min-entropy (PME). It is applicable to noisy real-time series from all classes of systems, deterministic and stochastic, without the need to require any assumption about the generating process. PME is defined as Equation (6). It is obvious that the bound of PME is [0, lnD!]:
PME = ln ( max i = 1 , , D ! p ( π i ) )
Taking into account the length of the RRI time series (300–400 data points), the embedding dimensions D = 3 and D = 4 and embedding delays τ between 1 and 10 have been considered in this paper.

2.4. Statistical Analysis

In this study, PE and PME were compared as discriminators of different emotional states. Normal distributions of the results were tested by the Kolmogorov–Smirnov test. Throughout all results, data were normally distributed. Hence, the differences between the neutral and the five non-neutral emotional states (happiness, sadness, anger, fear, disgust) were compared by paired Student’s t-test. All statistical analyses were performed using the SPSS software (Version 20, IBM, USA). The statistical significance was accepted at P < 0.05.

3. Results

3.1. Permutation and Permutation Min-Entropy Analysis

The value of permutation entropy is based on the distribution of ordinal patterns. The distribution of ordinal patterns during the six emotional states was observed. The distribution of ordinal patterns with D = 3 and τ = 1 of the neutral and happiness emotional states of 60 subjects are shown in Figure 2. Similar results were obtained for the other emotional states. We found that patterns {123} and {321}, {132} and {213}, {231} and {312}, respectively, had similar frequency percentages as Graff et al.’s study [28]. Subtle changes in the distribution of the ordinal patterns could reflect on the value of PE, and also on PME.
The effect of embedding delays on PE and PME was analyzed and it was found that PE and PME were near to theoretical maximum when τ was higher than 5. Furthermore, the difference between neutral and non-neutral emotional states with τ ∈ [6, 10] was similar to τ = 1. Mean and SD of PE and PME with embedding dimension D = 3 and D = 4 for τ ∈ [1, 5] are plotted in Figure 3 and Figure 4, respectively. The results show that there are significant differences in PE and PME when τ = 1 for four emotional states (happiness, sadness, anger, and disgust).
Hereafter, taking into account these findings, τ = 1 was chosen as the optimal one for distinguishing between neutral and non-neutral emotional states with embedding dimensions D = 3 and D = 4.

3.2. Comparisons between Neutral and Non-Neutral Emotional States

Figure 5 shows the values estimated for the PE (with (a) D = 3 and τ = 1 and (b) D = 4 and τ = 1) for each of the 60 subjects under the six emotional states. It was found that PE of the non-neutral emotions for the mass of subjects was higher than the neutral emotional state. The statistical results (mean and SD) are shown in the first plots of Figure 3a,b, respectively. Compared with the neutral emotional state (1.6834 ± 0.06), PE increased with embedding dimension D = 3 in all five non-neutral emotional states: happiness (1.7205 ± 0.05, P < 0.01), sadness (1.7065 ± 0.05, P < 0.01), anger (1.7098 ± 0.04, P < 0.01), fear (1.6937 ± 0.05, P = 0.10), and disgust (1.7186 ± 0.05, P < 0.01). With embedding dimension D = 4 and comparing with the neutral emotional state (2.8295 ± 0.16), PE also increased in all five non-neutral emotional states: happiness (2.9395 ± 0.14, P < 0.01), sadness (2.8762 ± 0.13, P < 0.01), anger (2.8701 ± 0.14, P < 0.01), fear (2.8413 ± 0.13, P = 0.32), and disgust (2.8876 ± 0.14, P < 0.01). The increases are significant for all except the fear emotional state.
Figure 6 shows the values estimated for the PME (with (a) D = 3 and τ = 1 and (b) D = 4 and τ = 1) for each of the 60 subjects under the six emotional states; the results are similar to PE. The first plots of Figure 4a,b show the results obtained for the PME (also with (a) D = 3 and τ = 1 and (b) D = 4 and τ = 1). Compared with the neutral emotional state (1.2099 ± 0.14), PME increased with embedding dimension D = 3 in all five non-neutral emotional states: happiness (1.3360 ± 0.14, P < 0.01), sadness (1.2699 ± 0.12, P < 0.01), anger (1.2781 ± 0.14, P < 0.01), fear (1.2282 ± 0.13, P = 0.27), and disgust (1.3114 ± 0.14, P < 0.01). With embedding dimension D = 4 and compared with the neutral emotional state (2.0870 ± 0.32), PME also increased in all five non-neutral emotional states: happiness (2.1778 ± 0.34, P < 0.01), sadness (2.1993 ± 0.26, P < 0.01), anger (2.1790 ± 0.26, P < 0.01), fear (2.1222 ± 0.27, P = 0.22), and disgust (2.2083 ± 0.26, P < 0.01). Similar to the PE, the change in PME is significant for all except the fear emotional state.
Since the number of females and males enrolled in this research were the same, possible differences due to gender were tested. As shown in Figure 7 and Figure 8, the results show that compared with the neutral emotional state, both PE and PME increased in all five non-neutral emotional states when embedding dimension D = 3 and D = 4.

3.3. Comparisons of Permutation and Permutation Min-Entropy

To analyze which of the two symbolic entropies is more sensitive for discriminating emotional states, Figure 9 shows the relative increment associated with the PE and PME for the non-neutral emotional states compared to the neutral emotion state with both embedding dimensions (D = 3 and D = 4). We observed that the relative increments of the PME are larger than those obtained for the PE for all non-neutral emotional states for the two values of D. Thus, PME seems to be more sensitive than PE for characterizing emotional states.

4. Discussion

PE and PME are quantifiers of nonlinear dynamic characteristics that could evaluate the inherent complexities of the RRI time series under different emotional states. In this study, we found that PE and PME increase among five non-neutral emotional states, indicating the higher randomness of their associated RRI time series. Moreover, PME appears to be more sensitive for discriminating non-neutral from neutral emotional states. Previous studies showing that other nonlinear measures change during non-neutral emotional states when comparing with neutral emotional states agree with our results that emotions affect cardiac activity [29,30].
Cardiac activity is under the control of the autonomic nervous system [31]. A previous study reported a significant reduction in the entropy related to the RRI time series during head-up tilting, and the results indicate that a sympathetic dominance increases regularity leading to a complexity reduction [32]. On the other hand, Porta et al. [33] found that a high dose of atropine, blocking the vagal system, significantly reduces the RRI time series complexity compared to baseline conditions, while sympathetic blocking with propranolol does not have any effect on the RR interval series complexity. Thus, the authors concluded that vagal efferent activity is the main contributor to the RRI time series complexity and complexity decreases due to a vagal dominance reduction. However, Weippert et al. [34] found that a reduction in complexity after dual autonomic blocking is not as strong as after vagal blocking alone. Therefore, they indicated that both vagal and sympathetic modulations contribute to RRI time series complexity. From all the above-mentioned studies, we hypothesize that subjects under happiness, sadness, fear, anger, and disgust emotional states have increased vagal activity and decreased sympathetic modulation, making the RRI time series more random.
This study has some limitations that should be considered in the future. First, all of the subjects enrolled were healthy young people. To further investigate this concern, a large number of subjects and patients with different emotional diseases will be studied. Second, the experiment lasted about an hour and a half for each subject. Long experiments may cause subjects to feel tired, which may affect emotion stimulation. Third, the length of the time series is short (around 300/400 values) and discrimination between neutral and non-neutral emotional states is only on average. A better discrimination may be reached with longer records.
In conclusion, this study is the first study to compare PE and PME from RRI time series during different emotional states. PE and PME have been found to have a significant difference for non-neutral emotional states when compared with neutral emotional states. This study could help to better understand the differences between the RRI time series of different emotional states.

Acknowledgments

This work was supported by the Neutral Science Foundation of Shandong Province in China (grant 2014ZRE27230), Key Research and Development Program of Shandong Province (grant 2016GSF120009), National Neutral Science Foundation of China (grants 61671275 and 61201049), Southeast-Lenovo Wearable Heart-Sleep-Emotion Intelligent Monitoring Lab, and Consejo Nacional de Investigaciones Científicasy Técnicas (CONICET), Argentina.

Author Contributions

Yirong Xia, Licai Yang, and Chengyu Liu designed the study. Yirong Xia, Hongyu Shi, Yuan Zhuang, and Chengyu Liu were responsible for data collection and analysis. Yirong Xia, Licai Yang, and Chengyu Liu reviewed relevant literature and interpreted the acquired data. Yirong Xia, Chengyu Liu and Luciano Zunino drafted the manuscript. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mehmood, R.M.; Lee, H.J. A novel feature extraction method based on late positive potential for emotion recognition in human brain signal patterns. Comput. Electr. Eng. 2016, 53, 444–457. [Google Scholar] [CrossRef]
  2. Chang, C.Y.; Chang, C.W.; Zheng, J.Y.; Chung, P.C. Physiological emotion analysis using support vector regression. Neurocomputing 2013, 122, 79–87. [Google Scholar] [CrossRef]
  3. Gouizi, K.; Reguig, F.B.; Maaoui, C. Emotion recognition from physiological signals. J. Med. Eng. Technol. 2011, 35, 300–307. [Google Scholar] [CrossRef] [PubMed]
  4. AlZoubi, O.; Sidney, K.D.; Calvo, R.A. Detecting Naturalistic Expressions of Nonbasic Affect Using Physiological Signals. IEEE Trans. Affect. Comput. 2012, 3, 298–310. [Google Scholar] [CrossRef]
  5. Guo, H.W.; Huang, Y.S.; Lin, C.H.; Chien, J.C.; Haraikawa, K.; Shieh, J.S. Heart Rate Variability Signal Features for Emotion Recognition by Using Principal Component Analysis and Support Vectors Machine. In Proceedings of the 16th International Conference on Bioinformatics and Bioengineering (BIBE), Taichung, Taiwan, 31 October–2 November 2016; pp. 274–277. [Google Scholar]
  6. Kim, K.H.; Bang, S.W.; Kim, S.R. Emotion recognition system using short-term monitoring of physiological signals. Med. Biol. Eng. Comput. 2004, 42, 419–427. [Google Scholar] [CrossRef] [PubMed]
  7. Jang, E.H.; Park, B.J.; Park, M.S.; Kim, S.H.; Sohn, J.H. Analysis of physiological signals for recognition of boredom, pain, and surprise emotions. J. Physiol. Anthropol. 2015, 34, 25. [Google Scholar] [CrossRef] [PubMed]
  8. Ma, C.W.; Liu, G.Y. Feature Extraction, Feature Selection and Classification from Electrocardiography to Emotions. In Proceedings of the International Conference on Computational Intelligence and Natural Computing, Wuhan, China, 6–7 June 2009; pp. 190–193. [Google Scholar]
  9. Valenza, G.; Lanata, A.; Scilingo, E.P. The Role of Nonlinear Dynamics in Affective Valence and Arousal Recognition. IEEE Trans. Affect. Comput. 2012, 3, 237–249. [Google Scholar] [CrossRef]
  10. Valenza, G.; Citi, L.; Lanata, A.; Scilingo, E.P.; Barbieri, R. Revealing Real-Time Emotional Responses: A Personalized Assessment based on Heartbeat Dynamics. Sci. Rep. 2014, 4, 4998. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Goshvarpour, A.; Abbasi, A.; Goshvarpour, A. Indices from lagged poincare plots of heart rate variability: An efficient nonlinear tool for emotion discrimination. Australas. Phys. Eng. Sci. Med. 2017, 40, 277–287. [Google Scholar] [CrossRef] [PubMed]
  12. Bandt, C.; Pompe, B. Permutation entropy: A natural complexity measure for time series. Phys. Rev. Lett. 2002, 88, 174102. [Google Scholar] [CrossRef] [PubMed]
  13. Gao, Y.; Villecco, F.; Li, M.; Song, W. Multi-Scale Permutation Entropy Based on Improved LMD and HMM for Rolling Bearing Diagnosis. Entropy 2017, 19, 176. [Google Scholar] [CrossRef]
  14. Villecco, F.; Pellegrino, A. Entropic measure of epistemic uncertainties in multibody system models by axiomatic design. Entropy 2017, 19, 291. [Google Scholar] [CrossRef]
  15. Sena, P.; D’Amore, M.; Pappalardo, M.; Pellegrino, A.; Fiorentino, A.; Villecco, F. Studying the influence of cognitive load on driver’s performances by a Fuzzy analysis of Lane Keeping in a drive simulation. IFAC Proc. Vol. 2013, 46, 151–156. [Google Scholar] [CrossRef]
  16. Villecco, F.; Pellegrino, A. Evaluation of Uncertainties in the Design Process of Complex Mechanical Systems. Entropy 2017, 19, 475. [Google Scholar] [CrossRef]
  17. Nicolaou, N.; Georgiou, J. The Use of Permutation Entropy to Characterize Sleep Electroencephalograms. Clin. EEG Neurosci. 2011, 42, 24–28. [Google Scholar] [CrossRef] [PubMed]
  18. Morabito, F.C.; Labate, D.; La Foresta, F.; Bramanti, A.; Morabito, G.; Palamara, I. Multivariate Multi-Scale Permutation Entropy for Complexity Analysis of Alzheimer’s Disease EEG. Entropy 2012, 14, 1186–1202. [Google Scholar] [CrossRef]
  19. Li, J.; Yan, J.; Liu, X.; Ouyang, G. Using Permutation Entropy to Measure the Changes in EEG Signals During Absence Seizures. Entropy 2014, 16, 3049–3061. [Google Scholar] [CrossRef]
  20. Martin, J.; Schneider, F.; Kowalewskij, A.; Jordan, D.; Hapfelmeier, A.; Kochs, E.; Wagner, K.; Schulz, C. Linear and non-linear heart rate metrics for the assessment of anaesthetists’ workload during general anaesthesia. Br. J. Anaesth. 2016, 117, 767–774. [Google Scholar] [CrossRef] [PubMed]
  21. Aziz, W.; Arif, M. Multiscale permutation entropy of physiological time series. In Proceedings of the 9th International Conference of Multitopic, Karachi, Pakistan, 24–25 December 2005; pp. 368–373. [Google Scholar]
  22. Graff, B.; Graff, G.; Makowiec, D.; Kaczkowska, A.; Wejer, D.; Budrejko, S.; Kozłowski, D.; Narkiewicz, K. Entropy measures in the assessment of heart rate variability in patients with cardiodepressive vasovagal syncope. Entropy 2015, 17, 1007–1022. [Google Scholar] [CrossRef]
  23. García-Martínez, B.; Martínez-Rodrigo, A.; Zangróniz, R.; Pastor, J.M.; Alcaraz, R. Symbolic analysis of brain dynamics detects negative stress. Entropy 2017, 19, 196. [Google Scholar] [CrossRef]
  24. Zunino, L.; Olivares, F.; Rosso, O.A. Permutation min-entropy: An improved quantifier for unveiling subtle temporal correlations. EPL 2015, 109, 10005. [Google Scholar] [CrossRef]
  25. Rioul, O.; Martin, V. Wavelets and signal processing. IEEE Signal Process. Mag. 1991, 8, 14–38. [Google Scholar] [CrossRef]
  26. Pan, J.; Tompkins, W.J. A real-time QRS detection algorithm. IEEE Trans. Biomed. Eng. 1985, 32, 230–236. [Google Scholar] [CrossRef] [PubMed]
  27. Zunino, L.; Olivares, F.; Scholkmann, F.; Rosso, O.A. Permutation entropy based time series analysis: Equalities in the input signal can lead to false conclusions. Phys. Lett. 2017, 381, 1883–1892. [Google Scholar] [CrossRef]
  28. Graff, G.; Graff, B.; Kaczkowska, A.; Makowiec, D.; Amigó, J.; Piskorski, J.; Narkiewicz, K.; Guzik, P. Ordinal pattern statistics for the assessment of heart rate variability. Eur. Phys. J. Spec. Top. 2013, 222, 525–534. [Google Scholar] [CrossRef]
  29. Chendi, W.; Feng, W. An emotional analysis method based on heart rate variability. In Proceedings of the 2012 IEEE-EMBS International Conference on Biomedical and Health Informatics, Hong Kong, China, 5–7 January 2012; pp. 104–107. [Google Scholar]
  30. Valenza, G.; Allegrini, P.; Lanata, A.; Scilingo, E.P. Dominant Lyapunov exponent and approximate entropy in heart rate variability during emotional visual elicitation. Front. Neuroeng. 2012, 5, 1–7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Catalano, J.T. Guide to ECG; J.B. Lippincott: Philadelphia, PA, USA, 1993. [Google Scholar]
  32. Porta, A.; Gnecchi-Ruscone, T.; Tobaldini, E.; Guzzetti, S.; Furlan, R.; Montano, N. Progressive decrease of heart period variability entropy-based complexity during graded head-up tilt. J. Appl. Physiol. 2007, 103, 1143–1149. [Google Scholar] [CrossRef] [PubMed]
  33. Porta, A.; Castiglioni, P.; Bari, V.; Bassani, T.; Marchi, A.; Cividjian, A.; Quintin, L.; Di Rienzo, M. K-nearest-neighbor conditional entropy approach for the assessment of the short-term complexity of cardiovascular control. Physiol. Meas. 2013, 34, 17–33. [Google Scholar] [CrossRef] [PubMed]
  34. Weippert, M.; Behrens, M.; Rieger, A.; Behrens, K. Sample Entropy and Traditional Measures of Heart Rate Dynamics Reveal Different Modes of Cardiovascular Control During Low Intensity Exercise. Entropy 2014, 16, 5698–5711. [Google Scholar] [CrossRef]
Figure 1. Still images from emotion-stimulating videos for inducing (a) neutral, (b) happiness, (c) sadness, (d) anger, (e) fear, and (f) disgust emotions.
Figure 1. Still images from emotion-stimulating videos for inducing (a) neutral, (b) happiness, (c) sadness, (d) anger, (e) fear, and (f) disgust emotions.
Entropy 20 00148 g001
Figure 2. Distributions of ordinal patterns with D = 3 and τ = 1 of the neutral and happiness emotional states.
Figure 2. Distributions of ordinal patterns with D = 3 and τ = 1 of the neutral and happiness emotional states.
Entropy 20 00148 g002
Figure 3. Mean and SD of permutation entropy (PE) with (a) D = 3, (b) D = 4 and embedding delays τ between 1 and 5 of the RRI time series during six emotional states. ** mean P < 0.01. * mean P < 0.05.
Figure 3. Mean and SD of permutation entropy (PE) with (a) D = 3, (b) D = 4 and embedding delays τ between 1 and 5 of the RRI time series during six emotional states. ** mean P < 0.01. * mean P < 0.05.
Entropy 20 00148 g003
Figure 4. Mean and SD of permutation min-entropy (PME) with (a) D = 3, (b) D = 4 and embedding delays τ between 1 and 5 of the RRI time series during six emotional states ** mean P < 0.01. * mean P < 0.05.
Figure 4. Mean and SD of permutation min-entropy (PME) with (a) D = 3, (b) D = 4 and embedding delays τ between 1 and 5 of the RRI time series during six emotional states ** mean P < 0.01. * mean P < 0.05.
Entropy 20 00148 g004
Figure 5. PE results with (a) D = 3, τ = 1 and (b) D = 4, τ = 1 for each of the 60 subjects under the six emotional states.
Figure 5. PE results with (a) D = 3, τ = 1 and (b) D = 4, τ = 1 for each of the 60 subjects under the six emotional states.
Entropy 20 00148 g005
Figure 6. PME results with (a) D = 3, τ = 1 and (b) D = 4, τ = 1 for each of the 60 subjects under the six emotional states.
Figure 6. PME results with (a) D = 3, τ = 1 and (b) D = 4, τ = 1 for each of the 60 subjects under the six emotional states.
Entropy 20 00148 g006
Figure 7. Mean and SD of PE and PME for males and females with (a) D = 3, τ = 1 and (b) D = 4, τ = 1 under the six emotional states. ** mean P < 0.01.* mean P < 0.05.
Figure 7. Mean and SD of PE and PME for males and females with (a) D = 3, τ = 1 and (b) D = 4, τ = 1 under the six emotional states. ** mean P < 0.01.* mean P < 0.05.
Entropy 20 00148 g007
Figure 8. Mean and SD of PME for males and females with (a) D = 3, τ = 1 and (b) D = 4, τ = 1 under the six emotional states. ** mean P < 0.01.* mean P < 0.05.
Figure 8. Mean and SD of PME for males and females with (a) D = 3, τ = 1 and (b) D = 4, τ = 1 under the six emotional states. ** mean P < 0.01.* mean P < 0.05.
Entropy 20 00148 g008
Figure 9. Relative increments of PE and PME with (a) D = 3 and (b) D = 4 during the non-neutral emotional states when comparing with the neutral emotional state (mean and standard error).
Figure 9. Relative increments of PE and PME with (a) D = 3 and (b) D = 4 during the non-neutral emotional states when comparing with the neutral emotional state (mean and standard error).
Entropy 20 00148 g009
Table 1. Basic clinical characteristics from all 60 subjects.
Table 1. Basic clinical characteristics from all 60 subjects.
VariablesMaleFemaleTotal
Number of subjects303060
Age (year)23 ± 123 ± 223 ± 2
Height (cm)175 ± 5163 ± 4169 ± 8
Weight (kg)69 ± 853 ± 661 ± 11
Body mass index (kg/m2)22 ± 320 ± 221 ± 3
Heart rate (beats/min)77 ± 1178 ± 1178 ± 11
Systolic blood pressure (mmHg)123 ± 10107 ± 7116 ± 12
Diastolic blood pressure (mmHg)76 ± 870 ± 673 ± 8
Note: data are expressed as mean +/− SD.

Share and Cite

MDPI and ACS Style

Xia, Y.; Yang, L.; Zunino, L.; Shi, H.; Zhuang, Y.; Liu, C. Application of Permutation Entropy and Permutation Min-Entropy in Multiple Emotional States Analysis of RRI Time Series. Entropy 2018, 20, 148. https://doi.org/10.3390/e20030148

AMA Style

Xia Y, Yang L, Zunino L, Shi H, Zhuang Y, Liu C. Application of Permutation Entropy and Permutation Min-Entropy in Multiple Emotional States Analysis of RRI Time Series. Entropy. 2018; 20(3):148. https://doi.org/10.3390/e20030148

Chicago/Turabian Style

Xia, Yirong, Licai Yang, Luciano Zunino, Hongyu Shi, Yuan Zhuang, and Chengyu Liu. 2018. "Application of Permutation Entropy and Permutation Min-Entropy in Multiple Emotional States Analysis of RRI Time Series" Entropy 20, no. 3: 148. https://doi.org/10.3390/e20030148

APA Style

Xia, Y., Yang, L., Zunino, L., Shi, H., Zhuang, Y., & Liu, C. (2018). Application of Permutation Entropy and Permutation Min-Entropy in Multiple Emotional States Analysis of RRI Time Series. Entropy, 20(3), 148. https://doi.org/10.3390/e20030148

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop