Researcher profile

Sara Sommariva

Sara Sommariva contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

5 published item(s)

preprint2026arXiv

A neighbour selection approach for identifying differential networks in conditional functional graphical models

Estimation of brain functional connectivity from EEG data is of great importance both for medical research and diagnosis. It involves quantifying the conditional dependencies among the activity of different brain areas from the time-varying electric field recorded by sensors placed outside the scalp. These dependencies may vary within and across individuals and be influenced by covariates such as age, mental status, or disease severity. Motivated by this problem, we propose a novel neighbour selection approach based on functional-on-functional regression for the characterization of conditional Gaussian functional graphical models. We provide a fully automated, data-driven procedure for inferring conditional dependence structures among observed functional variables. In particular, pairwise interactions are directly identified and allowed to vary as a function of covariates, enabling covariate-specific modulation of connectivity patterns. Our proposed method accommodates an arbitrary number of continuous and discrete covariates. Moreover, unlike existing methods for direct estimation of differential graphical models, the proposed approach yields directly interpretable coefficients, allowing discrimination between covariate-induced increases and decreases in interaction strength. The methodology is evaluated through extensive simulation studies and an application to experimental EEG data. The results demonstrate clear advantages over existing approaches, including higher estimation accuracy and substantially reduced computational cost, especially in high-dimensional settings.

preprint2022arXiv

A fast and convergent combined Newton and gradient descent method for computing steady states of chemical reaction networks

In this work we present a fast, globally convergent, iterative algorithm for computing the asymptotically stable states of nonlinear large--scale systems of quadratic autonomous Ordinary Differential Equations (ODEs) modeling, e.g., the dynamic of complex chemical reaction networks. Towards this aim, we reformulate the problem as a box--constrained optimization problem where the roots of a set of nonlinear equations need to be determined. Then, we propose to use a projected Newton's approach combined with a gradient descent algorithm so that every limit point of the sequence generated by the overall algorithm is a stationary point. More importantly, we suggest replacing the standard orthogonal projector with a novel operator that ensures the final solution to satisfy the box constraints while lowering the probability that the intermediate points reached at each iteration belong to the boundary of the box where the Jacobian of the objective function may be singular. The effectiveness of the proposed approach is shown in a practical scenario concerning a chemical reaction network modeling the signaling network of colorectal cancer cells. Specifically, in this scenario the proposed algorithm is proven to be faster and more accurate than a classical dynamical approach where the asymptotically stable states are computed as the limit points of the flux of the Cauchy problem associated with the ODEs system.

preprint2020arXiv

Gain and Loss of Function mutations in biological chemical reaction networks: a mathematical model with application to colorectal cancer cells

This paper studies a system of Ordinary Differential Equations modeling a chemical reaction network and derives from it a simulation tool mimicking Loss of Function and Gain of Function mutations found in cancer cells. More specifically, from a theoretical perspective, our approach focuses on the determination of moiety conservation laws for the system and their relation with the corresponding stoichiometric surfaces. Then we show that Loss of Function mutations can be implemented in the model via modification of the initial conditions in the system, while Gain of Function mutations can be implemented by eliminating specific reactions. Finally, the model is utilized to examine in detail the G1-S phase of a colorectal cancer cell.

preprint2020arXiv

Mathematical modelling of nuclear medicine data

Positron Emission Tomography using 2-[18F]-2deoxy-D-glucose as radiotracer (FDG-PET) is currently one of the most frequently applied functional imaging methods in clinical applications. The interpretation of FDG-PET data requires sophisticated mathematical approaches able to exploit the dynamical information contained in this kind of data. Most of these approaches are formulated within the framework of compartmental analysis, which connects the experimental nuclear data with unknown tracer coefficients measuring the effectiveness of the tracer metabolism by means of Cauchy systems of ordinary differential equations. This paper provides a coincise overview of linear compartmental methods, focusing on the analytical solution of the forward compartmental problem and on the specific issues concerning the corresponding compartmental inverse problem.

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.