| Literature DB >> 34988457 |
James Peters1, Mohammod Abdul Motin2,3, Laura Perju-Dumbrava1, Sheik Mohammed Ali2, Catherine Ding1, Michael Eller1, Sanjay Raghav1, Dinesh Kant Kumar2, Peter Kempster1,4.
Abstract
We investigated whether computerised analysis of writing and drawing could discriminate essential tremor (ET) phenotypes according to the 2018 Consensus Statement on the Classification of Tremors. The Consensus scheme emphasises soft additional findings, mainly motor, that do not suffice to diagnose another tremor syndrome. Ten men and nine women were classified by blinded assessors according to Consensus Axis 1 definitions of ET and ET plus. Blinded scoring of tremor severity and alternating limb movement was also conducted. Twenty healthy participants acted as controls. Four writing and three drawing tasks were performed on a Wacom Intuos Pro Large digital tablet with a pressure-sensor mounted ink pen. Sixty-seven computerised measurements were obtained, comprising static (dimensional and temporal), kinematic and pen pressure features. The mean age of ET participants was 67.2±13.0 years and mean tremor duration was 21.7±19.0 years. Six were classified as ET, five had one plus feature and eight had two plus features. The computerised analysis could predict the presence and number of ET plus features. Measures of acceleration and variation of pen pressure performed strongly to separate ET phenotypes (p<0.05). Plus features were associated with higher scores on the Fahn-Tolosa-Marin Tremor Rating Scale (p=0.001) and it appeared that ET groups were mainly being separated according to severity of tremor and by compensatory manoeuvres used by participants with more severe tremor. There were, in addition, a small number of negative kinematic correlations suggesting some slowness with ET plus. Abnormal repetitive limb movement was also correlated with tremor severity (R=0.57) by clinical grading. Critics of the Consensus Statement have drawn attention to weaknesses of the ET plus concept in relation to duration and severity of ET. This classification of ET may be too biased towards tremor severity to assist in distinguishing underlying biological differences by clinical measurement. © Author(s) (or their employer(s)) 2021. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.Entities:
Keywords: clinical neurology; movement disorders; tremor
Year: 2021 PMID: 34988457 PMCID: PMC8679070 DOI: 10.1136/bmjno-2021-000212
Source DB: PubMed Journal: BMJ Neurol Open ISSN: 2632-6140
Inventory of 67 computed features for writing and drawing
| Feature type | Feature abbreviation | Feature description | Feature dimension |
| Static features | |||
| Dimensional | avgStrLenAir, medStrLenAir, maxStrLenAir, minStrLenAir, stdStrLenAir, iqrStrLenAir, avgStrLen, medStrLen, maxStrLen, minStrLen, stdStrLen, iqrStrLen | Average, median, maximum, minimum, SD and IQR of pen stroke length in air and on tablet surface. | 6×2=12 |
| Temporal | inAirT, onSurT, totalT, ratioAirtoSurT, satPreDur, noPause | Time in air; time on tablet surface; total time; ratio of air-to-surface time; duration with saturated pressure; number of pauses. | 6 |
| Kinematic features | |||
| Velocity | avgVx, medVx, maxVx, stdVx, avgVy, medVy, maxVy, stdVy, ratioVxVy | Average, median, maximum and SD of pen tip velocity in | 4×2+1=9 |
| Acceleration | avgAx, medAx, maxAx, stdAx, avgAy, medAy, maxAy, stdAy, ratioAxAy | Average, median, maximum and SD of pen tip acceleration in | 4×2+1=9 |
| Jerkiness | avgJx, medJx, maxJx, stdJx, avgJy, medJy, maxJy, stdJy, ratioJxJy | Average, median, maximum and SD of pen tip jerk in | 4×2+1=9 |
| Linear speed | avgS, medS, maxS, stdS | Average, median, maximum and SD of pen tip speed. | 4 |
| Angular speed | avgAngS, medAngS, maxAngS, stdAngS | Average, median, maximum and SD of pen tip angular speed. | 4 |
| Direction changes | Dx, Dy, Dxy, diffDxy | Total number of direction changes in | 4 |
| Pen pressure features | |||
| Pressure | avgP, medP, stdP, skewP, kurtP | Average, median and SD of pen tip pressure applied on the tablet surface; skewness and kurtosis of pressure. | 5 |
| Pressure differences | avgDiffP, medDiffP, stdDiffP, skewDiffP, kurtDiffP | Average, median and SD of pen tip pressure difference applied on the tablet surface; skewness and kurtosis of pressure difference. | 5 |
Demographic and clinical features according to Consensus Axis 1 definitions
| Control | ET-0 | ET+1 | ET+2 | P value | |
| Number | 20 | 6 | 5 | 8 | – |
| M:F | 9:11 | 3:3 | 1:4 | 6:2 | – |
| Age | 64.3±12.0 | 62.8±16.5 | 61.4±11.6 | 74.0±8.5 | 0.191* |
| Age at tremor onset | – | 37.7±18.5 | 51.2±13.1 | 46.3±22.4 | 0.505* |
| Tremor duration | – | 25.2±23.2 | 10.2±5.8 | 26.4±19.7 | 0.300* |
| FTM score | – | 16.0±9.2 | 20.4±3.8 | 39.8±13.1 |
|
| MDS-UPDRS-III 3.4–3.8 | – | 1.0±1.6 | 2.5±4.8 | 4.4±4.0 | 0.063† |
| MoCA | – | 25.8±4.2 | 25.2±3.3 | 24.5±3.6 | 0.80* |
Values shown as mean±SD.
*By one-way analysis of variance.
†By Kruskal-Wallis.
FTM, Fahn-Tolosa-Marin Tremor Rating Scale; MDS-UPDRS-III, Movement Disorder Society-Unified Parkinson’s Disease Rating Scale Part III; MoCA, Montreal Cognitive Assessment.
Figure 1Clinical scale results by ET phenotype. Mean and SD shown by column height and error bars. ET, essential tremor; FTM, Fahn-Tolosa-Marin; MDS-UPDRS-III, Movement Disorder Society-Unified Parkinson’s Disease Rating Scale Part III, MoCA, Montreal Cognitive Assessment.
Discriminative features for each task by ET phenotype
| Task | Feature | Control | ET-0 | ET+1 | ET+2 | P value | Spearman correlation coefficient |
| Writing tasks | |||||||
| 1 | stdStrLen | 13.35±26.3 | 21.01±120 | 41.15±29.15 | 70.62±57.71 | 0.078* | 0.51 |
| avgStrLenAir | 83.77±35.47 | 99.18±63.96 | 76.75±18.41 | 128.08±37.9 | 0.077* | 0.35 | |
| 2 | totalT | 6.2±0.42 | 6.33±0.24 | 6.14±0.19 | 6.58±0.21 |
| 0.39 |
| ratioAxAy | 1±0.21 | 0.90±0.14 | 0.66±0.25 | 1.41±0.62 |
| 0.46 | |
| 3 | avgDiffP | 620.3±274.71 | 639.4±275.25 | 686±190.2 | 440±135.17 | 0.06* | −0.34 |
| ratioAirSurT | 0.67±0.37 | 0.58±0.38 | 0.78±0.25 | 0.51±0.17 | 0.07* | −0.25 | |
| 4 | medAx | 1±0.62 | 0.5±0.18 | 0.5±0.22 | 1±0.92 |
| 0.59 |
| avgAx | 2.5±4.3 | 1.39±0.6 | 1.23±3.7 | 2.87±7.44 | 0.053* | 0.55 | |
| Drawing tasks | |||||||
| 5a | avgDiffP | 93.92±40.71 | 88.22±35.65 | 159.32±25.56 | 171.05±155.51 |
| 0.60 |
| skewDiffP | 18.86±7.23 | 17.80±4.37 | 10.35±6.12 | 8.72±6.57 |
| −0.56 | |
| 5b | medAx | 0.25±0.07 | 0.25±0.20 | 0.50±0.11 | 0.875±1.92 |
| 0.58 |
| Dxy | 33.5±16.49 | 39±104.11 | 120±55.41 | 176.50±206.5 | 0.068* | 0.54 | |
| 6 | medVx | 4.75±1.24 | 3.25±3.80 | 6±4.07 | 16.25±15.31 |
| 0.66 |
| maxAngS | 21.17±28.0 | 21.15±9.44 | 22.01±8.73 | 16.76±59.76 |
| −0.54 | |
| minStrLenAir | 82.50±41.65 | 103.50±54.12 | 9±48.59 | 7.50±30.45 |
| −0.65 | |
| maxJx | 50.08±86.2 | 50.26±22.0 | 50.19±19.9 | 39.29±170.2 |
| −0.56 | |
All features with p<0.05 are shown, otherwise the two features with lowest p values for each task. The p values and Spearman correlation coefficients were calculated for ET-0, ET+1 and ET+2.
*By Kruskal-Wallis.
ET, essential tremor.
Most discriminative features by Spearman rank order correlation coefficient (in brackets) to differentiate ET from age-matched controls
| Rank of features | Task 1 | Task 2 | Task 3 | Task 4 | Task 5a | Task 5b | Task 6 |
| 1 | stdStrLen | Dx | kurtP | maxStrLenAir | Dxy | Dx | Dxy |
| 2 | iqrStrLen | maxStrLen | skewP | Dxy | Dx | medAx | medAx |
| 3 | Dxy | minStrLen | medAngS | Dx | medAx | Dxy | Dy |
Spearman correlation analysis of raw FTM tremor score: top five computerised features for each task
| Writing tasks | |||||||
| Features Task 1 | Spearman correlation coefficient | Features Task 2 | Spearman correlation coefficient | Features | Spearman correlation coefficient | Features | Spearman correlation coefficient |
| stdStrLen | 0.58 | satPreDur | 0.50 | kurtP | 0.46 | Dxy | 0.55 |
| iqrStrLen | 0.56 | maxStrLen | 0.49 | ratioAxAy | 0.45 | Dx | 0.53 |
| maxStrLen | 0.52 | minStrLen | 0.49 | medStrLen | 0.45 | Dy | 0.51 |
| minStrLen | 0.52 | Dx | 0.49 | skewP | −0.44 | maxStrLenAir | 0.50 |
| medStrLen | 0.49 | ratioJxJy | 0.49 | Dx | 0.43 | inAirTime | 0.49 |
Spearman correlation analysis of MDS-UPDRS-III limb bradykinesia score: top five computerised features for each task
| Writing tasks | |||||||
| Features | Spearman correlation coefficient | Features | Spearman correlation coefficient | Features | Spearman correlation coefficient | Features | Spearman correlation coefficient |
| iqrStrLen | 0.58 | satPreDur | 0.64 | ratioAxAy | 0.50 | Dxy | 0.56 |
| stdStrLen | 0.54 | ratioJxJy | 0.55 | Dx | 0.44 | Dx | 0.53 |
| Dxy | 0.49 | Dx | 0.49 | medAy | −0.40 | Dy | 0.52 |
| Dy | 0.47 | maxStrLen | 0.48 | medAngS | −0.40 | maxStrLenAir | 0.51 |
| Dx | 0.45 | minStrLen | 0.48 | Dxy | 0.39 | totalT | 0.49 |
Figure 2Tremor severity (FTM score) versus repetitive limb movement (MDS-UPDRS-III 3.4–3.8 score). FTM, Fahn-Tolosa-Marin Tremor Rating Scale; MDS-UPDRS-III, Movement Disorder Society-Unified Parkinson’s Disease Rating Scale Part III.