# Detecting position dependent tremor with the Empirical mode decomposition

- André Lee
^{1}Email author and - Eckart Altenmüller
^{1}

**2**:3

https://doi.org/10.1186/s40734-014-0014-z

© Lee and Altenmüller; licensee BioMed Central. 2015

**Received: **30 June 2014

**Accepted: **17 December 2014

**Published: **16 February 2015

## Abstract

### Background

Primary bowing tremor (PBT) occurs in violinists in the right bowing-arm and is a highly nonlinear and non-stationary signal. However, Fourier-transform based methods (FFT) make the a priori assumption of linearity and stationarity. We present an interesting case of a violinist with PBT and apply a novel method for nonlinear and non-stationary signals for tremor analysis: the empirical mode decomposition (EMD). We compare the results of FFT and EMD analyses.

### Methods

Tremor was measured and quantified in a 50-year-old professional violinist with an accelerometer. Data were analyzed using the EMD, the Hilbert transform, the Hilbert spectrum and the marginal Hilbert spectrum. Findings are compared to the FFT-spectrum and FFT-spectrogram.

### Results

We could show that the EMD yields intrinsic mode functions, which represent the tremor and IMFs, which are associated with voluntary movement. The instantaneous frequency and amplitude are obtained. In contrast the low time frequency resolution and the artifacts of voluntary movements are seen in the FFT results.

### Conclusions

PBT may present itself as a highly non-stationary and nonlinear phenomenon, which can be accurately analyzed with the EMD, since it gives the instantaneous amplitude and frequency and can identify voluntary from involuntary (tremor) movement.

## Keywords

## Background

Tremor is defined as an involuntary rhythmical oscillation of a body part [1]. Particularly pathological tremors are time-varying [2] and highly nonlinear and non-stationary in nature [1,3,4]. Task-specific tremors (TST) are pathological tremors that occur predominantly during certain tasks [1]. Primary bowing tremor (PBT) [5] occurs unilaterally in the right arm of bowed string-instrument players while playing the instrument. This is a highly disabling condition and may threaten the musician’s professional career. We describe a violinist in whom PBT occurred when he played a fast movement from the tip of the bow to the frog (the part of the bow held by the violinist), brought the movement to a sudden stop and tried to maintain the hand in a stable position (Additional file 1: Video). In that position tremor appeared and decreased in amplitude over the next 10–20 seconds, giving a highly non-stationary and non-linear signal.

The disadvantage of applying the Fourier transform (FFT) to these kinds of signals is the *a-priori* assumption of a linear and periodic or stationary signal i.e. a sine or cosine of constant amplitude and frequency spanning the whole signal. The FFT gives reliable results therefore only in case of linear and stationary signals [6]. However, periodicity cannot be assumed for tremors, since frequency not only changes with the waves in a dispersive system (interwave modulation) but likewise within one oscillation cycle or wave (intrawave modulation) [7]. Therefore the wave-profile cannot be considered a sine or cosine function. Furthermore FFT has a limited time-frequency resolution. Thus potentially meaningful local (in a temporal sense) oscillations may not be detected. Finally the FFT does not distinguish between noise (e.g. voluntary movement in this study) and the actual signal (e.g. tremor, as in this study), making the result less reliable. In recent studies a new method that takes into account the nonlinearity and non-stationarity of signals has been introduced [8] and has been applied in tremor research [2,9,10]. This approach combines two tools: Empirical mode decomposition (EMD) and the Hilbert transform.

*a-posteriori*approach, which does not need a-priori assumptions with regard to the signal [7,8]. The IMF can be regarded as a more general counterpart of the simple harmonic functions obtained by the FFT that may have a variable amplitude and frequency [7]. It has been shown that IMFs can be used to identify distinct frequency bands associated with physical or physiological phenomena, for example particular types of tremor, [9,11,12] however it has been noted that mainly tremor data from gyroscopes was analyzed [9,11,12], but not from accelerometers. Furthermore, the EMD may separate noise from the actual data [8].

The Hilbert transform yields the Hilbert spectrum, i.e. the instantaneous amplitude and frequency. It is thus a measure of the contribution of each frequency over time, from which the marginal Hilbert spectrum (MHS), a measure of the total amplitude contribution for each frequency value [7,8], can be derived. The advantage of this method over the FFT is that it is a windowing independent time-frequency representation with a high time-frequency resolution.

The aim of this paper was thus threefold: 1. to present an interesting case of a task-specific tremor in a violinist; 2. to investigate, whether the EMD and Hilbert transform can identify the tremor signal from the highly non-stationary and nonlinear signal obtained from the accelerometer and to separate artifacts from voluntary movements inherent to the task from the involuntary tremor (see Methods) 3. to demonstrate the advantages of the EMD over the FFT.

## Methods

The study was approved by the ethics committee of the Hanover medical school and written informed consent was obtained by the participant. Tremor was measured in a 50-year-old professional violinist who had played in a prestigious orchestra for more than 20 year. He was asked to play a fast up-bow-movement, which triggered tremor at the end of that movement when trying to hold the hand in a stable position. He then took back the bow from the frog of the bow to the tip to prepare another fast up-bow movement. The fast up-bow movement and retaking the bow are referred to as voluntary movement. This was repeated for five times. Measurement occurred with a 3D accelerometer (biovision, Wehrheim, Germany, 8×8×11 mm; 4 gram; DC–500 Hz; max 50 g), which was attached to the metacarpo-phalangeal (MCP) joint of the index finger of the right hand. Data were bandpass filtered using a 4th degree butterworth-filter (cutoff 1–50 Hz), applied back and forth to compensate for phase shift. With the accelerometer signal onset and end of the fast up-bow movement as well as retaking the bow could be identified.

### Empirical mode decomposition

The EMD consists of a sifting process of the original signal [X(t)], in which the intrinsic mode functions (IMF) are obtained. The EMD algorithm as described by Huang et al. is as follows [7,8].

_{1}) and subtract it from X(t):

_{1}(t) does not fulfill the criteria for an IMF, this procedure is repeated n-times until h

_{(1n)}(t) fulfills the criteria and is thus defined as the first IMF [c

_{1}(t)]:

_{1}(t) is subtracted from X(t):

_{1}(t) is called the residue which is substituted for X(t) in formula (1) and the first steps are repeated m-times until a residual r

_{m}is reached that is a monotonic function of which no more IMFs can be extracted [7,8]. Thus, the original signal can be obtained by

The frequency ranges of the IMF are ordered in such a way that IMF1 contains the highest and IMF_{m} the lowest frequencies. Being a data-driven approach, the frequency ranges depend on the original signal.

### Hilbert transform, Hilbert spectrum and marginal Hilbert spectrum

a_{j}(t) is the instantaneous amplitude, ω_{j}(t) the instantaneous frequency.

A detailed description can be found in Huang et al. [8] and Huang [7].

We expected to find an IMF, that would display the onset of tremor after stopping the up-bow movement with a decreasing amplitude without the artifact of the fast up-bow movement itself.

EMD was performed in Matlab using the EMD package by Rilling et al. [http://perso.ens-lyon.fr/patrick.flandrin/emd.html], applying the default stopping criterion [13].

### FFT for comparison

To compare the results of the EMD and Hilbert transform with the FFT we performed an FFT with a window of 5096 data points and an overlap of 1024 data points.

## Results

### Tremor detection

#### Empirical mode decomposition

Figure 1 displays the EMD with 5 IMFs and the original signal for the first two fast up-bow movements. In the original signal the artifact voluntary movements (i.e. at the onset of the fast up-bow movement and when taking back the bow) are visible, showing the non-stationarity of the signal. IMF3 best represents the tremor signal, however, part of signal becomes apparent in IMF2, known as mode-mixing. No tremor is detected in IMF1. Low-frequency tremor can be seen in IMFs 4 and 5. To investigate the mode mixing we chose the combination of IMF2 and IMF3 (IMF2 + 3) for further evaluation (see below).

#### Hilbert spectrum and marginal Hilbert spectrum

### Comparison with FFT

## Discussion

We describe a patient with an unusual form of task specific bowing tremor that appeared only when a sudden up-bow movement was brought to a sudden stop at the frog of the bow and held there as steadily as possible. The peak frequency was at 4.7 Hz (Figure 4) and thus was in the range of PBT described before [14,15]. Since we had a highly non-stationary signal with artifacts from voluntary movements (i.e. the fast up-bow movement as well as the movement for returning the bow to the tip) as seen in Figures 1 and 2, we applied a novel method for analyzing the signal, the EMD and Hilbert transform that do not require stationarity and linearity as a prerequisite. We could show that the EMD and the Hilbert transform are able to correctly identify the tremor signal. As expected the tremor signal was mainly contained in one IMF, namely IMF3. One interesting aspect of our analysis was the finding of a mode mixing, which has been described by Huang [7]. It describes the finding of part of the tremor signal being distributed between two (or more) IMFs, in our case between IMF2 and IMF3. According to formula (4), the original signal can be obtained by adding all IMF and the final residual r_{m}, which implies the possibility of adding two (or more) IMF. We therefore added IMF2 + IMF3. The Hilbert spectrum of this combined IMF2 + 3 (Figure 3, bottom) corroborates this finding, since the tremor amplitude of very precisely represents the course of the tremor observed during clinical examination at the instrument (Additional file 1: Video): It is highest after the up-bow movement and diminishes over time. Virtually no tremor is present when retaking the bow.

Importantly, our expectation that the EMD can reduce artifacts, e.g. from voluntary movement as has been shown before [9], was confirmed. As shown in Figure 2, the low frequency artifacts of the up-bow movement and of retaking the bow are removed in the IMFs. The up-bow movement lasted for 0.5 s, giving a low frequency artifact of 2 Hz. Even lower frequencies were to be expected from retaking the bow. Peaks at around 2 Hz and below are visible in IMF4 and IMF5. We therefore interpret the low frequency signal of these IMFs (Figure 1) to be caused by the voluntary movements. This is important, because it suggests that discontinuities in tremor signals and voluntary movements have a small effect on the precision of the method.

Given the algorithm of the EMD, fast voluntary movements in music making with a high frequency (e.g. tremolo) would have been represented in IMF 1. Our paradigm did not include those kinds of movements and thus, as expected, almost no signal was detected here.

The comparison with results of an FFT analysis revealed similar peak frequencies for both methods. However, the limited time-frequency resolution of the FFT becomes apparent when comparing Figures 3 and 5. The instantaneous frequency and amplitude allows a more precise course of the tremor, whereas for FFT-based methods, a compromise between time and frequency resolution has to be made. In the FFT spectrogram (Figure 6) the low-frequency artifact of the voluntary movement is overestimated and taken into account in the FFT power spectrum (Figure 5). In the IMFs, however, the signal is separated into frequency ranges that cover both, tremor (IMF2 + 3) and lower frequency voluntary movement (IMF 4, IMF5, see Figure 4). We therefore conclude that IMF2 + 3 represents the frequency range and thus the intrawave frequency modulation of task-specific tremor than the frequency range covered by the FFT power spectrum (Figure 5). This is important, because it has been shown that the wave-profile deformation seen in the FFT-power spectrum that are usually interpreted as harmonic distortions, are more likely due to the intrawave frequency modulation [8].

Since the EMD has been applied so far only for tremor data derived from a gyroscope [2,9,10] we finally conclude that it is applicable for accelerometer data, as well.

## Conclusion

We present an interesting and unusual case of a patient with a highly non-stationary and nonlinear PBT. We could show that the EMD can accurately detect and analyze the tremor, identify voluntary movement and is applicable for data obtained with an accelerometer. In comparison to the FFT, the EMD does not make a priori assumptions on the data, yields the instantaneous frequency and amplitude and is thus a more precise tool for analyzing tremor data.

## Declarations

### Acknowledgements

Funding for AL and EA: Hannover University of Music, Drama and Media.

## Authors’ Affiliations

## References

- Deuschl G, Bain P, Brin M. Consensus statement of the Movement Disorder Society on Tremor. Ad Hoc Scientific Committee. Mov Disord. 1998;3(13 Suppl):2–23.Google Scholar
- De Lima ER, Andrade AO, Pons JL, Kyberd P, Nasuto SJ. Empirical mode decomposition: a novel technique for the study of tremor time series. Med Biol Eng Comput. 2006;44:569–82.View ArticlePubMedGoogle Scholar
- Gantert C, Honerkamp J, Timmer J. Analyzing the dynamics of hand tremor time series. Biol Cybern. 1992;66:479–84.View ArticlePubMedGoogle Scholar
- Timmer J, Haussler S, Lauk M, Lucking C-H. Pathological tremors: deterministic chaos or nonlinear stochastic oscillators? Chaos. 2000;10:278–88.View ArticlePubMedGoogle Scholar
- Lee A, Altenmüller E. Primary task-specific tremor: an entity of its Own? Med Probl Perform Art. 2012;27:224–6.PubMedGoogle Scholar
- Oppenheim AV, Schafer RW. Discrete-Time Signal Processing. Prentice Hall: Englewood Cliffs, N.J; 1989.Google Scholar
- Huang NE. Introduction to the Hilbert–Huang Transform and Its Related Mathematical Problems. In: Hilbert-Huang Transform and its Applications. Singapore; Hackensack, NJ; London: World Scientific; 2005. p. 1–26.Google Scholar
- Huang NE, Shen Z, Long SR, Wu MC, Shih HH, Zheng Q, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc R Soc A: Math Phys Eng Sci. 1998;454:903–95.View ArticleGoogle Scholar
- Rocon E, Pons JL, Andrade AO, Nasuto SJ. Application of EMD as a novel technique for the study of tremor time series. Conf Proc IEEE Eng Med Biol Soc. 2006, Suppl:6533–6536.Google Scholar
- Gallego JA, Rocon E, Koutsou AD, Pons JL. Analysis of kinematic data in pathological tremor with the Hilbert-Huang transform. IEEE. 2011:80–83Google Scholar
- Silchenko AN, Adamchic I, Pawelczyk N, Hauptmann C, Maarouf M, Sturm V, et al. Data-driven approach to the estimation of connectivity and time delays in the coupling of interacting neuronal subsystems. J Neurosci Methods. 2010;191:32–44.View ArticlePubMedGoogle Scholar
- Li K, Hogrel J-Y, Duchêne J, Hewson DJ. Analysis of fatigue and tremor during sustained maximal grip contractions using Hilbert-Huang Transformation. Med Eng Phys. 2011;34:832–40.View ArticlePubMedGoogle Scholar
- Hogan N. Adaptive control of mechanical impedance by coactivation of antagonist muscles. IEEE Trans Automatic Control. 1984;29:681–90.View ArticleGoogle Scholar
- Lee A, Chadde M, Altenmüller E, Schoonderwaldt E. Characteristics of task-specific tremor in string instrument players. Tremor Other Hyperkinet Mov (NY). 2014;4:198.Google Scholar
- Lee A, Tominaga K, Furuya S, Miyazaki F, Altenmüller E. Coherence of coactivation and acceleration in task-specific primary bowing tremor. J Neural Transm. 2014;121:739–42.View ArticlePubMedGoogle Scholar

## Copyright

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.