Cholinergic-related pupil activity reflects level of emotionality during motor performance

Pupil size covaries with the diffusion rate of the cholinergic and noradrenergic neurons throughout the brain, which are essential to arousal. Recent findings suggest that slow pupil fluctuations during locomotion are an index of sustained activity in cholinergic axons, whereas phasic dilations are related to the activity of noradrenergic axons. Here, we investigated movement induced arousal (i.e., by singing and swaying to music), hypothesising that actively engaging in musical behaviour will provoke stronger emotional engagement in participants and lead to different qualitative patterns of tonic and phasic pupil activity. A challenge in the analysis of pupil data is the turbulent behaviour of pupil diameter due to exogenous ocular activity commonly encountered during motor tasks and the high variability typically found between individuals. To address this, we developed an algorithm that adaptively estimates and removes pupil responses to ocular events, as well as a functional data methodology, derived from Pfaffs' generalised arousal, that provides a new statistical dimension on how pupil data can be interpreted according to putative neuromodulatory signalling. We found that actively engaging in singing enhanced slow cholinergic-related pupil dilations and having the opportunity to move your body while performing amplified the effect of singing on pupil activity. Phasic pupil oscillations during motor execution attenuated in time, which is often interpreted as a measure of sense of agency over movement.

New approaches and technical considerations in detecting outlier measurements and trajectories in longitudinal children growth data

BACKGROUND

Growth studies rely on longitudinal measurements, typically represented as trajectories. However, anthropometry is prone to errors that can generate outliers. While various methods are available for detecting outlier measurements, a gold standard has yet to be identified, and there is no established method for outlying trajectories. Thus, outlier types and their effects on growth pattern detection still need to be investigated. This work aimed to assess the performance of six methods at detecting different types of outliers, propose two novel methods for outlier trajectory detection and evaluate how outliers affect growth pattern detection.

METHODS

We included 393 healthy infants from The Applied Research Group for Kids (TARGet Kids!) cohort and 1651 children with severe malnutrition from the co-trimoxazole prophylaxis clinical trial. We injected outliers of three types and six intensities and applied four outlier detection methods for measurements (model-based and World Health Organization cut-offs-based) and two for trajectories. We also assessed growth pattern detection before and after outlier injection using time series clustering and latent class mixed models. Error type, intensity, and population affected method performance.

RESULTS

Model-based outlier detection methods performed best for measurements with precision between 5.72-99.89%, especially for low and moderate error intensities. The clustering-based outlier trajectory method had high precision of 14.93-99.12%. Combining methods improved the detection rate to 21.82% in outlier measurements. Finally, when comparing growth groups with and without outliers, the outliers were shown to alter group membership by 57.9 -79.04%.

CONCLUSIONS

World Health Organization cut-off-based techniques were shown to perform well in few very particular cases (extreme errors of high intensity), while model-based techniques performed well, especially for moderate errors of low intensity. Clustering-based outlier trajectory detection performed exceptionally well across all types and intensities of errors, indicating a potential strategic change in how outliers in growth data are viewed. Finally, the importance of detecting outliers was shown, given its impact on children growth studies, as demonstrated by comparing results of growth group detection.

Central Posterior Envelopes for Bayesian Functional Principal Component Analysis

Bayesian methods provide direct inference in functional data analysis applications without reliance on bootstrap techniques. A major tool in functional data applications is the functional principal component analysis which decomposes the data around a common mean function and identifies leading directions of variation. Bayesian functional principal components analysis (BFPCA) provides uncertainty quantification on the estimated functional model components via the posterior samples obtained. We propose central posterior envelopes (CPEs) for BFPCA based on functional depth as a descriptive visualization tool to summarize variation in the posterior samples of the estimated functional model components, contributing to uncertainty quantification in BFPCA. The proposed BFPCA relies on a latent factor model and targets model parameters within a mixed effects modeling framework using modified multiplicative gamma process shrinkage priors on the variance components. Functional depth provides a center-outward order to a sample of functions. We utilize modified band depth and modified volume depth for ordering of a sample of functions and surfaces, respectively, to derive at CPEs of the mean and eigenfunctions within the BFPCA framework. The proposed CPEs are showcased in extensive simulations. Finally, the proposed CPEs are applied to the analysis of a sample of power spectral densities (PSD) from resting state electroencephalography (EEG) where they lead to novel insights on diagnostic group differences among children diagnosed with autism spectrum disorder and their typically developing peers across age.

How does obstructive sleep apnea alter cerebral hemodynamics?

STUDY OBJECTIVES

We aimed to characterize the cerebral hemodynamic response to obstructive sleep apnea/hypopnea events, and evaluate their association to polysomnographic parameters. The characterization of the cerebral hemodynamics in obstructive sleep apnea (OSA) may add complementary information to further the understanding of the severity of the syndrome beyond the conventional polysomnography.

METHODS

Severe OSA patients were studied during night sleep while monitored by polysomnography. Transcranial, bed-side diffuse correlation spectroscopy (DCS) and frequency-domain near-infrared diffuse correlation spectroscopy (NIRS-DOS) were used to follow microvascular cerebral hemodynamics in the frontal lobes of the cerebral cortex. Changes in cerebral blood flow (CBF), total hemoglobin concentration (THC), and cerebral blood oxygen saturation (StO2) were analyzed.

RESULTS

We considered 3283 obstructive apnea/hypopnea events from sixteen OSA patients (Age (median, interquartile range) 57 (52-64.5); females 25%; AHI (apnea-hypopnea index) 84.4 (76.1-93.7)). A biphasic response (maximum/minimum followed by a minimum/maximum) was observed for each cerebral hemodynamic variable (CBF, THC, StO2), heart rate and peripheral arterial oxygen saturation (SpO2). Changes of the StO2 followed the dynamics of the SpO2, and were out of phase from the THC and CBF. Longer events were associated with larger CBF changes, faster responses and slower recoveries. Moreover, the extrema of the response to obstructive hypopneas were lower compared to apneas (p < .001).

CONCLUSIONS

Obstructive apneas/hypopneas cause profound, periodic changes in cerebral hemodynamics, including periods of hyper- and hypo-perfusion and intermittent cerebral hypoxia. The duration of the events is a strong determinant of the cerebral hemodynamic response, which is more pronounced in apnea than hypopnea events.

Band depth based initialization of K-means for functional data clustering

The k-Means algorithm is one of the most popular choices for clustering data but is well-known to be sensitive to the initialization process. There is a substantial number of methods that aim at finding optimal initial seeds for k-Means, though none of them is universally valid. This paper presents an extension to longitudinal data of one of such methods, the BRIk algorithm, that relies on clustering a set of centroids derived from bootstrap replicates of the data and on the use of the versatile Modified Band Depth. In our approach we improve the BRIk method by adding a step where we fit appropriate B-splines to our observations and a resampling process that allows computational feasibility and handling issues such as noise or missing data. We have derived two techniques for providing suitable initial seeds, each of them stressing respectively the multivariate or the functional nature of the data. Our results with simulated and real data sets indicate that our Functional Data Approach to the BRIK method (FABRIk) and our Functional Data Extension of the BRIK method (FDEBRIk) are more effective than previous proposals at providing seeds to initialize k-Means in terms of clustering recovery.

Breathing patterns recognition: A functional data analysis approach

BACKGROUND AND OBJECTIVE

The ongoing pandemic proved fundamental is to assess a subject's respiratory functionality and breathing pattern measurement during quiet breathing is feasible in almost all patients, even those uncooperative. Breathing pattern consists of tidal volume and respiratory rate in an individual assessed by data tracks of lung or chest wall volume over time. State-of-art analysis of these data requires operator-dependent choices such as individuation of local minima in the track, elimination of anomalous breaths and individuation of breath clusters corresponding to different breathing patterns.

METHODS

A semi-automatic, robust and reproducible procedure was proposed to pre-process and analyse respiratory tracks, based on Functional Data Analysis (FDA) techniques, to identify representative breath curve and the corresponding breathing patterns. This was achieved through three steps: 1) breath separation through precise localization of the minima of the volume trace; 2) functional outlier breaths detection according to time-duration, magnitude and shape; 3) breath clustering to identify different pattern of interest, through K-medoids with Alignment. The method was firstly validated on simulated tracks and then applied to real data in conditions of clinical interest: operational volume change, exercise, mechanical ventilation, paradoxical breathing and age.

RESULTS

The total error in the accuracy of minima detection and in was less than 5%; with the artificial outliers being almost completely removed with an accuracy of 99%. During incremental exercise and independently on the bike resistance level, five clusters were identified (quiet breathing; recovery phase; onset of exercise; maximal and intermediate levels of exercise). During mechanical ventilation, the procedure was able to separate the non-ventilated from the ventilatory-supported breathing and to identify the worsening of paradoxical breathing due to the disease progression and the breathing pattern changes in healthy subjects due to age.

CONCLUSIONS

We proposed a robust validated automatic breathing patterns identification algorithm that extracted representative curves that could be implemented in clinical practice for objective comparison of the breathing patterns within and between subjects. In all case studies the identified patterns proved to be coherent with the clinical conditions and the physiopathology of the subjects, therefore enforcing the potential clinical translational value of the method.

Depthgram: Visualizing outliers in high-dimensional functional data with application to fMRI data exploration

Functional magnetic resonance imaging (fMRI) is a non‐invasive technique that facilitates the study of brain activity by measuring changes in blood flow. Brain activity signals can be recorded during the alternate performance of given tasks, that is, task fMRI (tfMRI), or during resting‐state, that is, resting‐state fMRI (rsfMRI), as a measure of baseline brain activity. This contributes to the understanding of how the human brain is organized in functionally distinct subdivisions. fMRI experiments from high‐resolution scans provide hundred of thousands of longitudinal signals for each individual, corresponding to brain activity measurements over each voxel of the brain along the duration of the experiment. In this context, we propose novel visualization techniques for high‐dimensional functional data relying on depth‐based notions that enable computationally efficient 2‐dim representations of fMRI data, which elucidate sample composition, outlier presence, and individual variability. We believe that this previous step is crucial to any inferential approach willing to identify neuroscientific patterns across individuals, tasks, and brain regions. We present the proposed technique via an extensive simulation study, and demonstrate its application on a motor and language tfMRI experiment.

A Geometric Perspective on Functional Outlier Detection

We consider functional outlier detection from a geometric perspective, specifically: for functional datasets drawn from a functional manifold, which is defined by the data’s modes of variation in shape, translation, and phase. Based on this manifold, we developed a conceptualization of functional outlier detection that is more widely applicable and realistic than previously proposed taxonomies. Our theoretical and experimental analyses demonstrated several important advantages of this perspective: it considerably improves theoretical understanding and allows describing and analyzing complex functional outlier scenarios consistently and in full generality, by differentiating between structurally anomalous outlier data that are off-manifold and distributionally outlying data that are on-manifold, but at its margins. This improves the practical feasibility of functional outlier detection: we show that simple manifold-learning methods can be used to reliably infer and visualize the geometric structure of functional datasets. We also show that standard outlier-detection methods requiring tabular data inputs can be applied to functional data very successfully by simply using their vector-valued representations learned from manifold learning methods as the input features. Our experiments on synthetic and real datasets demonstrated that this approach leads to outlier detection performances at least on par with existing functional-data-specific methods in a large variety of settings, without the highly specialized, complex methodology and narrow domain of application these methods often entail.

ExpanDrogram: Dynamic Visualization of Big Data Segmentation over Time

In dynamic and big data environments the visualization of a segmentation process over time often does not enable the user to simultaneously track entire pieces. The key points are sometimes incomparable, and the user is limited to a static visual presentation of a certain point. The proposed visualization concept, called ExpanDrogram, is designed to support dynamic classifiers that run in a big data environment subject to changes in data characteristics. It offers a wide range of features that seek to maximize the customization of a segmentation problem. The main goal of the ExpanDrogram visualization is to improve comprehensiveness by combining both the individual and segment levels, illustrating the dynamics of the segmentation process over time, providing “version control” that enables the user to observe the history of changes, and more. The method is illustrated using different datasets, with which we demonstrate multiple segmentation parameters, as well as multiple display layers, to highlight points such as new trend detection, outlier detection, tracking changes in original segments, and zoom in/out for more/less detail. The datasets vary in size from a small one to one of more than 12 million records.

A depth-based global envelope test for comparing two groups of functions with applications to biomedical data

Functional data are commonly observed in many emerging biomedical fields and their analysis is an exciting developing area in statistics. Numerous statistical methods, such as principal components, analysis of variance, and linear regression, have been extended to functional data. The statistical analysis of functions can be significantly improved using nonparametric and robust estimators. New ideas of depth for functional data have been proposed in recent years and can be extended to image data. They provide a way of ordering curves or images from center-outward, and of defining robust order statistics in a functional context. In this paper we develop depth-based global envelope tests for comparing two groups of functions or images. In addition to providing global P-values, the proposed envelope test can be displayed graphically and indicates the specific portion(s) of the functional data (eg, in pixels or in time) that may have led to rejection of the null hypothesis. We show in a simulation study the performance of the envelope test in terms of empirical power and size in different scenarios. The proposed depth-based global approach has good power even for small differences and is robust to outliers. The methodology introduced is applied to test whether children with normal and low birth weight have similar growth pattern. We also analyzed a brain image dataset consisting of positron emission tomography scans of severe depressed patients and healthy controls. The global envelope test was used to find and visualize differences between the two groups.

Is it time to stop sweeping data cleaning under the carpet? A novel algorithm for outlier management in growth data

All data are prone to error and require data cleaning prior to analysis. An important example is longitudinal growth data, for which there are no universally agreed standard methods for identifying and removing implausible values and many existing methods have limitations that restrict their usage across different domains. A decision-making algorithm that modified or deleted growth measurements based on a combination of pre-defined cut-offs and logic rules was designed. Five data cleaning methods for growth were tested with and without the addition of the algorithm and applied to five different longitudinal growth datasets: four uncleaned canine weight or height datasets and one pre-cleaned human weight dataset with randomly simulated errors. Prior to the addition of the algorithm, data cleaning based on non-linear mixed effects models was the most effective in all datasets and had on average a minimum of 26.00% higher sensitivity and 0.12% higher specificity than other methods. Data cleaning methods using the algorithm had improved data preservation and were capable of correcting simulated errors according to the gold standard; returning a value to its original state prior to error simulation. The algorithm improved the performance of all data cleaning methods and increased the average sensitivity and specificity of the non-linear mixed effects model method by 7.68% and 0.42% respectively. Using non-linear mixed effects models combined with the algorithm to clean data allows individual growth trajectories to vary from the population by using repeated longitudinal measurements, identifies consecutive errors or those within the first data entry, avoids the requirement for a minimum number of data entries, preserves data where possible by correcting errors rather than deleting them and removes duplications intelligently. This algorithm is broadly applicable to data cleaning anthropometric data in different mammalian species and could be adapted for use in a range of other domains.

Characterization of the microvascular cerebral blood flow response to obstructive apneic events during night sleep

Obstructive apnea causes periodic changes in cerebral and systemic hemodynamics, which may contribute to the increased risk of cerebrovascular disease of patients with obstructive sleep apnea (OSA) syndrome. The improved understanding of the consequences of an apneic event on the brain perfusion may improve our knowledge of these consequences and then allow for the development of preventive strategies. Our aim was to characterize the typical microvascular, cortical cerebral blood flow (CBF) changes in an OSA population during an apneic event. Sixteen patients (age 58 ± 8    years , 75% male) with a high risk of severe OSA were measured with a polysomnography device and with diffuse correlation spectroscopy (DCS) during one night of sleep with 1365 obstructive apneic events detected. All patients were later confirmed to suffer from severe OSA syndrome with a mean of 83 ± 15 apneas and hypopneas per hour. DCS has been shown to be able to characterize the microvascular CBF response to each event with a sufficient contrast-to-noise ratio to reveal its dynamics. It has also revealed that an apnea causes a peak increase of microvascular CBF ( 30 ± 17 % ) at the end of the event followed by a drop ( - 20 ± 12 % ) similar to what was observed in macrovascular CBF velocity of the middle cerebral artery. This study paves the way for the utilization of DCS for further studies on these populations.

Concurrent measurement of cerebral hemodynamics and electroencephalography during transcranial direct current stimulation

Transcranial direct current stimulation (tDCS) is currently being used for research and treatment of some neurological and neuropsychiatric disorders, as well as for improvement of cognitive functions. In order to better understand cerebral response to the stimulation and to redefine protocols and dosage, its effects must be monitored. To this end, we have used functional diffuse correlation spectroscopy (fDCS) and time-resolved functional near-infrared spectroscopy (TR-fNIRS) together with electroencephalography (EEG) during and after stimulation of the frontal cortex. Twenty subjects participated in two sessions of stimulation with two different polarity montages and twelve also underwent a sham session. Cerebral blood flow and oxyhemoglobin concentration increased during and after active stimulation in the region under the stimulation electrode while deoxyhemoglobin concentration decreased. The EEG spectrum displayed statistically significant power changes across different stimulation sessions in delta (2 to 4 Hz), theta (4 to 8 Hz), and beta (12 to 18 Hz) bands. Results suggest that fDCS and TR-fNIRS can be employed as neuromonitors of the effects of transcranial electrical stimulation and can be used together with EEG.

Error and anomaly detection for intra-participant time-series data

Identification of errors or anomalous values, collectively considered outliers, assists in exploring data or through removing outliers improves statistical analysis. In biomechanics, outlier detection methods have explored the ‘shape’ of the entire cycles, although exploring fewer points using a ‘moving-window’ may be advantageous. Hence, the aim was to develop a moving-window method for detecting trials with outliers in intra-participant time-series data. Outliers were detected through two stages for the strides (mean 38 cycles) from treadmill running. Cycles were removed in stage 1 for one-dimensional (spatial) outliers at each time point using the median absolute deviation, and in stage 2 for two-dimensional (spatial–temporal) outliers using a moving window standard deviation. Significance levels of the t-statistic were used for scaling. Fewer cycles were removed with smaller scaling and smaller window size, requiring more stringent scaling at stage 1 (mean 3.5 cycles removed for 0.0001 scaling) than at stage 2 (mean 2.6 cycles removed for 0.01 scaling with a window size of 1). Settings in the supplied Matlab code should be customised to each data set, and outliers assessed to justify whether to retain or remove those cycles. The method is effective in identifying trials with outliers in intra-participant time series data.

Dissimilarity for functional data clustering based on smoothing parameter commutation

Many studies measure the same type of information longitudinally on the same subject at multiple time points, and clustering of such functional data has many important applications. We propose a novel and easy method to implement dissimilarity measure for functional data clustering based on smoothing splines and smoothing parameter commutation. This method handles data observed at regular or irregular time points in the same way. We measure the dissimilarity between subjects based on varying curve estimates with pairwise commutation of smoothing parameters. The intuition is that smoothing parameters of smoothing splines reflect the inverse of the signal-to-noise ratios and that when applying an identical smoothing parameter the smoothed curves for two similar subjects are expected to be close. Our method takes into account the estimation uncertainty using smoothing parameter commutation and is not strongly affected by outliers. It can also be used for outlier detection. The effectiveness of our proposal is shown by simulations comparing it to other dissimilarity measures and by a real application to methadone dosage maintenance levels.

Data correction to determine a representative pattern of a set of 3D knee kinematic measurements

The purpose of this study is to determine a representative pattern of a set of three dimensional (3D) knee kinematic measurement curves recorded throughout several trials with a patient walking on a treadmill. The measurements are knee angles, (namely joint angles) with respect to the sagittal, frontal, and transverse planes, as a function of time during a gait cycle. Two serious difficulties met while extracting a representative pattern from the trials are that the curves possess phase variability and there are outliers. We propose a scheme which first removes outliers using the modified band depth index method, and follows with phase variability reduction by curve registration. This scheme leads to retaining the mean curve of the corrected set of curves, as the most representative.