Evaluating phase synchronization methods in fMRI: A comparison study and new approaches.
Ontology highlight
ABSTRACT:
SUBMITTER: Honari H
PROVIDER: S-EPMC8011682 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Synchronization or phase-locking between oscillating neuronal groups is considered to be important for coordination of information among cortical networks. Spectral coherence is a commonly used approach to quantify phase locking between neural signals. We systematically explored the validity of spectral coherence measures for quantifying synchronization among neural oscillators. To that aim, we simulated coupled oscillatory signals that exhibited synchronization dynamics using an abstract phase-oscillator model as well as interacting gamma-generating spiking neural networks. We found that, within a large parameter range, the spectral coherence measure deviated substantially from the expected phase-locking. Moreover, spectral coherence did not converge to the expected value with increasing signal-to-noise ratio. We found that spectral coherence particularly failed when oscillators were in the partially (intermittent) synchronized state, which we expect to be the most likely state for neural synchronization. The failure was due to the fast frequency and amplitude changes induced by synchronization forces. We then investigated whether spectral coherence reflected the information flow among networks measured by transfer entropy (TE) of spike trains. We found that spectral coherence failed to robustly reflect changes in synchrony-mediated information flow between neural networks in many instances. As an alternative approach we explored a phase-locking value (PLV) method based on the reconstruction of the instantaneous phase. As one approach for reconstructing instantaneous phase, we used the Hilbert Transform (HT) preceded by Singular Spectrum Decomposition (SSD) of the signal. PLV estimates have broad applicability as they do not rely on stationarity, and, unlike spectral coherence, they enable more accurate estimations of oscillatory synchronization across a wide range of different synchronization regimes, and better tracking of synchronization-mediated information flow among networks.
| S-EPMC4706353 | biostudies-literature
Project description:Detecting sight-threatening retinopathy of prematurity (ROP) relies on a diagnostic examination (DE) performed by an experienced ophthalmologist. An alternative may be a telemedicine system where retinal images of at-risk infants are graded by readers to determine features of ROP indicating the need for a DE.The multicenter Telemedicine Approaches to Evaluating Acute-phase ROP (e-ROP) Study is a cohort study of 2000 infants with birth weights <1251 g. At each visit, ophthalmologists perform DEs and non-physician imagers obtain iris and five retinal images with the disc positioned in the center, right, left, up and down. Images are uploaded to a secure server for grading by non-physician readers for the detection of plus disease, stage 3 ROP and/or zone I disease, any of which indicates "referral-warranted ROP" (RW-ROP). Images from all infants with RW-ROP and a random sample of infants without RW-ROP (based on DEs) are selected for grading. Gradings are compared to DEs to determine the validity and evaluate reliability, feasibility, safety, and cost-effectiveness of the telemedicine system.e-ROP is conducted in 12 Clinical Centers in the US and Canada with Study Headquarters, the Data Coordinating Center and the Image Reading Center in Philadelphia and the ROP Data Center in Oklahoma City. A total of 27 study center coordinators, 34 ophthalmologists, 26 imagers, and 4 readers have been certified. All study data are submitted using a secure web-based system.The design and findings of this study will be useful to conduct other ROP studies or evaluate telemedicine for other diseases.
| S-EPMC4861056 | biostudies-literature
Project description:Noise is a major challenge for the analysis of fMRI data in general and for connectivity analyses in particular. As researchers develop increasingly sophisticated tools to model statistical dependence between the fMRI signal in different brain regions, there is a risk that these models may increasingly capture artifactual relationships between regions, that are the result of noise. Thus, choosing optimal denoising methods is a crucial step to maximize the accuracy and reproducibility of connectivity models. Most comparisons between denoising methods require knowledge of the ground truth: of what is the 'real signal'. For this reason, they are usually based on simulated fMRI data. However, simulated data may not match the statistical properties of real data, limiting the generalizability of the conclusions. In this article, we propose an approach to evaluate denoising methods using real (non-simulated) fMRI data. First, we introduce an intersubject version of multivariate pattern dependence (iMVPD) that computes the statistical dependence between a brain region in one participant, and another brain region in a different participant. iMVPD has the following advantages: 1) it is multivariate, 2) it trains and tests models on independent partitions of the real fMRI data, and 3) it generates predictions that are both between subjects and between regions. Since whole-brain sources of noise are more strongly correlated within subject than between subjects, we can use the difference between standard MVPD and iMVPD as a 'discrepancy metric' to evaluate denoising techniques (where more effective techniques should yield smaller differences). As predicted, the difference is the greatest in the absence of denoising methods. Furthermore, a combination of removal of the global signal and CompCorr optimizes denoising (among the set of denoising options tested).
| S-EPMC6759145 | biostudies-literature
Project description:Sequencing technologies have enabled the investigation of whole genomes of many individuals in parallel. Studies have shown that the joint consideration of multiple rare variants may explain a relevant proportion of the genetic basis for disease so that grouping of rare variants, termed collapsing, can enrich the association signal. Following this assumption, we investigate the type I error and the power of two proposed collapsing methods (combined multivariate and collapsing method and the functional principal component analysis [FPCA]-based statistic) using the case-control data provided for the Genetic Analysis Workshop 18 with knowledge of the true model. Variants with a minor allele frequency (MAF) of 0.05 or less were collapsed per gene for combined multivariate and collapsing. Neither of the methods detected any of the truly associated genes reliably. Although combined multivariate and collapsing identified one gene with a power of 0.66, it had an unacceptably high false-positive rate of 75%. In contrast, FPCA covered the type I error level well but at the cost of low power. A strict filtering of variants by small MAF might lead to a better performance of the collapsing methods. Furthermore, the inclusion of information on functionality of the variants could be helpful.
| S-EPMC4143760 | biostudies-literature
Project description:Analysis methods in cognitive neuroscience have not always matched the richness of fMRI data. Early methods focused on estimating neural activity within individual voxels or regions, averaged over trials or blocks and modeled separately in each participant. This approach mostly neglected the distributed nature of neural representations over voxels, the continuous dynamics of neural activity during tasks, the statistical benefits of performing joint inference over multiple participants and the value of using predictive models to constrain analysis. Several recent exploratory and theory-driven methods have begun to pursue these opportunities. These methods highlight the importance of computational techniques in fMRI analysis, especially machine learning, algorithmic optimization and parallel computing. Adoption of these techniques is enabling a new generation of experiments and analyses that could transform our understanding of some of the most complex-and distinctly human-signals in the brain: acts of cognition such as thoughts, intentions and memories.
| S-EPMC5457304 | biostudies-other
Project description:A significant challenge for fMRI research is statistically controlling for false positives without omitting true effects. Although a number of traditional methods for multiple comparison correction exist, several alternative tools have been developed that do not rely on strict parametric assumptions, but instead implement alternative methods to correct for multiple comparisons. In this study, we evaluated three of these methods, Statistical non-Parametric Mapping (SnPM), 3DClustSim, and Threshold Free Cluster Enhancement (TFCE), by examining which method produced the most consistent outcomes even when spatially-autocorrelated noise was added to the original images. We assessed the false alarm rate and hit rate of each method after noise was applied to the original images.
| S-EPMC6721788 | biostudies-literature
Project description:Many cognitive and clinical neuroscience research studies seek to determine how contextual factors modulate cognitive processes. In fMRI, hypotheses about how context modulates distributed patterns of information processing are often tested by comparing functional connectivity between neural regions A and B as a function of task conditions X and Y, which is termed context-modulated functional connectivity (FC). There exist two exploratory statistical approaches to testing context-modulated FC: the beta-series method and psychophysiological interaction (PPI) analysis methods. While these approaches are commonly used, their relative power for detecting context-modulated FC is unknown, especially with respect to real-world experimental parameters (e.g., number of stimulus repetitions, inter-trial-interval, stimulus duration). Here, we use simulations to compare power for detecting context-modulated FC between the standard PPI formulation (sPPI), generalized PPI formulation (gPPI), and beta series methods. Simulation results demonstrate that gPPI and beta series methods are generally more powerful than sPPI. Whether gPPI or beta series methods performed more powerfully depended on experiment parameters: block designs favor the gPPI, whereas the beta series method was more powerful for designs with more trial repetitions and it also retained more power under conditions of hemodynamic response function variability. On a real dataset of adolescent girls, the PPI methods appeared to have greater sensitivity in detecting task-modulated FC when using a block design and the beta series method appeared to have greater sensitivity when using an event-related design with many trial repetitions. Implications of these performance results are discussed.
| S-EPMC4019671 | biostudies-literature
Project description:A large number of fMRI studies have shown that the temporal dynamics of evoked BOLD responses can be highly heterogeneous. Failing to model heterogeneous responses in statistical analysis can lead to significant errors in signal detection and characterization and alter the neurobiological interpretation. However, to date it is not clear that, out of a large number of options, which methods are robust against variability in the temporal dynamics of BOLD responses in block-design studies. Here, we used rodent optogenetic fMRI data with heterogeneous BOLD responses and simulations guided by experimental data as a means to investigate different analysis methods' performance against heterogeneous BOLD responses. Evaluations are carried out within the general linear model (GLM) framework and consist of standard basis sets as well as independent component analysis (ICA). Analyses show that, in the presence of heterogeneous BOLD responses, conventionally used GLM with a canonical basis set leads to considerable errors in the detection and characterization of BOLD responses. Our results suggest that the 3rd and 4th order gamma basis sets, the 7th to 9th order finite impulse response (FIR) basis sets, the 5th to 9th order B-spline basis sets, and the 2nd to 5th order Fourier basis sets are optimal for good balance between detection and characterization, while the 1st order Fourier basis set (coherence analysis) used in our earlier studies show good detection capability. ICA has mostly good detection and characterization capabilities, but detects a large volume of spurious activation with the control fMRI data.
| S-EPMC5472444 | biostudies-literature
Project description:The robustness of the phase retrieval methods is of critical importance for limiting and reducing radiation doses involved in x-ray phase contrast imaging. This work is to compare the robustness of two phase retrieval methods by analyzing the phase maps retrieved from the experimental images of a phantom.Two phase retrieval methods were compared. One method is based on the transport of intensity equation (TIE) for phase contrast projections, and the TIE-based method is the most commonly used method for phase retrieval in the literature. The other is the recently developed attenuation-partition based (AP-based) phase retrieval method. The authors applied these two methods to experimental projection images of an air-bubble wrap phantom for retrieving the phase map of the bubble wrap. The retrieved phase maps obtained by using the two methods are compared.In the wrap's phase map retrieved by using the TIE-based method, no bubble is recognizable, hence, this method failed completely for phase retrieval from these bubble wrap images. Even with the help of the Tikhonov regularization, the bubbles are still hardly visible and buried in the cluttered background in the retrieved phase map. The retrieved phase values with this method are grossly erroneous. In contrast, in the wrap's phase map retrieved by using the AP-based method, the bubbles are clearly recovered. The retrieved phase values with the AP-based method are reasonably close to the estimate based on the thickness-based measurement. The authors traced these stark performance differences of the two methods to their different techniques employed to deal with the singularity problem involved in the phase retrievals.This comparison shows that the conventional TIE-based phase retrieval method, regardless if Tikhonov regularization is used or not, is unstable against the noise in the wrap's projection images, while the AP-based phase retrieval method is shown in these experiments to be superior to the TIE-based method for the robustness in performing the phase retrieval.
| S-EPMC3170396 | biostudies-other
Project description:Absence seizures-generalized rhythmic spike-and-wave discharges (SWDs) are the defining property of childhood (CAE) and juvenile (JAE) absence epilepsies. Such seizures are the most compelling examples of pathological neuronal hypersynchrony. All the absence detection algorithms proposed so far have been derived from the properties of individual SWDs. In this work, we investigate EEG phase synchronization in patients with CAE/JAE and healthy subjects to explore the possibility of using the wavelet phase synchronization index to detect seizures and quantify their disorganization (fragmentation). The overlap of the ictal and interictal probability density functions was high enough to preclude effective seizure detection based solely on changes in EEG synchronization. We used a machine learning classifier with the phase synchronization index (calculated for 1 s data segments with 0.5 s overlap) and the normalized amplitude as features to detect generalized SWDs. Using 19 channels (10-20 setup), we identified 99.2% of absences. However, the overlap of the segments classified as ictal with seizures was only 83%. The analysis showed that seizures were disorganized in approximately half of the 65 subjects. On average, generalized SWDs lasted about 80% of the duration of abnormal EEG activity. The disruption of the ictal rhythm can manifest itself as the disappearance of epileptic spikes (with high-amplitude delta waves persisting), transient cessation of epileptic discharges, or loss of global synchronization. The detector can analyze a real-time data stream. Its performance is good for a six-channel setup (Fp1, Fp2, F7, F8, O1, O2), which can be implemented as an unobtrusive EEG headband. False detections are rare for controls and young adults (0.03% and 0.02%, respectively). In patients, they are more frequent (0.5%), but in approximately 82% cases, classification errors are caused by short epileptiform discharges. Most importantly, the proposed detector can be applied to parts of EEG with abnormal EEG activity to quantitatively determine seizure fragmentation. This property is important because a previous study reported that the probability of disorganized discharges is eight times higher in JAE than in CAE. Future research must establish whether seizure properties (frequency, length, fragmentation, etc.) and clinical characteristics can help distinguish CAE and JAE.
| S-EPMC10317177 | biostudies-literature