skip to main content


Title: Mining the Mind: Linear Discriminant Analysis of MEG Source Reconstruction Time Series Supports Dynamic Changes in Deep Brain Regions During Meditation Sessions
Abstract

Meditation practices have been claimed to have a positive effect on the regulation of mood and emotions for quite some time by practitioners, and in recent times there has been a sustained effort to provide a more precise description of the influence of meditation on the human brain. Longitudinal studies have reported morphological changes in cortical thickness and volume in selected brain regions due to meditation practice, which is interpreted as an evidence its effectiveness beyond the subjective self reporting. Using magnetoencephalography (MEG) or electroencephalography to quantify the changes in brain activity during meditation practice represents a challenge, as no clear hypothesis about the spatial or temporal pattern of such changes is available to date. In this article we consider MEG data collected during meditation sessions of experienced Buddhist monks practicing focused attention (Samatha) and open monitoring (Vipassana) meditation, contrasted by resting state with eyes closed. The MEG data are first mapped to time series of brain activity averaged over brain regions corresponding to a standard Destrieux brain atlas. Next, by bootstrapping and spectral analysis, the data are mapped to matrices representing random samples of power spectral densities in$$\alpha$$α,$$\beta$$β,$$\gamma$$γ, and$$\theta$$θfrequency bands. We use linear discriminant analysis to demonstrate that the samples corresponding to different meditative or resting states contain enough fingerprints of the brain state to allow a separation between different states, and we identify the brain regions that appear to contribute to the separation. Our findings suggest that the cingulate cortex, insular cortex and some of the internal structures, most notably the accumbens, the caudate and the putamen nuclei, the thalamus and the amygdalae stand out as separating regions, which seems to correlate well with earlier findings based on longitudinal studies.

 
more » « less
Award ID(s):
1951446
NSF-PAR ID:
10306307
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Brain Topography
Volume:
34
Issue:
6
ISSN:
0896-0267
Page Range / eLocation ID:
p. 840-862
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. A<sc>bstract</sc>

    The charge-parity (CP) structure of the Yukawa interaction between the Higgs (H) boson and the top quark is measured in a data sample enriched in the t$$ \overline{\textrm{t}} $$t¯H and tH associated production, using 138 fb1of data collected in proton-proton collisions at$$ \sqrt{s} $$s= 13 TeV by the CMS experiment at the CERN LHC. The study targets events where the H boson decays via H → WW or H →ττand the top quarks decay via t → Wb: the W bosons decay either leptonically or hadronically, and final states characterized by the presence of at least two leptons are studied. Machine learning techniques are applied to these final states to enhance the separation ofCP-even fromCP-odd scenarios. Two-dimensional confidence regions are set onκtand$$ \overset{\sim }{\kappa } $$κ~t, which are respectively defined as theCP-even andCP-odd top-Higgs Yukawa coupling modifiers. No significant fractionalCP-odd contributions, parameterized by the quantity|$$ {f}_{CP}^{\textrm{Htt}} $$fCPHtt|are observed; the parameter is determined to be|$$ {f}_{CP}^{\textrm{Htt}} $$fCPHtt|= 0.59 with an interval of (0.24,0.81) at 68% confidence level. The results are combined with previous results covering the HZZ and H→ γγdecay modes, yielding two- and one-dimensional confidence regions onκtand$$ \overset{\sim }{\kappa } $$κ~t, while|$$ {f}_{CP}^{\textrm{Htt}} $$fCPHtt|is determined to be|$$ {f}_{CP}^{\textrm{Htt}} $$fCPHtt|= 0.28 with an interval of|$$ {f}_{CP}^{\textrm{Htt}} $$fCPHtt| <0.55 at 68% confidence level, in agreement with the standard modelCP-even prediction of|$$ {f}_{CP}^{\textrm{Htt}} $$fCPHtt|= 0.

     
    more » « less
  2. Abstract

    Solar flares are explosions on the Sun. They happen when energy stored in magnetic fields around solar active regions (ARs) is suddenly released. Solar flares and accompanied coronal mass ejections are sources of space weather, which negatively affects a variety of technologies at or near Earth, ranging from blocking high-frequency radio waves used for radio communication to degrading power grid operations. Monitoring and providing early and accurate prediction of solar flares is therefore crucial for preparedness and disaster risk management. In this article, we present a transformer-based framework, named SolarFlareNet, for predicting whether an AR would produce a$$\gamma$$γ-class flare within the next 24 to 72 h. We consider three$$\gamma$$γclasses, namely the$$\ge$$M5.0 class, the$$\ge$$M class and the$$\ge$$C class, and build three transformers separately, each corresponding to a$$\gamma$$γclass. Each transformer is used to make predictions of its corresponding$$\gamma$$γ-class flares. The crux of our approach is to model data samples in an AR as time series and to use transformers to capture the temporal dynamics of the data samples. Each data sample consists of magnetic parameters taken from Space-weather HMI Active Region Patches (SHARP) and related data products. We survey flare events that occurred from May 2010 to December 2022 using the Geostationary Operational Environmental Satellite X-ray flare catalogs provided by the National Centers for Environmental Information (NCEI), and build a database of flares with identified ARs in the NCEI flare catalogs. This flare database is used to construct labels of the data samples suitable for machine learning. We further extend the deterministic approach to a calibration-based probabilistic forecasting method. The SolarFlareNet system is fully operational and is capable of making near real-time predictions of solar flares on the Web.

     
    more » « less
  3. Abstract

    There are many ways of measuring and modeling tail-dependence in random vectors: from the general framework of multivariate regular variation and the flexible class of max-stable vectors down to simple and concise summary measures like the matrix of bivariate tail-dependence coefficients. This paper starts by providing a review of existing results from a unifying perspective, which highlights connections between extreme value theory and the theory of cuts and metrics. Our approach leads to some new findings in both areas with some applications to current topics in risk management.

    We begin by using the framework of multivariate regular variation to show that extremal coefficients, or equivalently, the higher-order tail-dependence coefficients of a random vector can simply be understood in terms of random exceedance sets, which allows us to extend the notion of Bernoulli compatibility. In the special but important case of bivariate tail-dependence, we establish a correspondence between tail-dependence matrices and$$L^1$$L1- and$$\ell _1$$1-embeddable finite metric spaces via the spectral distance, which is a metric on the space of jointly 1-Fréchet random variables. Namely, the coefficients of the cut-decomposition of the spectral distance and of the Tawn-Molchanov max-stable model realizing the corresponding bivariate extremal dependence coincide. We show that line metrics are rigid and if the spectral distance corresponds to a line metric, the higher order tail-dependence is determined by the bivariate tail-dependence matrix.

    Finally, the correspondence between$$\ell _1$$1-embeddable metric spaces and tail-dependence matrices allows us to revisit the realizability problem, i.e. checking whether a given matrix is a valid tail-dependence matrix. We confirm a conjecture of Shyamalkumar and Tao (2020) that this problem is NP-complete.

     
    more » « less
  4. Abstract Background

    Many bioinformatics applications involve bucketing a set of sequences where each sequence is allowed to be assigned into multiple buckets. To achieve both high sensitivity and precision, bucketing methods are desired to assign similar sequences into the same bucket while assigning dissimilar sequences into distinct buckets. Existingk-mer-based bucketing methods have been efficient in processing sequencing data with low error rates, but encounter much reduced sensitivity on data with high error rates. Locality-sensitive hashing (LSH) schemes are able to mitigate this issue through tolerating the edits in similar sequences, but state-of-the-art methods still have large gaps.

    Results

    In this paper, we generalize the LSH function by allowing it to hash one sequence into multiple buckets. Formally, a bucketing function, which maps a sequence (of fixed length) into a subset of buckets, is defined to be$$(d_1, d_2)$$(d1,d2)-sensitive if any two sequences within an edit distance of$$d_1$$d1are mapped into at least one shared bucket, and any two sequences with distance at least$$d_2$$d2are mapped into disjoint subsets of buckets. We construct locality-sensitive bucketing (LSB) functions with a variety of values of$$(d_1,d_2)$$(d1,d2)and analyze their efficiency with respect to the total number of buckets needed as well as the number of buckets that a specific sequence is mapped to. We also prove lower bounds of these two parameters in different settings and show that some of our constructed LSB functions are optimal.

    Conclusion

    These results lay the theoretical foundations for their practical use in analyzing sequences with high error rates while also providing insights for the hardness of designing ungapped LSH functions.

     
    more » « less
  5. Abstract

    We study the behavior of solutions to the incompressible 2dEuler equations near two canonical shear flows with critical points, the Kolmogorov and Poiseuille flows, with consequences for the associated Navier–Stokes problems. We exhibit a large family of new, non-trivial stationary states that are arbitrarily close to the Kolmogorov flow on the square torus$$\mathbb {T}^2$$T2in analytic regularity. This situation contrasts strongly with the setting of some monotone shear flows, such as the Couette flow: there the linearized problem exhibits an “inviscid damping” mechanism that leads to relaxation of perturbations of the base flows back to nearby shear flows. Our results show that such a simple description of the long-time behavior is not possible for solutions near the Kolmogorov flow on$$\mathbb {T}^2$$T2. Our construction of the new stationary states builds on a degeneracy in the global structure of the Kolmogorov flow on$$\mathbb {T}^2$$T2, and we also show a lack of correspondence between the linearized description of the set of steady states and its true nonlinear structure. Both the Kolmogorov flow on a rectangular torus and the Poiseuille flow in a channel are very different. We show that the only stationary states near them must indeed be shears, even in relatively low regularity. In addition, we show that this behavior is mirrored closely in the related Navier–Stokes settings: the linearized problems near the Poiseuille and Kolmogorov flows both exhibit an enhanced rate of dissipation. Previous work by us and others shows that this effect survives in the full, nonlinear problem near the Poiseuille flow and near the Kolmogorov flow on rectangular tori, provided that the perturbations lie below a certain threshold. However, we show here that the corresponding result cannot hold near the Kolmogorov flow on$${\mathbb T}^2$$T2.

     
    more » « less