Source author record

Jorge P. Zubelli

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

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

11 published item(s)

preprint2022arXiv

Matrix Bispectrality and Noncommutative Algebras: beyond the prolate spheroidals

The bispectral problem is motivated by an effort to understand and extend a remarkable phenomenon in Fourier analysis on the real line: the operator of time-and-band limiting is an integral operator admitting a second-order differential operator with a simple spectrum in its commutator. In this article, we discuss a noncommutative version of the bispectral problem, obtained by allowing all objects in the original formulation to be matrix-valued. Deep attention is given to bispectral algebras and their presentations as a tool to get information about bispectral triples.

preprint2022arXiv

Noncommutative Bispectral Algebras and their Presentations

We prove a general result on presentations of finitely-generated algebras and apply it to obtain nice presentations for some noncommutative algebras arising in the matrix bispectral problem. By "nice presentation" we mean a presentation that has as few as possible defining relations. This in turn, has potential applications in computer algebra implementations and examples. Our results can be divided into three parts. In the first two, we consider bispectral algebras with the eigenvalue in the physical equation to be scalar-valued for $2\times 2$ and $3\times 3$ matrix-valued eigenfunctions. In the third part, we assume the eigenvalue in the physical equation to be matrix-valued and draw an important connection with spin Calogero-Moser systems. In all cases, we show that these algebras are finitely presented. As a byproduct, we answer positively a conjecture of F.~A.~Grünbaum about these algebras.

preprint2021arXiv

The Impact of COVID-19 Vaccination Delay: A Modelling Study for Chicago and NYC Data

Background: By the beginning of December 2020, some vaccines against COVID-19 already presented efficacy and security, which qualify them to be used in mass vaccination campaigns. Thus, setting up strategies of vaccination became crucial to control the COVID19 pandemic. Methods: We use daily COVID-19 reports from Chicago and NYC from 01-Mar2020 to 28- Nov-2020 to estimate the parameters of an SEIR-like epidemiological model that accounts for different severity levels. To achieve data adherent predictions, we let the model parameters to be time-dependent. The model is used to forecast different vaccination scenarios, where the campaign starts at different dates, from 01-Oct-2020 to 01-Apr-2021. To generate realistic scenarios, disease control strategies are implemented whenever the number of predicted daily hospitalizations reaches a preset threshold. Results: The model reproduces the empirical data with remarkable accuracy. Delaying the vaccination severely affects the mortality, hospitalization, and recovery projections. In Chicago, the disease spread was under control, reducing the mortality increment as the start of the vaccination was postponed. In NYC, the number of cases was increasing, thus, the estimated model predicted a much larger impact, despite the implementation of contention measures. The earlier the vaccination campaign begins, the larger is its potential impact in reducing the COVID-19 cases, as well as in the hospitalizations and deaths. Moreover, the rate at which cases, hospitalizations and deaths increase with the delay in the vaccination beginning strongly depends on the shape of the incidence of infection in each city.

preprint2019arXiv

A policy iteration algorithm for nonzero-sum stochastic impulse games

This work presents a novel policy iteration algorithm to tackle nonzero-sum stochastic impulse games arising naturally in many applications. Despite the obvious impact of solving such problems, there are no suitable numerical methods available, to the best of our knowledge. Our method relies on the recently introduced characterization of the value functions and Nash equilibrium via a system of quasi-variational inequalities. While our algorithm is heuristic and we do not provide a convergence analysis, numerical tests show that it performs convincingly in a wide range of situations, including the only analytically solvable example available in the literature at the time of writing.

preprint2016arXiv

Local Volatility Models in Commodity Markets and Online Calibration

We introduce a local volatility model for the valuation of options on commodity futures by using European vanilla option prices. The corresponding calibration problem is addressed within an online framework, allowing the use of multiple price surfaces. Since uncertainty in the observation of the underlying future prices translates to uncertainty in data locations, we propose a model-based adjustment of such prices that improves reconstructions and smile adherence. In order to tackle the ill-posedness of the calibration problem we incorporate a priori information through a judiciously designed Tikhonov-type regularization. Extensive empirical tests with market as well as synthetic data are used to demonstrate the effectiveness of the methodology and algorithms.

preprint2015arXiv

A Hedged Monte Carlo Approach to Real Option Pricing

In this work we are concerned with valuing optionalities associated to invest or to delay investment in a project when the available information provided to the manager comes from simulated data of cash flows under historical (or subjective) measure in a possibly incomplete market. Our approach is suitable also to incorporating subjective views from management or market experts and to stochastic investment costs. It is based on the Hedged Monte Carlo strategy proposed by Potters et al (2001) where options are priced simultaneously with the determination of the corresponding hedging. The approach is particularly well-suited to the evaluation of commodity related projects whereby the availability of pricing formulae is very rare, the scenario simulations are usually available only in the historical measure, and the cash flows can be highly nonlinear functions of the prices.

preprint2014arXiv

On the Choice of the Tikhonov Regularization Parameter and the Discretization Level: A Discrepancy-Based Strategy

We address the classical issue of appropriate choice of the regularization and discretization level for the Tikhonov regularization of an inverse problem with imperfectly measured data. We focus on the fact that the proper choice of the discretization level in the domain together with the regularization parameter is a key feature in adequate regularization. We propose a discrepancy-based choice for these quantities by applying a relaxed version of Morozov's discrepancy principle. Indeed, we prove the existence of the discretization level and the regularization parameter satisfying such discrepancy. We also prove associated regularizing properties concerning the Tikhonov minimizers.

preprint2010arXiv

A Structured Population Model of Cell Differentiation

We introduce and analyze several aspects of a new model for cell differentiation. It assumes that differentiation of progenitor cells is a continuous process. From the mathematical point of view, it is based on partial differential equations of transport type. Specifically, it consists of a structured population equation with a nonlinear feedback loop. This models the signaling process due to cytokines, which regulate the differentiation and proliferation process. We compare the continuous model to its discrete counterpart, a multi-compartmental model of a discrete collection of cell subpopulations recently proposed by Marciniak-Czochra et al. in 2009 to investigate the dynamics of the hematopoietic system. We obtain uniform bounds for the solutions, characterize steady state solutions, and analyze their linearized stability. We show how persistence or extinction might occur according to values of parameters that characterize the stem cells self-renewal. We also perform numerical simulations and discuss the qualitative behavior of the continuous model vis a vis the discrete one.

preprint2009arXiv

Global stability for a class of virus models with CTL immune response and antigenic variation

We study the global stability of a class of models for in-vivo virus dynamics, that take into account the CTL immune response and display antigenic variation. This class includes a number of models that have been extensively used to model HIV dynamics. We show that models in this class are globally asymptotically stable, under mild hypothesis, by using appropriate Lyapunov functions. We also characterise the stable equilibrium points for the entire biologically relevant parameter range. As a byproduct, we are able to determine what is the diversity of the persistent strains.

preprint2009arXiv

On the Calibration of a Size-Structured Population Model from Experimental Data

The aim of this work is twofold. First, we survey the techniques developed in (Perthame, Zubelli, 2007) and (Doumic, Perthame, Zubelli, 2008) to reconstruct the division (birth) rate from the cell volume distribution data in certain structured population models. Secondly, we implement such techniques on experimental cell volume distributions available in the literature so as to validate the theoretical and numerical results. As a proof of concept, we use the data reported in the classical work of Kubitschek [3] concerning Escherichia coli in vitro experiments measured by means of a Coulter transducer-multichannel analyzer system (Coulter Electronics, Inc., Hialeah, Fla, USA.) Despite the rather old measurement technology, the reconstructed division rates still display potentially useful biological features.

preprint2008arXiv

Numerical Solution of an Inverse Problem in Size-Structured Population Dynamics

We consider a size-structured model for cell division and address the question of determining the division (birth) rate from the measured stable size distribution of the population. We propose a new regularization technique based on a filtering approach. We prove convergence of the algorithm and validate the theoretical results by implementing numerical simulations, based on classical techniques. We compare the results for direct and inverse problems, for the filtering method and for the quasi-reversibility method proposed in [Perthame-Zubelli].