Next Article in Journal
Transfer Entropy
Next Article in Special Issue
Entropy Change of Biological Dynamics in Asthmatic Patients and Its Diagnostic Value in Individualized Treatment: A Systematic Review
Previous Article in Journal
A Co-Opetitive Automated Negotiation Model for Vertical Allied Enterprises Teams and Stakeholders
Previous Article in Special Issue
Coarse-Graining Approaches in Univariate Multiscale Sample and Dispersion Entropy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Centered and Averaged Fuzzy Entropy to Improve Fuzzy Entropy Precision

by
Jean-Marc Girault
1 and
Anne Humeau-Heurtier
2,*
1
Groupe ESEO, 49000 Angers, France
2
Laboratoire Angevin de Recherche en Ingénierie des Systèmes (LARIS), Univ Angers, 49000 Angers, France
*
Author to whom correspondence should be addressed.
Submission received: 16 March 2018 / Revised: 13 April 2018 / Accepted: 13 April 2018 / Published: 15 April 2018

Abstract

:
Several entropy measures are now widely used to analyze real-world time series. Among them, we can cite approximate entropy, sample entropy and fuzzy entropy (FuzzyEn), the latter one being probably the most efficient among the three. However, FuzzyEn precision depends on the number of samples in the data under study. The longer the signal, the better it is. Nevertheless, long signals are often difficult to obtain in real applications. This is why we herein propose a new FuzzyEn that presents better precision than the standard FuzzyEn. This is performed by increasing the number of samples used in the computation of the entropy measure, without changing the length of the time series. Thus, for the comparisons of the patterns, the mean value is no longer a constraint. Moreover, translated patterns are not the only ones considered: reflected, inversed, and glide-reflected patterns are also taken into account. The new measure (so-called centered and averaged FuzzyEn) is applied to synthetic and biomedical signals. The results show that the centered and averaged FuzzyEn leads to more precise results than the standard FuzzyEn: the relative percentile range is reduced compared to the standard sample entropy and fuzzy entropy measures. The centered and averaged FuzzyEn could now be used in other applications to compare its performances to those of other already-existing entropy measures.

1. Introduction

Approximate entropy (ApEn) and sample entropy (SampEn) algorithms are now widely used to quantify the irregularity of experimental time series [1,2]. They both rely on the evaluation of vectors’ similarity. However, in both ApEn and SampEn, the vectors’ similarity is based on the Heaviside function, a function that has rigid boundaries. Thus, the contributions of samples inside the boundary are treated equally, but the samples outside the boundary are left out. However, in the real world, boundaries between classes may be ambiguous: it is often difficult to determine if an input pattern belongs totally to a class. To overcome this lack of reality in ApEn and SampEn algorithms, Chen et al. proposed the fuzzy entropy (FuzzyEn) algorithm [3]. In the latter case, the vectors’ similarity is defined by the soft and continuous boundaries of a fuzzy function. Since its introduction, it has been reported that FuzzyEn leads to better performance than ApEn or SampEn [4,5,6]. FuzzyEn presents a stronger relative consistency and shows less dependence on data length than ApEn and SampEn [3].
Nevertheless, the number of samples in a signal still plays a role in the precision of FuzzyEn: the shorter the signal, the lower the number of vectors, and thus, the lower the precision of FuzzyEn (i.e., the larger the standard deviation). Therefore, to obtain more precise entropy values, the longer the signal, the better it is. In practical situations (real data), this may be a challenge. Indeed, it is often difficult to have long recordings, particularly in the biomedical field where patients may have difficulty to stay still or to cooperate.
This is why we herein propose a new fuzzy entropy measure that presents better precision than the traditional FuzzyEn measure. This is performed by increasing the number of samples used in the computation, without changing the length of the time series.
The paper is organized as follows. The original algorithm of FuzzyEn is first detailed in Section 2; then the new entropy measure is described. The synthetic and biomedical data (fetal heart rate time series) used in our work are introduced in Section 3. In Section 4, we first present, analyze, and discuss the results obtained with the synthetic data. We then describe and interpret the results obtained with the biomedical time series. We finally end with the conclusion.

2. Standard Fuzzy Entropy and the New Entropy Measure

In this section, we recall the FuzzyEn concept based on the use of a membership function. For this purpose, the generalized Gaussian membership function is used since it allows the derivation of both the rectangular function used in the calculation of SampEn and the standard Gaussian function used in the calculation of FuzzyEn.

2.1. Fuzzy Entropy Algorithm

For a given discrete time series X = { x ( 1 ) , x ( 2 ) , , x ( N ) } of length N, the algorithm to compute FuzzyEn relies on the following steps [1]:
  • Split X into a series of subsequences X m ( i ) of length m starting at x ( i ) : X m ( i ) = { x ( i ) , x ( i + 1 ) , , x ( i + m 1 ) } , 1 i N m + 1 .
  • For each vector X m ( i ) , compute the similarity degree D i j m of its neighboring vector X m ( j ) using a similarity function as:
    D i j m = μ p ( d [ X m ( i ) , X m ( j ) ] , r ) ,
    where the membership function μ p reported in Figure 1 is defined d 0 as:
    μ p ( d , r ) = exp ( ( d / r ) p ) ,
    and where the distance function d is the maximum absolute difference d [ X m ( i ) , X m ( j ) ] = max 0 k m 1 ( | x ( i + k ) x ( j + k ) | ) . For p = 2 , we have the Gaussian function, and for p = , we have the rectangular function.
  • For each i ( 1 i N m + 1 ), compute ϕ i m as:
    ϕ i m ( r ) = 1 N m 1 j = 1 , j i N m D i j m .
  • Construct φ m and φ m + 1 as:
    φ m ( r ) = 1 N m i = 1 N m ϕ i m ( r ) ,
    φ m + 1 ( r ) = 1 N m i = 1 N m φ i m + 1 ( r ) .
  • Fuzzy entropy is then calculated as:
    F u z z y E n ( m , r ) = lim N ln φ m ( r ) φ m + 1 ( r ) ,
    which, for finite datasets, can be estimated by the statistic:
    F u z z y E n ( m , r , N ) = ln φ m ( r ) φ m + 1 ( r ) .
As shown in Figure 2, the 2-pattern ‘1’ has only one similar 2-pattern among the 27 possible 2-patterns in the time series. From the time series reported in Figure 2, the total number of similar 2-patterns is 12: (‘1’,‘15’), (‘5’,‘21’), (‘7’,‘19’), (‘8’,‘20’), (‘13’,‘24’), (‘14’,‘25’).
As for ApEn and SampEn, the statistical stability of the FuzzyEn estimation depends on the length N of the time series as reported in Equation (7). To decrease this length-dependency, several strategies can be proposed.

2.2. New Approaches

As mentioned above, from a fixed number of samples N in the time series, a way to improve the statistical stability of the entropy measurement consists in artificially increasing the number of similar m-patterns taken into account in the entropy calculation. To do so, three different ways are proposed:
  • The first approach is inspired by [3,7]. In the latter studies, the interest in centering each m-pattern has been shown. In this case, instead of limiting the search of m-patterns with the same mean value, any pattern can be taken into account. Therefore, the number of similar patterns drastically increases.
    Therefore, in the first approach, a centered m-pattern X c m ( j ) is compared to a reference centered m-pattern X c m ( i ) . The similarity degree is calculated with X c m ( i ) = { x ( i ) , x ( i + 1 ) , , x ( i + m 1 ) } x 0 ( i ) , where 1 i N m + 1 and x 0 ( i ) = 1 m j = 0 j = m 1 x ( i + j ) , through a similarity function:
    D c i j m = μ p ( d [ X c m ( i ) , X c m ( j ) ] , r ) ,
    with the same membership function as the one reported in Equation (2). The centered fuzzy entropy FuzzyEn c is thus defined as:
    F u z z y E n c ( m , r , N ) = ln φ c m ( r ) φ c m + 1 ( r ) ,
    with φ c m ( r ) = 1 N m i = 1 N m ϕ c i m ( r ) and with ϕ c i m ( r ) = 1 N m 1 j = 1 , j i N m D c i j m .
    As shown in Figure 3a, removing the mean value of 2-patterns increases the number of centered similar 2-patterns since the number of centered 2-patterns similar to ‘1’ is six compared to one when the centering approach is not used. From Figure 3b, the total number of centered similar 2-patterns is 25: (‘1’,‘9’,‘13’,‘15’,‘17’,‘24’), (‘2’,‘14’,‘25’), (‘3’,‘8’,‘20’), (‘4’,‘23’), (‘5’,‘7’,‘10’,‘19’,‘21’), (‘11’,‘18’), (‘12’,‘16’), (‘22’,‘26’). The total number of similar centered 2-patterns is much larger than no-centered 2-patterns.
  • The second approach is inspired by [8], where transformed patterns are compared to reference patterns. Thus, in the second approach, a transformed m-pattern Γ k [ X m ( j ) ] (see below) is compared to a reference m-pattern X m ( i ) . The similarity degree is calculated with the same membership function as the one reported in Equation (2):
    k D i j m = μ p ( d [ X m ( i ) , Γ k [ X m ( j ) ] ] , r ) .
    Four types of Γ k [ X m ( j ) ] operations with k = { T , R , I , G } are evaluated:
    • Γ T [ X m ( j ) ] = X m ( j + n ) corresponds to a translation of n samples, k = T ;
    • Γ R [ X m ( j ) ] = X m ( j + n ) corresponds to a reflection at the position n, k = R ;
    • Γ I [ X m ( j ) ] = X m ( j + n ) corresponds to an inversion at the position n, k = I ;
    • Γ G [ X m ( j ) ] = X m ( j + n ) corresponds to a glide reflection of n samples, k = G .
    At first sight, any type of operation could be used. However, from our point of view, only isometries (translation T, reflection R, inversion I and glide reflection G) are suitable. This statement is supported by the recent work reported in [8] where the concept of symmetry was placed back on stage in the study of time series. Indeed, in [8], it was shown that the concept of recurrences could be generalized by taking into account the symmetry properties of m-patterns. As entropy can be derived from the recurrence concept (the recurrence plot [9] is defined as R P = ( N m + 1 ) D i j with μ ( d , r ) ), from [8], four new kinds of entropy ( A p E n T , A p E n R , A p E n I , A p E n G or S a m p E n T , S a m p E n R , S a m p E n I , S a m p E n G or F u z z y E n T , F u z z y E n R , F u z z y E n I , F u z z y E n G ) can be proposed. Finally, as our ultimate goal is to increase the precision of FuzzyEn, it is more appropriate here to calculate the mean value of the four new fuzzy entropies. In this case, the averaged fuzzy entropy FuzzyEn a is defined as:
    F u z z y E n a ( m , r , N ) = ( F u z z y E n T + F u z z y E n R + F u z z y E n I + F u z z y E n G ) 4 ,
    with:
    F u z z y E n k ( m , r , N ) = ln φ k m ( r ) φ k m + 1 ( r ) ,
    with k = { T , R , I , G } for φ k m ( r ) = 1 N m i = 1 N m ϕ k i m ( r ) and ϕ k i m ( r ) = 1 N m 1 j = 1 , j i N m k D i j m . F u z z y E n T corresponds to the standard FuzzyEn measure when m > 1 .
    As shown in Figure 3b, the transformation of the 2-patterns increases the number of similar 2-patterns. From Figure 3b, for the 2-pattern (‘1’), four kinds of 2-patterns can be obtained: 2-patterns with translation (‘T’) in black (‘1’,‘15’), 2-patterns with vertical reflection (‘R’) in red (‘7’,‘19’), 2-patterns with inversion (‘I’) in green (‘13’,‘24’) and 2-patterns with glide reflection (‘G’) in blue (‘5’,‘21’). By considering all 2-patterns ranging from ‘1’–‘27’, the mean total number of symmetrical 2-patterns is N s y m = 92 with N s y m T = 12 , N s y m R = 30 , N s y m I = 24 , N s y m G = 26 .
  • The last approach compares a centered m-pattern X c m ( i ) to a transformed centered m-pattern Γ k [ X c m ( j ) ] . In this case, the centered and averaged fuzzy entropy FuzzyEn ca is defined as:
    F u z z y E n c a ( m , r , N ) = ( F u z z y E n c T + F u z z y E n c R + F u z z y E n c I + F u z z y E n c G ) 4 ,
    with:
    F u z z y E n c k ( m , r , N ) = ln φ c k m ( r ) φ c k m + 1 ( r ) ,
    with k = { T , R , I , G } for φ c k m ( r ) = 1 N m i = 1 N m ϕ c k i m ( r ) and ϕ c k i m ( r ) = 1 N m 1 j = 1 , j i N m k D c i j m . k D c i j m is defined as k D c i j m = μ p ( d [ X c m ( i ) , Γ k [ X c m ( j ) ] ] , r ) .
    As shown in Figure 3, one can observe that the combination of the centering and averaging operations globally increases the number of m-patterns taken into account in the calculation of the entropy measure. Furthermore, a centered m-pattern transformed by an inversion (‘I’) is similar to a centered m-pattern transformed by a translation (’T’). The same remark applies for glide and vertical reflection transformations of centered m-patterns.
    From Figure 3c, regarding the 2-pattern (‘1’), two kinds of centered 2-patterns can be obtained: 2-patterns (‘T’,‘I’) in black (‘1’,‘9’,‘13’,‘15’,‘17’,‘24’) and 2-patterns (‘R’,‘G’) in blue (‘5’,‘7’,‘10’,‘19’,‘21’). By considering all 2-patterns ranging from ‘1’–‘27’, the mean total number of symmetrical 2-patterns is N c s y m = 312 with N c s y m T = 86 , N c s y m R = 70 , N c s y m I = 86 and N c s y m G = 70 .
The novelty of our method therefore relies on two main points: (i) the mean value of the patterns is no longer a constraint in the computation as the patterns are centered; (ii) translated patterns, but also reflected, inversed, and glide-reflected patterns are taken into account (in the standard sample and fuzzy entropy measures, only translated patterns are considered). Therefore, for a given number of samples N in the time series, we managed to increase the number of similar m-patterns taken into account in the entropy calculation. In what follows, the new entropy measure will be applied to synthetic 1 / f β time series and biomedical datasets. Its precision will be compared to the one of the standard FuzzyEn.

3. Data Processed

3.1. Synthetic Signals

In order to analyze the new fuzzy entropy measures and to compare their performances with the ones of the standard FuzzyEn, we used 1 / f β time series, with different β values: β varied from −1 to 2 in steps of 0.2. For β > 0 , the 1 / f β signals are persistent processes with long-term correlations [10]. However, for β < 0 , the 1 / f β signals are anti-persistent processes with short-term anti-correlations [10]. From a theoretical point of view, the higher the value of β , the larger the number of correlations in the time series and, therefore, the larger the number of similar samples used in the computation of FuzzyEn. For each β value, 50 time series were simulated.

3.2. Biomedical Data

The new descriptors mentioned above were also applied to biomedical data and more precisely to fetal heart rate (FHR) time series. The latter were acquired using a homemade pulse Doppler system co-developed with Altaïs Technologies (Tours, France). This Doppler fetal monitor transmits ultrasound waves of 2.25 MHz for an acoustic power limited to 1 mW/cm 2 (for more details, see [11]). It was developed to measure both the FHR and fetal movements (pseudo-breathing, limb movements).
The study was approved by the Ethics Committee of the Clinical Investigation Centre for Innovative Technology of Tours (CIC-IT 806 CHRUof Tours). Before acquisition, the consent of each parent was obtained. All parents were over eighteen years of age, and pregnancies were single. After locating the fetal heart with an echographic scanner, 18 Doppler recordings of 30 min each were acquired at CHRU Bretonneau Tours, France. This corresponds to approximately 3600 heart beats for each recording. In order to constitute homogeneous groups without spurious data, gestations complicated by other kinds of disorders (hypertension, diabetes) were discarded. Two groups of fetuses were selected: normal and those with severe intra-uterine growth retardation (IUGR). The severe IUGR group included nine fetuses delivered prematurely by cesarean section. The normal group included nine fetuses without disorders, delivered at term by spontaneous labor. For this clinical protocol, the gestational ages of fetuses ranged from 30–34 weeks.
In what follows, the 30 min of data were processed, but also segments of 10 min and 20 min. Our goal was thus to compare the results obtained as the data length decreases. Moreover, in order to compare the results obtained between normal and IUGR groups, a Mann–Whitney test was used. A p-value strictly less than 0.05 was considered to define statistical significance.

4. Results and Discussion

In all that follows, the value of r is set at 0.1 × the standard deviation of the time series.

4.1. Results for the Synthetic Signals

In order to validate our hypothesis (that is, the greater the number of similar m-patterns taken into account in the computation, the more precise the entropy measure), we started by counting the number of similar m-patterns from 50 synthetic time series.
From 1 / f β noises generated with N = 5000 samples with β ranging from 1 to 2, the median of the mean number M N of similar 3-patterns and the median of the mean number M N c a of centered and averaged similar 3-patterns were evaluated and are reported in Table 1. As expected, the higher the sample correlation in the time series, the higher the value of β and the higher the number of similar 3-patterns. Indeed, from Table 1, when β increases from 0 to 2, M N goes from 1 to 162. When symmetrical properties and the centering operation are taken into account, M N c a goes from 21 to 9278 for β ranging from 0 to 2. From this, it can be claimed that the averaging and the centering operations increase the number of similar patterns. Furthermore, whatever the m-value, we obtain rising trends as β increases (data not shown).
In order to evaluate the performance of our new approaches, for a fixed m-value and for 50 1 / f β time series with different β values, different measures have been computed: the medians M F u z z y E n , M F u z z y E n c , M F u z z y E n a , M F u z z y E n c a and the percentiles at 75% and 25% P F u z z y E n ( 75 ) , P F u z z y E n ( 25 ) , P F u z z y E n c ( 75 ) , P F u z z y E n c ( 25 ) , P F u z z y E n a ( 75 ) , P F u z z y E n a ( 25 ) , P F u z z y E n c a ( 75 ) , P F u z z y E n c a ( 25 ) have been compared.
To quantitatively evaluate the gain brought by our new approaches in comparison with FuzzyEn, two kinds of statistics have been evaluated: percentile ranges and relative percentile ranges. The following percentile ranges have thus been computed:
  • R F = P F u z z y E n ( 75 ) P F u z z y E n ( 25 ) ;
  • R F c = P F u z z y E n c ( 75 ) P F u z z y E n c ( 25 ) ;
  • R F a = P F u z z y E n a ( 75 ) P F u z z y E n a ( 25 ) ;
  • R F c a = P F u z z y E n c a ( 75 ) P F u z z y E n c a ( 25 ) .
Finally, from the percentile ranges, the following relative percentile ranges have been evaluated:
  • ( R F R F c ) / R F c ;
  • ( R F R F a ) / R F a ;
  • ( R F R F c a ) / R F c a .
The global results are presented in Table A1, Table A2 and Table A3 reported in the Appendix and are shown in Figure 4. We observe from the tables that SampEn leads to worse results than FuzzyEn, as already shown by others. Moreover, we observe that the new approach leads to results that show a reduced percentile range compared to the standard fuzzy entropy measure. Its precision is therefore better than the other entropy measures. However, our work also has some drawbacks: the gain provided by the method depends on the signal properties. The gain differs with β values.

4.2. Results for the Fetal Heart Rate Time Series

The results obtained from FHR time series for m = 2 are presented in Figure 5 for data lengths of 10 min, 20 min, and 30 min. For the three data lengths, we observe that the normal fetuses show a significantly higher entropy value than the pathological fetuses. This is true for the two entropy measures: F u z z y E n c a and the standard F u z z y E n . This means that FHR time series are more irregular for the normal fetuses than for the pathological ones. We also observe that the p-value between the two groups decreases as the data length increases. Therefore, the longer the data, the better the separation between the two groups. However, we note that, whatever the length studied, the p-value is lower for F u z z y E n c a than for the standard F u z z y E n . Our new entropy measure is therefore more interesting for this classification purpose than the standard F u z z y E n . Other data may now be processed; see, e.g., [12,13,14].

5. Conclusions

A new entropy measure, F u z z y E n ca , is proposed to improve the precision of the standard FuzzyEn. The new measure relies on centering and averaging approaches that lead to a larger number of similar patterns used in the computation of the entropy algorithm. This is performed by removing the constraint of the mean value in the comparison of the patterns. Moreover, translated patterns are not the only ones considered: reflected, inversed, and glide-reflected patterns are also taken into account. The results obtained on 1 / f β time series reveal that F u z z y E n ca shows a greater precision than FuzzyEn. Moreover, when applied to FHR time series acquired from normal and pathological fetuses, F u z z y E n ca leads to a better discrimination between the two groups than the standard FuzzyEn. These findings could allow one to obtain entropy-based relevant information by processing shorter datasets (we could obtain the same precision as the standard FuzzyEn, but with less data). This is particularly interesting for the biomedical field. F u z z y E n ca now has to be applied to other datasets, and its performance has to be compared to those of other already-existing entropy measures.

Acknowledgments

This study was supported financially by the “Agence Nationale de la Recherche” (Project ANR-07-TECSAN-023, Surfoetus), France. Furthermore, the authors would like to thank the Clinical Investigation Centre for Innovative Technology of Tours (CIC-IT 1415 CHRU of Tours) and F. Perrotin’s team in the Obstetrics Department for their support in recording the signals.

Author Contributions

Jean-Marc Girault and Anne Humeau-Heurtier conceived and designed the study. Jean-Marc Girault and Anne Humeau-Heurtier analyzed the synthetic data. Jean-Marc Girault analyzed the FHR time series. Jean-Marc Girault and Anne Humeau-Heurtier wrote the paper. Both authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Results reported in Table A1, Table A2 and Table A3 show performances that differ with the β values. We observe that the higher the β value, the lower the gain obtained in terms of relative percentile range. This is probably due to the level of correlation between samples in the time series.
Table A1. Results obtained for 1 / f β time series, for m = 2 , N = 5000 samples and for different β values.
Table A1. Results obtained for 1 / f β time series, for m = 2 , N = 5000 samples and for different β values.
β value−1−0.8−0.6−0.4−0.200.20.40.60.811.21.41.61.82
M S a m p E n 3.463.53.533.543.573.593.583.543.453.293.042.682.211.681.130.64
M F u z z y E n 3.173.213.253.263.283.283.283.253.163.002.762.401.931.420.930.54
M F u z z y E n c 3.583.583.583.583.563.533.493.423.303.112.842.461.991.460.940.51
M F u z z y E n a 3.143.173.23.223.243.243.233.203.122.972.732.371.901.400.930.57
M F u z z y E n c a 3.573.583.583.573.563.533.483.413.303.112.832.441.961.440.940.53
R S a m p E n 0.050.080.070.080.090.080.080.100.080.050.080.110.160.240.250.23
R F u z z y E n 0.030.040.030.040.040.040.040.050.040.040.080.110.170.220.220.18
R F u z z y E n c 0.010.010.010.010.010.020.020.010.020.030.070.110.180.230.230.20
R F u z z y E n a 0.020.030.020.020.020.020.020.030.030.020.040.080.110.130.130.10
R F u z z y E n c a 0.010.010.010.010.010.010.010.010.010.020.050.080.110.140.140.11
( R S a m p E n R F u z z y E n ) / R F u z z y E n 0.630.751.241.091.171.041.100.990.900.430.120.00-0.050.100.160.29
( R S a m p E n R F u z z y E n c ) / R F u z z y E n c 3.085.775.67.066.113.804.275.853.170.680.21-0.05-0.120.080.090.14
( R S a m p E n R F u z z y E n a ) / R F u z z y E n a 1.952.093.263.344.382.532.782.722.031.330.890.400.500.891.031.39
( R S a m p E n R F u z z y E n c a ) / R F u z z y E n c a 5.478.679.19.6911.312.208.968.575.551.590.770.360.430.740.861.07
( R F u z z y E n R F u z z y E n c ) / R F u z z y E n c 1.492.861.952.872.271.351.512.441.190.180.080.050.070.020.070.12
( R F u z z y E n R F u z z y E n a ) / R F u z z y E n a 0.80.760.911.081.480.730.800.870.590.630.690.400.580.720.750.86
( R F u z z y E n R F u z z y E n c a ) / R F u z z y E n c a 2.964.513.524.134.665.463.743.812.440.810.590.360.510.590.600.61
Table A2. Same as Table A1, but for m = 3 .
Table A2. Same as Table A1, but for m = 3 .
β value−1−0.8−0.6−0.4−0.200.20.40.60.811.21.41.61.82
M S a m p E n 3.473.513.483.543.563.573.543.583.433.273.012.652.181.661.130.64
M F u z z y E n 3.023.13.123.153.23.183.163.153.042.882.612.261.801.300.840.49
M F u z z y E n c 3.283.273.283.293.273.263.223.163.052.872.592.221.771.260.800.44
M F u z z y E n a 2.993.023.043.073.113.103.103.052.972.822.582.221.771.290.840.51
M F u z z y E n c a 3.213.223.233.243.233.223.183.123.022.852.592.211.741.250.800.47
R S a m p E n 0.510.490.430.50.660.580.590.390.400.350.190.180.180.210.250.23
R F u z z y E n 0.170.160.130.170.230.170.180.140.170.110.090.110.170.190.200.16
R F u z z y E n c 0.050.040.030.030.050.030.040.040.030.050.080.110.170.200.200.16
R F u z z y E n a 0.090.070.080.070.060.070.080.070.060.060.050.070.100.130.120.09
R F u z z y E n c a 0.020.020.020.020.020.020.020.020.020.020.030.080.110.130.120.09
( R S a m p E n R F u z z y E n ) / R F u z z y E n 1.952.132.461.861.942.482.231.781.382.211.150.540.080.080.250.42
( R S a m p E n R F u z z y E n c ) / R F u z z y E n c 10.0811.611.814.5412.616.1914.518.3311.035.971.490.580.050.030.250.41
( R S a m p E n R F u z z y E n a ) / R F u z z y E n a 4.685.554.515.749.747.206.064.936.254.903.121.600.740.641.171.62
( R S a m p E n R F u z z y E n c a ) / R F u z z y E n c a 26.8822.7417.4128.4333.8130.3330.3921.0823.7816.414.911.340.590.571.091.55
( R F u z z y E n R F u z z y E n c ) / R F u z z y E n c 2.753.032.74.433.633.933.802.364.051.170.160.030.030.050.000.00
( R F u z z y E n R F u z z y E n a ) / R F u z z y E n a 0.921.10.591.362.651.351.181.132.050.840.910.690.610.520.740.85
( R F u z z y E n R F u z z y E n c a ) / R F u z z y E n c a 8.446.594.329.2810.848.008.716.949.414.421.740.520.470.450.680.79
Table A3. Same as Table A1, but for m = 4 . “-” means that an undefined value is obtained due the absence of similar m-patterns in the time series.
Table A3. Same as Table A1, but for m = 4 . “-” means that an undefined value is obtained due the absence of similar m-patterns in the time series.
β value−1−0.8−0.6−0.4−0.200.20.40.60.811.21.41.61.82
M S a m p E n ----------2.952.712.191.591.130.64
M F u z z y E n 3.093.022.953.083.013.263.173.132.992.752.502.131.701.170.780.45
M F u z z y E n c 3.153.153.143.153.193.143.133.082.952.792.512.111.701.160.760.43
M F u z z y E n a 2.732.732.842.772.812.922.842.872.722.552.352.051.651.190.780.47
M F u z z y E n c a 3.083.113.123.113.113.123.093.042.932.762.512.141.671.190.760.44
R S a m p E n ----------1.090.490.260.280.260.22
R F u z z y E n 0.610.390.710.590.580.800.630.580.290.400.220.180.160.210.190.15
R F u z z y E n c 0.210.160.160.150.120.140.120.110.080.060.090.140.160.230.190.15
R F u z z y E n a 0.270.190.290.280.290.270.290.260.240.190.090.100.100.100.110.08
R F u z z y E n c a 0.090.10.10.080.070.070.050.060.040.030.030.060.110.100.120.09
( R S a m p E n R F u z z y E n ) / R F u z z y E n ----------4.031.750.640.310.380.49
( R S a m p E n R F u z z y E n c ) / R F u z z y E n c ----------11.292.600.590.220.400.44
( R S a m p E n R F u z z y E n a ) / R F u z z y E n a ----------11.064.131.541.931.441.70
( R S a m p E n R F u z z y E n c a ) / R F u z z y E n c a ----------34.586.911.331.921.271.54
( R F u z z y E n R F u z z y E n c ) / R F u z z y E n c 1.891.373.382.973.944.774.264.102.766.281.440.310.030.070.010.03
( R F u z z y E n R F u z z y E n a ) / R F u z z y E n a 1.281.031.491.1312.021.161.260.231.071.400.860.551.240.770.82
( R F u z z y E n R F u z z y E n c a ) / R F u z z y E n c a 5.52.796.166.226.8910.1211.258.056.2512.676.071.870.421.230.640.71

References

  1. Pincus, S.M. Approximate entropy as a measure of system complexity. Proc. Natl. Acad. Sci. USA 1991, 88, 2297–2301. [Google Scholar] [CrossRef] [PubMed]
  2. Richman, J.S.; Moorman, J.R. Physiological time-series analysis using approximate entropy and sample entropy. Am. J. Physiol.-Heart Circ. Physiol. 2000, 278, H2039–H2049. [Google Scholar] [CrossRef] [PubMed]
  3. Chen, W.; Zhuang, J.; Yu, W.; Wang, Z. Measuring complexity using FuzzyEn, ApEn, and SampEn. Med. Eng. Phys. 2009, 31, 61–68. [Google Scholar] [CrossRef] [PubMed]
  4. Hu, J. An approach to EEG-based gender recognition using entropy measurement methods. Knowl.-Based Syst. 2018, 140, 134–141. [Google Scholar] [CrossRef]
  5. Tibdewal, M.N.; Dey, H.R.; Mahadevappa, M.; Ray, A.; Malokar, M. Multiple entropies performance measure for detection and localization of multi-channel epileptic EEG. Biomed. Signal Process. Control 2017, 38, 158–167. [Google Scholar] [CrossRef]
  6. Hu, J.; Wang, P. Noise robustness analysis of performance for EEG-based driver fatigue detection using different entropy feature sets. Entropy 2017, 19, 385. [Google Scholar]
  7. Liu, C.; Li, K.; Zhao, L.; Liu, F.; Zheng, D.; Liu, C.; Liu, S. Analysis of heart rate variability using Fuzzy measure entropy. Comput. Biol. Med. 2013, 43, 100–108. [Google Scholar]
  8. Girault, J.-M. Recurrence and symmetry of time series: Application to transition detection. Chaos Solitons Fractals 2015, 77, 11–28. [Google Scholar] [CrossRef] [Green Version]
  9. Eckmann, J.P.; Oliffson Kamphorts, S.; Ruelle, D. Recurrence plots of dynamical systems. Europhys. Lett. 1987, 4, 973–977. [Google Scholar] [CrossRef]
  10. Tarnopolski, M. On the relationship between the Hurst exponent, the ratio of the mean square successive difference to the variance, and the number of turning points. Phys. A 2016, 461, 662–673. [Google Scholar] [CrossRef]
  11. Voicu, I.; Menigot, S.; Kouamé, D.; Girault, J.-M. New estimators and guidelines for better use of fetal heart rate estimators with Doppler ultrasound devices. Comput. Math. Methods Med. 2014, 2014, 784862. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Fang, Y.; Zhou, D.; Li, K.; Liu, H. Interface Prostheses With Classifier-Feedback-Based User Training. IEEE Trans. Biomed. Eng. 2017, 64, 2575–2583. [Google Scholar] [CrossRef] [PubMed]
  13. Zhou, D.; Fang, Y.; Botzheim, J.; Kubota, N.; Liu, H. Bacterial memetic algorithm based feature selection for surface EMG based hand motion recognition in long-term use. In 2016 IEEE Symposium Series on Computational Intelligence (SSCI); IEEE: Piscataway Township, NJ, USA, 2016; pp. 1–7. [Google Scholar]
  14. Humeau-Heurtier, A.; Mahé, G.; Durand, S.; Abraham, P. Multiscale entropy study of medical laser speckle contrast images. IEEE Trans. Biomed. Eng. 2013, 60, 872–879. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Membership functions μ p ( d , r ) = exp ( ( d / r ) p ) with r = 0.1 . Gaussian function (blue) with p = 2 ; rectangular function (magenta) with p = , for d 0 .
Figure 1. Membership functions μ p ( d , r ) = exp ( ( d / r ) p ) with r = 0.1 . Gaussian function (blue) with p = 2 ; rectangular function (magenta) with p = , for d 0 .
Entropy 20 00287 g001
Figure 2. Stochastic time series where 2-patterns are pointed out. Each number corresponds to the place of the corresponding segment. No-centered 2-patterns are considered. The two 2-patterns ‘1’ and ‘15’ (black bullets) have the same mean value; they are similar. The total number of similar 2-patterns is 12: (‘1’,‘15’), (‘5’,‘21’), (‘7’,‘19’), (‘8’,‘20’), (‘13’,‘24’), (‘14’,‘25’).
Figure 2. Stochastic time series where 2-patterns are pointed out. Each number corresponds to the place of the corresponding segment. No-centered 2-patterns are considered. The two 2-patterns ‘1’ and ‘15’ (black bullets) have the same mean value; they are similar. The total number of similar 2-patterns is 12: (‘1’,‘15’), (‘5’,‘21’), (‘7’,‘19’), (‘8’,‘20’), (‘13’,‘24’), (‘14’,‘25’).
Entropy 20 00287 g002
Figure 3. Stochastic time series with different types of 2-patterns. (a) Centered 2-patterns are considered. Centered 2-patterns similar to ‘1’ are represented with magenta bullets; there are six patterns similar to ‘1’. The total number of centered similar 2-patterns is 25: (‘1’,‘9’,‘13’,‘15’,‘17’,‘24’), (‘2’,‘14’,‘25’), (‘3’,‘8’,‘20’), (‘4’,‘23’), (‘5’,‘7’,‘10’,‘19’,‘21’), (‘11’,‘18’), (‘12’,‘16’), (‘22’,‘26’). The total number of similar centered 2-patterns is much larger than that of no-centered 2-patterns. (b) Regarding the 2-pattern (‘1’), four kinds of 2-patterns can be obtained: 2-patterns with translation (‘T’) in black (‘1’,‘15’), 2-patterns with vertical reflection (‘R’) in red (‘7’, ’19’), 2-patterns with inversion (‘I’) in green (‘13’,‘24’), 2-patterns with glide reflection (‘G’) in blue (‘5’,‘21’). By considering all 2-patterns ranging from ‘1’–‘27’, the mean total number of symmetrical 2-patterns is N s y m = 92 with N s y m T = 12 , N s y m R = 30 , N s y m I = 24 , N s y m G = 26 . (c) Regarding the 2-pattern (‘1’), two kinds of centered 2-patterns can be obtained: 2-patterns (‘T’,‘I’) in black (‘1’,‘9’,‘13’,‘15’,‘17’,‘24’), 2-patterns (‘R’,‘G’) in blue (‘5’,‘7’,‘10’,‘19’,‘21’). By considering all 2-patterns ranging from ‘1’–‘27’, the mean total number of symmetrical 2-patterns is N c s y m = 312 with N c s y m T = 86 , N c s y m R = 70 , N c s y m I = 86 and N c s y m G = 70 .
Figure 3. Stochastic time series with different types of 2-patterns. (a) Centered 2-patterns are considered. Centered 2-patterns similar to ‘1’ are represented with magenta bullets; there are six patterns similar to ‘1’. The total number of centered similar 2-patterns is 25: (‘1’,‘9’,‘13’,‘15’,‘17’,‘24’), (‘2’,‘14’,‘25’), (‘3’,‘8’,‘20’), (‘4’,‘23’), (‘5’,‘7’,‘10’,‘19’,‘21’), (‘11’,‘18’), (‘12’,‘16’), (‘22’,‘26’). The total number of similar centered 2-patterns is much larger than that of no-centered 2-patterns. (b) Regarding the 2-pattern (‘1’), four kinds of 2-patterns can be obtained: 2-patterns with translation (‘T’) in black (‘1’,‘15’), 2-patterns with vertical reflection (‘R’) in red (‘7’, ’19’), 2-patterns with inversion (‘I’) in green (‘13’,‘24’), 2-patterns with glide reflection (‘G’) in blue (‘5’,‘21’). By considering all 2-patterns ranging from ‘1’–‘27’, the mean total number of symmetrical 2-patterns is N s y m = 92 with N s y m T = 12 , N s y m R = 30 , N s y m I = 24 , N s y m G = 26 . (c) Regarding the 2-pattern (‘1’), two kinds of centered 2-patterns can be obtained: 2-patterns (‘T’,‘I’) in black (‘1’,‘9’,‘13’,‘15’,‘17’,‘24’), 2-patterns (‘R’,‘G’) in blue (‘5’,‘7’,‘10’,‘19’,‘21’). By considering all 2-patterns ranging from ‘1’–‘27’, the mean total number of symmetrical 2-patterns is N c s y m = 312 with N c s y m T = 86 , N c s y m R = 70 , N c s y m I = 86 and N c s y m G = 70 .
Entropy 20 00287 g003
Figure 4. Relative percentile ranges derived from Table A1, Table A2 and Table A3 reported in the Appendix. (a) For β = 1 , relative percentile range values obtained for different m-values: for the centered fuzzy entropy compared to the fuzzy entropy ( ( R F R F c ) / R F c ), for the averaged fuzzy entropy compared to the fuzzy entropy ( ( R F R F a ) / R F a ) and for the centered and averaged fuzzy entropy compared to the fuzzy entropy ( ( R F R F c a ) / R F c a ); (bd) similar to (a), but for β = 0 , β = 1 and β = 2 , respectively.
Figure 4. Relative percentile ranges derived from Table A1, Table A2 and Table A3 reported in the Appendix. (a) For β = 1 , relative percentile range values obtained for different m-values: for the centered fuzzy entropy compared to the fuzzy entropy ( ( R F R F c ) / R F c ), for the averaged fuzzy entropy compared to the fuzzy entropy ( ( R F R F a ) / R F a ) and for the centered and averaged fuzzy entropy compared to the fuzzy entropy ( ( R F R F c a ) / R F c a ); (bd) similar to (a), but for β = 0 , β = 1 and β = 2 , respectively.
Entropy 20 00287 g004
Figure 5. Centered and averaged fuzzy entropy ( F u z z y E n c a ) and standard fuzzy entropy ( F u z z y E n ) for normal (N) in blue and pathological fetuses (P) in green with m = 2 . The results for three data lengths are shown. means statistically significant between the two groups.
Figure 5. Centered and averaged fuzzy entropy ( F u z z y E n c a ) and standard fuzzy entropy ( F u z z y E n ) for normal (N) in blue and pathological fetuses (P) in green with m = 2 . The results for three data lengths are shown. means statistically significant between the two groups.
Entropy 20 00287 g005
Table 1. For the calculation of F u z z y E n and F u z z y E n c a , the median of the mean number M N of similar 3-patterns and the median of the mean number of centered and averaged M N c a of similar 3-patterns obtained from 1 / f β noises ( N = 5000 samples) with β ranging from 1 to 2. M N c a = ( M N c a T + M N c a R + M N c a I + M N c a G ) , where M N c a k is the median of the number of centered symmetric similar 3-patterns obtained in the calculation of F u z z y E n c a , k = { T , R , I , G } . For the computation, m = 3 and r = 0.1 × standard deviation of the time series.
Table 1. For the calculation of F u z z y E n and F u z z y E n c a , the median of the mean number M N of similar 3-patterns and the median of the mean number of centered and averaged M N c a of similar 3-patterns obtained from 1 / f β noises ( N = 5000 samples) with β ranging from 1 to 2. M N c a = ( M N c a T + M N c a R + M N c a I + M N c a G ) , where M N c a k is the median of the number of centered symmetric similar 3-patterns obtained in the calculation of F u z z y E n c a , k = { T , R , I , G } . For the computation, m = 3 and r = 0.1 × standard deviation of the time series.
β −1−0.8−0.6−0.4−0.200.20.40.60.81.01.21.41.61.82.0
M N 0.730.690.650.620.60.630.630.670.801.071.763.649.2426.4771.68162.38
M N c a 16.7117.0317.4818.1319.1820.7523.1927.3535.7153.1593.73206.09540.001580.674317.409277.86

Share and Cite

MDPI and ACS Style

Girault, J.-M.; Humeau-Heurtier, A. Centered and Averaged Fuzzy Entropy to Improve Fuzzy Entropy Precision. Entropy 2018, 20, 287. https://0-doi-org.brum.beds.ac.uk/10.3390/e20040287

AMA Style

Girault J-M, Humeau-Heurtier A. Centered and Averaged Fuzzy Entropy to Improve Fuzzy Entropy Precision. Entropy. 2018; 20(4):287. https://0-doi-org.brum.beds.ac.uk/10.3390/e20040287

Chicago/Turabian Style

Girault, Jean-Marc, and Anne Humeau-Heurtier. 2018. "Centered and Averaged Fuzzy Entropy to Improve Fuzzy Entropy Precision" Entropy 20, no. 4: 287. https://0-doi-org.brum.beds.ac.uk/10.3390/e20040287

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