Evaluation of a spring-loaded upper-limb exoskeleton in cleaning activities.

In the past few years, companies have started considering the adoption of upper-limb occupational exoskeletons as a solution to reduce the health and cost issues associated with work-related shoulder overuse injuries. Most of the previous research studies have evaluated the efficacy of these devices in laboratories by measuring the reduction in muscle exertion resulting from device use in stereotyped tasks and controlled conditions. However, to date, uncertainties exist about generalizing laboratory results to more realistic conditions of use. The current study aims to investigate the in-field efficacy (through electromyography and perceived exertion), usability, and acceptance of a commercial spring-loaded upper-limb exoskeleton in cleaning job activities. The operators were required to maintain prolonged overhead postures while holding and moving a pole equipped with tools for window and ceiling cleaning. Compared to the normal working condition, the exoskeleton significantly reduced the total shoulder muscle activity (∼17%), the activity of the anterior deltoid (∼26%), medial deltoid (∼28%), and upper trapezius (∼24%). With the exoskeleton, the operators perceived reduced global effort (∼17%) as well as a reduced local effort in the shoulder (∼18%), arm (∼22%), upper back (∼14%), and lower back (∼16%). The beneficial effect of the exoskeleton and its suitability in cleaning settings are corroborated by the acceptance and usability scores assigned by operators, which averaged ∼5.5 out of 7 points. To the authors' knowledge, this study is the first to present an experience of exoskeleton use in cleaning contexts. The outcomes of this research invite further studies to test occupational exoskeletons in various realistic applications to foster scientific-grounded ergonomic evaluations and encourage the informed adoption of the technology.

Exoskeletons for workers: A case series study in an enclosures production line.

This case-series study aims to investigate the effects of a passive shoulder support exoskeleton on experienced workers during their regular work shifts in an enclosures production site. Experimental activities included three sessions, two of which were conducted in-field (namely, at two workstations of the painting line, where panels were mounted and dismounted from the line; each session involved three participants), and one session was carried out in a realistic simulated environment (namely, the workstations were recreated in a laboratory; this session involved four participants). The effect of the exoskeleton was evaluated through electromyographic activity and perceived effort. After in-field sessions, device usability and user acceptance were also assessed. Data were reported individually for each participant. Results showed that the use of the exoskeleton reduced the total shoulder muscular activity compared to normal working conditions, in all subjects and experimental sessions. Similarly, the use of the exoskeleton resulted in reductions of the perceived effort in the shoulder, arm, and lower back. Overall, participants indicated high usability and acceptance of the device. This case series invites larger validation studies, also in diverse operational contexts.

Extension of the Rigid-Constraint Method for the Heuristic Suboptimal Parameter Tuning to Ten Sensor Fusion Algorithms Using Inertial and Magnetic Sensing.

The orientation of a magneto-inertial measurement unit can be estimated using a sensor fusion algorithm (SFA). However, orientation accuracy is greatly affected by the choice of the SFA parameter values which represents one of the most critical steps. A commonly adopted approach is to fine-tune parameter values to minimize the difference between estimated and true orientation. However, this can only be implemented within the laboratory setting by requiring the use of a concurrent gold-standard technology. To overcome this limitation, a Rigid-Constraint Method (RCM) was proposed to estimate suboptimal parameter values without relying on any orientation reference. The RCM method effectiveness was successfully tested on a single-parameter SFA, with an average error increase with respect to the optimal of 1.5 deg. In this work, the applicability of the RCM was evaluated on 10 popular SFAs with multiple parameters under different experimental scenarios. The average residual between the optimal and suboptimal errors amounted to 0.6 deg with a maximum of 3.7 deg. These encouraging results suggest the possibility to properly tune a generic SFA on different scenarios without using any reference. The synchronized dataset also including the optical data and the SFA codes are available online.

Analysis of the Accuracy of Ten Algorithms for Orientation Estimation Using Inertial and Magnetic Sensing under Optimal Conditions: One Size Does Not Fit All.

The orientation of a magneto and inertial measurement unit (MIMU) is estimated by means of sensor fusion algorithms (SFAs) thus enabling human motion tracking. However, despite several SFAs implementations proposed over the last decades, there is still a lack of consensus about the best performing SFAs and their accuracy. As suggested by recent literature, the filter parameters play a central role in determining the orientation errors. The aim of this work is to analyze the accuracy of ten SFAs while running under the best possible conditions (i.e., their parameter values are set using the orientation reference) in nine experimental scenarios including three rotation rates and three commercial products. The main finding is that parameter values must be specific for each SFA according to the experimental scenario to avoid errors comparable to those obtained when the default parameter values are used. Overall, when optimally tuned, no statistically significant differences are observed among the different SFAs in all tested experimental scenarios and the absolute errors are included between 3.8 deg and 7.1 deg. Increasing the rotation rate generally leads to a significant performance worsening. Errors are also influenced by the MIMU commercial model. SFA MATLAB implementations have been made available online.

Locomotory behaviour of the intertidal marble crab (Pachygrapsus marmoratus) supports the underwater spring-loaded inverted pendulum as a fundamental model for punting in animals.

In aquatic pedestrian locomotion the dynamics of terrestrial and aquatic environments are coupled. Here we study terrestrial running and aquatic punting locomotion of the marine-living crab Pachygrapsus marmoratus. We detected both active and passive phases of running and punting through the observation of crab locomotory behaviour in standardized settings and by three-dimensional kinematic analysis of its dynamic gaits using high-speed video cameras. Variations in different stride parameters were studied and compared. The comparison was done based on the dimensionless parameter the Froude number (Fr) to account for the effect of buoyancy and size variability among the crabs. The underwater spring-loaded inverted pendulum (USLIP) model better fitted the dynamics of aquatic punting. USLIP takes account of the damping effect of the aquatic environment, a variable not considered by the spring-loaded inverted pendulum (SLIP) model in reduced gravity. Our results highlight the underlying principles of aquatic terrestrial locomotion by comparing it with terrestrial locomotion. Comparing punting with running, we show and increased stride period, decreased duty cycle and orientation of the carapace more inclined with the horizontal plane, indicating the significance of fluid forces on the dynamics due to the aquatic environment. Moreover, we discovered periodicity in punting locomotion of crabs and two different gaits, namely, long-flight punting and short-flight punting, distinguished by both footfall patterns and kinematic parameters. The generic fundamental model which belongs to all animals performing both terrestrial and aquatic legged locomotion has implications for control strategies, evolution and translation to robotic artefacts.

Grasp force estimation from the transient EMG using high-density surface recordings.

OBJECTIVE: Understanding the neurophysiological signals underlying voluntary motor control and decoding them for prosthesis control are among the major challenges in applied neuroscience and bioengineering. Usually, information from the electrical activity of residual forearm muscles (i.e. the electromyogram, EMG) is used to control different functions of a prosthesis. Noteworthy, forearm EMG patterns at the onset of a contraction (transient phase) have shown to contain predictive information about upcoming grasps. However, decoding this information for the estimation of grasp force (GF) was so far overlooked. APPROACH: High density-EMG signals (192 channels) were recorded from twelve participants performing a pick-and-lift task. The final GF was estimated offline using linear regressors, with four subsets of channels and ten features obtained using three channels-features selection methods. Two different evaluation metrics (absolute error and R 2), complemented with statistical analysis, were used to select the optimal configuration of the parameters. Different windows of data starting at the GF onset were compared to determine the time at which the GF can be ascertained from the EMG signals. MAIN RESULTS: The prediction accuracy improved by increasing the window length from the moment of the onset and kept improving until the steady state at which a plateau of performances was reached. With our methodology, estimations of the GF through 16 EMG channels reached an absolute error of 2.52% the maximum voluntary force using only transient information and 1.99% with the first 500 ms of data following the onset. SIGNIFICANCE: The final GF estimation from transient EMG was comparable to the one obtained using steady state data, confirming our hypothesis that the transient phase contains information about the final GF. This result paves the way to fast online myoelectric controllers capable of decoding grasp strength from the very early portion of the EMG signal.

Tongue Rehabilitation Device for Dysphagic Patients.

Dysphagia refers to difficulty in swallowing often associated with syndromic disorders. In dysphagic patients' rehabilitation, tongue motility is usually treated and monitored via simple exercises, in which the tongue is pushed against a depressor held by the speech therapist in different directions. In this study, we developed and tested a simple pressure/force sensor device, named "Tonic Tongue (ToTo)", intended to support training and monitoring tasks for the rehabilitation of tongue musculature. It consists of a metallic frame holding a ball bearing support equipped with a sterile disposable depressor, whose angular displacements are counterbalanced by extensional springs. The conversion from angular displacement to force is managed using a simple mechanical model of ToTo operation. Since the force exerted by the tongue in various directions can be estimated, quantitative assessment of the outcome of a given training program is possible. A first prototype of ToTo was tested on 26 healthy adults, who were trained for one month. After the treatment, we observed a statistically significant improvement with a force up to 2.2 N (median value) in all tested directions of pushing, except in the downward direction, in which the improvement was slightly higher than 5 N (median value). ToTo promises to be an innovative and reliable device that can be used for the rehabilitation of dysphagic patients. Moreover, since it is a self-standing device, it could be used as a point-of-care solution for in-home rehabilitation management of dysphasia.

Ambulatory Assessment of the Dynamic Margin of Stability Using an Inertial Sensor Network.

Loss of stability is a precursor to falling and therefore represents a leading cause of injury, especially in fragile people. Thus, dynamic stability during activities of daily living (ADLs) needs to be considered to assess balance control and fall risk. The dynamic margin of stability (MOS) is often used as an indicator of how the body center of mass is located and moves relative to the base of support. In this work, we propose a magneto-inertial measurement unit (MIMU)-based method to assess the MOS of a gait. Six young healthy subjects were asked to walk on a treadmill at different velocities while wearing MIMUs on their lower limbs and pelvis. We then assessed the MOS by computing the lower body displacement with respect to the leading inverse kinematics approach. The results were compared with those obtained using a camera-based system in terms of root mean square deviation (RMSD) and correlation coefficient (ρ). We obtained a RMSD of ≤1.80 cm and ρ ≥ 0.85 for each walking velocity. The findings revealed that our method is comparable to camera-based systems in terms of accuracy, suggesting that it may represent a strategy to assess stability during ADLs in unstructured environments.

Combining Evolutionary and Adaptive Control Strategies for Quadruped Robotic Locomotion.

In traditional robotics, model-based controllers are usually needed in order to bring a robotic plant to the next desired state, but they present critical issues when the dimensionality of the control problem increases and disturbances from the external environment affect the system behavior, in particular during locomotion tasks. It is generally accepted that the motion control of quadruped animals is performed by neural circuits located in the spinal cord that act as a Central Pattern Generator and can generate appropriate locomotion patterns. This is thought to be the result of evolutionary processes that have optimized this network. On top of this, fine motor control is learned during the lifetime of the animal thanks to the plastic connections of the cerebellum that provide descending corrective inputs. This research aims at understanding and identifying the possible advantages of using learning during an evolution-inspired optimization for finding the best locomotion patterns in a robotic locomotion task. Accordingly, we propose a comparative study between two bio-inspired control architectures for quadruped legged robots where learning takes place either during the evolutionary search or only after that. The evolutionary process is carried out in a simulated environment, on a quadruped legged robot. To verify the possibility of overcoming the reality gap, the performance of both systems has been analyzed by changing the robot dynamics and its interaction with the external environment. Results show better performance metrics for the robotic agent whose locomotion method has been discovered by applying the adaptive module during the evolutionary exploration for the locomotion trajectories. Even when the motion dynamics and the interaction with the environment is altered, the locomotion patterns found on the learning robotic system are more stable, both in the joint and in the task space.

Reducing Fall Risk with Combined Motor and Cognitive Training in Elderly Fallers.

BACKGROUND: Falling is a major clinical problem in elderly people, demanding effective solutions. At present, the only effective intervention is motor training of balance and strength. Executive function-based training (EFt) might be effective at preventing falls according to evidence showing a relationship between executive functions and gait abnormalities. The aim was to assess the effectiveness of a motor and a cognitive treatment developed within the EU co-funded project I-DONT-FALL. METHODS: In a sample of 481 elderly people at risk of falls recruited in this multicenter randomised controlled trial, the effectiveness of a motor treatment (pure motor or mixed with EFt) of 24 one-hour sessions delivered through an i-Walker with a non-motor treatment (pure EFt or control condition) was evaluated. Similarly, a 24 one-hour session cognitive treatment (pure EFt or mixed with motor training), delivered through a touch-screen computer was compared with a non-cognitive treatment (pure motor or control condition). RESULTS: Motor treatment, particularly when mixed with EFt, reduced significantly fear of falling (F(1,478) = 6.786, p = 0.009) although to a limited extent (ES -0.25) restricted to the period after intervention. CONCLUSIONS: This study suggests the effectiveness of motor treatment empowered by EFt in reducing fear of falling.

Automatic classification of gait in children with early-onset ataxia or developmental coordination disorder and controls using inertial sensors.

Early-Onset Ataxia (EOA) and Developmental Coordination Disorder (DCD) are two conditions that affect coordination in children. Phenotypic identification of impaired coordination plays an important role in their diagnosis. Gait is one of the tests included in rating scales that can be used to assess motor coordination. A practical problem is that the resemblance between EOA and DCD symptoms can hamper their diagnosis. In this study we employed inertial sensors and a supervised classifier to obtain an automatic classification of the condition of participants. Data from shank and waist mounted inertial measurement units were used to extract features during gait in children diagnosed with EOA or DCD and age-matched controls. We defined a set of features from the recorded signals and we obtained the optimal features for classification using a backward sequential approach. We correctly classified 80.0%, 85.7%, and 70.0% of the control, DCD and EOA children, respectively. Overall, the automatic classifier correctly classified 78.4% of the participants, which is slightly better than the phenotypic assessment of gait by two pediatric neurologists (73.0%). These results demonstrate that automatic classification employing signals from inertial sensors obtained during gait maybe used as a support tool in the differential diagnosis of EOA and DCD. Furthermore, future extension of the classifier's test domains may help to further improve the diagnostic accuracy of pediatric coordination impairment. In this sense, this study may provide a first step towards incorporating a clinically objective and viable biomarker for identification of EOA and DCD.

Ambulatory Assessment of Instantaneous Velocity during Walking Using Inertial Sensor Measurements.

A novel approach for estimating the instantaneous velocity of the pelvis during walking was developed based on Inertial Measurement Units (IMUs). The instantaneous velocity was modeled by the sum of a cyclical component, decomposed in the Medio-Lateral (ML), VerTical (VT) and Antero-Posterior (AP) directions, and the Average Progression Velocity (APV) over each gait cycle. The proposed method required the availability of two IMUs, attached to the pelvis and one shank. Gait cycles were identified from the shank angular velocity; for each cycle, the Fourier series coefficients of the pelvis and shank acceleration signals were computed. The cyclical component was estimated by Fourier-based time-integration of the pelvis acceleration. A Bayesian Linear Regression (BLR) with Automatic Relevance Determination (ARD) predicted the APV from the stride time, the stance duration, and the Fourier series coefficients of the shank acceleration. Healthy subjects performed tasks of Treadmill Walking (TW) and Overground Walking (OW), and an optical motion capture system (OMCS) was used as reference for algorithm performance assessment. The widths of the limits of agreements (±1.96 standard deviation) were computed between the proposed method and the reference OMCS, yielding, for the cyclical component in the different directions: ML: ±0.07 m/s (±0.10 m/s); VT: ±0.03 m/s (±0.05 m/s); AP: ±0.06 m/s (±0.10 m/s), in TW (OW) conditions. The ARD-BLR achieved an APV root mean square error of 0.06 m/s (0.07 m/s) in the same conditions.

Assessing the Performance of Sensor Fusion Methods: Application to Magnetic-Inertial-Based Human Body Tracking.

Information from complementary and redundant sensors are often combined within sensor fusion algorithms to obtain a single accurate observation of the system at hand. However, measurements from each sensor are characterized by uncertainties. When multiple data are fused, it is often unclear how all these uncertainties interact and influence the overall performance of the sensor fusion algorithm. To address this issue, a benchmarking procedure is presented, where simulated and real data are combined in different scenarios in order to quantify how each sensor's uncertainties influence the accuracy of the final result. The proposed procedure was applied to the estimation of the pelvis orientation using a waist-worn magnetic-inertial measurement unit. Ground-truth data were obtained from a stereophotogrammetric system and used to obtain simulated data. Two Kalman-based sensor fusion algorithms were submitted to the proposed benchmarking procedure. For the considered application, gyroscope uncertainties proved to be the main error source in orientation estimation accuracy for both tested algorithms. Moreover, although different performances were obtained using simulated data, these differences became negligible when real data were considered. The outcome of this evaluation may be useful both to improve the design of new sensor fusion methods and to drive the algorithm tuning process.

Prior-to- and Post-Impact Fall Detection Using Inertial and Barometric Altimeter Measurements.

This paper investigates a fall detection system based on the integration of an inertial measurement unit with a barometric altimeter (BIMU). The vertical motion of the body part the BIMU was attached to was monitored on-line using a method that delivered drift-free estimates of the vertical velocity and estimates of the height change from the floor. The experimental study included activities of daily living of seven types and falls of five types, simulated by a cohort of 25 young healthy adults. The downward vertical velocity was thresholded at 1.38 m/s, yielding 80% sensitivity (SE), 100% specificity (SP) and a mean prior-to-impact time of 157 ms (range 40-300 ms). The soft falls, i.e., those with downward vertical velocity above 0.55 m/s and below 1.38 m/s were analyzed post-impact. Six fall detection methods, tuned to achieve 100% SE, were considered to include features of impact, change of posture and height, singularly or in association with one another. No single feature allowed for 100% SP. The detection accuracy marginally improved when the height change was considered in association with either the impact or the change of posture; the post-impact fall detection method that analyzed the impact and the change of posture together achieved 100% SP.

Dealing with Magnetic Disturbances in Human Motion Capture: A Survey of Techniques.

Magnetic-Inertial Measurement Units (MIMUs) based on microelectromechanical (MEMS) technologies are widespread in contexts such as human motion tracking. Although they present several advantages (lightweight, size, cost), their orientation estimation accuracy might be poor. Indoor magnetic disturbances represent one of the limiting factors for their accuracy, and, therefore, a variety of work was done to characterize and compensate them. In this paper, the main compensation strategies included within Kalman-based orientation estimators are surveyed and classified according to which degrees of freedom are affected by the magnetic data and to the magnetic disturbance rejection methods implemented. By selecting a representative method from each category, four algorithms were obtained and compared in two different magnetic environments: (1) small workspace with an active magnetic source; (2) large workspace without active magnetic sources. A wrist-worn MIMU was used to acquire data from a healthy subject, whereas a stereophotogrammetric system was adopted to obtain ground-truth data. The results suggested that the model-based approaches represent the best compromise between the two testbeds. This is particularly true when the magnetic data are prevented to affect the estimation of the angles with respect to the vertical direction.

A Simulation Environment for Benchmarking Sensor Fusion-Based Pose Estimators.

In-depth analysis and performance evaluation of sensor fusion-based estimators may be critical when performed using real-world sensor data. For this reason, simulation is widely recognized as one of the most powerful tools for algorithm benchmarking. In this paper, we present a simulation framework suitable for assessing the performance of sensor fusion-based pose estimators. The systems used for implementing the framework were magnetic/inertial measurement units (MIMUs) and a camera, although the addition of further sensing modalities is straightforward. Typical nuisance factors were also included for each sensor. The proposed simulation environment was validated using real-life sensor data employed for motion tracking. The higher mismatch between real and simulated sensors was about 5% of the measured quantity (for the camera simulation), whereas a lower correlation was found for an axis of the gyroscope (0.90). In addition, a real benchmarking example of an extended Kalman filter for pose estimation from MIMU and camera data is presented.

Fourier-based integration of quasi-periodic gait accelerations for drift-free displacement estimation using inertial sensors.

BACKGROUND: In biomechanical studies Optical Motion Capture Systems (OMCS) are considered the gold standard for determining the orientation and the position (pose) of an object in a global reference frame. However, the use of OMCS can be difficult, which has prompted research on alternative sensing technologies, such as body-worn inertial sensors. METHODS: We developed a drift-free method to estimate the three-dimensional (3D) displacement of a body part during cyclical motions using body-worn inertial sensors. We performed the Fourier analysis of the stride-by-stride estimates of the linear acceleration, which were obtained by transposing the specific forces measured by the tri-axial accelerometer into the global frame using a quaternion-based orientation estimation algorithm and detecting when each stride began using a gait-segmentation algorithm. The time integration was performed analytically using the Fourier series coefficients; the inverse Fourier series was then taken for reconstructing the displacement over each single stride. The displacement traces were concatenated and spline-interpolated to obtain the entire trace. RESULTS: The method was applied to estimate the motion of the lower trunk of healthy subjects that walked on a treadmill and it was validated using OMCS reference 3D displacement data; different approaches were tested for transposing the measured specific force into the global frame, segmenting the gait and performing time integration (numerically and analytically). The width of the limits of agreements were computed between each tested method and the OMCS reference method for each anatomical direction: Medio-Lateral (ML), VerTical (VT) and Antero-Posterior (AP); using the proposed method, it was observed that the vertical component of displacement (VT) was within ±4 mm (±1.96 standard deviation) of OMCS data and each component of horizontal displacement (ML and AP) was within ±9 mm of OMCS data. CONCLUSIONS: Fourier harmonic analysis was applied to model stride-by-stride linear accelerations during walking and to perform their analytical integration. Our results showed that analytical integration based on Fourier series coefficients was a useful approach to accurately estimate 3D displacement from noisy acceleration data.

How Angular Velocity Features and Different Gyroscope Noise Types Interact and Determine Orientation Estimation Accuracy.

In human movement analysis, 3D body segment orientation can be obtained through the numerical integration of gyroscope signals. These signals, however, are affected by errors that, for the case of micro-electro-mechanical systems, are mainly due to: constant bias, scale factor, white noise, and bias instability. The aim of this study is to assess how the orientation estimation accuracy is affected by each of these disturbances, and whether it is influenced by the angular velocity magnitude and 3D distribution across the gyroscope axes. Reference angular velocity signals, either constant or representative of human walking, were corrupted with each of the four noise types within a simulation framework. The magnitude of the angular velocity affected the error in the orientation estimation due to each noise type, except for the white noise. Additionally, the error caused by the constant bias was also influenced by the angular velocity 3D distribution. As the orientation error depends not only on the noise itself but also on the signal it is applied to, different sensor placements could enhance or mitigate the error due to each disturbance, and special attention must be paid in providing and interpreting measures of accuracy for orientation estimation algorithms.

Estimating orientation using magnetic and inertial sensors and different sensor fusion approaches: accuracy assessment in manual and locomotion tasks.

Magnetic and inertial measurement units are an emerging technology to obtain 3D orientation of body segments in human movement analysis. In this respect, sensor fusion is used to limit the drift errors resulting from the gyroscope data integration by exploiting accelerometer and magnetic aiding sensors. The present study aims at investigating the effectiveness of sensor fusion methods under different experimental conditions. Manual and locomotion tasks, differing in time duration, measurement volume, presence/absence of static phases, and out-of-plane movements, were performed by six subjects, and recorded by one unit located on the forearm or the lower trunk, respectively. Two sensor fusion methods, representative of the stochastic (Extended Kalman Filter) and complementary (Non-linear observer) filtering, were selected, and their accuracy was assessed in terms of attitude (pitch and roll angles) and heading (yaw angle) errors using stereophotogrammetric data as a reference. The sensor fusion approaches provided significantly more accurate results than gyroscope data integration. Accuracy improved mostly for heading and when the movement exhibited stationary phases, evenly distributed 3D rotations, it occurred in a small volume, and its duration was greater than approximately 20 s. These results were independent from the specific sensor fusion method used. Practice guidelines for improving the outcome accuracy are provided.

Walking speed estimation using foot-mounted inertial sensors: comparing machine learning and strap-down integration methods.

In this paper we implemented machine learning (ML) and strap-down integration (SDI) methods and analyzed them for their capability of estimating stride-by-stride walking speed. Walking speed was computed by dividing estimated stride length by stride time using data from a foot mounted inertial measurement unit. In SDI methods stride-by-stride walking speed estimation was driven by detecting gait events using a hidden Markov model (HMM) based method (HMM-based SDI); alternatively, a threshold-based gait event detector was investigated (threshold-based SDI). In the ML method a linear regression model was developed for stride length estimation. Whereas the gait event detectors were a priori fixed without training, the regression model was validated with leave-one-subject-out cross-validation. A subject-specific regression model calibration was also implemented to personalize the ML method. Healthy adults performed over-ground walking trials at natural, slower-than-natural and faster-than-natural speeds. The ML method achieved a root mean square estimation error of 2.0% and 4.2%, with and without personalization, against 2.0% and 3.1% by HMM-based SDI and threshold-based SDI. In spite that the results achieved by the two approaches were similar, the ML method, as compared with SDI methods, presented lower intra-subject variability and higher inter-subject variability, which was reduced by personalization.