Sample Entropy

Sample entropy (SampEn) is a modification of approximate entropy (ApEn), used for assessing the complexity of physiological time-series signals, diagnosing diseased states. SampEn has two advantages over ApEn: data length independence and a relatively trouble-free implementation. Also, there is a small computational difference: In ApEn, the comparison between the template vector (see below) and the rest of the vectors also includes comparison with itself. This guarantees that probabilities $C_'^(r)$ are never zero. Consequently, it is always possible to take a logarithm of probabilities. Because template comparisons with itself lower ApEn values, the signals are interpreted to be more regular than they actually are. These self-matches are not included in SampEn. However, since SampEn makes direct use of the correlation integrals, it is not a real measure of information but an approximation. The foundations and differences with ApEn, as well as a step-by-step tutorial for its application is available at.

There is a multiscale version of SampEn as well, suggested by Costa and others. SampEn can be used in biomedical and biomechanical research, for example to evaluate postural control.

Definition

Like approximate entropy (ApEn), Sample entropy (SampEn) is a measure of complexity. But it does not include self-similar patterns as ApEn does. For a given embedding dimension $m$, tolerance $r$ and number of data points $N$, SampEn is the negative natural logarithm of the probability that if two sets of simultaneous data points of length $m$ have distance $< r$ then two sets of simultaneous data points of length $m+1$ also have distance $< r$. And we represent it by $SampEn(m,r,N)$ (or by $SampEn(m,r,\tau,N)$ including sampling time $\tau$).

Now assume we have a time-series data set of length $N =$ with a constant time interval $\tau$. We define a template vector of length $m$, such that $X_m (i)=$ and the distance function d _m(i),X_m(j) (i≠j) is to be the Chebyshev distance (but it could be any distance function, including Euclidean distance). We define the sample entropy to be

:$SampEn=-\ln$

Where

$A$ = number of template vector pairs having d _(i),X_(j)< r

$B$ = number of template vector pairs having d _m(i),X_m(j)< r

It is clear from the definition that $A$ will always have a value smaller or equal to $B$. Therefore, $SampEn(m,r,\tau)$ will be always either be zero or positive value. A smaller value of $SampEn$ also indicates more self-similarity in data set or less noise.

Generally we take the value of $m$ to be $2$ and the value of $r$ to be $0.2 \times std$.
Where std stands for standard deviation which should be taken over a very large dataset. For instance, the r value of 6 ms is appropriate for sample entropy calculations of heart rate intervals, since this corresponds to $0.2 \times std$ for a very large population.

Multiscale SampEn

The definition mentioned above is a special case of multi scale sampEn with $\delta=1$, where $\delta$ is called skipping parameter. In multiscale SampEn template vectors are defined with a certain interval between its elements, specified by the value of $\delta$. And modified template vector is defined as
$X_(i)=$
and sampEn can be written as
$SampEn \left ( m,r,\delta \right )=-\ln$
And we calculate $A_\delta$ and $B_\delta$ like before.

Implementation

Sample entropy can be implemented easily in many different programming languages. Below lies a vectorized example written in Python.

import numpy as np

def sampen(L, m, r):
"""Sample entropy."""
N = len(L)
B = 0.0
A = 0.0

# Split time series and save all templates of length m
xmi = np.array( [i : i + mfor i in range(N - m)"> : i + m">[i : i + mfor i in range(N - m)
xmj = np.array( [i : i + mfor i in range(N - m + 1)"> : i + m">[i : i + mfor i in range(N - m + 1)

# Save all matches minus the self-match, compute B
B = np.sum([np.sum(np.abs(xmii - xmj).max(axis=1) <= r) - 1 for xmii in xmi])

# Similar for computing A
m += 1
xm = np.array( : i + m">[i : i + mfor i in range(N - m + 1)

A = np.sum([np.sum(np.abs(xmi - xm).max(axis=1) <= r) - 1 for xmi in xm"> : i + mfor i in range(N - m + 1)"> : i + m">[i : i + mfor i in range(N - m + 1)

A = np.sum([np.sum(np.abs(xmi - xm).max(axis=1) <= r) - 1 for xmi in xm

# Return SampEn
return -np.log(A / B)

An example written in other languages can be found:

Matlab

R

See also

* Kolmogorov complexity
* Approximate entropy

References

{{reflist

Statistical signal processing
Entropy
Articles with example Python (programming language) code