analyzed time series can be used to present
evidences of nonlinear components in the signal
(Kantz & Schreiber, 1997).
Nonlinear approaches to calculate biological
time series complexity can be obtained using
nonlinear entropy measures such as Approximate
Entropy (Pincus, 1991); Lempel-Ziv complexity
(Doganaksoy & Göloglu, 2009) and Sample Entropy
(Richman & Moorman, 2000).
The Sample Entropy is a metric based on
Approximate Entropy that has been already applied
to evaluate the complexity in cardiovascular time
series (Richman & Moorman, 2000).
This work presents a nonlinear analysis of FHR,
which is the most significant signal monitored in a
CTG, based on Sample Entropy (SampEn)
behaviour in time.
2 MATERIALS AND METHODS
2.1 Sample Entropy
The SampEn of a time series is the negative natural
logarithm of the probability that two similar
sequences of m points remain similar at the next
point. The polarization of self-matches is not
considered. Low values of SampEn indicate lower
complexity and more regularity in the time series.
For a better understanding of these definitions as
a measure of system complexity, the mathematical
foundations of its calculation is provided.
Given S
i
as a time series with N samples S
i
= {S
1
,
S
2, ….,
S
N
} the first step to calculate SampEn(S
i
,m,r),
is the determination of two input parameters m and
r, where m is the length of a subset of S
i
and r is the
similarity criteria. Consider also that p
m
(i) is the
subsequence or pattern of S
N
beginning in sample i
and with m samples of length.
Consider two patterns p
m
(i) and p
m
(j), beginning,
at index i and j respectively. These patterns can be
considered similar if the scalar distance between
them, i.e., the module of the difference between
sample pairs is less than r. (1).
rSS
kjki
<−
++
(1)
for 0 < k < m.
Consider P
m
as the set of all patterns from S
N
with length m. The relation C
i,m
(r)
can now be
defined (2).
1
)(
)(
,
,
+−
=
mN
rn
rC
mi
mi
(2)
where n
i,m
(r) is the number of patterns similar to
p
m
(i) in P
m
. The parameter C
m
(r) must be calculated
as the average of all C
i,m
(r) for the entire P
m
set.
Finally, SampEn(S
i
,m,r) can be found (3).
)(
)(
ln),,(
1
rC
rC
rmSSampEn
m
m
i
+
=
(3)
2.2 Development Environment
and CTG Database
The development environment was based on the
Matlab software, version 7.6.0.324 R2008a
(Mathworks, 2009).
The results were obtained using a database from
Trium GmBH, a project partner from Munich,
Germany. This database was labelled as the CTG-I
and contains 22 intrapartum exams (during labour
and delivery), all classified by medical staff.
This database was chosen because of the high
level of dynamics found during labour, especially
when the influence of uterine contractions can be
found in FHR signals.
2.3 Entropy Calculation Parameters
Sample entropy can be calculated for the whole FHR
signal, providing a long term index or, alternatively,
windows of samples can be used for short term
nonlinear characteristics evaluation.
Data were submitted to many tests with different
Δt
e
window sizes and also changing m and r
parameters.
In this paper, the FHR signal entropy
calculations consider a subset of 1200
samples,
which represents a 5 minute-long window. The aim
is the monitoring of the signal entropy based on its
time evolution. For example, the entropy behaviour
during pathological FHR events, such as prolonged
decelerations, could be a predictive tool for
electronic fetal monitoring.
The SampEn input parameters used were m=5
and r=0.2σ[FHR(t)], where FHR(t) is the FHR
signal and σ[FHR(t)] is the standard deviation of the
time series (Kaplan & Staffin, 2008).
3 RESULTS AND DISCUSSION
In this section, the general results obtained with the
CTG-I database are presented with some illustrative
examples of high and low values of SampEn and the
correspondent visualization of FHR trace in time.
NON-LINEAR ANALYSIS OF FETAL HEART RATE IN CARDIOTOCOGRAPHY USING SAMPLE ENTROPY
245