A Dynamic Fitting Strategy for Physiological Models: A Case Study of a Cardiorespiratory Model for the Simulation of Incremental Aerobic Exercise.

Using mathematical models of physiological systems in medicine has allowed for the development of diagnostic, treatment, and medical educational tools. However, their complexity restricts, in most cases, their application for predictive, preventive, and personalized purposes. Although there are strategies that reduce the complexity of applying models based on fitting techniques, most of them are focused on a single instant of time, neglecting the effect of the system's temporal evolution. The objective of this research was to introduce a dynamic fitting strategy for physiological models with an extensive array of parameters and a constrained amount of experimental data. The proposed strategy focused on obtaining better predictions based on the temporal trends in the system's parameters and being capable of predicting future states. The study utilized a cardiorespiratory model as a case study. Experimental data from a longitudinal study of healthy adult subjects undergoing aerobic exercise were used for fitting and validation. The model predictions obtained in a steady state using the proposed strategy and the traditional single-fit approach were compared. The most successful outcomes were primarily linked to the proposed strategy, exhibiting better overall results regarding accuracy and behavior than the traditional population fitting approach at a single instant in time. The results evidenced the usefulness of the dynamic fitting strategy, highlighting its use for predictive, preventive, and personalized applications.

Electroencephalographic assessment in patients with Rett syndrome during cognitive stimulation by means of eye tracking technology and alternative and augmentative communication systems.

Rett syndrome (RTT) is considered a rare disease despite being the leading genetic disorder to cause severe intellectual disability in women. There is no cure for RTT, so the treatment is symptomatic and supporting, requiring a multidisciplinary approach. Occupational therapy can help girls and their families to improve communication, being one of the main concerns when verbal language and intentional hand movement are impaired or lost. This paper presents a pilot study of cognitive training through the combined use of eye-tracking technology (ETT) and augmentative and alternative communication (AAC) using the Peabody Picture Vocabulary Test (PPVT-IV). The objective was to evaluate brain activation by means of electroencephalography (EEG) during the stimulation of non-verbal communication. EEG data were recorded during an eyes-open resting state (EO-RS) period and during cognitive stimulation via AAC activity. To assess their effect, both signals were compared at the spectral level, focusing on frequency, brain symmetry and connectivity. During the task, a redistribution of power towards fast frequency bands was observed, as well as an improvement in the brain symmetry index (BSI) and functional synchronicity through increased coherence. Therefore, the results of the spectral analysis showed a possible deviation from the pathological pattern, manifesting a positive effect in the use of non-verbal cognitive stimulation activities. In conclusion, it was observed that it is possible to establish a cognitive training system that produces brain activation and favors communication and learning despite intentional language loss.Clinical Relevance- This manifests a method of cognitive training that would induce brain activation in RTT patients with absence of intentional communication. The evaluation system through spectral analysis could complement the standardized protocols to asses communication that are based on verbal and motor production.

Spinal Cord Injury Patients Exhibit Changes in Motor-Related Activity and Topographic Distribution.

Spinal Cord Injury (SCI) is a common disease that usually limits the patient's independence by affecting their motor function. SCI patients usually present neuroplasticity, which allows brain signals transmission through spread pathways. Some innovative rehabilitation therapies, such as functional electrical stimulation (FES) or Brain-computer interfaces (BCIs) jointly with motor neuroprostheses, provide hope for functional restoration. BCIs require the analysis of event-related EEG potentials (ERPs). Movement-related cortical potentials (MRCPs) and event-related desynchroni-zation and synchronization (ERD/ERS) are the most commonly studied ERPs during motor activity. ERPs of healthy subjects may vary from SCI patients. Thus, this study aimed to compare ERPs between healthy subjects and SCI patients during upper-limb movements (forearm supination and pronation, and hand open). Differences between controls and SCI patients were shown in terms of ERPs' amplitude as well as in topographic maps. Changes in amplitude were more substantial in ERD potentials than in MRCPs, while topographic maps showed better localization of all features in healthy patients. The level of SCI injury determines the patients' mobility. A comparison between complete, partial and no motor function subjects showed lower values of feature's amplitudes in the latter group.Clinical Relevance- This demonstrates the existence of significant statistical differences between healthy and SCI subjects, and might be helpful when performing SCI rehabilitation techniques such as designing BCI and neuroprostheses, or analyzing and understanding the brain plasticity process.

A Novel Strategy to Fit and Validate Physiological Models: A Case Study of a Cardiorespiratory Model for Simulation of Incremental Aerobic Exercise.

Applying complex mathematical models of physiological systems is challenging due to the large number of parameters. Identifying these parameters through experimentation is difficult, and although procedures for fitting and validating models are reported, no integrated strategy exists. Additionally, the complexity of optimization is generally neglected when the number of experimental observations is restricted, obtaining multiple solutions or results without physiological justification. This work proposes a fitting and validation strategy for physiological models with many parameters under various populations, stimuli, and experimental conditions. A cardiorespiratory system model is used as a case study, and the strategy, model, computational implementation, and data analysis are described. Using optimized parameter values, model simulations are compared to those obtained using nominal values, with experimental data as a reference. Overall, a reduction in prediction error is achieved compared to that reported for model building. Furthermore, the behavior and accuracy of all the predictions in the steady state were improved. The results validate the fitted model and provide evidence of the proposed strategy's usefulness.

An Accurate and Real-Time Method for Resolving Superimposed Action Potentials in MultiUnit Recordings.

OBJECTIVE: Spike sorting of muscular and neural recordings requires separating action potentials that overlap in time (superimposed action potentials (APs)). We propose a new algorithm for resolving superimposed action potentials, and we test it on intramuscular EMG (iEMG) and intracortical recordings. METHODS: Discrete-time shifts of the involved APs are first selected based on a heuristic extension of the peel-off algorithm. Then, the time shifts that provide the minimal residual Euclidean norm are identified (Discrete Brute force Correlation (DBC)). The optimal continuous-time shifts are then estimated (High-Resolution BC (HRBC)). In Fusion HRBC (FHRBC), two other cost functions are used. A parallel implementation of the DBC and HRBC algorithms was developed. The performance of the algorithms was assessed on 11,000 simulated iEMG and 14,000 neural recording superpositions, including two to eight APs, and eight experimental iEMG signals containing four to eleven active motor units. The performance of the proposed algorithms was compared with that of the Branch-and-Bound (BB) algorithm using the Rank-Product (RP) method in terms of accuracy and efficiency. RESULTS: The average accuracy of the DBC, HRBC and FHRBC methods on the entire simulated datasets was 92.16±17.70, 93.65±16.89, and 94.90±15.15 (%). The DBC algorithm outperformed the other algorithms based on the RP method. The average accuracy and running time of the DBC algorithm on 10.5 ms superimposed spikes of the experimental signals were 92.1±21.7 (%) and 2.3±15.3 (ms). CONCLUSION AND SIGNIFICANCE: The proposed algorithm is promising for real-time neural decoding, a central problem in neural and muscular decoding and interfacing.

Influence of the number of trials on evoked motor cortical activity in EEG recordings.

Objective. Improvements in electroencephalography enable the study of the localization of active brain regions during motor tasks. Movement-related cortical potentials (MRCPs), and event-related desynchronization (ERD) and synchronization are the main motor-related cortical phenomena/neural correlates observed when a movement is elicited. When assessing neurological diseases, averaging techniques are commonly applied to characterize motor related processes better. In this case, a large number of trials is required to obtain a motor potential that is representative enough of the subject's condition. This study aimed to assess the effect of a limited number of trials on motor-related activity corresponding to different upper limb movements (elbow flexion/extension, pronation/supination and hand open/close).Approach. An open dataset consisting on 15 healthy subjects was used for the analysis. A Monte Carlo simulation approach was applied to analyse, in a robust way, different typical time- and frequency-domain features, topography, and low-resolution electromagnetic tomography.Main results. Grand average potentials, and topographic and tomographic maps showed few differences when using fewer trials, but shifts in the localization of motor-related activity were found for several individuals. MRCP and beta ERD features were more robust to a limited number of trials, yielding differences lower than 20% for cases with 50 trials or more. Strong correlations between features were obtained for subsets above 50 trials. However, the inter-subject variability increased as the number of trials decreased. The elbow flexion/extension movement showed a more robust performance for a limited number of trials, both in population and in individual-based analysis.Significance. Our findings suggested that 50 trials can be an appropriate number to obtain stable motor-related features in terms of differences in the averaged motor features, correlation, and changes in topography and tomography.

Improving the ripple classification in focal pediatric epilepsy: identifying pathological high-frequency oscillations by Gaussian mixture model clustering.

Objective. High-frequency oscillations (HFOs) have emerged as a promising clinical biomarker for presurgical evaluation in childhood epilepsy. HFOs are commonly classified in stereo-encephalography as ripples (80-200 Hz) and fast ripples (200-500 Hz). Ripples are less specific and not so directly associated with epileptogenic activity because of their physiological and pathological origin. The aim of this paper is to distinguish HFOs in the ripple band and to improve the evaluation of the epileptogenic zone (EZ).Approach. This study constitutes a novel modeling approach evaluated in ten patients from Sant Joan de Deu Pediatric Hospital (Barcelona, Spain), with clearly-defined seizure onset zones (SOZ) during presurgical evaluation. A subject-by-subject basis analysis is proposed: a probabilistic Gaussian mixture model (GMM) based on the combination of specific ripple features is applied for estimating physiological and pathological ripple subpopulations.Main Results. Clear pathological and physiological ripples are identified. Features differ considerably among patients showing within-subject variability, suggesting that individual models are more appropriate than a traditional whole-population approach. The difference in rates inside and outside the SOZ for pathological ripples is significantly higher than when considering all the ripples. These significant differences also appear in signal segments without epileptiform activity. Pathological ripple rates show a sharp decline from SOZ to non-SOZ contacts and a gradual decrease with distance.Significance. This novel individual GMM approach improves ripple classification and helps to refine the delineation of the EZ, as well as being appropriate to investigate the interaction of epileptogenic and propagation networks.

Choosing Strategies to Deal with Artifactual EEG Data in Children with Cognitive Impairment.

Rett syndrome is a disease that involves acute cognitive impairment and, consequently, a complex and varied symptomatology. This study evaluates the EEG signals of twenty-nine patients and classify them according to the level of movement artifact. The main goal is to achieve an artifact rejection strategy that performs well in all signals, regardless of the artifact level. Two different methods have been studied: one based on the data distribution and the other based on the energy function, with entropy as its main component. The method based on the data distribution shows poor performance with signals containing high amplitude outliers. On the contrary, the method based on the energy function is more robust to outliers. As it does not depend on the data distribution, it is not affected by artifactual events. A double rejection strategy has been chosen, first on a motion signal (accelerometer or EEG low-pass filtered between 1 and 10 Hz) and then on the EEG signal. The results showed a higher performance when working combining both artifact rejection methods. The energy-based method, to isolate motion artifacts, and the data-distribution-based method, to eliminate the remaining lower amplitude artifacts were used. In conclusion, a new method that proves to be robust for all types of signals is designed.

SGM: a novel time-frequency algorithm based on unsupervised learning improves high-frequency oscillation detection in epilepsy.

OBJECTIVE: We propose a novel automated method called the S-Transform Gaussian Mixture detection algorithm (SGM) to detect high-frequency oscillations (HFO) combining the strengths of different families of previously published detectors. APPROACH: This algorithm does not depend on parameter tuning on a subject (or database) basis, uses time-frequency characteristics, and relies on non-supervised classification to determine if the events standing out from the baseline activity are HFO or not. SGM consists of three steps: the first stage computes the signal baseline using the entropy of the autocorrelation; the second uses the S-Transform to obtain several time-frequency features (area, entropy, and time and frequency widths); and in the third stage Gaussian mixture models cluster time-frequency features to decide if events correspond to HFO-like activity. To validate the SGM algorithm we tested its performance in simulated and real environments. MAIN RESULTS: We assessed the algorithm on a publicly available simulated stereoelectroencephalographic (SEEG) database with varying signal-to-noise ratios (SNR), obtaining very good results for medium and high SNR signals. We further tested the SGM algorithm on real signals from patients with focal epilepsy, in which HFO detection was performed visually by experts, yielding a high agreement between experts and SGM. SIGNIFICANCE: The SGM algorithm displayed proper performance in simulated and real environments and therefore can be used for non-supervised detection of HFO. This non-supervised algorithm does not require previous labelling by experts or parameter adjustment depending on the subject or database considered. SGM is not a computationally intensive algorithm, making it suitable to detect and characterize HFO in long-term SEEG recordings.

Alterations in EEG connectivity in healthy young adults provide an indicator of sleep depth.

Current sleep analyses have used electroencephalography (EEG) to establish sleep intensity through linear and nonlinear measures. Slow wave activity (SWA) and entropy are the most commonly used markers of sleep depth. The purpose of this study is to evaluate changes in brain EEG connectivity during sleep in healthy subjects and compare them with SWA and entropy. Four different connectivity metrics: coherence (MSC), synchronization likelihood (SL), cross mutual information function (CMIF), and phase locking value (PLV), were computed focusing on their correlation with sleep depth. These measures provide different information and perspectives about functional connectivity. All connectivity measures revealed to have functional changes between the different sleep stages. The averaged CMIF seemed to be a more robust connectivity metric to measure sleep depth (correlations of 0.78 and 0.84 with SWA and entropy, respectively), translating greater linear and nonlinear interdependences between brain regions especially during slow wave sleep. Potential changes of brain connectivity were also assessed throughout the night. Connectivity measures indicated a reduction of functional connectivity in N2 as sleep progresses. The validation of connectivity indexes is necessary because they can reveal the interaction between different brain regions in physiological and pathological conditions and help understand the different functions of deep sleep in humans.

An Improved Dynamic Model for the Respiratory Response to Exercise.

Respiratory system modeling has been extensively studied in steady-state conditions to simulate sleep disorders, to predict its behavior under ventilatory diseases or stimuli and to simulate its interaction with mechanical ventilation. Nevertheless, the studies focused on the instantaneous response are limited, which restricts its application in clinical practice. The aim of this study is double: firstly, to analyze both dynamic and static responses of two known respiratory models under exercise stimuli by using an incremental exercise stimulus sequence (to analyze the model responses when step inputs are applied) and experimental data (to assess prediction capability of each model). Secondly, to propose changes in the models' structures to improve their transient and stationary responses. The versatility of the resulting model vs. the other two is shown according to the ability to simulate ventilatory stimuli, like exercise, with a proper regulation of the arterial blood gases, suitable constant times and a better adjustment to experimental data. The proposed model adjusts the breathing pattern every respiratory cycle using an optimization criterion based on minimization of work of breathing through regulation of respiratory frequency.

Automated detection of epileptic ripples in MEG using beamformer-based virtual sensors.

OBJECTIVE: In epilepsy, high-frequency oscillations (HFOs) are expressively linked to the seizure onset zone (SOZ). The detection of HFOs in the noninvasive signals from scalp electroencephalography (EEG) and magnetoencephalography (MEG) is still a challenging task. The aim of this study was to automate the detection of ripples in MEG signals by reducing the high-frequency noise using beamformer-based virtual sensors (VSs) and applying an automatic procedure for exploring the time-frequency content of the detected events. APPROACH: Two-hundred seconds of MEG signal and simultaneous iEEG were selected from nine patients with refractory epilepsy. A two-stage algorithm was implemented. Firstly, beamforming was applied to the whole head to delimitate the region of interest (ROI) within a coarse grid of MEG-VS. Secondly, a beamformer using a finer grid in the ROI was computed. The automatic detection of ripples was performed using the time-frequency response provided by the Stockwell transform. Performance was evaluated through comparisons with simultaneous iEEG signals. MAIN RESULTS: ROIs were located within the seizure-generating lobes in the nine subjects. Precision and sensitivity values were 79.18% and 68.88%, respectively, by considering iEEG-detected events as benchmarks. A higher number of ripples were detected inside the ROI compared to the same region in the contralateral lobe. SIGNIFICANCE: The evaluation of interictal ripples using non-invasive techniques can help in the delimitation of the epileptogenic zone and guide placement of intracranial electrodes. This is the first study that automatically detects ripples in MEG in the time domain located within the clinically expected epileptic area taking into account the time-frequency characteristics of the events through the whole signal spectrum. The algorithm was tested against intracranial recordings, the current gold standard. Further studies should explore this approach to enable the localization of noninvasively recorded HFOs to help during pre-surgical planning and to reduce the need for invasive diagnostics.

Acute Sleep Deprivation Induces a Local Brain Transfer Information Increase in the Frontal Cortex in a Widespread Decrease Context.

Sleep deprivation (SD) has adverse effects on mental and physical health, affecting the cognitive abilities and emotional states. Specifically, cognitive functions and alertness are known to decrease after SD. The aim of this work was to identify the directional information transfer after SD on scalp EEG signals using transfer entropy (TE). Using a robust methodology based on EEG recordings of 18 volunteers deprived from sleep for 36 h, TE and spectral analysis were performed to characterize EEG data acquired every 2 h. Correlation between connectivity measures and subjective somnolence was assessed. In general, TE showed medium- and long-range significant decreases originated at the occipital areas and directed towards different regions, which could be interpreted as the transfer of predictive information from parieto-occipital activity to the rest of the head. Simultaneously, short-range increases were obtained for the frontal areas, following a consistent and robust time course with significant maps after 20 h of sleep deprivation. Changes during sleep deprivation in brain network were measured effectively by TE, which showed increased local connectivity and diminished global integration. TE is an objective measure that could be used as a potential measure of sleep pressure and somnolence with the additional property of directed relationships.

Influence of metallic artifact filtering on MEG signals for source localization during interictal epileptiform activity.

OBJECTIVE: Medical intractable epilepsy is a common condition that affects 40% of epileptic patients that generally have to undergo resective surgery. Magnetoencephalography (MEG) has been increasingly used to identify the epileptogenic foci through equivalent current dipole (ECD) modeling, one of the most accepted methods to obtain an accurate localization of interictal epileptiform discharges (IEDs). Modeling requires that MEG signals are adequately preprocessed to reduce interferences, a task that has been greatly improved by the use of blind source separation (BSS) methods. MEG recordings are highly sensitive to metallic interferences originated inside the head by implanted intracranial electrodes, dental prosthesis, etc and also coming from external sources such as pacemakers or vagal stimulators. To reduce these artifacts, a BSS-based fully automatic procedure was recently developed and validated, showing an effective reduction of metallic artifacts in simulated and real signals (Migliorelli et al 2015 J. Neural Eng. 12 046001). The main objective of this study was to evaluate its effects in the detection of IEDs and ECD modeling of patients with focal epilepsy and metallic interference. APPROACH: A comparison between the resulting positions of ECDs was performed: without removing metallic interference; rejecting only channels with large metallic artifacts; and after BSS-based reduction. Measures of dispersion and distance of ECDs were defined to analyze the results. MAIN RESULTS: The relationship between the artifact-to-signal ratio and ECD fitting showed that higher values of metallic interference produced highly scattered dipoles. Results revealed a significant reduction on dispersion using the BSS-based reduction procedure, yielding feasible locations of ECDs in contrast to the other two approaches. SIGNIFICANCE: The automatic BSS-based method can be applied to MEG datasets affected by metallic artifacts as a processing step to improve the localization of epileptic foci.

Auditory P3a and P3b neural generators in schizophrenia: An adaptive sLORETA P300 localization approach.

The present study investigates the neural substrates underlying cognitive processing in schizophrenia (Sz) patients. To this end, an auditory 3-stimulus oddball paradigm was used to identify P3a and P3b components, elicited by rare-distractor and rare-target tones, respectively. Event-related potentials (ERP) were recorded from 31 Sz patients and 38 healthy controls. The P3a and P3b brain-source generators were identified by time-averaging of low-resolution brain electromagnetic tomography (LORETA) current density images. In contrast with the commonly used fixed window of interest (WOI), we proposed to apply an adaptive WOI, which takes into account subjects' P300 latency variability. Our results showed different P3a and P3b source activation patterns in both groups. P3b sources included frontal, parietal and limbic lobes, whereas P3a response generators were localized over bilateral frontal and superior temporal regions. These areas have been related to the discrimination of auditory stimulus and to the inhibition (P3a) or the initiation (P3b) of motor response in a cognitive task. In addition, differences in source localization between Sz and control groups were observed. Sz patients showed lower P3b source activity in bilateral frontal structures and the cingulate. P3a generators were less widespread for Sz patients than for controls in right superior, medial and middle frontal gyrus. Our findings suggest that target and distractor processing involves distinct attentional subsystems, both being altered in Sz. Hence, the study of neuroelectric brain information can provide further insights to understand cognitive processes and underlying mechanisms in Sz.

Automatic BSS-based filtering of metallic interference in MEG recordings: definition and validation using simulated signals.

OBJECTIVE: One of the principal drawbacks of magnetoencephalography (MEG) is its high sensitivity to metallic artifacts, which come from implanted intracranial electrodes and dental ferromagnetic prosthesis and produce a high distortion that masks cerebral activity. The aim of this study was to develop an automatic algorithm based on blind source separation (BSS) techniques to remove metallic artifacts from MEG signals. APPROACH: Three methods were evaluated: AMUSE, a second-order technique; and INFOMAX and FastICA, both based on high-order statistics. Simulated signals consisting of real artifact-free data mixed with real metallic artifacts were generated to objectively evaluate the effectiveness of BSS and the subsequent interference reduction. A completely automatic detection of metallic-related components was proposed, exploiting the known characteristics of the metallic interference: regularity and low frequency content. MAIN RESULTS: The automatic procedure was applied to the simulated datasets and the three methods exhibited different performances. Results indicated that AMUSE preserved and consequently recovered more brain activity than INFOMAX and FastICA. Normalized mean squared error for AMUSE decomposition remained below 2%, allowing an effective removal of artifactual components. SIGNIFICANCE: To date, the performance of automatic artifact reduction has not been evaluated in MEG recordings. The proposed methodology is based on an automatic algorithm that provides an effective interference removal. This approach can be applied to any MEG dataset affected by metallic artifacts as a processing step, allowing further analysis of unusable or poor quality data.

Neurophysiological correlates in Mild Cognitive Impairment detected using group Independent Component Analysis.

Alzheimer's disease is the most prevalent cause of dementia. Mild Cognitive Impairment (MCI) is defined as a grey area between intact cognitive functioning and clinical dementia. Electroencephalography (EEG) has been used to identify biomarkers in dementia. Currently, there is a great interest in translating the study from raw signals to signal generators, trying to keep the relationship with neurophysiology. In the current study, EEG recordings during an encoding task were acquired in MCI subjects and healthy controls. Data was decomposed using group Independent Component Analysis (gICA) and the most neuronal components were analyzed using Phase Intertrial Coherence (PIC) and Phase shift Intertrial Coherence (PsIC). MCI subjects exhibited an increase of PIC in the theta band, while controls showed increase in PsIC in the alpha band. Correlation between PIC and PsIC and clinical scales were also found. Those findings indicate that the methodology proposed based in gICA can help to extract information from EEG recordings with neurophysiological meaning.

High-density surface EMG maps from upper-arm and forearm muscles.

BACKGROUND: sEMG signal has been widely used in different applications in kinesiology and rehabilitation as well as in the control of human-machine interfaces. In general, the signals are recorded with bipolar electrodes located in different muscles. However, such configuration may disregard some aspects of the spatial distribution of the potentials like location of innervation zones and the manifestation of inhomogineties in the control of the muscular fibers. On the other hand, the spatial distribution of motor unit action potentials has recently been assessed with activation maps obtained from High Density EMG signals (HD-EMG), these lasts recorded with arrays of closely spaced electrodes. The main objective of this work is to analyze patterns in the activation maps, associating them with four movement directions at the elbow joint and with different strengths of those tasks. Although the activation pattern can be assessed with bipolar electrodes, HD-EMG maps could enable the extraction of features that depend on the spatial distribution of the potentials and on the load-sharing between muscles, in order to have a better differentiation between tasks and effort levels. METHODS: An experimental protocol consisting of isometric contractions at three levels of effort during flexion, extension, supination and pronation at the elbow joint was designed and HD-EMG signals were recorded with 2D electrode arrays on different upper-limb muscles. Techniques for the identification and interpolation of artifacts are explained, as well as a method for the segmentation of the activation areas. In addition, variables related to the intensity and spatial distribution of the maps were obtained, as well as variables associated to signal power of traditional single bipolar recordings. Finally, statistical tests were applied in order to assess differences between information extracted from single bipolar signals or from HD-EMG maps and to analyze differences due to type of task and effort level. RESULTS: Significant differences were observed between EMG signal power obtained from single bipolar configuration and HD-EMG and better results regarding the identification of tasks and effort levels were obtained with the latter. Additionally, average maps for a population of 12 subjects were obtained and differences in the co-activation pattern of muscles were found not only from variables related to the intensity of the maps but also to their spatial distribution. CONCLUSIONS: Intensity and spatial distribution of HD-EMG maps could be useful in applications where the identification of movement intention and its strength is needed, for example in robotic-aided therapies or for devices like powered- prostheses or orthoses. Finally, additional data transformations or other features are necessary in order to improve the performance of tasks identification.

A synchronization system for the analysis of biomedical signals recorded with different devices from mechanically ventilated patients.

Conducting research associated with mechanically ventilated patients often requires the recording of several biomedical signals to dispose of multiple sources of information to perform a robust analysis. This is especially important in the analysis of the relationship between pressure, volume and flow, signals available from mechanical ventilators, and other biopotentials such as the electromyogram of respiratory muscles, intrinsically related with the ventilatory process, but not commonly recorded in the clinical practice. Despite the usefulness of recording signals from multiple sources, few medical devices include the possibility of synchronizing its data with other provided by different biomedical equipment and some may use inaccurate sampling frequencies. Even thought a variant or inaccurate sampling rate does not affect the monitoring of critical patients, it restricts the study of simultaneous related events useful in research of respiratory system activity. In this article a device for temporal synchronization of signals recorded from multiple biomedical devices is described as well as its application in the study of patients undergoing mechanical ventilation with research purposes.