BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Alberto Sorrentino

Alberto Sorrentino appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
Applicationsmath.NAMethodologyQuantitative MethodsNumerical Analysisphysics.ao-phphysics.optics

Catalog footprint

What is connected

7works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2020arXiv

Where Bayes tweaks Gauss: Conditionally Gaussian priors for stable multi-dipole estimation

We present a very simple yet powerful generalization of a previously described model and algorithm for estimation of multiple dipoles from magneto/electro-encephalographic data. Specifically, the generalization consists in the introduction of a log-uniform hyperprior on the standard deviation of a set of conditionally linear/Gaussian variables. We use numerical simulations and an experimental dataset to show that the approximation to the posterior distribution remains extremely stable under a wide range of values of the hyperparameter, virtually removing the dependence on the hyperparameter.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On the two-step estimation of the cross--power spectrum for dynamical inverse problems

We consider the problem of reconstructing the cross--power spectrum of an unobservable multivariate stochatic process from indirect measurements of a second multivariate stochastic process, related to the first one through a linear operator. In the two--step approach, one would first compute a regularized reconstruction of the unobservable signal, and then compute an estimate of its cross--power spectrum from the regularized solution. We investigate whether the optimal regularization parameter for reconstruction of the signal also gives the best estimate of the cross--power spectrum. We show that the answer depends on the regularization method, and specifically we prove that, under a white Gaussian assumption: (i) when regularizing with truncated SVD the optimal parameter is the same; (ii) when regularizing with the Tikhonov method, the optimal parameter for the cross--power spectrum is lower than half the optimal parameter for the signal. We also provide evidence that a one--step approach would likely have better mathematical properties of the two--step approach. Our results apply particularly to the brain connectivity estimation from magneto/electro-encephalographic recordings and provide a formal interpretation of recent empirical results.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bayesian smoothing of dipoles in Magneto-/Electro-encephalography

We describe a novel method for dynamic estimation of multi-dipole states from Magneto/Electro-encephalography (M/EEG) time series. The new approach builds on the recent development of particle filters for M/EEG; these algorithms approximate, with samples and weights, the posterior distribution of the neural sources at time t given the data up to time t. However, for off-line inference purposes it is preferable to work with the smoothing distribution, i.e. the distribution for the neural sources at time t conditioned on the whole time series. In this study, we use a Monte Carlo algorithm to approximate the smoothing distribution for a time-varying set of current dipoles. We show, using numerical simulations, that the estimates provided by the smoothing distribution are more accurate than those provided by the filtering distribution, particularly at the appearance of the source. We validate the proposed algorithm using an experimental dataset recorded from an epileptic patient. Improved localization of the source onset can be particularly relevant in source modeling of epileptic patients, where the source onset brings information on the epileptogenic zone.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Expectation Maximization and the retrieval of the atmospheric extinction coefficients by inversion of Raman lidar data

We consider the problem of retrieving the aerosol extinction coefficient from Raman lidar measurements. This is an ill--posed inverse problem that needs regularization, and we propose to use the Expectation--Maximization (EM) algorithm to provide stable solutions. Indeed, EM is an iterative algorithm that imposes a positivity constraint on the solution, and provides regularization if iterations are stopped early enough. We describe the algorithm and propose a stopping criterion inspired by a statistical principle. We then discuss its properties concerning the spatial resolution. Finally, we validate the proposed approach by using both synthetic data and experimental measurements; we compare the reconstructions obtained by EM with those obtained by the Tikhonov method, by the Levenberg-Marquardt method, as well as those obtained by combining data smoothing and numerical derivation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Bayesian Multi--Dipole Modeling of a Single Topography in MEG by Adaptive Sequential Monte Carlo Samplers

In the present paper, we develop a novel Bayesian approach to the problem of estimating neural currents in the brain from a fixed distribution of magnetic field (called \emph{topography}), measured by magnetoencephalography. Differently from recent studies that describe inversion techniques, such as spatio-temporal regularization/filtering, in which neural dynamics always plays a role, we face here a purely static inverse problem. Neural currents are modelled as an unknown number of current dipoles, whose state space is described in terms of a variable--dimension model. Within the resulting Bayesian framework, we set up a sequential Monte Carlo sampler to explore the posterior distribution. An adaptation technique is employed in order to effectively balance the computational cost and the quality of the sample approximation. Then, both the number and the parameters of the unknown current dipoles are simultaneously estimated. The performance of the method is assessed by means of synthetic data, generated by source configurations containing up to four dipoles. Eventually, we describe the results obtained by analyzing data from a real experiment, involving somatosensory evoked fields, and compare them to those provided by three other methods.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sequential Monte Carlo samplers for semilinear inverse problems and application to magnetoencephalography

We discuss the use of a recent class of sequential Monte Carlo methods for solving inverse problems characterized by a semi-linear structure, i.e. where the data depend linearly on a subset of variables and nonlinearly on the remaining ones. In this type of problems, under proper Gaussian assumptions one can marginalize the linear variables. This means that the Monte Carlo procedure needs only to be applied to the nonlinear variables, while the linear ones can be treated analytically; as a result, the Monte Carlo variance and/or the computational cost decrease. We use this approach to solve the inverse problem of magnetoencephalography, with a multi-dipole model for the sources. Here, data depend nonlinearly on the number of sources and their locations, and depend linearly on their current vectors. The semi-analytic approach enables us to estimate the number of dipoles and their location from a whole time-series, rather than a single time point, while keeping a low computational cost.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamic filtering of static dipoles in magnetoencephalography

We consider the problem of estimating neural activity from measurements of the magnetic fields recorded by magnetoencephalography. We exploit the temporal structure of the problem and model the neural current as a collection of evolving current dipoles, which appear and disappear, but whose locations are constant throughout their lifetime. This fully reflects the physiological interpretation of the model. In order to conduct inference under this proposed model, it was necessary to develop an algorithm based around state-of-the-art sequential Monte Carlo methods employing carefully designed importance distributions. Previous work employed a bootstrap filter and an artificial dynamic structure where dipoles performed a random walk in space, yielding nonphysical artefacts in the reconstructions; such artefacts are not observed when using the proposed model. The algorithm is validated with simulated data, in which it provided an average localisation error which is approximately half that of the bootstrap filter. An application to complex real data derived from a somatosensory experiment is presented. Assessment of model fit via marginal likelihood showed a clear preference for the proposed model and the associated reconstructions show better localisation.

0 citations0 reviewsNo code linkedOpen access