The present disclosure relates to a system and device for monitoring movement for objectively quantifying a motor control disorder in a subject, as well as an automated method for doing the same. It relates particularly, but not exclusively, to a system, device and automated method for objectively quantifying a motor control disorder using movement detection device.
Cerebellar Ataxia (CA or Ataxia) is a complex clinical term which carries two broad meanings. First, it is used to describe the abnormal movement that arises from dysfunction of the cerebellum or of its afferent paths. Thus, ataxia can result from strokes or Multiple sclerosis (MS) affecting the cerebellum, from neurodegenerations, from infections (e.g. post chicken pox), from toxins (alcohol), from many genetic disorders and from disorders of the afferent pathways (e.g. injury to the peripheral nerves or vestibular apparatus). There are also numerous episodic ataxias and there is an uncertain relationship with migraine. Second, ataxia appears in the specific name of specific disorders described in the 19th century —such as Friedreich's Ataxia and Spinocerebellar Ataxia.
Friedreich Ataxia (FA) is a life-shortening, progressive, autosomal recessive disorder. It is a neurodegenerative disease affecting the central and peripheral nervous systems, the heart, the pancreas and the musculoskeletal system. The neurological manifestations of FA include gait ataxia, loss of limb reflexes, dysarthria, reduction in proprioception and loss of vibrational sense. Other manifestations include cardiomyopathy, scoliosis, foot deformity and diabetes. Clinical features typically appear during adolescence (between 5 and 15 years), usually manifesting as gait ataxia and loss of coordination. A person with FA is wheelchair bound on average 15 years after symptom onset. Limb ataxia affects coordination and dexterity making it increasingly difficult to engage in basic daily activities such as dressing, drinking, eating and toileting known collectively as Activities of Daily Life (ADL). Difficulties with ADL may lead to embarrassment, social isolation and loss of independence and is an important index of the quality of life experienced by an individual with FA.
The overall prevalence of ataxia is not known. Prevalence of inherited ataxias is around 0.004-0.008%. The prevalence of other forms of ataxia (excluding alcohol) is not known but together are much more common, possibly 5 times more common.
Treatment is currently limited. Physiotherapy is commonly used, but with the exception of vestibular hypofunction (impaired function of the inner ear balance system) the evidence of benefit is limited. On the other hand, there is active research, with the possibility of imminent molecular therapies for the inherited ataxias.
The clinical characteristics of ataxia were described over 100 years ago by Gordon Holmes and still form the basis of current clinical examination which is the main means of assessing ataxia. Assessments involve observing a subject performing several specific tasks chosen because they emphasise ataxic movements. Accurate assessment requires experience and, as most doctors (even neurologists) have limited exposure to ataxia, many people with ataxia are initially misdiagnosed. This clinical examination has been codified by rating scales such as the Friedreich Ataxia Rating Scale (FARS), and the International Cooperative Scale for Rating Ataxia (ICARS). These rating scales have been criticised for their subjectivity. The Nine Hole Peg test (9HPT) and the Box and Block test (BBT) are two tests that measure upper limb function in individuals with FA. The 9HPT measures the time taken to place pegs in a pegboard and the BBT measures the number of blocks placed in sections of a box in 60 seconds. These tests address the subjectivity of the rating scales but are limited by ‘floor’ and ‘ceiling’ effects. Rehabilitation centres (mainly physiotherapists) also use devices such as the Balance Master and Gaitrite to assess balance and gait however these are costly and of limited value for diagnosis.
There have been studies that use sensors to measure movement with an objective of measuring and objectively quantifying movement disorders. Some have involved use of video cameras to capture upper limb functionality of children, and game based rehabilitation exercise have been developed to measure upper limb displacements in individuals with FA. However, none to date have provided a reliable alternative to clinical assessment methods.
It would be desirable to provide a means for providing a measure or assessment of ataxia which is objective and ideally, measures severity as well as existence of ataxia.
The discussion of the background is included herein including reference to documents, acts, materials, devices, articles and the like is included to explain the context of the present disclosure. This is not to be taken as an admission or a suggestion that any of the material referred to was published, known or part of the common general knowledge in Australia or in any other country as at the priority date of any of the claims.
A criticism of current clinical tests for determining presence or severity of ataxia is that they are subjective and so clinical assessments can vary between clinicians, even experienced clinicians and neurologists. Thus, an aspect of the disclosure is an instrumented method for assessment that involves a movement detection device for objective monitoring of the subject's movement.
Viewed from one aspect, the present disclosure provides a movement monitoring system for objectively quantifying a motor control disorder in a subject, the system comprising: (a) a movement detection device generating movement data representing movement of a limb of the subject, wherein the movement detection device comprises sensors measuring at least motion of the device and pressure applied to the device by the subject; and (b) an analyser for analysing the movement data. The analyser comprises a processor and a memory containing code which, when executed by the processor: (i) receives the movement data generated by the movement detection device; (ii) applies the received movement data to an algorithmic model stored in the memory and identifies one or more features from the movement data that represent disordered movement by the subject; and (ii) calculates from the one or more identified features a score corresponding to the existence of the motor control disorder in the subject.
Preferably the movement detection device generates movement data representing movement of an upper limb of the subject. The movement detection device may be e.g. body worn or may be attached to or comprise part of an object of daily living and ideally is compressed during movement a task performed by a subject.
In some embodiments, the analyser applies the received movement data to one or more of: (a) a first algorithmic model to identify a first set of features used by the processor to calculate a selection score which is indicative of presence or absence of the motor control disorder in the subject; (b) a second algorithmic model to identify a second set of features used by the processor to calculate a severity score which is indicative of severity of the motor control disorder in the subject; and (c) a third algorithmic model to identify a third set of features used by the processor to calculate a progression score which is indicative of progression of the motor control disorder in the subject. Typically, the first set of features and the second set of features are not identical, although some features may be common to both sets and one feature set may be a subset of the other feature set.
In some embodiments, the severity score calculated by the processor corresponds to a score obtained according to a clinical scale. For example, the severity score may correspond to a score that would otherwise be determined by a clinician performing assessment of the subject using assessment scales such as modified-FARS (mFARS), FARS, 9HPT, BBT, Scale for the Assessment and Rating of Ataxia (SARA), ICARS, ADL and the like.
In some embodiments, e.g. where the motor control disorder is ataxia, the first set of features used to calculate a selection score may comprise a single feature, ideally a pressure feature such as resonant frequency of pressure (PrRF). In other embodiments, the first set of features comprises some or all of the features PrRF, resonant frequency of acceleration in the x-axis (ACC
In some embodiments, e.g. where the motor control disorder is ataxia, the second set of features used to calculate a severity score may comprise a single feature, or a plurality of or all the features selected from the group comprising: PrRFPrM, At and θRFc. In some embodiments, the feature set comprising PrM, At and θRFc is preferred.
In some embodiments, e.g. where the motor control disorder is ataxia, the third set of features used to calculate a progression score may comprise a single feature or a plurality of or all the features selected from a group comprising: MRpr, SRFgyr, MSEIMF
In some embodiments, e.g. when the motor control disorder is spasticity, the set of features used to indicate presence of spasticity may comprise one or more or all of Standard deviation value of Pressure (PrSD), RMS value of Pressure (Prrms), magnitude at resonance of pressure (PrMR) and PrRF.
In some embodiments, the analyser categorises movement dysfunction in the subject by the processor calculating a contribution made by each of the first, second or third set of features to each of a plurality of movement characteristics that are attributable to movement dysfunction in the subject. The set of movement characteristics may correlate to clinically accepted descriptions of movement disorder. Such descriptions may relate, in some embodiments, to e.g. stability, timing, accuracy and rhythmicity of the movement. In some embodiments, the analyser sums the contribution made by each of the features to each of the plurality of movement characteristics to determine a collective contribution to each of the plurality of movement characteristics.
In some embodiments, the movement detection device comprises sensors measuring one or more of (a) position of the limb; (b) acceleration of the limb; and (c) angular position of the limb. In some embodiments, the movement detection device RECTIFIED SHEET (RULE 91) ISA/AU simulates or is incorporated into an object of daily living such as a cup, spoon, brush or the like and comprises one or more of: (a) a pressure sensor, (b) an accelerometer, and (c) a gyroscope.
Another criticism of current clinical tests is that they lack functional relevance to the daily activities which are perceived as important for people with ataxia or measures salient for rehabilitation. An important component of daily life for a person with FA is the capacity to drink a beverage independently. This requires the capacity to reach, grasp, transport and release a cup in a continuum of motion in preparation for drinking the contents of a canister or container. These subcomponents of the task are also reflected in other daily living tasks. Ensuring individuals with FA can continue to independently administer liquids may be an important achievement of therapeutic interventions.
Thus in some embodiments, the movement detection device comprises a canister with a grasping portion and a pressure sensor for measuring pressure applied to the grasping portion by the subject.
Preferably the movement monitoring system is configured to measure the motor control disorder during a movement task.
The movement monitoring system may be utilised to quantify a range of motor control disorders such as, but not limited to ataxia (in its various forms) and spasticity.
Viewed from another aspect, the present disclosure provides a movement detection device for use with a system for objectively quantifying motor control disorder in a subject, the movement detection device comprising: (a) a grasping portion; and (b) a movement sensor comprising at least a pressure sensor generating pressure data representing pressure applied to the grasping portion; wherein the movement detection device simulates or is incorporated into an object of daily living.
The object of daily living may be selected from a group including but not limited to a cup or drinking vessel, a spoon or eating utensil, brush, comb or the like.
The movement sensor typically comprises one or both of an accelerometer and a gyroscope generating motion data representing movement of the device in multiple axes.
In some embodiments, the movement detection device comprises a microcontroller receiving movement data from the motion sensor and optionally, a wireless communication module for wireless transmission of the movement data from the microcontroller to a receiving device. The receiving device may be any suitable device such as, for example, a smart phone, tablet, laptop, desktop computer, or dedicated device receiving the movement data from the motion sensor.
In some embodiments, the movement detection device comprises a canister simulating a cup or drinking vessel, the canister comprising a flexible body portion forming a fluid filled chamber and defining the grasping portion. In some embodiments, the pressure sensor is a differential pressure sensor with a first input in fluid communication with the chamber and a second input in fluid communication with atmospheric pressure. A one-way valve may be provided for releasable coupling with a fluid source to restore fluid pressure in the chamber.
In some embodiments, the canister comprises a rigid base containing a microcontroller and ideally, a wireless communication module.
The movement detection device is typically provided use with the movement monitoring system according to embodiments disclosed herein.
Viewed from another aspect, the present disclosure provides an automated method for objectively quantifying a motor control disorder in a subject, comprising the steps of: (a) receiving at a processor movement data corresponding to movements of a limb of the subject, the movement data comprising at least pressure data and motion data; (b) the processor applying the received movement data to an algorithmic model and identifying one or more features from the movement data that represent disordered movement in the subject; (c) the processor calculating, from the one or more identified features, a score quantifying the motor control disorder in the subject; and (d) the processor generating a display signal causing the calculated score to be presented on a display device.
In some embodiments, the processor applies the received movement data to one or more of: (a) a first algorithmic model to identify a first set of features used by the processor to calculate a selection score which is indicative of presence or absence of the motor control disorder in the subject; (b) a second algorithmic model to identify a second set of features used by the processor to calculate a severity score which is indicative of severity of the motor control disorder in the subject; and (c) a third algorithmic model to identify a third set of features used by the processor to calculate a progression score which is indicative of progression of the motor control disorder in the subject.
It is to be understood that the one or more features identified by the automated method correspond with the features as disclosed previously in relation to the earlier described aspect relating to a system for objectively quantifying a motor control disorder in a subject.
In some embodiments, the method further comprises categorising movement dysfunction in the subject, by the processor calculating a contribution made by each of the first or second set of features to each of a plurality of movement characteristics that are attributable to movement dysfunction in the subject. In some embodiments, the plurality of movement characteristics correlate to clinically accepted descriptions of movement disorder or aspects of movement disorder such as stability, timing, accuracy and rhythmicity.
In some embodiments, the received movement data is obtained from a movement detection device and comprises at least one or both of: pressure data corresponding to pressure applied to the device by the subject; and motion data comprising one or more of position of the limb, acceleration of the limb and angular position of the limb. Ideally, the received movement data is collected while the subject performs a movement task and preferably wherein the movement task is or simulates an activity of daily living. In some embodiments, the motor control disorder may include e.g. ataxia (in its various forms) and/or spasticity.
Viewed from another aspect, the present disclosure provides an automated method for determining an algorithmic model for use in a system or method for objectively quantifying a motor control disorder in a subject comprising the steps of an algorithm designing processor: (a) receiving movement data from a plurality of training subjects who had undergone clinical assessment for the motor control disorder; the movement data comprising pressure data, acceleration data and angular velocity data obtained from a movement detection device while each of the training subjects performed a movement exercise; (b) receiving for each of the plurality of training subjects clinical scores corresponding to the clinical assessment; (c) extracting from the received movement data a plurality of features representing aspects of movement during the movement exercise; (d) selecting from the plurality of features a subset of features using a feature selection algorithm; and (e) using a machine learning approach to build the algorithmic model.
In some embodiments, the feature selection algorithm is Neighbourhood Component Analysis with regularization (NCA-R) although other feature selection algorithms may be utilised.
In one embodiment, the algorithmic model is a first algorithmic model used to calculate a selection score which is indicative of presence or absence of the motor control disorder in the subject, and the machine learning approach is a k-Nearest Neighbour (k-NN) classification model. In some embodiments, the subset of features identified by the feature selection algorithm for the first algorithmic model is selected from the group of subsets comprising: (a) PrRF; (b) PrRF, ACC
In another embodiment, the algorithmic model is a second algorithmic model used to calculate a severity score which is indicative of severity of the motor control disorder in the subject, and the machine learning approach is a random-forest regression model (RFR). In some embodiments, the calculated selection and severity scores may be scaled to correspond to an accepted clinical scale such as FARS, mFARS, 9HPT, BBT, SARA, ICARS, ADL or the like, or some other clinical measure. In some embodiments, the subset of features identified by the feature selection algorithm for the second algorithmic model is selected from a group of comprising PrRF PrM, At and θRFc and preferably comprises the feature set PrM, At and θRFc.
In some embodiments, the pressure data is processed by the algorithm designing processor to identify a compressed phase and an uncompressed phase in each instance of the movement exercise.
In another embodiment, the algorithmic model is a third algorithmic model used to calculate a progression score which is indicative of progression of the motor control disorder in the subject between an initial and a subsequent clinical assessment during which movement data was received, and the feature selection algorithm comprises (i) extracting features of the movement data in one or more of time domain, frequency domain and time-frequency domain; (ii) identifying in the extracted features, relevant features showing a change in value at T2 relative to T1; and (iii) selecting from the relevant features, significant features which show a statistically significant change in value from T1 to T2. The statistical significance in the change in value may be assessed by a Wilcoxon signed-rank test for matched pairs. In some embodiments, the machine learning approach comprises Principal Component Analysis followed by linear regression, optionally with an assumption that the motor control disorder increases in severity between the initial and subsequent assessments. In some embodiments, the subset of features identified by the feature selection algorithm for the third algorithmic model is selected from a group of features comprising MRpr, SRFgyr, MSEIMF
Viewed from another aspect, the present disclosure provides a movement monitoring system for objectively quantifying a motor control disorder in a subject, the system comprising: (a) a movement detection device generating movement data representing movement of a limb of the subject; and (b) an analyser for analysing the movement data. The analyser comprises a processor and a memory containing code which, when executed by the processor: (i) receives the movement data generated by the movement detection device; (ii) applies the received movement data to an algorithmic model stored in the memory and identifies one or more features from the movement data that represent disordered movement by the subject; and (ii) calculates from the one or more identified features a score corresponding to the existence of the motor control disorder in the subject.
In another aspect, the present disclosure provides a movement detection device for use in objectively quantifying a motor control disorder, the movement detection device comprising a plurality of sensors measuring at least motion of the device and pressure applied to the device by the subject, wherein the device represents an object of daily living such as a cup or drinking vessel, a spoon or eating utensil, a brush or comb or the like. Preferably the device has a grasping portion and a pressure sensor of the device generates pressure data representing pressure applied to the gasping portion. Ideally, the device includes one or both of an accelerometer and a gyroscope generating motion data representing movement of the device in multiple axes. The device may also include a wireless communication module for wireless transmission of the movement data from the microcontroller to a receiving device which may be used to process the movement data.
In another aspect, the present disclosure provides for use of the systems, methods, and devices disclosed herein in evaluation, e.g. in a clinical trial, of a therapy for treating a motor control disorder, such as ataxia. Therapies may include medicines or physical therapies, and treatment may involve one or more of mitigating symptoms or reducing or eliminating disease.
It is to be noted that any one of the aspects mentioned above may include any of the features of any of the other aspects mentioned above and may include any of the features of any of the embodiments described below, as appropriate.
It is to be understood each of the various aspects described herein may incorporate features, modifications and alternatives described in the context of one or more other aspects, such as but not limited to the features of the movement monitoring system, the movement detection device, and the automated method for objectively quantifying a motor control disorder in a subject. For efficiency, such features, modifications and alternatives have not been repetitiously disclosed for each and every aspect although one of skill in the art will appreciate that such combinations of features, modifications and alternatives disclosed for some aspects apply similarly for other aspects and are within the scope of and form part of the subject matter of this disclosure.
The present disclosure will now be described in greater detail with reference to the accompanying drawings. It is to be understood that the embodiments shown are examples only and are not to be taken as limiting the scope of the disclosure as defined in the claims appended hereto.
Referring firstly to
The embodiment illustrated in
In some embodiments, movement detection device 200 contains a communications interface, such as a Wi-Fi interface, which provides for wireless real time transmission of movement data collected by movement detection device 200 to a receiving device 320 such as a smart phone or tablet operated by an operator such as the subject or a clinician responsible for the subject's assessment. Receiving device 320 operates as a gateway for transmission of the movement data to analyser 300 for application of the received motion data to a first algorithmic model to identify a first set of features that may be used by a processor of the analyser 300 to calculate a selection score which is indicative of presence or absence of the motor control disorder in the subject. The selection score may be used e.g. in a clinical setting to diagnose a subject presenting with motor disorder symptoms. In some embodiments, the selection score confirms the presence or absence of ataxia or spasticity in the subject.
Alternatively/additionally, movement data received by analyser 300 may be applied to a second algorithmic model to identify a second set of features that may be used by a processor of the analyser to calculate a severity score which is indicative of severity of the motor control disorder in the subject. The severity score may be used in clinical settings to objectively assess the severity of motor disorder symptoms. In preferred embodiments, kinematic (pressure) and kinetic (motion) features extracted from the motion data are used, together with the algorithmic models, to produce scores, such as severity scores, that correlate with clinical scales such as modified-FARS, the Activities of Daily Living (ADL) scale, BBT and 9HPT. Severity scores from objective assessments repeated according to embodiments of the disclosure over time may be plotted or compared to evaluate progression or improvement in motor control disorder symptoms and be used to guide or evaluate therapies and rehabilitation.
Alternatively/additionally, movement data received by analyser 300 may be applied to a third algorithmic model to identify a third set of features that may be used by a processor of the analyser to calculate a progression score which is indicative of extent or rate of progression of disease and may be positive (i.e. advancement of disease) or negative (i.e. improvement e.g. with therapy). The progression score may be used in clinical settings to objectively assess progression of disease and/or effectiveness of therapy and/or to titrate therapy such that e.g. dosage increases are matched to a corresponding rate of progression. The progression score may alternatively/additionally be used in research settings or e.g. in clinical trials to evaluate effectiveness of novel therapies.
In some embodiments, base unit 260 contains most electronic components of the device 200 including the motion sensor which may comprise a 9-axis inertial movement unit (IMU) 210 (e.g. Invensense MPU9250) interfaced with a micro-controller unit 220 (e.g. 32-bit ARM R Cortex-M3) having a Wi-Fi interface (e.g. GS2011), and a rechargeable battery (not shown) equipped with a power management circuit (not shown).
IMU 210 may comprise a triaxial accelerometer and a triaxial gyroscope with a range of +/−2000 degrees per second and optionally, a triaxial geomagnetic sensor. The accelerometer provides information about accelerations in all three directions and the gyroscope provides information about rotations around each axis. To mitigate small errors that build up in each axis over time causing drift in the absolute direction, algorithms utilised by a microprocessor in IMU 210 use the extra magnetic field information from the geomagnetic sensor to compensate for small drifts. Using IMU 210, inertial data is captured as changes in acceleration (±8 g) and angular velocity (±2000/s).
Connector 240A provides for releasable coupling of base unit 260 with a corresponding connector 240B on main body 250 using a friction fit, screw fit or other means. It is to be understood connectors 240A,B as shown are merely examples and that other means for connecting main body 250 and base unit 260 would be known to one of skill in the art.
Main body 250 has rigid cap 251, flexible body 255 defining chamber 256 and rigid wall portion 252 providing on/off switch 253. Rigid cap 251 has sealing lid 254 with air inlet 270, consisting of a one-way valve which permits air to be injected into the chamber 256 to restore the chamber's air pressure should there be leakage over time. Support structure 257 extends through chamber 256 to impart structural integrity to the canister device 200 although it is to be understood that a plurality of support structures or other means could be provided to impart the requisite structural integrity of the canister whilst also achieving the required flexibility and compressibility of flexible body 255.
Flexible body 255 provides a grasping portion for grasping by the subject and defines an air filled chamber 256. Pressure sensor 230 (e.g. MS4525DO DS5AI001DP, TE connectivity) is incorporated into main body 250 with one input port in fluid communication with air chamber 256 and the other port exposed to atmospheric pressure. Pressure sensor 230 transmits the 16 bit differential pressure data as an electrical signal to the micro-controller 220 e.g. via an I2C communication interface or similar. The differential pressure data, being the difference between atmospheric pressure and the varying pressure inside the flexible chamber 256 is transmitted to micro-controller 220. Both kinetic and kinematic data are transferred to receiving device 320 by Wi-Fi transmission.
Movement detection device 200 in this example is provided in the form of a canister of compact size and light weight (e.g. 185 grams) which is comparable to a cup or other drinking vessel typically used by the subject in daily living. In one example, canister device 200 may be manufactured using a multi-material additive manufacturing technique (Objet350 Connex3 3D printer, Stratasys Ltd., USA) wherein inkjet deposited droplets of photopolymer controlled at the minimal length scales 40 μm×80 μm×30 μm×35 μm are UV cured. The soft and hard surfaces of the composite may be printed respectively using VeroCyan (RGD841, shore hardness 85) and Agilus (transparent, shore hardness 30, flexible). The volume of the chamber 256 defined by flexible body 255 is, in the illustrated embodiment, 2.27×107 mm3 (13.8 cubic in) and the pressure of air inside is maintained at −0.05 psi.
For assessment, a subject 300 is asked to sit at a table with forearms pronated and wrists in a natural position resting in front of them with hips and knees flexed to 90° and feet flat on the floor. The device is placed on the table in front of the subject and the subject is instructed to mimic the task of preparing to drink a glass of water using their dominant arm. The phases of movement involved in this preparatory task are illustrated in
Data collected by movement detection device 200 is transmitted wirelessly to receiving device 320 and forwards it on to analyser 300 where it is processed according to methods disclosed herein. It is to be understood, however, that transmission via receiving device 320 need not be wireless and the electronics unit of base unit 260 may be connected by a cable or other connector to the receiving device. In other embodiments, movement detection device 200 contains a SIM card or other means for communicating directly with analyser 300, omitting the intervening step of movement data transmission via receiving device 320. Analyser 300, processes received movement data for the subject according to embodiments disclosed herein, to determine a score such as a selection score or a severity score. In some embodiments the analyser 300 processes received movement data for the subject to determine scores for categories of movement dysfunction in the subject by a processor associated with the analyser calculating a contribution made by each of the first or second set of features to each of a plurality of movement characteristics that are attributable to movement dysfunction in the subject.
42 participants consisting of 20 controls (11 males, 9 females) with an average age of 35.8±13.43 years and twenty two individuals with FA of varying severities (12 males and 10 females) with an average age of 37.05±12.23 years took part in this study. All subjects performed a trial in which upper limb movement data was recorded using the movement detection device 200 described as an Ataxia Instrumented Measure-Canister (AIM-C).
All the participants were assessed using 9HPT and BBT and using clinical scales. All subjects performed a trial in which upper limb movement data was recorded using the movement detection device 200. Individuals with FA also underwent assessment using known clinical scales. The presence of spasticity in the wrist and long finger flexors was identified by clinical testing.
The rotational, translational and pressure data were sampled at 100 Hz and subsequently low pass filtered with a cut-off frequency of 20 Hz. A median filter was used to suppress noise transients in the data. A 6th order Savitzky-Golay filter smoothed data and contributed to minimising the drift effects with averaging. The functional architecture of the overall methodology is given in
IMU 210 measures acceleration and angular velocity from which the velocity and position were derived using known techniques (
A) Complimentary filtering to estimate roll and pitch angles: The Euler angles related to angular measurements uncover distinct movement characteristics intrinsically linked to ataxia. Euler angles were obtained through complimentary filtering techniques while negating the contribution from the gravitational force. The sensor orientation was observed in terms of Euler angle (Roll φ, Pitch θ) variations in the preliminary analysis. The complimentary filter designed to estimate roll {circumflex over (ϕ)}a
With the initial conditions for roll {circumflex over (ϕ)}=ϕgco and pitch {circumflex over (θ)}=θd
{circumflex over (ϕ)}g
{circumflex over (ϕ)}g
the complimentary filter can be stated as,
{circumflex over (ϕ)}c=(1−ω)*{circumflex over (ϕ)}g
{circumflex over (θ)}c=(1−ω)*{circumflex over (θ)}g
where p, q and rare the gyro readings in X, Y and Z axis respectively with the time constant ω. The complimentary filtering technique with ω=0.15 performed better than Madgwick's algorithm and Kalman filtering techniques. Table 1 describes features extracted in this Example.
B) Dimensionless jerk: The stability of the performance of the task was assessed by finding the smoothness Sm, which is log dimensionless jerk where the jerk is defined as the rate of change of acceleration in time t and is given as:
where (t) is the jerk, t1 and t2 are initial and final times, Lm denotes the maximum peak velocity where Lm=max∀t∈[t1,t2] {dot over (x)}(t) and Tp is the time period.
C) Phase of movement-based features: Four distinct phases of movement were identified during the movement task, namely: ‘grasping’, ‘transporting’, ‘releasing’ and ‘stabilising’ phases (
Euclidean distance and time duration of ‘accommodating’ and ‘stabilising’ phases (
D) Pressure data measures: The pressure data revealed fluctuations in the pressure applied to the flexible body 255 of the device 200, which can be considered as two distinct phases within a cycle:
The change in phase was detected using a known change-point detection algorithm (R. Killick, P. Fearnhead, and I. A. Eckley, “Optimal detection of changepoints with a linear computational cost,” Journal of the American Statistical Association, vol. 107, no. 500, pp. 1590-1598, 2012). Accordingly Mab is the mean of a given time series of data S(i)i=ai=b, the residual or cost functional (J) associated with the changing point algorithm for the data series S(i) is defined as:
J(k)=Σi=1k-1(s(i)−M1k)2+Σj=kn(s(j)−Mkn)2 (5)
The standard deviation in the amplitude of the pks (PrSD) was a feature pertaining to the CP phase. Asymmetry ratio feature (PrCP/UP) is the ratio of the duration of compressed CP and un-compressed pressure UP phases. Also, the average time duration of the complete cycle (PrM), i.e. time duration between cp1 and cp3 was considered as a feature of interest.
Of the features described in Table 1, only those that (i) effectively ‘separated’ the individuals with FA from controls and (i) showed correlation with clinical scales indicating ‘severity’ are described.
Separation: The feature separation for the two cohorts was quantified using Area Under the Curve AUC measure from the Receiver Operating Characteristics (ROC) which discriminates between data from the two cohorts using trapezoidal approximation. The AUC value also indicates the extent to which the two cohorts are classified within the estimated confident bounds.
Correlation: The movement and pressure data features were correlated with the commonly used ataxia scales and with scores for spasticity using three measures: Distance correlation Dc, Pearson correlation Pc, and Spearman rank correlation Sc. This was done so that the features with best correlation (combined/uncombined) could be identified and to establish the type of correlation.
Statistical analysis: Wilcoxon-Mann-Whitney signed-rank test was used to determine if the features extracted for the cohorts were statistically significant. G*Power (version 3.1.9.4) software with significance level (p-value) 0.01 and AUC>0.75 was used to determine whether the features held sufficient statistical-power (α) for the study (min. threshold: α≥80%, the effect-size for our sample size (42) is 1.06).
Feature selection techniques: Most of the extracted features could separate controls from individuals with FA with good statistical power (Table 1). To avoid the risk of over estimations when a large number of features are combined, feature selection algorithms were used to select the feature-sets that best estimated severity while separating between controls and individuals with FA. Four feature section algorithms; Sequential Forward Selection (SFS), Random Subset Feature Selection (RSFS), Statistical Dependency (SDe) and Neighbourhood Component Analysis with regularization (NCA-R) were compared. NCA-R was selected because its low cross-validation error resulted in significantly fewer features. NCA-R learn the weights of features by minimising an objective function which measures the average leave-one-out (LOO) classification or regression loss over the training data and fine-tuning a regularisation parameter ‘λ’ to 0.0329 while evaluating the weights of the features such that the significance is indicative of the rank. The features were selected using NCA-R for: (i) selecting subjects as ataxic or not (NCARC) (ii) regression with clinical scores as an estimation of the severity of ataxia (NCARR).
The most significant feature sets (NCARC, NCARR) selected by the NCA-R algorithm were then used in machine learning approaches to build models that could select (i.e. identify subjects as controls or individuals with FA) or compare against clinical scales (i.e. predict severity of ataxia according to clinical scales).
Model prediction for selection: A k-Nearest Neighbour (KNN) classification model was trained using (NCARC) features to classify subjects as individuals with FA or controls based on parameters such as accuracy and AUC from ROC.
Model prediction for comparison to clinical scores: A random-forest regression model (RFR) which uses fine-tree learning approach, was employed for enhancing prediction of clinical scores (mFARS) using NCARR features. Other models such as Random Subset regression, Multiple linear polynomial regression and the like were tested however using regression model parameters such as the goodness of fit measure R-squared (R2), mean absolute error (MAE) and correlation co-efficient Pc, the RFR performance was evaluated as the best performing. The classifier and regression model performance was validated using a 10 fold cross-validation technique in order to minimize misclassification and potential overfitting issues. Table 2 compares performance of the RFR based model prediction of mFARS with 10-fold cross validation (AIM-C) against 9HBT and BBT test.
To categorise movement dysfunction, it is convenient to use known descriptions of the motor deficit of ataxia. Gordon Holmes's widely accepted descriptions can be conveniently categorised into four dimensions as follows: Stability (S), Timing (T), Accuracy (A) and Rhythmicity (R). In the present example, four propositions P1, P2, P3 and P4 described below were designed to enable each of the features identified in the movement data to be assigned one of the S, T, A and R dimensions. The primary axis is the vertical (X-axis) direction of motion. Resultant motion (X, Y, Z combined) of acceleration, and angular velocity were also considered as the primary motion and adheres to either proposition P1 or P4. Since uni-directional temporal variations of motion data were considered, principle P3 was not applicable for the pressure data.
S.T.A.R. dimension assignment: Each selected feature was examined according to the above P1, P2, P3 and P4 propositions and assigned one of the dimensions S, T, A, or R. According to P1, the strength of the repetitiveness or magnitude at resonant frequency along the primary motion of the task movement were identified for the Rhythmicity dimension. According to P2, the features demonstrating deviations in the efficient path required to execute the task were identified for the Accuracy dimension. According to P3, features manifesting in deviations in directions other than the direction of dominant motion were identified as being in the Stability dimension. According to P4, the features which are associated with delay or timing deficit while carrying out the device accommodation or transport phase of the task were identified as being in the Timing dimension.
S.T.A.R contribution of the selected features: The S.T.A.R dimension assignment was applied to the feature subset obtained for both selection (NCARC) and severity estimation (NCARR). Thus, subsequent to the feature selection, features obtained from the NCA-R algorithm were assigned to one of the S.T.A.R. dimensions using propositions P1 to P4 and summed (according to the feature weightings from the relevant model) to determine the collective contribution to each dimension of the S.T.A.R classification.
A preliminary examination reveals that direct measures of movement such as acceleration (Acc), angular velocity (Gyr), pressure (Pr) as well as derived measures velocity (Vel) and Euler angle (Euler) can separate controls from individuals with FA. The mean values for acceleration, angular velocity, Euler angles and velocity were lower in individuals with FA than controls and mean pressure value was higher. The maximum pressure value was also higher when individuals with FA (−0.1 psi) performed the task than control participants (−0.05 psi) as shown in
Diagnostic discrimination based on separation: Initially phase-of-movement based features were assessed as measures for the separating the two cohorts.
The average accommodation time (AT) and the average number of movement units (MU) were greater when individuals with FA performed the task than when controls did (AT is 1 s and MU is 1.23 units respectively for individuals with FA compared to 0.3 s and 0.58 units for controls) (
Frequency characteristics based features such as ACC
Smoothness (Sm) is considered to be a measure of the quality of motor performance, and is measured by dimensionless jerk. It was higher in movements made by controls (mean Sm=−19.7) than those made by people with ataxia (mean Sm=−26.6) from Table 3 and
As noted above, individuals with FA compress the wall of the AIM-C device 200, more than controls (
Individuals with FA could almost be completely distinguished from controls using PrRF, ACC
Correlation assessment with clinical scores: Extracted pressure and motion features correlated significantly with clinical scores (Table 4). The individual feature with the highest correlation with mFARS was PrRF (Sc=0.87) and θRFc (Sc=0.83). A notable finding of this study was that the features PrSD, Prrms, PrMR, PrRF correlated with the presence of spasticity (Pc=0.83) while maintaining significant cohort separation. It is hypothesised that these features may be important in
describing the inability in ‘accommodating’ the AIM-C device 200 in the hand during the movement task.
Features ACC
Since most features separated individuals with FA from controls and provided significant p-values and high correlations (Table 4), the NCA-R feature selection scheme with 10-fold cross validation was used to find the best features in terms of diagnostic and severity prediction as discussed below.
Feature Selection for Classification and Severity Estimation Using NCA-R with Predictive Modelling
Using NCARC subset trained with a KNN classifier, the information and data obtained from the AIM-C device 200 correctly classified control participants (specificity≈1) from individuals with FA (sensitivity≈1) with high diagnostic accuracy (99%). Table 5 shows a comparison of KNN with other classification models with NCARC subset after 10 fold cross-validation, where SVM is Support Vector Machine classifier, RFC is Ensemble Tree classifier, LDA is Linear discriminant classifier. ACC is classification accuracy, TPR (sensitivity) represents: the proportion of individuals with FA identified as having the condition, TNR (specificity) represents the number of controls correctly classified as having no condition, PPV: represents precision, No. represents number of features.
The feature subset (NCARR) selected for severity assessment included PrM, At and θRFc. The subset NCARR when trained with RFR, predicted mFARS (Table 2) with a high goodness of fit parameter (R2=0.92) and superior correlation (Pc=0.96) with mFARS scores. A severity prediction model of 9HPT and BBT provided a lower R2 and correlation value (Table 2).
The power of NCARC and NCARR feature sets were 99.99% with an effect size of 1.182 (Table 1) which was greater than the required statistical power of 80% for the sample size.
Each feature in either the NCARC or NCARR subset can be allocated to one of the four propositions P1, P2, P3 and P4. The feature-weights from the relevant model (KNN or RFR) can then be used to determine the percentage contribution of each feature to the dimensions of the S.T.A.R. descriptions which in turn provide clinical context for describing the manifestation of the movement disorder. Table 6 shows the S.T.A.R dimension assignment for the NCARC and NCARR feature subsets selected using NCA-R algorithm.
For the NCARC subset, the Timing (T) dimension contributed the most (75%) (Table 6 and
For the NCARR subset, the Timing (T) dimension contributed the most again (76%) (Table 6 and
The S.T.A.R. assessment approach has provided objective indications that domains such as timing and stability provide significant contributions to both diagnostic accuracy and severity estimation, and point to differences between ataxic and non-ataxic movement when maintaining the stability of the execution platform.
Thus, embodiments of the disclosure provide an automated method as shown in the schematic illustration of
In preferred embodiments, the received movement data is collected while the subject performs a movement task which is or simulates an activity of daily living such as eating or drinking. The movement data comprises pressure data and motion data obtained from a movement detection device. Typically, the movement detection device comprises a pressure sensor for monitoring pressure applied to the device during a movement task, and motion sensors such as an accelerometer and gyroscope monitoring kinematic motion of the device in multiple axes. In one embodiment, the motion sensors are provided in the form of an IMU.
In a step 903, the processor applies the received movement data to a first algorithmic model to identify a first set of features used by the processor to calculate a selection score which is indicative of presence or absence of the motor control disorder in the subject. Where the motor control disorder is Ataxia, the first set of features may comprise a single feature, ideally a pressure feature such as PrRF. In other embodiments, the first set of features comprises some or all of the features PrRF, ACC
In another step 904, the processor applies the received movement data to a second algorithmic model to identify a second set of features used by the processor to calculate a severity score which is indicative of severity of the motor control disorder in the subject. Where the motor control disorder is Ataxia, the second set of features may comprise a single feature, or a plurality of or all of the features selected from the group comprising: PrRF PrM, At and θRFc. In one embodiment, the feature set comprising PrM, At and θRFc is preferred.
In a step 910 the processor generates a display signal causing the calculated selection and/or severity scores to be presented on a display device such as a screen of a computer, smart phone, tablet or the like, or presented in a report transmitted to an operator by email or other suitable communication means.
In a step 905, the processor categorises movement dysfunction in the subject, by calculating a contribution made by each of the first, second or third (see Example 2) set of features to each of a plurality of movement characteristics that are attributable to movement dysfunction in the subject. In one embodiment, the plurality of movement characteristics are determined according to propositions P1 to P4 disclosed herein and the contribution made by the features to each of the propositions may be aggregated and presented on the display device. In another embodiment, propositions P1 to P4 may mapped to clinically accepted descriptions or categories of movement disorder which, in the case of Ataxia, can be conveniently categorised into four dimensions as follows: Stability (S), Timing (T), Accuracy (A) and Rhythmicity (R). Thus, the contribution of each of the S, T, A and R dimensions in the movement data may be presented in a step 910 on the display device e.g. as percentages and/or rankings or other measures that represent the dominance of each dimension on the movement.
While the example herein demonstrates utility in objectively detecting the presence of and scoring the severity of ataxia, it is to be understood that the embodiments disclosed also have utility in assessing other aspects of motor control disorder such as spasticity. In some embodiments, one or more or all of the features PrSD, PrRMS, PrMR, PrRF as described in Table 1 can be used to identify the presence/absence of spasticity.
Ten individuals diagnosed with FRDA (mean standard deviation (SD) age=37.05±12.23, 4 males and 6 females) were assessed at two time points (T1 and T2) that were on average 24 weeks apart. At T1, the average Functional Staging of Ataxia score of participants was 4.8 (SD 0.42), indicating moderate disability (walking requires an assistive device) and 68.3 (SD 10.01) on the mFARS. Table 7 shows clinical features of the cohort at T1 and T2 (GC: Group characteristics, M: mean, SD: standard deviation, R: Range, n: number of participants (n=10), age is given in years, no. of males: 4, females: 6; right hand dominant: 10, * is age of onset in years, + is disease duration in years).
4-5.5
Additionally, all participants performed a task that simulated drinking from the AIM-C device 200 while movement data was recorded by sensors in the canister as described previously. Recordings were made at T1 (first visit) and T2 (second visit) approximately 24 weeks later. At each visit, clinical assessments including Functional Staging of Ataxia score (Staging), mFARS, 9HPT and BBT were administered specifically. The upper limb subscale from the mFARS (NeuroUL) was used in the analysis. Other clinical parameters were collected, included GAA1 repeat size (GAA repeat size on the smallest allele) and GAA2 repeat size (GAA repeat size on the largest allele), age of disease onset (the age that symptoms related to FRDA were first observed), and disease duration (difference between current age and age of disease onset).
Using Wilcoxin-signed rank statistical test (with 1−αs=0.01) on a sample size of 10 subjects, the effect-size (d) of 1.59 was obtained after feature selection with a power Pα≥80%. Hence, the sample size of 10 was considered adequate for this study.
The kinematic measures obtained from the IMU were acceleration (acc) and angular velocity (gyr) while the kinetic measure obtained from the pressure sensor was the differential pressure (pr) representing the grip pressure exerted on the AIM-C device 200. Both measures were sampled at 100 Hz and low pass filtered with a cut-off frequency of 20 Hz. A median filter was used to suppress noise transients in the data. A 6th order Savitzky-Golay filter smoothed data and contributed to minimising the drift effects with averaging.
Spatio-temporal features of kinematic measurements were obtained using a complimentary filtering technique based on Madgwick's algorithm. This technique extracted Euler angles (roll φ), pitch (θ) and yaw (Φ) as ‘deduced’ measures with the gravitational contributions removed while converting the sensor frame measurements to the global frame. The velocity of movement was a ‘deduced’ measure obtained from acceleration data after the bias removal.
The features were considered in the time, frequency, and time-frequency domains. Twenty eight features of interest were extracted using the following described signal processing techniques to obtain time domain features, frequency domain features, and time-frequency domain features.
A) Time Domain Features: Time domain features from the kinematic and kinetic data included statistical parameters such as mean and standard deviation. Range of Motion (ROM) of the AIM-C device 200 as it was transported from the table to the lip was obtained as the difference between highest and lowest position of the pitch (elevation) angle for the time duration TN over five cycles.
B) Frequency Domain Features: In the frequency domain analysis, the spectral leakages were reduced using Blackman-Harris windowing technique to improve the algorithmic performance when capturing features such as resonant frequency (RF) and magnitude at resonance (MR). The frequency domain features such as RF and MR were extracted for both kinematic (in all three orthogonal X, Y and Z axes combined) and kinetic measures.
C) Time-frequency domain features: Only the linear and periodic behaviour of the data was extracted by frequency domain analysis. However it is possible that disease progression is related to nonlinear and non-stationary behaviour of the kinetic and kinematic measures. Hence, the Hilbert-Huang Transform (HHT) which characterises non-linear behaviour by ‘time-frequency’ domain variations was used to probe the movement data.
The HHT combines two methods: First, an Empirical Mode Decomposition (EMD) to decompose the time series into Intrinsic Mode Functions (IMF). EMD consisted of a sifting process of the original signal, in which ‘n’ IMFs were obtained for each of the measures. Second, the Hilbert transform was used to compute instantaneous amplitude and frequency and present the IMF's in an energy-time-frequency domain representation. The application of HHT algorithm on EMD as described by Huang et al. is given as follows.
Consider n=6 IMFs (Ck(t), . . . , Cn(t)∀k∈[1, 2, . . . n]) and a residue signal rn(t). The decomposed signal of each of the measures (acc, gyr, pr) denoted by E(t) is given by:
The Hilbert transform Hk(t) of each IMF Component Ck(t) is:
The analytic signal is then defined as:
R
k(t)=Ck(t)+jHk(t)=αk(t)ejØ(t) (9)
where the amplitude αk(t)=√{square root over (Ck(t)2+Hk(t)2)} and the phase
Corresponding to this, the instantaneous frequency can be given as:
The Hilbert spectrum representing amplitude on the time-frequency domain is given by:
H(ω,t)=Re[Σk=1nαk(t)ej∫ω
Features extracted in time-frequency domain. The Mean Square Energy (MSE) of the kth IMF is given by,
where T is the sampling time, N is the length of the series and k is the IMF index. MSEk∀k∈[1, 2, . . . , 6] is considered to be a feature that might represent changes in disease progression.
A z×q matrix was obtained as the HHT spectrum, where z corresponds to frequency sampling points, and q the time sampling points. In combination with HHT, a Singular Value Decomposition (SVD) technique was further engaged for extracting features (this approach is hereinafter referred to as SvHT). The singular values of HHT spectrum were extracted using SVD, where SVD decomposes the matrix into three: one diagonal matrix and two orthogonal matrices. The SVD of the z×q matrix A is given as:
A=U√{square root over (λ)}VT, (13)
where U and V are orthogonal matrices of size z×z and q×q respectively, and √{square root over (λ)} is a diagonal matrix of singular values. The first ten singular values i.e. √{square root over (λ)}=[Sv1, Sv2, Sv3 . . . Sv10]T of the HHT matrix were extracted as features where the initial values contained most of the matrix information and can be related to the energy of the measure. For clarity, notational abbreviations of extracted features using the techniques described (FFT, EMD, HHT) are given in Table 8, noting the three direct measures utilised in this Example were acc, gyr and pr.
A) Feature selection and combination for AIM-C score for progression. The most important features contributing to disease progression were identified based on two premises. Firstly, the feature indicating change in value at T2 is given by
where FiT1, FiT2 are feature values at T1 and T2 respectively ∀i∈[1, 2, . . . Q] (Q is the number of features extracted). Secondly, the statistical significance of the change from T1 to T2 is assessed by the Wilcoxon signed-rank test for matched pairs since not all data sets were normally distributed. The Wilcoxon signed rank test found the difference between each set of matched pairs (T1 and T2) in terms of error rate as and power Pα.
The features adhering to these two premises were selected as the best features. Amalgamation of these features was used to generate a score that represents the overall disease progression. This was achieved through Principal Component Analysis (PCA) followed by linear regression, collectively denoted as Principal Component Regression [43]. The dimensionality reduction of all combined features after normalisation was accomplished through the Principal Component Analysis (PCA), using visual observation of data distributions. Using linear orthogonal transformation, the data is transformed into a new coordinate system where the diagonal covariance matrix maximizes feature variance and then projects the data in accordance with the variance of the distributions. The Principal Component Regression algorithm consisted of passing the first principal component (representing the maximum variance of the feature contribution) to a linear regression to generate the AIM-C score for progression.
B) Comparison of AIM-C score for progression and clinical parameters used to measure performance overtime (mFARS, BBT 9HPT). The sensitivity of AIM-C score for progression to detect changes in scores at T2 compared to T1 was statistically quantified using the effect-size (d) (using G*Power software). The capacity of data from the AIM-C device 200 to capture change over time was compared to the mFARS, BBT and 9HPT. The differences between multiple groups (mFARS, BBT, 9HPT and AIM-C scores) were compared using the Holm-Bonferroni multi comparison post-hoc method to compare the (T1 and T2) scores that exhibited significant difference to each other. This method was also robust against type-I errors by controlling the family-wise error rate (FWER) at 0.05 for Holm-correction.
C) Feature correlation with clinical parameters. All the kinetic and kinematic features extracted were also correlated with the clinical measures (mFARS score, the NeuroUL score, the Functional Staging of Ataxia score, ADL score, 9HPT and BBT) and clinical parameters (GAA1, GAA2, disease duration and age at disease onset) at T1 and T2 using Spearman's rank correlation (PcT1 or T2). The functional architecture of the overall methodology for Example 2 is provided in
Features Extracted to Generate AIM-C Score and Comparison with Clinical Scores
A) Extracted features: The MRpr feature value obtained from the frequency analysis of all participants increased from timepoint T1 to T2. This is shown in
In the time-frequency domain, the mean square energy (MSE) of all the 6 IMFs was calculated using the empirical mode decomposition of the three measures (acc, gyr and pr). The MSE of IMF2 of each subject obtained from acceleration (MSEIMF
The initial singular values of the HHT matrix contained most of the matrix information and were considered features related to the energy. The first three singular values for Hilbert spectra of absolute acceleration were significantly higher at T2 than at T1 (Table 11). The remaining singular values did not demonstrate noticeable change from T1 to T2. In particular, the first singular value of the Hilbert spectra of acceleration (Sv1−HTacc) for all participants was higher at T2 than at T1. This is illustrated in
B) Comparison of clinical score versus AIM-C score for progression: The features whose values at T2 had significantly changed from their values at T1 are shown in
that the AIM-C score for progression is more sensitive to disease progression than the standard clinical tests.
C) Performance of AIM-C score for progression compared to clinical scores over time: The sensitivity of the AIM-C score for progression from T1 to T2 was expressed in terms of effect size (d) as shown in Table 12, where d is a dependent p-value obtained using Wilcoxon signed rank test for matched pairs (*d0.8 and d0.9 refer to the effect size for achieving 80% and 90% power (Pa) respectively). The effect-size for each test was such that AIM-C (1.59)>NHPT (1.23)>BBT (0.633)>mFARS (0.00) when then power (Pα) was set at 80%.
The scores were further compared using multivariate analysis by post-hoc tests using Holm-Bonferonni method for non-normal data distributions (FWER+0.05, alphaSi=0.01, alphaBonf=0.008). Table 13 presents the comparison of scores and shows significant difference between AIM-C scores for progression and other clinical scales (p-value=0.0 indicates highly significant values). The performance of AIM-C score for progression was significantly different (p=0.000) from the other clinical scores in T1 and T2. The BBT score was also significantly different from the other scores. The 9HPT score was not significantly different (p=0.104) to the mFARS score in capturing the change in upper limb function over time but noticeably different to the AIM-C and BBT scores.
D) Correlation between T1 and T2 values of extracted features and clinical scores. Clinical measures (mFARS score, Neuro UL, Functional Staging of Ataxia score and the ADL score) and fixed clinical parameters (GAA1 and GAA2) significantly correlated with the kinematic and kinetic features as shown in
The time domain feature ROMθ, highly correlated (PT1=0.99, PT2=0.7) with the Functional Staging of Ataxia. The mean angular velocity (gyrm) also correlated with the Functional Staging of Ataxia Score at T1.
The frequency domain feature MRvel correlated highest with GAA1 (PT1=0.84, PT2=0.70), the Functional Staging of Ataxia (PT1=0.91, PT2=0.82) and the mFARS score (PT1=0.85, PT2=0.72). RFpr highly correlated (PT1>0.75) with the clinical scores (GAA1, GAA2, ADL) at T1, but not at T2. The kinematic feature SRFQcc correlated strongest (PT1=0.78, PT2=0.80) at T1 and T2 with the ADL scores. Most of the features correlated (PCT1>0.8) with the Functional Staging of Ataxia.
The highest correlation (PCT1=0.76, PCT2=0.82) obtained for (NeuroUL) upper limb scores was with the Sv1-HTpr feature. Sv2-HTpr also correlated the highest (PCT1=0.83, PCT2=0.95) with traditional measures of upper limb function (the 9HPT and BBT). The feature correlation obtained for disease duration was not significant (Pc<0.5) at T2. For detailed reporting, the statistical significance (p<0.01) of features that significantly correlated with the clinical scores are provided in Table 14.
Example 2 demonstrates that there were extractable features whose values significantly changed (α<0.05) over the duration of this study for most participants and that when these features were represented in the AIM-C score for progression, their statistical power was greater than conventional clinical tests. This finding implies that there is information in the kinetic and kinematic data that is relevant to FRDA and in particularly to its progression which cannot be determined using clinical tests alone, pointing to the utility of instrumented assessments as disclosed herein.
Relevant features identified in Example 2 include the magnitude of pressure (MRpr) applied to the wall of the AIM-C device 200 while it was being held. People with ataxia have difficulty maintaining a constant grip force either as an impairment in predicting appropriate forces required for the grip and/or as a strategy for increasing mean grip force to avoid dropping objects during the execution of the task. As disease progresses, it is plausible that variability in grip force and the compensatory mean grip force may both increase, explaining why MRpr was useful as a feature for monitoring disease progression.
Another feature identified in Example 2 as providing information about disease progression included the second peak of angular acceleration after application of FFT (SRFgyr). This feature reveals rotational motion in axes not related to the primary axis of motion that may have arisen from dysmetria as well as instability at more proximal joints. Again it is plausible that this may have been due to increasing ataxia, making this feature a useful indicator of progression.
The Mean Square Energy of the second Intrinsic Mode Function obtained from acceleration (MSEIMF
The first singular values of the HHT spectrum of acceleration (Sv1−HTgyr) was chosen as the fourth feature that showed significant difference between T1 and T2. Sv1−HTgyr contained most of the time-frequency information, highlighting the inherent non-linear characteristics of the motion of the vessel during the task which are also likely to be more prominent as ataxia worsens.
Although the four features (MRpr, SRFgyr, MSEIMF
It is to be noted, however, that a larger sample of subjects with FRDA could yield a slightly different set of candidate features for the AIM-C score for progression. It is also possible that some of the four candidate features identified in Example 2 carry shared information (e.g. SRFgyr and energy (Sv1-HTgyr)) although it is likely that such overlap has been addressed to at least some extent by PCA. Nevertheless, a larger population would permit a more detailed analysis with possibility of a more refined feature set to avoid overfitting. One of skill in the art could readily utilise the teachings of this disclosure in evaluating a larger dataset to determine such a refined feature set which is explicitly stated to be within the scope of this disclosure.
Irrespective of the dataset size utilised to identify the features of interest in Example 2, it is to be noted that the AIM-C score for progression has been demonstrated to have a greater effect size than the three clinical scores (mFARS, 9HPT and BBT). Therefore, the AIM-C score for progression disclosed herein is likely a more sensitive test for measuring disease progression than currently available clinical alternatives. A particular advantage of an instrumented form of assessment of disease progression with enhanced sensitivity as disclosed herein is that disease modifying therapies with small effects are more likely to be detected by the instrumented methods and systems disclosed herein, than by other measures, and therapies with larger effects could be identified in studies with a smaller sample size or shorter duration. This offers the potential to provide accurate, objective evaluation of novel therapies for ataxia, and to provide a tool to determine response to therapies which may be utilised in dosage titration.
The present disclosure provides an instrumented approach to measuring ataxia in the form of a device, such as a pressure-measuring canister, capable of sensing certain kinetic and kinematic parameters of interest to quantify the impairment levels of participants particularly when engaged in an activity that is closely associated with daily living. In particular, the functional task of simulated preparation for drinking can be utilised to capture characteristic features of disability manifestation in terms of diagnosis (separation of individuals with FA and controls) and severity assessment of individuals diagnosed with the debilitating condition of FA. Time and frequency domain analysis of these features enables the classification of individuals with FA and control subjects to reach an accuracy of 99% and a correlation level reaching 96% with the clinical scores. The instrumented approach to measuring ataxia also permits scoring of disease progression in an objective, repeatable manner that gives greater sensitivity than currently available clinical assessments of disease progression. Measuring disease progression is a useful proxy for testing the ability to measure disease modification. Because existing scales may not necessarily map onto functional capacity of daily activities, they may also neglect to capture the effect of emerging disease modifying therapies on factors impairing daily life. To date, longitudinal progression of FRDA using kinetic and kinematic sensing while simulating self-drinking has not been considered.
Unlike the 9HPT, which cannot be performed by severely impaired participants, assessments according to embodiments of the present disclosure can be performed by people who have significant motor deficits demonstrating the potential to add value to existing assessments available for use in clinical trials. In addition to the separation, severity and progression scoring capabilities, characterisation of movement dysfunction based on Holmesian dimensions (Stability, Timing delay, Accuracy, and Rhythmicity) provides clinically meaningful descriptions to the practitioner. For example, for the movement task of preparing for drinking, S.T.A.R. categorisation demonstrated ataxia manifestation predominantly in timing dimension. Indeed, using techniques disclosed herein, the ability to unravel and quantify the disability information is evident. Given the significant correlation with the clinical assessment scales, the present disclosure provides the basis for an objective assessment tool with simple and unique physical attributes vital for the use in the motor impairment management arena. Additionally, the device component is robust, compact and light weight and in embodiments where it resembles an object of daily living, it is readily taken up by subjects for assessment. Indeed these attributes and the relevance to ADL, suggest assessment in non-clinical settings allowing more regular and accurate testing to facilitate use in clinical trials as well as monitoring in homes and in community based rehabilitation programs.
Where the terms “comprise”, “comprises”, “comprised” or “comprising” are used in this specification (including the claims) they are to be interpreted as specifying the presence of the stated features, integers, steps or components, but not precluding the presence of one or more other features, integers, steps or components or group thereof.
It is to be understood that various modifications, additions and/or alterations may be made to the parts previously described without departing from the ambit of the present disclosure as defined in the claims appended hereto.
It is to be understood that the following claims are provided by way of example only, and are not intended to limit the scope of what may be claimed in the current or any future application. Features may be added to or omitted from the claims at a later date so as to further define or re-define the invention or inventions.
Number | Date | Country | Kind |
---|---|---|---|
2020902051 | Jun 2020 | AU | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/AU2021/050637 | 6/18/2021 | WO |