Sensors 2010, 10, 1399-1422; doi:10.

3390/s100201399
OPEN ACCESS

sensors
ISSN 1424-8220 www.mdpi.com/journal/sensors Review

Neurological Tremor: Sensors, Signal Processing and Emerging Applications

Giuliana Grimaldi and Mario Manto * FNRS, Neurologie ULB-Erasme, 808 Route de Lennik, 1070 Bruxelles, Belgium; E-Mail: [email protected] * Author to whom correspondence should be addressed; E-Mail: [email protected]; Tel.: +32-2-555-39-92; Fax: +32-2-555-39-92. Received: 23 December 2009; in revised form: 22 January 2010 / Accepted: 20 February 2010 / Published: 24 February 2010

Abstract: Neurological tremor is the most common movement disorder, affecting more than 4% of elderly people. Tremor is a non linear and non stationary phenomenon, which is increasingly recognized. The issue of selection of sensors is central in the characterization of tremor. This paper reviews the state-of-the-art instrumentation and methods of signal processing for tremor occurring in humans. We describe the advantages and disadvantages of the most commonly used sensors, as well as the emerging wearable sensors being developed to assess tremor instantaneously. We discuss the current limitations and the future applications such as the integration of tremor sensors in BCIs (brain-computer interfaces) and the need for sensor fusion approaches for wearable solutions. Keywords: tremor; sensors; signal analysis; brain-computer interface (BCI)

1. Introduction Tremor is the most common movement disorder. Its incidence and prevalence increase with ageing, affecting more than 4% of the patients older than 65 years [1]. More than two-thirds of the population with upper limb tremor face serious difficulties in daily life. Although tremor is not life-threatening, it causes functional disability and social inconvenience, contaminating daily life activities, such as writing, pouring, eating, and so on [2]. Figure 1 illustrates an example of the effect of tremor in the spiral copying test (Archimedes spiral).

Sensors 2010, 10 Figure 1. Superimposition of spirals drawn on a digitized tablet. Left panel: control subject, right panel: patient with tremor - the spirals are irregular with swerves.

1400

Tremor can be defined as a rhythmic shaking of a body part [3]. It occurs in healthy individuals as the so-called physiological tremor, which is composed of two distinct oscillations (mechanical reflex and central neurogenic) superimposed on a background of irregular fluctuations in muscle forces and displacements [4,5]. The mechanical reflex component is governed by the inertial and elastic properties of the body, whereas the central neurogenic component is associated with the modulation of motor unit activity under the control of active generators in the brain. Participating motor units usually discharge at about 813 Hz [6]. The central component is inertia-insensitive. Physiological tremor may be enhanced by anxiety, stress, fatigue and medications. In neurological patients, tremor is clinically described as rest, postural and/or kinetic tremor [7] according to its mode of clinical presentation. Rest tremor appears during resting. Postural tremor is triggered by maintenance of a posture or a position against gravity. Kinetic tremor is evoked by a voluntary movement and is maximal near the target. In a clinical setting, tremor is characterized by its dominant frequency and its power spectral density. Rest tremor frequency is typically in the 36 Hz frequency range [7,8] and may increase with mental stress (i.e., counting backwards) or contralateral voluntary motion [8]. The most common cause of rest tremor is idiopathic Parkinsons Disease. The frequency of postural tremor is usually between 4 and 12 Hz. Many disorders are associated with postural tremor in upper limbs. Essential tremor is the commonest cause [9]. Kinetic tremor has a frequency between 2 and 7 Hz in the large majority of cases [7]. Tremor occurring in cerebellar diseases is the typical example of a kinetic tremor. Cerebellar disorders are either sporadic or genetic [10]. The cerebellum is considered to be a major site for tremorgenesis. Detailed lesionmapping studies have now delineated the regions of the cerebellum that cause deficits in limbs or trunk [11]. Concerning imaging, brain MRI has become a routine procedure in the evaluation and follow-up of neurological patients [12]. Tremor may be associated with a single or multiple generators in the brain. Figure 2 illustrates the role of brain imaging in the identification of anatomical lesions participating in the generators of tremor.

Sensors 2010, 10 Figure 2. Brain MRI of a patient exhibiting a disabling neurological tremor in the upper limbs. Note the left inferior olivary hypertrophy (Guillain-Mollaret triangle; arrow in left panel, (A) and the lesion in putamen (arrow in right panel, (B) . Axial T2-weighted images. R: right, L: left.

1401

Although tremor can be estimated clinically, the non stationary feature and the difficulties related to a pure clinical evaluation (with inherent subjectivity) make the use of sensitive, reliable and stable sensors mandatory. Repercussions on daily living activities (ADL) can be evaluated using specific questionnaires such as the ADL-T24, which has a good inter-session reproducibility (Table 1). Clinical scales include the Tolosa-Fahn-Marin scale, as well as other specifically-designed and detailed composite scales combining clinical and functional evaluation such as the composite CNF-TES (clinical neurophysiological functional tremor evaluation scale; Table 2) which compares clinical, neurophysiological and functional results. Differential diagnoses of tremor includes dystonia (prolonged muscle contractions that lead to abnormal postures), dyskinesias (acute or chronic involuntary movements similar to tics, chorea, spasm, or myoclonus), chorea (irregular movements with a dance-like aspect), athetosis (continuous slow hyperkinesia of distal segments), tics (repetitive hyperkinetic movements usually with a head topography) and myoclonus (sudden and short movements). In essence, tremor is stereotyped and repetitive.

Sensors 2010, 10

1402

Table 1. Questionnaire of tremor-induced difficulties encountered during daily life (ADL-T24 score).
-to move a glass full of water on a table No problem 0 Slight difficulties 1 Important difficulties 2 Impossible 3 -to drink No problem 0 Slight difficulties 1 Important difficulties 2 Impossible 3 -to eat (use of forks and knives) No problem 0 Slight difficulties 1 Important difficulties 2 Impossible 3 -to shave No problem 0 Slight difficulties 1 Important difficulties 2 Impossible 3 -to write words on a sheet of paper or to sign No problem 0 Slight difficulties 1 Important difficulties 2 Impossible 3 -to read a book No problem 0 Slight difficulties 1 Important difficulties 2 Impossible 3 -to drive a car No problem 0 Slight difficulties 1 Important difficulties 2 Impossible 3 -to dress one-self No problem 0 Slight difficulties 1 Important difficulties 2 Impossible 3 Total Score:./24

Sensors 2010, 10

1403

Table 2. The composite clinical/neurophysiological/functional tremor evaluation scale (CNF-TES). C-TES (Clinical-TES) Anamnesis Assessment of disability (Activities of Daily Living scales / ADL-T24) Physical and Neurological examination Tremor evaluation - Body segments involved (head, trunk, upper limbs, lower limbs) - Enhancing/reducing effect (effect of mental calculation/contralateral contractions) - Distribution (symmetry/asymmetry) - Grade/Amplitude - Frequency Brain imaging (CT-scan - MRI - SPECT- PET) Blood studies N-TES (Neurophysiological-TES) EMG and EEG recordings (time-frequency analysis, coherence, ERS/ERD) Evaluation of pinch force Analysis of writing (digitizing tablet) F-TES (Functional-TES; see [13-14]) Mechanical counters Box and block Test 9-Hole-Peg Test Coin test Several sensors for the monitoring and analysis of tremor have been developed these last decades (Table 3). The advantages and disadvantages of techniques for tremor evaluation are summarized (see also references [8], 15-17). Table 3. Comparison of the most commonly used sensors for quantification and monitoring of tremor *.
Assessment of kinematics** EMG Surface EMG (SEMG) Needle electrodes Fine-wire electrodes Long-term recordings Force transducers and force-feedback devices (haptic devices)

Accelerometer Gravity effect influence Accuracy of frequency information Signal-to-noise ratio Electrical contacts with subjects

Gyroscope

Video

yes good

no good

no may be low variable no

no good

no good

low to high no

high no

high yes

high yes

Sensors 2010, 10 Table 3. Cont.

Size Painful Cost Easy to use Data processing required Measurement of tremor amplitude

1404

small no cheap yes yes

small no cheap yes yes

relatively large no cheap to expensive yes yes

small yes (needle EMG) variable variable yes

Large no expensive relatively difficult yes

calculation from time/acceleration

measurement of inertial angular rate

from calibrated video frames

no

from force/mass or position encoder

*Emerging sensors include wearable textiles with integrated sensors, sensors based on the fusion approach and hybrid sensors. **Includes also flexible angular sensors/goniometers and laser-based devices.

They provide information about which are complementary to the clinical evaluation and often guide therapies, hence their critical role in this field. The combination of electromyography (EMG) with kinematic sensors is widely used. Each technique will be described below. General recommendations As a prerequisite, any biomedical instrumentation, independently from its technical characteristics, requires the appropriate placement of the sensor on the body segment and a quiet recording environment free of electrical interferences. Regular calibration of the instrumentation should not be underestimated. Standardized recording procedures are essential for intra- and inter-patients comparisons [8]. Both the clinical examination and quantification of tremor should be performed by the same examiner if possible, in order to reduce the anxiety effect, which may cause noticeable differences in results [18]. Tremor recordings should be performed immediately after clinical examination. Functional tests provide useful additional information for a better understanding of patients disabilities. Indeed, functional evaluation is now part of the evaluation of neurological tremor. The need for blind or even double blind studies when possible, as well as future studies with control conditions should be underlined. 2. Sensors and Instrumentation for the Characterization of Tremor The issue of the localization of tremor generators remains central in the field of tremor research [19]. The selection of the sensor helps in improving (1) our understanding of the respective contributions of supra-spinal and spinal structures, and (2) how the state of the neuromuscular periphery influences tremor. Using sensors that do not significantly impede wrist movements (so-

Sensors 2010, 10

1405

called soft tremor sensors) and more rigid sensors (approximating the condition to isometric), it has been shown that the fixation of a wrist joint decreases the amplitude of tremor bilaterally with a relatively low impact on tremor frequency [19], suggesting an important role for the neuromuscular and spinal cord mechanisms in determining tremor amplitude. This is an example of how mechanical properties of a sensor influence results and interpretation of the mechanisms of tremor. 2.1. Electromyography (EMG) Surface EMG (SEMG), needle EMG and recordings with fine wires record electrical potentials generated by the muscle fibers [20]. For SEMG, sensors are fixed on the skin at the level of muscle(s) of interest, using either self-adhesive disposable electrodes or differential electrodes with built-in preamplifiers. Several electrodes commercially available do not require conductive gel or skin preparation. Surface electrodes are arranged approximately 2 cm apart over the muscle. The use of anatomical landmarks is recommended to avoid variations due to electrode misplacement. Regarding intra-muscular studies with conventional needles or fine wire electrodes, electrodes are implanted in a given muscle to extract the properties of motor unit firings (firing times, firing rates, recruitment threshold, cross-correlation of firing rates of concomitant motor units). Analysis of the duration of bursts of EMG activities can be useful to distinguish the various disorders (Figure 3). Patterns of muscle contractions (co-contraction/alternating) are identified with EMG recordings (Figure 4). Figure 3. Duration of EMG bursts in forearm muscles in various neurological disorders associated with tremor. Abbreviations: PD: Parkinsons disease, ET: essential tremor, PN: peripheral neuropathy, OT: orthostatic tremor. Adapted from Grimaldi and Manto, 2008.

Patterns of motoneuronal discharges might constitute motoneuronal signatures which takes into account the size of the recruited motor units, the timing of their respective firings, their firing rates and their modulation [21]. Interestingly, tremor recordings during a tapping exercise with the controlateral arm provide a quantified electrophysiological entrainment test which is a sensitive and specific procedure to distinguish psychogenic tremor from dystonic and other organic tremors [22]. Long-term EMG has proven to be a valid and reliable method for the quantification of pathological tremors. Recent works indicate that these recordings allow a rater-independent classification of parkinsonian versus essential tremor [23]. Arrays of SEMG are being adopted by many laboratories.

Sensors 2010, 10

1406

The arrays can be integrated in textiles for wearable solutions (Figure 5), with the aim of combining spatio-temporal accuracy with aesthetic requirements and comfort for daily use. Figure 4. Recordings of action tremor in a neurological patient. Note the oscillations in channel 2 (accelerometry, right side; accelerometer fixed on the extremity of the index finger) and the alternating pattern of EMG bursts in channels 5 (Flexor carpi radialis, right side small arrows) and 6 (Extensor carpi radialis, right side large arrows). Channel 1, channel 3 and channel 4 correspond to accelerometry and EMG recordings on the left side.

Figure 5. Illustration of a textile with a matrix of surface EMG electrodes used to record upper limb tremor. The textile is designed to fit with comfort issues. Courtesy of Smartex, Italy.

SMARTE

Sensors 2010, 10 2.2. Accelerometers

1407

Accelerometer measures acceleration along the sensitive axis of the sensor based on Newtons second law (Force = Mass x Acceleration). The loading force drives a second-order damped harmonic oscillator in a mass-spring-damper system [15]. There are three main categories of accelerometers: piezoelectric, piezoresistive and capacitive types. Impedance, signal level, effects of gravitational component and cost are some important parameters related to this technique [15]. Table 4 compares the three types of accelerometers. The natural high-pass filtering of piezoresistive accelerometers reduces significantly slow drifts [24]. The need for a second integration to obtain displacement values and the need for regular calibration are considered as negative aspects. Another disadvantage is the fact that accelerometers measure linear acceleration, while articular motions are mainly based on rotations of joints. Data from accelerometers are composed of a combination of linear acceleration, gravity, and additive noise [16]. There is still no validated analytic model to distinguish between data due to acceleration and gravity. Low-pass filtering is commonly used [25]. Table 4. Comparison of accelerometers. Parameter Gravitational component Bandwidth Impedance Signal level Ruggedness Cost Piezoelectric No Wide High High Good High Piezoresistive Yes Low to moderate Low Low Moderate Low Capacitive Yes Wide Very high Moderate Good High

Adapted from Wong et al., 2007 [15].

Accelerometry is simple, relatively reliable and remains a convenient technique to measure frequency and amplitude of oscillations of body segments [7]. Sensors are fixed on the skin at given anatomical landmarks. An example of upper limb tremor is shown in Figure 6, recorded with a combination of accelerometry and surface EMG recordings. Various levels of support of the limbs have been proposed [26]. The feasibility of using data recorded with accelerometers to estimate the severity of symptoms and motor complications in patients with Parkinson's disease has been shown recently [27]. Accelerometers are also used for intraoperative assessment of the best position to implant electrodes for deep brain stimulation (DBS) and neurophysiological monitoring of stereotactic intervention of movement disorders [28]. Accelerometers allow the estimation of orientation of body segments, acceleration of trunk, acceleration/velocity/translations of limbs [8]. These are clinically useful parameters for the management of neurological patients. There are also applications outside human disorders, for instance in animal models of tremor [29]. The model of harmaline-induced tremor in rodents is an example.

Sensors 2010, 10 Figure 6. Example of tremor in a patient with multiple brain lesions. One single-axis accelerometer is affixed on the extremity of right index (Acc. R.) while the patient is performing the finger-to-finger test (the index fingers are maintained horizontally at a distance of about 1 cm in front of the patient). Surface EMG recordings of the flexor carpi radialis (FCR) and extensor carpi radialis (ECR) muscles are shown for both sides (L: left, R: right). Note the oscillations in the accelerometry and the regular bursts of EMG activities. Peak tremor frequency: 4.7 Hz.

1408

2.3. Gyroscopes A gyroscope can be defined as an angular velocity sensor. The angular orientation can be obtained from integration of the velocity signal [15]. Gyroscopes are based on the measurement of the Coriolis force of vibrating devices. Coriolis force is an apparent force arising in a rotating reference frame and proportional to the angular rate of rotation. There are three basic types of gyroscopes: spinning rotor gyroscope, ring laser gyroscope and vibrating mass gyroscope. The last one presents many advantages for portable applications due to its size, weight, lower power consumption and cost. An internal mass in constantly vibrating inside the sensor. Gyroscopes are being used to measure velocity and stride length, joint angle of lower limbs, angular velocity of trunk rotation, and angular displacement of trunk motions [8]. They are often selected for implementation in wearable exoskeletons. Gyroscopes are considered as presenting long term stability, eliminating the need for periodic recalibration. However, a disadvantage is the presence of a low frequency bias, mainly due to temperature effects. Gyroscope stability and behavior over variations in temperature is essential for the long-term performances. Lowcost strap down gyroscopes are now internally compensated in temperature and therefore the bias vector slowly oscillates around a constant average value [30].

Sensors 2010, 10 2.4. Flexible Angular Sensors and Goniometers

1409

Strain gauges are usually inserted in goniometers, in order to measure a joint angle dynamically. The electrical output is proportional to the angle. Flexible angular sensors are commonly used in rehabilitation and sports, but have not been disseminated at a large scale for studies of tremor [8,31,32]. Some goniometers include flexible springs to compensate for joint migration related to movement. Double differentiation of angular position is performed to compute joint torque T associated with oscillations of limbs, according to the following equation: T=J Where J is the limb moment of inertia and is the angular acceleration. 2.5. Videos Videotaping of patients and computerized video motion detecting systems are qualitative and quantitative methods useful of analysis of motion disorders [33]. Measuring devices based on video imaging are effective in quantifying motor impairments in clinical settings [34]. For instance, videos can be used to evaluate patients presenting Parkinsons disease. However, the videotaped motor examination may present some limitations when patients have mild symptoms [35]. Experimental procedures for assessment of tremor as well as validation of new methods benefit from videotaping [36-37]. Videotaping is also a valuable teaching tool to improve the uniform application of tremor rating scales by raters having different levels of experience in movement disorders [38]. The use of at least two cameras is recommended to compute three-dimensional estimations. A calibration method, as well as a quiet room for recording are mandatory. It should be pointed out that analysis of videos may be time-consuming. This may become an obstacle in the evaluation of neurological disorders. 2.6. Optoelectronic Devices Optoelectronic devices are popular in the field of motion analysis. Most optoelectronic devices use reflective (passive) or active markers (such as infra-red light emitting diodes LEDs) fixed on anatomical landmarks of the body, such as bone landmarks [39,40]. Although they provide accurate measurements and allow for multi-site recordings, the procedure is often time-consuming. Most of them require a calibration procedure. The estimation of three-dimensional positioning from digitized data is based on the geometrical properties of a central projection [41]. Recognition of passive markers can be performed either via recognition software or with a dedicated hardware circuit [42-44]. Active markers are detected thanks to detection of sequential pulses. The absence of wires, batteries and pulsing circuitry are advantages of passive markers, but accuracy of active markers is higher. An important drawback is due to marker visibility constraints. Also, instrumental errors can arise [42]. Calibration inaccuracies, image distortion, electrical noise, errors in positioning of markers and marker flickering (imprecision related to image conversion), as well as experimental mistakes related to the selection of the environment may occur [42,45]. Compensation methods that have been developed, such as least-squares techniques applied to coordinates of markers allow improvement of the bone pose estimate [42]. Skin deformation and displacement causes marker movement with respect to the

Sensors 2010, 10

1410

underlying bone [46]. This soft tissue artefact has a frequency content similar to the actual bone movement, is task dependent and not reproducible among subjects [46]. For lower limbs, motion about the flexion/extension axis of hips, knees and ankles can be determined reliably, but motion about other axes at those joints should be interpreted cautiously, as the soft tissue artefact produces spurious effects with amplitudes comparable to the amount of motion actually occurring in those joints [46]. 2.7. Force Sensors Detailed instantaneous monitoring of force during functional tasks is promising both for the diagnosis and management of tremor. In the last decade, robotic devices have become more popular. They combine the haptic technology with conventional neurophysiological assessment including EMG (myohaptic devices). The haptic tools are also proposed as future therapeutical methods for rehabilitation [47]. These instruments record angular motion, torques and corresponding EMG activities (in case of combination with EMG sensors), and can be used to obtain real-time spectral analysis, to apply loads and mechanical oscillations as well as characterize the effects of damping on voluntary motion in neurological patients [48]. Comparison with conventional accelerometers has demonstrated a high accuracy in terms of reproducibility of frequency and amplitudes of imposed oscillations [49]. However, several haptic devices are non portable and expensive, two factors which limit their current dissemination. Most clinical and instrumental methods evaluate tremor amplitude and frequency, without considering joint movement in the affected limb. Estimation of joint movement facilitates the identification of the muscle groups which directly contribute to tremor. Rozman et al. have developed a method based on force transducers to evaluate rest tremor of the wrist and metacarpophalangeal joints in two degrees of freedom for each joint [50]. This technique of direct measurements of joint movement in the hand could be particularly useful when therapies acting on selective muscles are considered, especially local botulinum toxin injections. 2.8. Evaluation of Handwriting and Drawing Quantification of drawings of Archimedes spiral is either based on the off-line analysis of the digitized picture obtained from a commercial scanner [51], or from data collected from one or two accelerometers fixed on the pen. Spirals can decomposed and reconstituted, for subsequent analysis. The technique is relatively cheap and easy to assess two-dimensional tremor. 2.9. Wearable Orthosis The field of wearable sensors in bioengineering applications and medicine is growing rapidly [52]. The impact of wearable sensors is anticipated to be high in tremor research, given the need to obtain accurate estimations of tremor parameters during daily life activities. Wearable sensing can be used to monitor motor difficulties, to estimate fluctuations of deficits with time and response to therapies. Parkinsons disease is a typical disorder under investigation [53]. Fine adjustment of doses and comparison of modes of administration of drugs are examples of applications.

Sensors 2010, 10

1411

Quantification of tremor can be achieved using a wearable orthosis affixed on one upper limb [2]. An exoskeleton acting in parallel to the limb makes real-time estimation of upper limb tremor. In addition to measuring oscillations, this technique allows to investigate the effects of adaptative viscous control. Gyroscopes or accelerometers can be integrated in the orthosis for tremor monitoring. The main disadvantages of the wearable orthosis are the issues of aesthetics, cosmetics and the difficulty to make them affordable for each individual. They require multiple shapes and sizes. Future studies are required for the integration of sensors and actuators in comfortable textiles designed in order to be accepted by users. Tremor suppression using smart textiles takes the advantage of material flexibility, not hindering the bearers motor functions [54]. Active and semi-active control devices have been proposed to suppress oscillations [54]. Wearable hybrid sensors are being developed, in particular sensors combining SEMG and gyroscopes. Studies on their reliability and stability are required. 2.10. Other Sensors Sensors based upon magnetic fields (search-coil technique) are occasionally used for the characterization of tremor. Their main clinical application remains the assessment or upper limb kinematics [55]. 2.11. Choice of the Sensor and Future Developments in the Assessment of Tremor The choice of the sensors should be based upon their characteristics in terms of sensitivity to gravity, signal-to-noise ratio, bandwidth, size and weight (see Table 4) [8,30]. In addition, reliability and stability of the sensors used should be considered. For historical and practical reasons (size, price, ease of use), conventional accelerometers tend to be the most commonly used sensors [15]. They are disseminated in numerous laboratories worldwide. New advances in microelectromechanical sensors (MEMS are based on the integration of mechanical elements, sensors, actuators, and electronics on a common silicon substrate using microfabrication technology) are changing the field of sensors for tremor evaluation. The most used types of MEMS sensors are solid state accelerometers, gyroscopes, and also magnetometers. Since MEMS devices are being manufactured with batch fabrication techniques (similar to the techniques used for integrated circuits), high levels of functionality and reliability can be reached at a relatively low cost. However, because most of MEMS devices contain movable parts, material fatigue and aging under long-term repeated cycling load may cause potential device failure, which in turn decreases the device reliability [56]. Simulation results show that the fatigue lifetime of MEMS accelerometers made by poly-silicon material is relatively good. For MEMS gyroscopes, fabrication imperfections result in cross stiffness and cross damping effects that may interfere with the measurement of angular velocity [57]. Minimization of the cross coupling between two axes are challenging issues in vibrating gyroscopes. Several control algorithms have been proposed, including a sliding mode control with a force balancing control strategy and an adaptative controller for tuning the natural frequency of the drive axis of a vibratory gyroscope [58-60]. The proportional-integral sliding mode adaptative control proposed by Fei and Batur takes into account both matched and unmatched uncertainties and shows consistent estimation of gyroscope parameters including angular velocity and large robustness to parameter variations and external disturbances [57].

Sensors 2010, 10

1412

Overall, inertial sensing appears promising for tremor. Wireless inertial sensor with the combination of different types of sensors (for instance three linear accelerometers, three gyroscopes and three magnetometers [61]) is now being applied. Figure 7. Example of change in alpha rhythm recorded with needle EEG on the scalp before a distal movement in contralateral upper limb. Preparation of movement is associated with a desynchronization of the EEG rhythm, with a decrease in the alpha/beta ratio. Upper panels correspond to power spectrum of EEG segments. Acc: accelerometer fixed on the hand; EMG: surface EMG in forearm muscle. Dotted vertical line: movement onset. 11 Hz 9 Hz

/
25 Hz

25 Hz

EEG

Movement

Acc.
EMG
0.5 sec
A different approach to the direct measurement of trembling legs or arms is to use the EEG signal in the framework of a Brain Computer Interface (BCI) [62]. BCIs aim at providing a communication channel between brain and the external environment. BCIs can be either invasive or non-invasive. Invasive recordings relate to the recordings of brain electrical activity at the surface of the cortex (electrocorticography) or directly in the cortex (recordings of action potentials or local field potentials), whereas non-invasive recordings are obtained from the scalp by analysing electrical activity (EEG) or magnetic activity (magnetoencephalography) [63]. Most of the current BCIs are

Sensors 2010, 10

1413

EEG-based, since EEG recordings have a high temporal resolution and can be portable, avoiding the need for implants. Control signals include movement-related cortical potentials (MRCP), slow cortical potentials, sensorimotor rhythms (SMR) and steady-state visual evoked potentials (SSVEP) [63]. For MRCP, a surface negativity occurs about two seconds before the movement onset (BP: Bereitschaftspotential). The amplitude of MRCP is influenced by the characteristics of the movement performed. Changes in SMR at the alpha band/beta band correspond to event-related desynchronization (ERD) and event-related synchronization (ERS), respectively for a power decrease and a power increase. An example is shown in Figure 7. Research is going on to extract the best parameter(s) from the BCI to detect tremor onset and differentiate tremor from voluntary movements initiated by patients. Combinations of ERD/ERS, with parameters extracted from instantaneous assessment of muscle discharges and changes in corticomotor coherence are being considered. The ERD begins up to two seconds before movement onset. MRCP and ERD/ERS take their origin from different brain sources [64]. The early BP might reflect slowly increasing cortical excitability and subconscious readiness for the forthcoming movement [64]. Pattern recognition techniques have been proposed for the classification of single-trial MRCPs, based upon optimization of wavelets and support vector machine (SVM) [63]. SVM has shown good generalization property in various applications, including BCI [63]. Recent inter-disciplinary studies aim to achieve the early detection and monitoring of tremor through a multimodal brain computer interface (BCI) combining neurophysiological (EEG, EMG) and biomechanical data (MEMS inertial sensors) in a sensor fusion approach [65]. Estimation of upperlimb orientation using miniature accelerometers and gyroscopes is accurate and renders the measurement system less obtrusive [66]. The use of low-mass sensors is important to decrease the lowpass filtering effect due to mass addition. Both the sensor fusion approach and hybrid sensors offer the possibility to extract and combine the most reliable data from the respective instruments (EEG and EMG, EMG and accelerometers, EEG and gyroscopes,...). Sensors are intercommunicative, usually with the wireless technology. It is anticipated that BCIs will be used in the near future to distinguish voluntary motion from involuntary oscillations of the body in order to trigger an external device reducing tremor (namely exoskeleton applying viscosity, myohaptic device applying damping, neuromuscular functional electrical stimulation FES). 3. Methods of Signal Processing for Tremor Tremorous activity is composed of deterministic (non random) and stochastic components. Signal processing is required to interpret time-series data of nonlinear systems and instances in which the frequency content of a signal provides more information than the original waveform [8]. 3.1. Editing Editing remains a critical step in digital processing of tremor series [8,67,68]. Data editing can be compared to a pre-analysis procedure allowing the detection and removal of low quality signals (DC offset removal,). For tremor recordings, noise may result from a problem which occurred during acquisition (instrumentation error, recording in a noisy environment, low frequency noise due to movement artefacts related to tremor,). Many experts edit through visual inspection of recordings

Sensors 2010, 10

1414

presented graphically. Automatic detection of tremor is under development. A non-linear analysis technique based on a running second order moment function (SOMF) has been suggested [68]. Alternatives include the use of extended Kalman filter (EKF) to estimate the instantaneous tremor frequency, as applied in spike trains detected from microelectrode recordings [69]. The use of normalized variance of intensity of the signal adapted from studies of pulses in radio waves could be an alternative to identify the oscillations of tremor. 3.2. Noise Minimization and Wavelets Noise minimization is generally required for tremor time series [8]. Frequency-selective filters or adaptative filters are used. Wavelet (time-scale distribution) denoising may be applied also to bandpass filter a given signal. Wavelet transforms co-localize in both frequency domain and time domain and are effective for non stationary signals such as neurological tremor [70]. The basic idea behind wavelets is to express a signal as a linear combination of given sets of functions (wavelet transform). These are obtained by shifting and dilating a function called mother wavelet. Because the signal in tremor has often many transient components which are interesting to isolate and analyze and due to noisy components, wavelet-based denoising can efficiently isolate activities of interest such as EMG discharges. Wavelet denoising allows a separation of the signal from noise [8]. 3.3. Spectral Estimations Most spectral estimation algorithms are devised for complete data sets. However, missing samples may occur. Non parametric adaptative filtering-based tools are now available to deal with missing-data samples. The Gibbs phenomenon relates to errors at points of discontinuities (a value of about 9% occurs in any region of discontinuity) [8]. A discrete Fourier transform (FT) is often used for signals composed of data sampled at given spaced intervals and a continuous signal may be reconstructed without information loss if the sampling frequency is greater than twice the highest frequency component in the signal (Nyquist critical frequency), to avoid aliasing. Direct application of the FT is typically modified by windowing. Spectra are averaged in order to reduce the variance of spectral estimates [8]. The spectral estimation methods allowing the calculation of the power spectra are the autocorrelation function, the Fast Fourier Transform (FFT) and the autoregression. Units for power spectra are in power per frequency band (e.g., (m/sec2)2/Hz). Results of spectral densities may be gathered in a time-frequency representation, which provides interesting and meaningful information for daily clinical practice [8]. Fourier transforms break down time-domain signals into constituent sinusoids of different frequencies. For periodic functions, the original waveform may be reconstructed from the sinusoidal components by application of the FT. The FFT is the most popular method to perform transformation in the frequency domain, being computationally simple and fast [24]. 3.4. Parameters Commonly Extracted for Spectral Estimations For human tremor applications, the following parameters are extracted [71-73]: -spectrum shape with identification of single or multiple peaks

Sensors 2010, 10

1415

-peak frequency -median frequency -power spectral density (PSD), in particular in the band 120 Hz (PSD1-20Hz) and 133 Hz (PSD1-33Hz), with identification of the peak intensity (PI) -power in specific frequency bands ( band: 813 Hz, 1 band: 1320 Hz, 2 band 2026 Hz; 3 band: 2633 Hz) and power ratios for each frequency band as compared to total PSD -crest factor: ratio of the PI divided by PSD1-20Hz or PSD1-33Hz -center frequency (F50): frequency of the power spectrum dividing the area under the spectrum in two equal parts -harmonic index (HI): index comparing the tremor frequency pattern with the pattern of single harmonic oscillation (HI = 1 in case of single harmonic) -frequency dispersion: corresponds to center frequency SD. Frequency dispersion is low for regular tremors. The PSD is usually averaged over several epochs or smoothed over frequencies for reducing the variance of estimation, as explained above. The spectral content of the signal is presented in either linear or logarithmic scale. The FFT takes advantage of the symmetrical properties of sinusoidal periodic waveforms. Using FFT analysis the signal is linearly decomposed as combination of sines and cosines, but most tremors are non periodic. The compromise between time and frequency resolution of these methods may not underline the presence of local oscillations in the signal, which might bring key information for the understanding of tremor. Two additional techniques have been shown to provide useful information for human tremor. First, analysis of short data segments can be performed with an autoregressive method of spectral analysis. This is interesting for the analysis of individual movements. Second, the empirical mode decomposition (EMD) can be used to isolate tremor movements from non oscillatory movements [74]. EMD decomposes time-series into intrinsic mode functions (IMFs), whose energy is quantified. Cross-spectral studies are widely used to study the relation between signals which are recorded simultaneously from distinct sources [8]. In particular, the relations between EMG, EEG and MEG have been investigated in details. The estimation of the power spectrum of a zero mean process X(t) is defined as the Fourier Transform (FT) of the auto-covariance function and is performed by a direct spectral estimation based on the FT of the measured data [8]. The cross-spectrum CS() of two zero mean process X(t) and Y(t), similarly to the univariate case, is defined as the FT of the cross-correlation function: CCF(t)=( X(t)Y(t-t)) The modulus of the cross-spectrum CS() normalized by the respective auto-spectra Sx() and Sy() gives the coherence-spectrum: Coh()= |CS( |/

Sx() Sy()

The coherence can be seen as a measure of linear predictability. Coherence value is 1 whenever X(t) is obtained from Y(t) by a linear operator. Coherence value is equal to zero whether there is no relation or the relation between the processes is a quadratic one; interpretation of coherence does not rely on the linearity of the processes X(t) and Y(t). Simultaneously recorded signals (X(t) and Y(t)) may be

Sensors 2010, 10

1416

uncorrelated or present a coherence unequal 1 due to a non linear relationship, additional influences on X(t) apart from Y(t), estimation bias or noise. 4. The Need for Sensors Providing Tremor Parameters Correlated with Clinical Scores and Functional Evaluation In order to render sensors popular in the clinical community, several criteria need to be met. First, the sensors must be designed to fit with the clinical environment and the needs of users. This takes into account aesthetics and cosmetic criteria. Second, data should be correlated as much as possible with current clinical standards in terms of neurological evaluation. Figure 8 illustrates an example of a good linear correlation for upper limb tremor between the crest factor and the clinical score in a widely used clinical manoeuver. Figure 8. Correlation between the clinical grade of postural tremor and crest factor obtained with accelerometry during a postural task in a group of eight patients presenting a neurological tremor. Accelerometer fixed on the extremity of the index finger. Sampling rate: 512 Hz, duration of acquisition: 30 seconds.

Third, data obtained should be correlated with functional scores. Figure 9 shows an example. The assessment of neurological patients is greatly improved when clinical scores, functional evaluation and parameters extracted using sensitive sensors are combined, especially for the monitoring of patients. The need for reliable sensors will increase in the future, because new therapeutic strategies for disorders previously considered as not treatable are emerging [75,76].

Sensors 2010, 10 Figure 9. Correlation between the crest factor (CF) during maintenance of a postural task and the score obtained using mechanical counters (repeated tapping between 2 counters separated by 30 cm during 30 seconds) in a group of six patients presenting a neurological tremor. Accelerometer fixed on the extremity of the index finger. Sampling rate: 512 Hz, duration of acquisition: 30 seconds.

1417

It should be underlined that the fact that sensor data should be correlated with clinical scores does not necessarily mean that a linear dependence is necessary. However, the pure clinical evaluation of tremor is often considered as a good clinical correlate of tremor severity, despite the variability of tremor in a given patient and amongst different patients [8,67]. Frequency dispersion is an example of a parameter which is not correlated with clinical staging of tremor. Patients may exhibit a strong variability in this parameter, despite a relatively stable crest factor and a high inter-observer reliability for the clinical scores, especially in postural tasks. 5. Conclusions Neurological tremor is a disabling and common symptom related to various disorders, such as Parkinsons disease, essential tremor or cerebellar ataxias. Its assessment, mainly based on popular sensors such as accelerometers and EMG, is critical for (1) the evaluation of patients, (2) the monitoring of progression of the disease and (3) the adjustment of therapies. Effort should be devoted to the selection of the sensors of interest for a given experiment in order to provide data in agreement with the clinical-functional assessment of tremor severity. Discrepancies are expected, but they should be as low as possible. Sensitive methods are required. Clinically-validated sensors small in size, light in weight and reconfigurable have the highest chances of becoming successful.

Sensors 2010, 10

1418

In addition to conventional sensors, new avenues of research are emerging to characterize and reduce tremor during daily activities. The so-called sensor fusion approach combines information from various sensors and might provide new ways to characterize oscillatory phenomena with direct clinical implications at middle term, including for BCIs. The field of wearable sensors is expanding quickly and there is a hope that body area network will greatly help the management of disabling tremor soon, using clusters of conventional or hybrid sensors. Acknowledgments The work has been supported by the EC (ICT-2007-224051). MM is supported by the FNRS Belgium. References and Notes 1. Louis, E.D.; Marder, K.; Cote, L. Differences in the prevalence of essential tremor among elderly African Americans, whites, and Hispanics in northern Manhattan, NY. Arch. Neurol. 1995, 52, 1201-1205. Rocon, E.; Manto, M.; Pons, J.; Camut, S.; Belda, J.M. Mechanical suppression of essential tremor. Cerebellum 2007, 6, 73-78. Deuschl, G.; Bain, P.; Brin M. Consensus statement of the Movement Disorder Society on. Tremor. Ad Hoc Scientific Committee. Mov. Disord. 1998, 13, 2-23. Elble, R.J. Characteristics of physiologic tremor in young and elderly adults. Clin. Neurophysiol. 2003, 114, 624-35. Elble, R.J. Central mechanisms of tremor. J. Clin. Neurophysiol. 1996, 13, 133-144. Dietz, V.; Hillesheimer, W.; Freund, H.J. Correlation between tremor, voluntary contraction, and firing pattern of motor units in Parkinson's disease. J. Neurol Neurosurg. Psychiatry 1974, 37, 927-937. Findley, L.J.; Koller, W.C. Handbook of Tremor Disorders; Marcel Dekker: New York, NY, USA, 1995. Grimaldi, G.; Manto, M. Tremor: From Pathogenesis to Treatment; Morgan & Claypool: San Rafael, CA, USA, 2008. Kronenbuerger, M.; Konczak, J.; Ziegler, W.; Buderath, P.; Frank, B.; Coenen, V.A.; Kiening, K.; Reinacher, P.; Noth, J.; Timmann, D. Balance and motor speech impairment in essential tremor. Cerebellum 2009, 8, 389-398. Tsuji, S.; Onodera, O.; Goto, J.; Nishizawa, M. Study Group on Ataxic Diseases. Sporadic ataxias in Japan--a population-based epidemiological study. Cerebellum 2008, 7, 189-197. Timmann, D.; Brandauer, B.; Hermsdrfer, J.; Ilg, W.; Konczak, J.; Gerwig, M.; Gizewski, E.R.; Schoch, B. Lesion-symptom mapping of the human cerebellum. Cerebellum 2008, 7, 602-606. Dhlinger, S.; Hauser, T.K.; Borkert, J.; Luft, A.R.; Schulz, J.B. Magnetic resonance imaging in spinocerebellar ataxias. Cerebellum 2008, 7, 204-214. Mathiowetz, V.;Weber, K.; Kashman, N.;Volland, G. Adult norms for the nine hole peg test of finger dexterity. Occup. Ther. J. Res. 1985, 5, 24-37.

2. 3. 4. 5. 6.

7. 8. 9.

10. 11. 12. 13.

Sensors 2010, 10

1419

14. Mathiowetz, V.; Volland, G.; Kashman, N.; Weber, K. Adult norms for the box and block test of manual dexterity. Am. J. Occup. Ther. 1985, 39, 386-391. 15. Wong, W.Y.; Wong, M.S.; Lo, K.H. Clinical applications of sensors for human posture and movement analysis: a review. Prosthet. Orthot. Int. 2007, 31, 62-75. 16. Luinge, H.J.; Veltink, P.H. Measuring orientation of human body segments using miniature gyroscopes and accelerometers. Med. Biol. Eng. Comput. 2005, 43, 273-282. 17. Raethjen, J.; Lauk, M.; Kster, B.; Fietzek, U.; Friege, L.; Timmer, J.; Lcking, C.H.; Deuschl, G. Tremor analysis in two normal cohorts. Clin. Neurophysiol. 2004, 9, 2151-2156. 18. Wastensson, G.; Lamoureux, D.; Sllsten, G.; Beuter, A.; Barregrd, L. Quantitative assessment of neuromotor function in workers with current low exposure to mercury vapor. Neurotoxicology 2008, 29, 596-604. 19. Safronov, V.A.; Shevelev, I.N.; Shabalov, V.A.; Panin, V.V. Localization of a Tremor Generator. Human Physiol. 2001, 4, 454-457. 20. Daube, J.R.; Rubin, D.I. Needle electromyography. Muscle Nerve 2009, 39, 244-270. 21. Manto, M.; Sauvage, C.; Roark, R.M. Unifying hypothesis for the motoneuronal code in neurological disorders. Biosci. Hypoth. 2008, 1, 93-99. 22. McAuley, J.; Rothwell, J. Identification of psychogenic, dystonic, and other organic tremors by a coherence entrainment test. Mov. Disord. 2004, 19, 253-267. 23. Breit, S.; Spieker, S.; Schulz, J.B., Gasser, T. Long-term EMG recordings differentiate between parkinsonian and essential tremor. J. Neurol. 2008, 255, 103-11. 24. Engin, M. A recording and analysis system for human tremor. Measurement 2007, 40, 288-293. 25. Roetenberg, D.; Luinge, H.J.; Baten, C.T.M.; Veltink, P.H. Compensation of magnetic disturbances improves inertial and magnetic sensing of human body segment orientation. IEEE Trans. Neural. Syst. Rehab. Eng. 2005, 13, 395-405. 26. Morrison, S.; Newell, K.M. Postural and resting tremor in the upper limb. Clin. Neurophysiol. 2000, 111, 651-663. 27. Patel, S.; Lorincz, K.; Hughes, R.; Huggins, N.; Growdon, J.; Standaert, D.; Akay, M.; Dy, J.; Welsh, M.; Bonato, P. Monitoring motor fluctuations in patients with Parkinson's disease using wearable sensors. IEEE Trans. Inf. Technol. Biomed. 2009, 13, 864-873. 28. Journee, H.L.; Postma, A.A.; Staal, M.J. Intraoperative neurophysiological assessment of disabling symptoms in DBS surgery. Neurophysiol. Clin. 2007, 37, 467-475. 29. Manto, M.; Marmolino, D. Animal models of human cerebellar ataxias: a cornerstone for the therapies of the twenty-first century. Cerebellum 2009, 8, 137-154. 30. Metni, N.; Pflimlin, J.M.; Hamel, T.; Soures, P. Attitude and gyro bias estimation for a VTOL UAV. Contr. Eng. Pract. 2006, 14, 1511-1520. 31. Salarian, A.; Russmann, H.; Wider, C.; Burkhard, PR.; Vingerhoets, FJG.; Aminian, K. Quantification of tremor and bradykinesia in Parkinsons disease using a novel ambulatory monitoring system. IEEE Trans. Biomed. Eng. 2007, 2, 313-322. 32. Zhou, H.; Hu, H.; Tao, Y. Inertial measurements of upper limb motion. Med. Bio. Eng. Comput. 2006, 44, 479-487 33. Zhang, T.; Wei, G.; Yan, Z.; Ding, M.; Li, C.; Ding, H.; Xu, S. Quantitative assessment of Parkinson's disease deficits. Chin. Med. J. (English) 1999, 112, 812-815.

Sensors 2010, 10

1420

34. Swider, M. The application of video image processing to quantitative analysis of extremity tremor in humans. J. Neurosci. Methods 1998, 84, 167-172. 35. Louis, E.D.; Levy, G.; Cte, L.J.; Mejia, H.; Fahn, S.; Marder, K. Diagnosing Parkinson's disease using videotaped neurological examinations: validity and factors that contribute to incorrect diagnoses. Mov. Disord. 2002, 17, 513-517. 36. Louis, E.D.; Ford, B.; Bismuth, B. Reliability between two observers using a protocol for diagnosing essential tremor. Mov. Disord. 1998, 13, 287-293. 37. Stacy, M.A.; Elble, R.J.; Ondo, W.G.; Wu, S.C.; Hulihan, J.; TRS study group. Assessment of interrater and intrarater reliability of the Fahn-Tolosa-Marin Tremor Rating Scale in essential tremor. Mov. Disord. 2007, 22, 833-838. 38. Louis, E.D.; Barnes, L.; Wendt, K.J.; Ford, B.; Sangiorgio, M.; Tabbal, S.; Lewis, L.; Kaufmann, P.; Moskowitz, C.; Comella, C.L.; Goetz, C.C.; Lang, A.E. A teaching videotape for the assessment of essential tremor. Mov. Disord. 2001, 16, 89-93. 39. Klebe, S.; Stolze, H.;Grensing, K.;Volkmann, J.; Wenzelburger, R.; Deuschl, G. Influence of alcohol on gait in patients with essential tremor. Neurology 2005, 65, 96-101. 40. Pellegrini, B.; Faes, L.; Nollo, G.; Schena, F. Quantifying the contribution of arm postural tremor to the outcome of goal-directed pointing task by displacement measures. J. Neurosci. Meth. 2004, 139, 185-193. 41. Braune, W.; Fischer, O. The human gait; Springer: Berlin, Germany, 1987. 42. Chiari, L.; Della Croce, U.; Leardini, A.; Cappozzo, A. Human movement analysis using stereophotogrammetry. Part 2: instrumental errors. Gait Posture 2005, 21, 197-211. 43. Taylor, K.D.; Mottier, F.M.; Simmons, D.W.; Cohen, W.; Pavlak, R. Jr.; Cornell, D.P.; Hankins, G.B. An automated motion measurement system for clinical gait analysis. J. Biomech. 1982, 15, 505-516. 44. Ferrigno, G., Borghese, N.A.; Pedotti, A. Pattern-recognition in 3-D automatic human motion analysis. ISPRS J. Photogramm. 1990, 45, 227-246. 45. Gazzani, F. Comparative assessment of two algorithms for calibrating stereophotogrammetric systems. J. Biomech. 1993, 26, 1449-1454. 46. Leardini, A.; Chiari, L.; Della Croce, U.; Cappozzo, A. Human movement analysis using stereophotogrammetry. Part 3. Soft tissue artifact assessment and compensation. Gait Posture 2005, 21, 212-25. 47. Avizzano, C. A.; Bergamasco, M. Technological aids for the treatment of the tremor. International Conference on Rehabilitation Robotics (ICORR), Stanford University, Stanford, CA, USA, 1999. 48. Grimaldi, G.;Lammertse, P.; Van Den Braber, N.; Meuleman, J.; Manto, M. Effects of inertia and wrist oscillations on contralateral neurological postural tremor using the wristalyzer, a new myohaptic device. IEEE Trans. Biomed. Circuits Syst. 2008, 4, 269-279 49. Grimaldi, G.; Lammertse, P.; Manto, M. Effects of wrist oscillations on contralateral neurological postural tremor using a new myohaptic device (wristalyzer). In Proceedings of the 4th IEEE/EMBS International Summer School and Symposium on Medical Devices and Biosensors, Cambridge, UK, August 2007; pp. 44-48. 50. Rozman, J.; Bartoli, A., Ribaric, S. A new method for selective measurement of joint movement in hand tremor in Parkinson's disease patients. J. Med. Eng. Technol. 2007, 31, 305-311.

Sensors 2010, 10

1421

51. Miralles, F.; Tarong, S.; Espino, A. Quantification of the drawing of an Archimedes spiral through the analysis of its digitized picture. J. Neurosci. Meth. 2006, 152, 18-31. 52. Bonato, P. Wearable sensors/systems and their impact on biomedical engineering. IEEE Eng. Med. Biol. Magazine 2003, 3, 18-20. 53. Patel, S.; Mancinelli, C.; Hughes, R.; Dalton, A.; Shih, L.; Bonato, P. Optimizing deep brain stimulation settings using wearable sensing technology. Neural. Eng. 2009, 6-9. 54. Swallow, L.; Siores, E. Tremor Suppression Using Smart Textile Fibre Systems. J.F.B.I. 2009, 4, 261-266. 55. Timmann, D.; Lee, P.; Watts, S.; Hore, J. Kinematics of arm joint rotations in cerebellar and unskilled subjects associated with the inability to throw fast. Cerebellum 2008, 7, 366-378. 56. Xiong, X.; Wu, Y.L.; Jone, W.B. Material Fatigue and Reliability of MEMS Accelerometers. In Proceedings of the 2008 IEEE International Symposium on Defect and Fault Tolerance of VLSI Systems, Cambridge, MA, USA, October 1-3, 2008; pp. 314-322. 57. Fei, J.; Batur, C. A novel adaptive sliding mode control with application to MEMS gyroscope. ISA Trans. 2009, 48, 73-78 58. Batur, C.; Sreeramreddy, T.; Khasawneh, Q. Sliding mode control of a simulated MEMS gyroscope. ISA Trans. 2006, 45, 99-108. 59. Leland, R. Adaptive mode tuning for vibrational gyroscopes. IEEE Trans. Control Syst. Technol. 2003, 11, 242-247. 60. Park, S.; Horowitz, R. New adaptive mode of operation for MEMS gyroscopes. ASME Trans. Dyn. Syst. Meas. Control 2004, 126, 800-810. 61. Brunetti, F.J.; Rocon, E.; Pons, J.L.; Manto, M. The tremor coherence analyzer (TCA): a portable tool to assess instantaneous inter-muscle coupling in tremor. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2004, 1, 61-64. 62. Lin, C.T.; Ko, L.W.; Chang, M.H.; Duann, J.R.; Chen, J.Y.; Su, T.P.; Jung, T.P. Review of Wireless and Wearable Electroencephalogram Systems and Brain-Computer InterfacesA MiniReview. Gerontology 25 July 2009 (Epub ahead of print). 63. Gu, Y.; do Nascimento, O.F.; Lucas, M.F.; Farina, D. Identification of task parameters from movement-related cortical potentials. Med. Biol. Eng. Comput. 4 September 2009 (Epub ahead of print). 64. Shibasaki, H.; Hallett, M. What is the Bereitschaftspotential? Clin. Neurophysiol. 2006, 117, 2341-2356. 65. The tremor project (ICT-2007-224051). Available online: http://www.iai.csic.es/tremor/ (accessed on 11 February 2010). 66. Hyde, R.A.; Ketteringham, L.P.; Neild, S.A.; Jones, R.S. Estimation of upper-limb orientation based on accelerometer and gyroscope measurements. IEEE Trans. Biomed. Eng. 2008, 55, 746-754. 67. Manto, M.; Grimaldi, G.; Lorivel, T.; Farina, D.; Popovic, L.; Conforto, S.; D'Alessio, T.; BeldaLois, J.M.; Pons, J.L.; Rocon, E. Bioinformatic Approaches Used in Modelling Human Tremor. Curr. Bioinf. 2009, 2, 154-172.

Sensors 2010, 10

1422

68. Journe, H.L.; Postmab, A.A.; Sunc, M.; Staala, M.J. Detection of tremor bursts by a running second order moment function and analysis using interburst histograms. Med. Eng. Physics 2008, 1, 75-83. 69. Kim, S.; McNames, J. Tracking Tremor Frequency in Spike Trains Using the Extended Kalman Filter Engineering in Medicine and Biology Society. In Proceedings of 27th Annual Conference of the IEEE EMBS, Shanghai, China, September 1-4, 2005; pp. 7576-7579. 70. Strambi, S.K.; Rossi, B.; De Michele, G.; Sello, S. Effect of medication in Parkinsons disease: a wavelet analysis of EMG signals. Med. Eng. Physics 2004, 26, 279-290. 71. McAuley, J.H.; Rothwell, J.C.; Marsden, C.D. Frequency peaks of tremor, muscle vibration and electromyographic activity at 10 Hz, 20 Hz and 40 Hz during human finger muscle contraction may reflect rhythmicities of central neural firing. Exp. Brain Res. 1997, 114, 525-541. 72. Duval, C. Rest and postural tremors in patients with Parkinsons disease. Brain Res. Bull. 2006, 70, 44-48. 73. Edwards, R.; Beuter, A. Indexes for identification of abnormal tremor using computer tremor evaluation systems. IEEE Trans. Biomed. Eng. 1999, 7, 895-898. 74. Rocon, E.; Pons, J.L.; Andrade, A.O.; Nasuto, S.J. 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. 75. Matilla-Dueas, A. The highly heterogeneous spinocerebellar ataxias: from genes to targets for therapeutic intervention. Cerebellum 2008, 7, 97-100. 76. Underwood, B.R.; Rubinsztein, D.C. Spinocerebellar ataxias caused by polyglutamine expansions: a review of therapeutic strategies. Cerebellum 2008, 7, 215-221. 2010 by the authors; licensee Molecular Diversity Preservation International, Basel, Switzerland. This article is an open-access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/).