 Research
 Open Access
 Published:
Estimating upperextremity function from kinematics in stroke patients following goaloriented computerbased training
Journal of NeuroEngineering and Rehabilitation volume 18, Article number: 186 (2021)
Abstract
Introduction
After a stroke, a wide range of deficits can occur with varying onset latencies. As a result, assessing impairment and recovery are enormous challenges in neurorehabilitation. Although several clinical scales are generally accepted, they are timeconsuming, show high interrater variability, have low ecological validity, and are vulnerable to biases introduced by compensatory movements and action modifications. Alternative methods need to be developed for efficient and objective assessment. In this study, we explore the potential of computerbased body tracking systems and classification tools to estimate the motor impairment of the more affected arm in stroke patients.
Methods
We present a method for estimating clinical scores from movement parameters that are extracted from kinematic data recorded during unsupervised computerbased rehabilitation sessions. We identify a number of kinematic descriptors that characterise the patients’ hemiparesis (e.g., movement smoothness, work area), we implement a doublenoise model and perform a multivariate regression using clinical data from 98 stroke patients who completed a total of 191 sessions with RGS.
Results
Our results reveal a new digital biomarker of arm function, the Total GoalDirected Movement (TGDM), which relates to the patients work area during the execution of goaloriented reaching movements. The model’s performance to estimate FMUE scores reaches an accuracy of \(R^2\): 0.38 with an error (\(\sigma\): 12.8). Next, we evaluate its reliability (\(r=0.89\) for testretest), longitudinal external validity (\(95\%\) true positive rate), sensitivity, and generalisation to other tasks that involve planar reaching movements (\(R^2\): 0.39). The model achieves comparable accuracy also for the Chedoke Arm and Hand Activity Inventory (\(R^2\): 0.40) and Barthel Index (\(R^2\): 0.35).
Conclusions
Our results highlight the clinical value of kinematic data collected during unsupervised goaloriented motor training with the RGS combined with data science techniques, and provide new insight into factors underlying recovery and its biomarkers.
Introduction
Stroke is the second major cause of death and disability worldwide, with about 15 million new cases every year [1]. Onethird of these cases lead to persistent cognitive and motor disabilities [2]. About 80% of stroke survivors present weakness and partial loss of voluntary control in the upperextremities [3], or hemiparesis, which is often associated with other sensorimotor alterations, such as hypertonia or tremor.
Although hemiparesis is a highly prevalent symptom and severely limits the independence of affected patients, its causes and recovery dynamics are not fully understood [4]. Recent literature converges on the idea that recovery is mainly due to a combination of residual corticospinal tract capacity and an upregulation of the reticulospinal tract [5, 6]. Further, recovery seems to follow a temporal structure where most of the improvement occurs during the first months poststroke [7, 8]. So far the assessment of the hemiparesis phenotype and its progression, however, are based on assessment methods with known limitations (e.g. FuglMeyer Assessment [9, 10], Action Research Arm Test [11, 12]) and there is a need for more sensitive, objective, and reliable alternatives that are also compatible with contemporary digital health technologies.
A recent systematic review [13] of a total of 225 studies (N = 6197) using 151 different kinematic metrics found that kinematic assessments of upper limb sensorimotor function are poorly standardised and rarely measure clinimetrics in an unbiased manner. Specifically, using descriptors of accuracy, efficacy, efficiency, movement planning, precision, spatial posture, speed, temporal posture, and range of movement together with clinimetric properties of these descriptors (i.e., reliability, measurement error, convergent validity, and external validity), the authors showed that the studies analysed exclusively focused on finding correlations between measures of impairment, and only two of the studies reported correlations in change. Overall, there is very limited information regarding testretest reliability and the external validity of the change of kinematic outcome measures of reaching performance [14]. Exceptionally, Murphy et al. [15] explored external validity of the change in a number of kinematic descriptors and found a significant covariation of the Action Research Arm Test (ARAT) scores with movement time (\(R^2\) = 0.36), smoothness (\(R^2\) = 0.31), and trunk displacement (\(R^2\) = 0.35). Although the results are promising, this study involves a limited number of subjects (N = 24) from a highly homogeneous sample (i.e., acute patients only). Further, the ARAT clinical scale presents poor robustness to compensation and is especially vulnerable to the use of explicit strategies to improve performance. Majeed et al. [16] explored the application of models based on LASSO regression to predict changes in motor ability (FMUE) and motor function (Wolf Motor Function Test, WMFT). These models proposed that recovery in both scales can be approximated by the patient’s age, the patient’s motor control during the execution of fast movements, and other demographic and clinical features, altogether accounting for 65% and 86% of the variability for the FMUE and WMFT scales respectively. Although these models reached exceptional accuracy, their utility is limited because they make use of kinematic data obtained during the execution of very specific pointing movements supervised by clinicians and/or researchers, and are based on generic unbounded linear models, with the consequence that their estimated values could be largely outside the meaningful range of the scale.
We propose a new approach towards using kinematic data obtained in unsupervised rehabilitation sessions to estimate the level of impairment and functional recovery. Data is obtained from patients engaging with goaloriented embodied individualised training with the Rehabilitation Gaming System (RGS) [17, 18]. The RGS combines the paradigm of action execution with that of observation of the corresponding movement in Virtual Reality (VR), this goal is achieved by having the patients perform tasks from a firstperson perspective, where the movement of their limbs are captured by a camera or a depth sensor (i.e. Microsoft Kinect) and mapped to an analogous virtual representation on a computer screen. RGS includes individualisation mechanisms to adjust the difficulty of the task to the capabilities of the patient, contextual restrictions, and explicit and implicit feedback.
We first explore the potential of hand movements collected during unsupervised RGS sessions to characterise hemiparesis in stroke patients. Secondly, we build and analyse the performance (i.e., testretest reliability, validity, sensitivity, and generalisation) [13] of a model for estimating impairment and recovery scores captured by standardised clinical scales.
Methods
Subjects
Our retrospective analysis uses data of 191 RGS sessions from 98 individuals with hemiparesis (age in [23, 87], mean 63; days poststroke in [5, 3045], mean 400; Table 1) who were recruited between 2010 and 2015 to participate in studies conducted in Barcelona and Tarragona, Spain [18,19,20]. Participants met the following inclusion criteria: (1) ischemic strokes (middle cerebral artery territory) or hemorrhagic strokes (intracerebral); (2) mildtomoderate upper limb hemiparesis (Medical Research Council scale for proximal muscles \(>2\)) after a firstever stroke; (3) age between 20 and 90 years old; and (4) the absence of any significant cognitive impairment (MiniMental State Evaluation \(>22\)).
Protocol
Participants followed a rehabilitation protocol including 3–5 weekly sessions of 30 minutes each for 3–12 weeks using the RGS in Fig. 1. The RGS consists of a PC, a 17 inch LCD touch display, an imagebased motion capture device positioned on top of the screen [18]. The virtual tasks logic and graphics are implemented using the Torque 3D and Unity 3D game engines. The joint movements of the user’s head, shoulders and elbows are tracked and mapped onto an avatar through a biomechanical model using a customdeveloped visionbased motion capture system. Arm movements are displayed on a screen from a firstperson perspective, realising a rehabilitation paradigm that combines goaloriented embodied action execution and observation.
For the RGS sessions in the main dataset (Table 1) the participants are instructed to intercept virtual spheres by executing horizontal bimanual movements over the surface of a table (‘Spheroids’ protocol, cf. inset in Fig. 1). The task parameters (the frequency of sphere appearance, their speed, their range, and size) are combined in a single parameter (‘difficulty level’) and automatically adjusted during the session in order to maintain the user’s performance between 70 and 80% success rate [17, 21]. The system allows for the storage and extraction of performance parameters as well as hand path trajectories derived from joints’ positions and rotations recorded at about 100 Hz.
During the rehabilitation patients are evaluated using standard clinical scales: FuglMeyer Assessment for the upper extremity (FMUE), Chedoke Arm and Hand Activity Inventory (CAHAI) and Barthel Index (BI). When collecting the 191 samples (Table 1), the following measures are taken to improve data quality:

The clinical score measurements (FMUE, CAHAI, BI) are coupled to the RGS session closest in time, with a maximum time separation of 4 days between the measurement and the RGS session;

The first two RGS sessions of a patient at the start of the rehabilitation trajectory are excluded to ensure that patients are familiar with the RGS environment for all collected samples.
Outcomes and analysis
To analyse the potential of RGSderived movement descriptors for capturing both impairment and recovery in standardised clinical scales, we first extract a set of variables that are known to correlate with the severity of hemiparesis [13, 16]. Next, we define a model that combines the information of several variables to estimate the patient’s score on the FMUE clinical scale. The model includes an error estimate. We use repeated crossvalidation to avoid overfitting, allocating 50% of the samples to the training set and 50% to the validation set. We study the model evaluating its convergent validity, testretest reliability, external validity, and generalisation to other tasks and clinical scales.
Analysis of convergent validity: estimating impairment
To explore the convergent validity of RGSderived kinematic descriptors in comparison to standardised clinical assessments, we compute Pearson correlations between all the variables (the RGSderived descriptors and the baseline characteristics) and the clinical scores. By comparing these to a randomised outcome distribution we identify a threshold of \(r \simeq 0.081\) for all the Pearson correlations variablescores (Fig. 12 in Appendix). Subsequently, we adopt a nonlinear parametric model to combine information from several variables for estimating the patients’ impairment level (i.e., FMUE scale). We use only the patient’s baseline characteristics and RGSderived movement descriptors (extracted from a single session logfile) as predictors. In particular, we rule out the two variables sessions completed so far and time since stroke and functions of them. In this way, the estimate is intended to assess the clinical status of the patient at a given moment without knowledge of the rehabilitation history. To avoid overfitting, we perform repeated crossvalidation with 50% of samples for training and 50% for validation, obtaining the optimal active set of variables possible for our dataset. In the crossvalidation procedure, we define the accuracy of the model as the percentage of samples that are correctly estimated above or below the median score (e.g. 47 for FMUE). Additionally, we report the values of the Leaveoneout crossvalidation (LOOCV) for the optimal active set. To quantify the performance of the model, we compare true and estimated FMUE scores reporting the average error, the Pearson correlation, and the coefficient of determination. Finally, we report the performance of the estimation of FMUE obtained by a linear model with the same active set of variables.
Analysis of test–retest reliability
To evaluate the testretest reliability, we consider two unseen datasets each composed of 921 RGS sessions, for a total of 1842 unseen RGS sessions. Each session in the first dataset ‘test’ is associated with a session in the second dataset ‘retest’ and obtained from the same patient within less than 48 h. The small time frame makes it plausible that the clinical state of the patient is unchanged between the two test–retest sessions, and so they can be used to assess the reliability of the regression models. These data were collected in the same trials as the main dataset, but they correspond to rehabilitation sessions for which we do not have an associated measurement of a clinical scale (so they cannot be used for training the model). To quantify the reliability in the estimation of the FMUE scores, we evaluate the intraclass correlation coefficient (ICC) between the estimations obtained by the model for test and retest sessions.
Analysis of external validity: estimating the change in impairment
Starting from the original dataset (Spheroids, Table 1), we design a new dataset (recovery dataset, Table 2) composed of 54 samples where each sample represents a couple of sessions of the same patient for which the time lapsed between the two is larger than 16 days. Next, we compute the correlations between the change in the movement descriptors and clinical improvements. We utilise this dataset to analyse the external validity of the previous model (obtained for estimating impairment) in detecting changes in the clinical status of the same patient. Therefore, we adopt the same set of variables used for the estimate of impairment. We evaluate the performance of the model by comparing the true and predicted change of FMUE score. Next, we identify 38 out of the 54 samples from the recovery dataset (Table 2) for which the associated \(\Delta {\text {FMUE}}\) values exceed an MDC of 4 points. We then determine the sensitivity of the model evaluating the percentage of \(\Delta {\text {FMUE}}\) that are correctly predicted above the MDC.
Analysis of task generalisation
The Spheroids protocol requires predominantly movements in the lateral (left/right) direction. To explore the potential of the model to generalise to other tasks that involve bimanual 2D planar reaching movements, we identify a second dataset of 37 samples from a previous study in which 19 subacute stroke patients with hemiparesis trained with a different RGSbased activity which is a variation on the wellknown arcade game ‘WhacAMole’ [18, 22]. The observed FMUE scores in this dataset are in the range [5, 60], with a mean of 37 and an SD of 14. The gameplay of WhacAMole requires movements on the full 2d plane so, besides the movement descriptors used for Spheroids, we define independent descriptors for the movements in the front/back direction. Subsequently, we adopt the nonlinear parametric model to combine information from different variables to estimate the impairment level as measured by the FMUE scale. To quantify the performance of the model, we compare true and estimated FMUE scores reporting the average error, the Pearson correlation, and the coefficient of determination.
Analysis of generalisation to other clinical scales
To explore the potential of the model to generalise for the estimation of impairment captured by other standardised clinical approaches, we considered the estimation of the CAHAI and BI scales that are available in the Spheroids dataset (Table 1). We, therefore, follow the same steps as above to analyse the convergent validity of the model in estimating CAHAI and BI scores. It is useful to anticipate here that the results for CAHAI are close to those for FMUE. This similarity is expected as the two scales have a high relative correlations and comparable correlations to most of the variables, cf. Table 1. To be quantitative, we can consider the case in which the FMUE score \(S_i^{\text {FM}}\) of the ‘Spheroids’ dataset is estimated by simply rescaling the corresponding CAHAI value \(S_i^{\text {CAHAI}}\) by 66/91; this leads to the standard error
and an \(R^2\) value of \(1  \sigma _{\text {FM,CAHAI}}^2 / \sigma ^2_{\text {FM}} \simeq 0.62\). Estimating BI scores from kinematic descriptors is instead generally harder. This can be explained by the lower correlation of BI scores with most of the variables (Table 1). If we estimate the FMUE score by a simple rescaling of the BI value by 66/100 we get a standard error of
and a \(R^2\) value of \(1  \sigma _{\text {FM,BI}}^2 / \sigma ^2_{\text {FM}} \simeq 0.12\). So we see again that BI carries different information than FMUE (or CAHAI). The above values give us a natural benchmark to assess the performance of the model estimation of the clinical scales.
We conducted these analyses using inhouse software (SaddlePoint Signature v2.9.3).
Identification of variables
Given the technical constraints of the setup (e.g., sampling rate) and the nature of the training protocol (i.e., execution of horizontal reaching movements), we extracted all taskrelevant kinematic measures of arm use that we could identify in the literature [13] . In total we obtained 31 variables, 12 firstorder variables and 19 secondorder variables obtained as functions of the firstorder ones.
We identify two groups of firstorder variables: (1) demographic and physiological data at recruitment, variables 1–5 in Table 1, and (2) kinematic descriptors extracted for the more affected limb during training sessions, variables 6–12 in Table 1. For the evaluation of all kinematic descriptors, the first and last two minutes of each training session are discarded to avoid the interference of behaviours or events related to the start and the ending of the training session (i.e., revision of instructions, postural adjustments, exposure to the final score screen, etc.).
Secondorder variables include chronicity (i.e. acute, subacute, and chronic categories) and the difference between the less and the more affected upperlimb in each of the quantitative firstorder variables, as well as their logarithmic transformations. The descriptive statistics of secondorder variables are given in Table 5 in the Appendix.
Among the abovementioned variables, most are wellknown [13]. The work area is computed as the area of the complex hull of the hand movements using standard methods (Jarvis’ Algorithm [23]). The distance covered refers to the total length of the hand paths. The performance success rate is defined as the ratio of spheres intercepted over the total. However, we introduce also the new descriptors ‘Smoothness’ and ‘Totalgoal directed movement’ (TGDM). Specifically, to extract information on UE motor function we introduce an original kinematic descriptor, \(J(\sigma )\), to assess the patient’s movements at a specific temporal resolution, \(\sigma\). This metric allows us to isolate goaloriented movements from the hand trajectory in a certain direction, assumed to be stored over time as a function f(t). For the main dataset, we consider the left/right direction, as it is the principal axes in the movement dynamics of the ‘Spheroids’ protocol. We assume measurements are taken at discrete time points \(t_i = i\Delta\) for \(i = 0,1,\ldots ,T  1\), with \(\Delta\) being the timestep (for the ‘Spheroids’ dataset we have \(\Delta \simeq 0.01\) s). We define the total hand displacement during goaloriented movements \(J(\sigma )\) as the difference between the actual movements and a smoothed version of the discrete movements. The smoothed hand path \(f_\sigma (t)\) is obtained using a Gaussian smoothing process
where the parameter \(\sigma\) defines how smooth the new trajectory will be, see an example in Fig. 2. Therefore \(J(\sigma )\) is obtained as
Following this analysis method, we derive the two new variables corresponding to the value of \(J(\sigma )\) at the two peaks in the \(\sigma\)dependent Pearson correlation with the clinical scales, Fig. 3: ‘Smoothness’ in correspondence to the highfrequency peak, and ’TGDM’ in correspondence to the lowfrequency peak. The location of the two peaks is weakly dependent on the clinical scale considered, yet it appears to be related to the data structure: the highfrequency peak is linked to the time resolution of the data (\(\Delta \simeq 0.01\) s), while the lowfrequency peak is related to the typical timescale of the Spheroids protocol (i.e. a set of spheres is launched every \(\sim 10\) s).
Models description
To combine variables for estimating clinical scores of impairment and recovery, we introduce a model that allows for the presence of noise on both the variables \({\varvec{Z}}\) and the score S, and we hence name it a doublenoise parametric model. Its generative functional form is
where \(\varvec{\theta } = \{\varvec{\beta }, \beta _0, A, B, \sigma _1, \sigma _2 \}\) are the \(p+5\) model hyperparameters to be inferred: \(\varvec{\beta }\) (association parameters of p active variables), \(\beta _0\) (parametric offset), A and B (range offsets), \(\sigma _1\) and \(\sigma _2\) (noise strengths). The sources of noise u (the covariate noise) and v (the score noise) are both assumed to be standard normally distributed. Since \(\tanh (x) =  \tanh (x)\), we remove the resulting parameter sign ambiguity by enforcing \(B \ge A\). Note that the saturation of the sigmoidal function captures the limits of the clinical scores so that the average over the noise of the estimated score S is constrained in the interval [A, B].
The probability of a particular score S, given the variables and the hyperparameters, is given by
with the shorthands \(D v = (2 \pi )^{\frac{1}{2}} \mathrm{e}^{\frac{1}{2} v^2 } dv\), \(a=(BA)/2\), \(b=(B+A)/2\). We use the following (improper) prior distribution over the parameters \(\varvec{\theta }\):
We take \(p(\sigma _1) \sim \mathrm{e}^{q/\sigma 1}\) and similar for \(\sigma _2\) with q being a very small number, typically of the order of the accuracy of the numerical work, e.g. \(q \sim 10^{10}\). These priors guarantee that the loglikelihood function is bounded from below.
We adopt Maximum A Posteriori (MAP) inference: given the dataset \({\mathcal {D}}=\{({\varvec{Z}}_1,S_1), \ldots , ({\varvec{Z}}_n,S_n)\}\), the optimal parameters \(\varvec{\theta }\) correspond to the maximum of the posterior probability or, alternatively, to the minimum of the regularised loglikelihood function:
The errors on the inferred parameters \(\varvec{\theta }\) are estimated from the curvature of the regularised loglikelihood function at the minimum.
Two simpler models derived from Eq. 8 can be considered corresponding to having either score noise only or covariate noise only:

Score noise model, \(\sigma _1 = 0\). Taking the limit \(\sigma _1\rightarrow 0\) in Eq. 8 gives
$$\begin{aligned} \Omega _{\text {sco}}(\varvec{\theta })&=\frac{1}{2\sigma _2^2} \sum _{i = 1}^n \left[ S_i  b  a \tanh (\varvec{\beta } \cdot {\varvec{Z}}_i + \beta _0) \right] ^2 \nonumber \\&\quad + n \log (\sigma _2 ) + \frac{1}{2} d \varvec{\beta }^2 +\frac{q}{\sigma _2} , \end{aligned}$$(9) 
Covariate noise model, \(\sigma _2 = 0\) Taking the limit \(\sigma _2\rightarrow 0\) in Eq. 8 gives
$$\begin{aligned} \Omega _{\text {cov}}(\varvec{\theta })&=\frac{1}{2\sigma _1^2}\sum _{i=1}^n \Big [ \tanh ^{1}\Big (\frac{S_i  b}{a}\Big )  \varvec{\beta } \cdot {\varvec{Z}}_i  \beta _0 \Big ] ^2 \nonumber \\&\quad + \sum _{i=1}^n \log \big [a^2  (S_i  b)^2\big ] +n\log \Big (\frac{\sigma _1}{a}\Big ) \nonumber \\&\quad +\frac{1}{2} d \varvec{\beta }^2 +\frac{q}{\sigma _1} \end{aligned}$$(10)provided \(S_i  b<a\) for all i, otherwise \(\Omega _{\text {cov}}(\varvec{\theta }) = \infty\). This implies that for each b the minimisation over a is to be carried out strictly over the open interval \(a>\mathrm{max}_i S_i  b\).
Results
Identification of features
Several of the kinematic descriptors, in particular TGDM and the difficulty level reached (difficulty) correlate highly with all clinical scales, cf. Table 1 and Fig. 4. Generally, there is a high level of consistency between correlations of the variables with the different clinical scales examined. The main exceptions are age and the time since stroke (dayss). The former does not display a relevant correlation with FMUE and CAHAI scores but it is negatively correlated with BI (\(r=0.30\), \(p < 0.0001\)), while the latter correlates negatively with the FMUE (\(r=0.25\), \(p=0.00049\)) and CAHAI (\(r=0.17\), \(p=0.019\)) scores but positively with BI (\(r=0.22\), \(p=0.0033\)). The correlation between FMUE and CAHAI scores is high (\(r=0.89\), \(p<0.0001\)), while the consistency with BI is significantly lower: (\(r=0.34\), \(p<0.0001\) with FMUE; \(r=0.5\), \(p<0.0001\) with CAHAI). Generally, the correlations of the kinematic descriptors are higher with FMUE and CAHAI than with BI.
All kinematic variables show consistent intervariables correlation, cf. Fig. 4. In particular, the maximum speed achieved during reaching (maxsp) correlates highly with the work area (warea) (\(r=0.76\), \(p<0.0001\)). The variable age correlates negatively with TGDM (\(r=0.31\), \(p < 0.0001\)) and difficulty (\(r=0.24\), \(p = .00098\)), while it is uncorrelated with FMUE and CAHAI scores, but not to BI (\(r=0.30\), \(p < 0.0001\)). This suggests that agerelated technology proficiency may affect the kinematic descriptors. Nevertheless, this effect is relatively weak, i.e., the correlations of difficulty and TGDM with the clinical scores are significantly higher than the ones with age.
To better understand the meaning of the two variables obtained with the smoothing techniques (smoothness and TGDM), we also extract finite timewindowed variables for comparison (64 s for the range variable; 14 s for the speed). Specifically, we compute the maximum range of movement (left/right direction) within overlapping time windows of size \(\sigma\), where \(\sigma\) is again fixed by the condition of maximum correlation with the clinical scale of interest, and we average the measurements across all possible windows along the whole RGS session. The resulting values show a very high correlation with the TGDM variable (\(r=0.98\), \(p<0.01\)) and show very similar Pearson coefficients with the FMUE (\(r=0.54\), \(p<0.01\)), CAHAI (\(r=0.57\), \(p<0.01\)) and BI (\(r=0.44\), \(p<0.01\)) clinical scales. These results suggest that the TGDM is capturing information about the typical range of movement associated with the scenario events occurring within relevant time windows (i.e. about 10 seconds in the Spheroids scenario). Following the same method, we extract timewindowed maximum reaching speed. The resulting values show a very high correlation with the smoothness variable (\(r=0.87\), \(p<0.01\)) together with very similar Pearson coefficients with the FMUE (\(r=0.42\), \(p<0.01\)), CAHAI (\(r=0.39\), \(p<0.01\)) and BI (\(r=0.21\), \(p<0.01\)) clinical scales. These results suggest that smoothness is linked to the ability of the patient to perform fast movements to complete the RGS tasks.
Convergent validity: estimating impairment
We found that the three models implemented (double noise model, score noise model, covariate noise model) offer comparable performance in estimating the FMUE scale on the dataset in Table 1. The typical error in the final estimate is of order \(\sim 10\), while the inferred noise score error \(\sigma _2\) is typically \(\sim 0.1\). This result supports that the score noise has a low impact, therefore we select the covariate noise model to decrease the number of parameters to be inferred.
The most relevant variables for the estimate of the FMUE score are difficulty and Diff. TGDM (difference in TGDM between the nonparetic arm and the paretic arm). The total number of active variables is 6 and they are shown in Table 3. We have enforced the presence of Diff. distance covered in the active variable set since it is relevant in the estimate of clinical change, Tables 2 and 6. Note that the resulting active variable set does not contain the patient’s baseline variables and the estimation of scores on clinical scales are solely obtained from RGSderived data. This also means that estimates made by this model for different RGS sessions of the same patient are considered independent measurements. In total, we use 10 parameters (6 association parameters and 4 hyperparameters) inferred from the 191 sessions.
The hyperparameters of the model that predicts FMUE are given by the following mean values and standard deviations: \(a=32.72 (0.85)\), \(b=34.55 (0.73)\), \(\sigma _1=0.551 (0.043)\) and \(\beta _0 = 0.370 (0.051)\). Using LOOCV, we obtain that for the active variable set Table 3 the accuracy on the training set is 0.755 while the accuracy on the validation set is 0.746.
Eventually, we compare the true and predicted score values, shown in Fig. 5. The model predicts the FMUE score with an average error of \(E_\text {FM} \sim 12.8\), Pearson r truepredicted of 0.63, and a coefficient of determination \(R^2 = 0.38\). These values are intermediate between an estimation based on simple rescaling of CAHAI and BI scores, Eqs. 1,2.
To further exemplify the features of the above model, we compare with the simple linear model \(S = \varvec{\beta } \cdot {\varvec{Z}} + \beta _0 + \sigma u\) with the same covariates as in Table 3. The inferred values of \(\beta\) for normalised variables are in this case (same order as in Table 3) 2.32(0.71), 1.86(0.81), \(0.84 (0.72)\), \(2.58 (0.74)\), 1.27(0.75), 1.46(0.79) with \(\beta _0 = 42.98 (0.66)\) and \(\sigma = 12.90 (0.93)\). The linear model predicts the FMUE scores with an average error of \(E_\text {FM} \sim 17.3\), Pearson r truepredicted of 0.60, and a coefficient of determination \(R^2 = 0.14\). Interestingly, the range of estimated FM scores in the dataset is now [16, 84], with 46 scores (over 191) estimated above 66. This result exemplifies the problem of mapping unbounded variables on a finite range and illustrates the need for a nonlinear function as in Eq. 5.
Robustness: test–retest reliability
In Fig. 6 we show the FMUE estimations obtained by the model for the test and retest sets. The ICC between the ‘test’ and ‘retest’ is 0.89, which is at the high end of the ‘Good’ agreement measure according to the Koo and Li guidelines [24]. In Fig. 6 we also show the interval defined by the standard error of the regression \(E \simeq 12.7\). We measure an average retest error \(\sqrt{\sum _i (S_i^\text {test}S_i^\text {retest})^2/N}\) equal to 5.9. Finally, we estimate a minimally detectable change (MDC) [25] of 11.6 points.
These results support the internal consistency of this assessment method to estimate the clinical scores of the patients.
External validity: estimating the change in impairment
We observe that the Pearson correlation between change in FMUE (\(\Delta {\text {FMUE}}\)) and change in CAHAI (\(\Delta \text {CAHAI}\)) is \(r=0.68\), Pearson between \(\Delta {\text {FMUE}}\) and change in BI (\(\Delta \text {BI}\)) is \(r=0.67\), while the Pearson correlation between \(\Delta \text {CAHAI}\) and \(\Delta \text {BI}\) is \(r=0.72\).
For each sample, we consider now as variables the change (between the two sessions) of the original variables. By comparing it to a randomised outcome distribution, we identify a threshold of \(r \sim 0.14\) for the Pearson correlations variablescore. Several of the variables correlate highly with the change in all three clinical scores, cf. Tables 2 and 6 in the Appendix. The highest correlated variable is Change in TGDM (\(r=0.48\), \(p=0.00024\) with \(\Delta {\text {FMUE}}\); \(r=0.55\), \(p<0.0001\) with \(\Delta \text {CAHAI}\); \(r=0.49\), \(p=0.0011\) with \(\Delta \text {BI}\)). Note that, in comparison to estimating a single session’s score, the variables age and chronic correlate more with the outcome when predicting the score change. The initial scores (the clinical scores at the first session) have a high correlation with change because of ceiling effects. In Fig. 7 we show the correlogram of all the variables and clinical scales. We observe that generally the correlations between kinematic descriptors in a single session (Fig. 4) are preserved also when considering the change between sessions; for example, the highest intervariables correlation is for Change in warea and Change in maxsp at (\(r=0.87\), \(p<0.0001\)).
The association parameters of the model that estimates \(\Delta {{\text {FMUE}}}\) are shown in Table 3. The corresponding model hyperparameters are \(a=23.3 (3.3)\), \(b=20.9 (3.2)\), \(\sigma _1=0.430 (0.060)\), and \(\beta _0 = 0.68 (0.12)\).
We compare the true and predicted \(\Delta {\text {FMUE}}\) in Fig. 8. The Pearson correlation between true \(\Delta {\text {FMUE}}\) and predicted \(\Delta {\text {FMUE}}\) is 0.76. The value of the coefficient of determination \(R^2\) is 0.57. These results show that the model has a good external validity to clinical change with a precision comparable to the one obtained for the crosssectional FMUE score.
Sensitivity: estimating recovery
We obtain that the model correctly predicts recovery in 36 over 38 of the cases with \(\Delta {\text {FMUE}} \ge 4\), indicating a TPR of 95%.
Task generalisation
For the task generalisation analysis, we consider a dataset composed of 37 RGS sessions from 19 hemiparetic participants that trained in a rehabilitation protocol derived from WhacAMole [18].
Unlike the Spheroids protocol, the gameplay of WhacAMole requires movements on the full 2d plane. In response, we utilise the smoothing technique in both cardinal axes of the task. i.e. front/back and left/right directions. Pearson correlations between the clinical scales and the variable \(J(\sigma )\) reveal a similar pattern to the one observed in the Spheroids scenario, with a peak of the Pearson coefficients at about 1s corresponding to the variable TGDM in each direction (Fig. 13 in Appendix). The location of the main peak is again close to the typical timescale of the protocol. For the FMUE score, the highest Pearson coefficient is observed in the frontal direction (\(r=0.54\) for \(\sigma = 1.3 s\)); the lateral hand displacement peak is (\(r=0.50\) at \(\sigma = 1.1 s\)).
When predicting clinical scales, we use now only 2 active variables in order to limit overfitting: the variables TGDMfb (Totalgoal directed movement for front/back direction) and TGDMlr (Totalgoal directed movement for left/right direction). We then infer 6 parameters (2 association parameters + 4 hyperparameters) from the 37 RGS sessions. The two association parameters are (for normalised variables) \(\beta _\text {TGDMfb} = 0.15(0.13)\) and \(\beta _\text {TGDMlr} = 0.18(0.14)\). The hyperparameters of the model that predicts FMUE for the ‘WhacAMole’ scenario are given by \(a=31.9 (3.5)\), \(b=31.0 (2.5)\), \(\sigma _1=0.49 (0.11)\), and \(\beta _0 = 0.21 (0.13)\).
The FMUE estimates are shown against the true values in Fig. 9. The Pearson correlation between true FMUE and predicted FMUE is 0.63, the average error is \(E \sim 11.2\), and the value of the coefficient of determination \(R^2\) is 0.39.
Generalisation to CAHAI and BI
Overall, the CAHAI scale has similar properties to FMUE in relation to the kinematic descriptors, Fig. 4. To stress the generalisation potential of the model, we can then adopt the same model introduced in Table 3 for estimating the FMUE and CAHAI scores. The association parameters for CAHAI are reported in Table 4. The most important variables are difficulty and TGDM. The hyperparameters of the covariate noise model that predicts CAHAI scores are \(a=39.158 (0.099)\), \(b=51.962 (0.079)\), \(\sigma _1=0.953 (0.064)\), and \(\beta _0 = 0.0319 (0.071)\). The predicted scores are plotted against the true CAHAI values in Fig. 10. The model predicts the CAHAI score with an average error of \(E_\text {CAHAI} \sim 20.1\), Pearson r truepredicted of 0.66, and a coefficient of determination \(R^2 = 0.40\). This accuracy is close to what we obtained for estimating the FMUE score, Fig. 5.
In Fig. 11 we compare FMUE and CAHAI, both for the true scores and the predicted scores. We note that the relationship between FM and CAHAI is generally well preserved in the estimates; for example, the Pearson between FMUE and CAHAI scores is \(r=0.89\) for true values and \(r=0.88\) for estimates. The fact that the variability in the true FMUE vs true CAHAI is seemingly comparable to the one in the model, reinforces the idea that the precision we achieve is similar to the one of estimating FMUE directly from CAHAI, as we estimated using Eq. 1.
Finally, we observe that the model that predicts FMUE and CAHAI scores does not work well for the BI. Most kinematics variables have a significantly smaller correlation with the BI (in particular work area, TGDM, smoothness) while baseline information and clinical history of the patient are comparatively more relevant (for example the patient’s age, Table 1 and Fig. 4). The active variable set for BI is composed of 5 variables (\(\beta\) values for normalised variables): age (\(\beta =0.21(0.15)\)), sessions completed so far (\(\beta =0.24(0.19)\)), difficulty (\(\beta =1.21(0.72)\)), Log. time since stroke (\(\beta =0.14(0.14)\)), Log. Difficulty (\(\beta =0.91(0.70)\)). Using a double noise model, Eq. 8, we infer the hyperparameters \(a=51.56 (0.10)\), \(b=49.22 (0.12)\), \(\sigma _1=0.5948 (0.0072)\), \(\sigma _2=0.178 (0.027)\), and \(\beta _0 = 1.025 (0.011)\). Comparing true and predicted BI scores, we measure an average standard error of \(E_{\text {BI}} \sim 16.8\), a Pearson correlation trueprediction of 0.62, and a coefficient of determination \(R^2\) of 0.35. This accuracy is comparable to the one achieved by the models for FMUE and CAHAI scores. Nevertheless, the dataset Table 1 is very unbalanced towards high BI scores (mean score 80, with only 3 samples with a score below 25), so may not generalise well to homogeneous unseen BI data (i.e., the precision for low scores may be relatively poor).
Discussion and conclusion
Our understanding of poststroke motor recovery depends on our capacity to evaluate and characterise impairment and disability. Current standardised assessment methods are mostly subjective and present relevant unsystematic variability due to differences in the evaluators’ training, lack of systematicity in the administration of the assessments, and often are excessively focused on one single aspect of the impairment and/or disability.
Different rehabilitation approaches show a preference for using (and even targeting) specific assessment methods for the evaluation of their therapeutic efficacy, and often these methods have been developed by the same team of authors. For example, the effectiveness of ConstraintInduced Movement Therapy [26] is usually evaluated using the Wolf Motor Function Test [27] and the Motor Analog Scale [28], while the effectivity of occupational therapy has been frequently assessed using the Barthel Index [29] and the Functional Independence Measure [30]. There is an urgent need to establish alternative methods for a common evaluation protocol and characterisation of the hemiparesis phenotype, thus allowing us to identify specific impairment features that could advance our understanding of the recovery dynamics and guide the design of effective rehabilitation therapies. In pursuing this objective, we have conducted a careful analysis of the kinematic data from the upperextremities of 191 individuals with poststroke hemiparesis, and we have constructed a model for estimating impairment and recovery. Our results reveal a new digital biomarker of upperlimbs motor impairment, the Total GoalDirected Movement (TGDM), which relates to the patients’ range of motion during the execution of meaningful goaloriented reaching movements. The TGDM strongly correlates with the level of impairment captured by the FMUE and the level of disability captured by the CAHAI, and also carries relevant information about the patients’ progress, showing a high correlations with the magnitudes of improvement and deterioration estimated by both scales. The model presents high external validity of impairment estimates (\(R^2\): 0.38), robustness (test–retest reliability) (ICC: 0.89), external validity of recovery estimates (\(R^2\): 0.57), sensitivity (TPR: 95%) and task generalisation (\(R^2\): 0.39). Despite the high heterogeneity of the sample and the high level of noise of the selected kinematic parameters, the model’s accuracy to estimate the FMUE is comparable to other standardised clinical scales, such as CAHAI. These results are especially interesting given the currently limited evidence about the external validity of the change of kinematic outcome measures of reaching performance in people with hemiparesis after stroke [31]. According to a recent systematic review on the clinimetric properties of kinematic upper limb assessments [13], only two papers captured external validity of change (i.e., ability to capture longitudinal changes in the measured construct), and just nine parameters showed enough evidence to estimate recovery (i.e., number of velocity peaks, trunk displacement, task/movement time). The quality of evidence however was very low for all metrics. Further, our results outperform previous methods to estimate the level of impairment [15, 16] in two fundamental aspects: (1) it shows robustness to compensation and is resistant to using explicit strategies to boost performance, and (2) it exclusively relies on metrics collected under unsupervised rehabilitation sessions. Although current recommendations point out that wearables with integrated Inertial Measurement Units and visionbased tracking systems are insufficient to measure the quality of movement and improvement in motor function, our findings, together with the growing evidence supporting distance travelled as an accurate and responsive digital biomarker of recovery [32], suggest the opposite and advocate for further studies to clarify the limitations of these estimations, in particular regarding the responsiveness.
Although these results support that our method is able to capture improved sensorimotor control during gameplay (i.e., motor synergies), we cannot discard the contribution of taskrelated learning processes to the overall change measured. In addition, it is important to notice that the RGS setup did not provide antigravity support, thus the different capabilities of the participants to elevate the elbow and wrist may explain part of the unsystematic variability we observe in the model’s performance.
The relevance of our results is emphasised by their consistency across clinimetric properties and by their generalisation potential, relying on a large and heterogeneous dataset of patients at different stages poststroke. Notice however that its generalisation potential has been evaluated in similar setups using visionbased systems (camera or depthsensorbased) tracking the joints of the upper limbs during the performance of planar reaching arm movements. We believe that the applicability of the TGDM and its derived models to evaluate impairment and motor recovery is promising for a number of reasons: (1) it can be derived from unimanual displacements executed in the horizontal plane, (2) it does generalise to other tasks involving twodimensional horizontal reaching movements towards targets, and (3) it can be estimated during unsupervised motor training. Our results provide an example of how digital biomarkers of motor deficits derived from behavioural data could guide the design of automated assessment platforms for a continuous and remote monitoring of impairment.
Availability of data and materials
The datasets used and analysed in this study are available from the corresponding author on reasonable request.
Abbreviations
 RGS:

Rehabilitation gaming system
 FMUE:

FuglMeyer Assessment for Upper Extremities
 CAHAI:

Chedoke Arm and Hand Activity Inventory
 BI:

Barthel Index
 SD:

Standard deviation
 ICC:

Intraclass correlation coefficient
 UE:

Upper extremities
 TGDM:

Totalgoal directed movement
 maxsp:
 perchit:
 distcov:
 warea:
 numsess:
 dayss:
 pareticdom:
 MDC:

Minimal detectable change
 TPR:

True positive rate
References
Thrift AG, Thayabaranathan T, Howard G, Howard VJ, Rothwell PM, Feigin VL, Norrving B, Donnan GA, Cadilhac DA. Global stroke statistics. Int J Stroke. 2017;12(1):13–32.
Shah S, Vanclay F, Cooper B. Improving the sensitivity of the barthel index for stroke rehabilitation. J Clin Epidemiol. 1989;42(8):703–9.
Lai SM, Studenski S, Duncan PW, Perera S. Persisting consequences of stroke measured by the stroke impact scale. Stroke. 2002;33(7):1840–4.
Krakauer JW, Cortés JC. A nontaskoriented approach based on highdose playful movement exploration for rehabilitation of the upper limb early after stroke: a proposal. NeuroRehabilitation. 2018;43(1):31–40.
Asboth L, Friedli L, Beauparlant J, MartinezGonzalez C, Anil S, Rey E, Baud L, Pidpruzhnykova G, Anderson MA, Shkorbatova P, et al. Corticoreticulospinal circuit reorganization enables functional recovery after severe spinal cord contusion. Nat Neurosci. 2018;21(4):576–88.
Makin T, Diedrichsen J, Krakauer J. Reorganisation in adult primate sensorimotor cortex: does it really happen. The cognitive neurosciences. 6th ed. Cambridge: MIT Press; 2020.
McDonnell M, Hillier S. Vestibular rehabilitation for unilateral peripheral vestibular dysfunction. Cochrane Database Syst Rev. 2015. https://doi.org/10.1002/14651858.CD005397.pub4.
Ballester BR, Duff A, Maier M, Cameirao M, Bermudez S, Duarte E, Cuxart A, Rodríguez S, Verschure PF. Revealing an extended critical window of recovery poststroke. bioRxiv. 2018.
FuglMeyer AR, Jääskö L, Norlin V. The poststroke hemiplegic patient. II. Incidence, mortality, and vocational return in göteborg, Sweden with a review of the literature. Scand J Rehabil Med. 1975;7(2):73–83.
Gladstone DJ, Danells CJ, Black SE. The FuglMeyer assessment of motor recovery after stroke: a critical review of its measurement properties. Neurorehabil Neural Repair. 2002;16(3):232–40.
Carroll, D.: A quantitative test of upper extremity function. Journal of Chronic Diseases 18(5), 479–491 (1965). DOI: 10.1016/00219681(65)900305.
Lang CE, Wagner JM, Dromerick AW, Edwards DF. Measurement of upperextremity function early after stroke: properties of the action research arm test. Arch Phys Med Rehabil. 2006;87(12):1605–10.
Schwarz, A., Kanzler, C.M., Lambercy, O., Luft, A.R., Veerbeek, J.M.: Systematic review on kinematic assessments of upper limb movements after stroke. Stroke 50(3), 718–727 (2019).
Wegner, C., Filippi, M., Korteweg, T., Beckmann, C., Ciccarelli, O., De Stefano, N., Enzinger, C., Fazekas, F., Agosta, F., Gass, A., et al.: Relating functional changes during hand movement to clinical parameters in patients with multiple sclerosis in a multicentre fmri study. European journal of neurology 15(2), 113–122 (2008).
Murphy, M.A., Willén, C., Sunnerhagen, K.S.: Kinematic variables quantifying upperextremity performance after stroke during reaching and drinking from a glass. Neurorehabilitation and neural repair 25(1), 71–80 (2011).
Abdel Majeed, Y., Awadalla, S.S., Patton, J.L.: Regression techniques employing feature selection to predict clinical outcomes in stroke. PloS one 13(10), 0205639 (2018).
Cameirão, M.S., Badia, S.B., Oller, E.D., Verschure, P.F.: Neurorehabilitation using the virtual reality based rehabilitation gaming system: methodology, design, psychometrics, usability and validation. Journal of neuroengineering and rehabilitation 7(1), 1–14 (2010)
Ballester, B.R., Maier, M., Mozo, R.M.S.S., Castañeda, V., Duff, A., Verschure, P.F.: Counteracting learned nonuse in chronic stroke patients with reinforcementinduced movement therapy. Journal of neuroengineering and rehabilitation 13(1), 1–15 (2016).
da Silva Cameirão M, Bermudez i Badia S, Duarte E, Verschure PF. Virtual reality based rehabilitation speeds up functional recovery of the upper extremities after stroke: a randomized controlled pilot study in the acute phase of stroke using the rehabilitation gaming system. Restorative neurology and neuroscience. 2011;29(5):287–98.
Ballester BR, Maier M, Domingo DA, Aguilar A, Mura A, Pareja LT, Esteve MFG, Verschure PFMJ. Adaptive vrbased rehabilitation to prevent deterioration in adults with cerebral palsy. In: 2019 International Conference on Virtual Rehabilitation (ICVR); 2019. pp. 1–7. https://doi.org/10.1109/ICVR46560.2019.8994754
Nirme J, Duff A, Verschure PF. Adaptive rehabilitation gaming system: online individualization of stroke rehabilitation. In: 2011 Annual International Conference of the IEEE Engineering in Medicine and Biology Society; 2011. pp. 6749–6752. IEEE
Han, C.E., Kim, S., Chen, S., Lai, Y.H., Lee, J.Y., Osu, R., Schweighofer, N.: Quantifying arm nonuse in individuals poststroke. Neurorehabilitation and neural repair 27(5), 439–447 (2003).
Jarvis, R.A.: On the identification of the convex hull of a finite set of points in the plane. Information Processing Letters 2(1), 18–21 (1973). DOI: 10.1016/00200190(73)900203.
Koo, T.K., Li, M.Y.: A guideline of selecting and reporting intraclass correlation coefficients for reliability research. Journal of chiropractic medicine 15(2), 155–163 (2016).
de Vet, H.C., Terwee, C.B., Knol, D.L., Bouter, L.M.: When to use agreement versus reliability measures. Journal of clinical epidemiology 59(10), 1033–1039 (2006).
Wolf SL, Winstein CJ, Miller JP, Taub E, Uswatte G, Morris D, Giuliani C, Light KE, NicholsLarsen D, et al. Effect of constraintinduced movement therapy on upper extremity function 3 to 9 months after stroke: the excite randomized clinical trial. JAMA. 2006;296(17):2095–104.
Wolf, S.L., Catlin, P.A., Ellis, M., Archer, A.L., Morgan, B., Piacentino, A.: Assessing wolf motor function test as outcome measure for research in patients after stroke. Stroke 32(7), 1635–1639 (2001).
Taub E, Morris DM, Crago J, King DK, Bowman M, Bryson C, Bishop S, Pearson S, Shaw SE. Wolf motor function test (wmft) manual. Birmingham: University of Alabama, CI Therapy Research Group; 2011.
Mahoney FI, Barthel DW. Functional evaluation: the barthel index: a simple index of independence useful in scoring improvement in the rehabilitation of the chronically ill. Maryland State Med J. 1965.
Ottenbacher KJ, Hsu Y, Granger CV, Fiedler RC. The reliability of the functional independence measure: a quantitative review. Arch Phys Med Rehabil. 1996;77(12):1226–32.
Veerbeek, J.M., LangbroekAmersfoort, A.C., Van Wegen, E.E., Meskers, C.G., Kwakkel, G.: Effects of robotassisted therapy for the upper limb after stroke: a systematic review and metaanalysis. Neurorehabilitation and neural repair 31(2), 107–121 (2017).
Wagner JM, Rhodes JA, Patten C. Reproducibility and minimal detectable change of threedimensional kinematic analysis of reaching tasks in people with hemiparesis after stroke. Phys Ther. 2008;88(5):652–63.
Funding
This study was supported by the European Commission (EC), the European Research Council under Grant agreement 341196 (cDAC), EC H2020 project socSMCs (H2020EU.1.2.2. 641321), and by the RGS@home, EIT Health project 19277. EIT Health is supported by EIT, a body of the European Union.
Author information
Authors and Affiliations
Contributions
Contributors BR, FA, and AC planned the analysis. BR, MM, and PV contributed to data acquisition. FA and AC analysed the data. All authors wrote the draft paper. PV submitted the study and is responsible for the overall content as guarantor. All authors have read and approve the final version of the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
All research on human subjects reported in this manuscript was prospectively approved by the Institutional Review Board of Hospital Vall d’Hebron, Hospital del Mar i l’Esperança from Barcelona, and Hospital Joan XXIII from Tarragona and all participants provided written informed consent.
Consent for publication
Not applicable.
Competing interests
P. F. M. J. Verschure leads the research group SPECS that developed RGS, and is the CEO/founder of the spinoff company Eodyne Systems, SL, which commercialises RGS with the goal to achieve a largescale distribution of lowcost sciencebased rehabilitation technologies.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix
Appendix
In Table 5 we report the descriptive statistics of the secondary variables (functions of primary variables in Table 1) for dataset Spheroids.
The distribution of the Pearson correlation coefficients of the FMUE score with all the variables in Tables 1 and 5 is shown in Fig. 12. This is compared with the distribution obtained with the randomised outcome, whose standard deviation is \(r \simeq 0.081\).
In Table 6 we report the descriptive statistics of the secondary variables (functions of primary variables in Table 2) for the recovery dataset (change of variables and clinical scales between two RGS sessions of the same patient with a delay of at least 16 days).
In Fig. 13 we show the Pearson’s r between \(J(\sigma )\) Eq. 4 and the clinical scales as a function of the timescale \(\sigma\) for the database WhacAMole, Sec. Task Generalisation. This is a dataset composed of 37 RGS sessions with a limited range of the clinical scales: FMUE range of observations [5, 60], mean 37, SD 14; CAHAI range of observations [7, 49], mean 32, SD 15; BI range of observations [48, 100], mean 85, SD 13. The WhacAMole scenario requires movements in all 2D place (unlike Spheroids scenario, Table 1), so we define two independent \(J(\sigma )\): one associated with the front/back trajectory and one associated with the lateral (left/right) trajectory. In Fig. 13 we show that in both directions there is a peak in the Pearson’s r at about \(\sigma = 1\) s for all three clinical scales. In correspondence with these two peaks, we define two variables: TGDMfb (Totalgoal directed movement in front/back direction) and TGDMlr (Totalgoal directed movement in the leftright direction). We stress that the timescale of the peak is roughly equivalent to the timescale of the gameplay of the WhacAMole scenario. Finally, we note that for the same analysis done in the Spheroids scenario we observe two clear peaks in the lateral direction (Fig. 3). In WhacAMole instead, the highfrequency peak is not clearly visible. One factor that may affect this difference is that the gameplay of WhacAMole is faster than Spheroids (roughly 1 s instead of 10 s) so that the separation between the two potential peaks is smaller. Other factors that may influence the absence of the second peak are the fact that the gameplay is 2D (so that the speed in one direction is less informative) or the inadequate timeresolution of the camera (different from the one used for Spheroids).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Ballester, B.R., Antenucci, F., Maier, M. et al. Estimating upperextremity function from kinematics in stroke patients following goaloriented computerbased training. J NeuroEngineering Rehabil 18, 186 (2021). https://doi.org/10.1186/s12984021009718
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12984021009718
Keywords
 Rehabilitation
 Stroke
 Interactive feedback
 Upper extremities
 Posture monitoring
 Motion sensing
 Motion classification
 Multivariate regression