-
Transformations preserving the effective spectral radius of a matrix
Authors:
Jean-François Delmas,
Dylan Dronnier,
Pierre-André Zitt
Abstract:
We discuss transformations on matrices that preserve the effective spectrum and/or the effective spectral radius.
We discuss transformations on matrices that preserve the effective spectrum and/or the effective spectral radius.
△ Less
Submitted 27 February, 2024;
originally announced June 2024.
-
Atoms and associated spectral properties for positive operators on L^p
Authors:
Jean-François Delmas,
Kacem Lefki,
Pierre-André Zitt
Abstract:
Inspired by Schwartz, Jang-Lewis and Victory, who study in particular generalizations of triangularizations of matrices to operators, we shall give for positive operators on Lebesgue spaces equivalent definitions of atoms (maximal irreducible sets). We also characterize positive power compact operators having a unique non-zero atom which appears as a natural generalization of irreducible operators…
▽ More
Inspired by Schwartz, Jang-Lewis and Victory, who study in particular generalizations of triangularizations of matrices to operators, we shall give for positive operators on Lebesgue spaces equivalent definitions of atoms (maximal irreducible sets). We also characterize positive power compact operators having a unique non-zero atom which appears as a natural generalization of irreducible operators and are also considered in epidemiological models. Using the different characterizations of atoms, we also provide a short proof for the representation of the ascent of a positive power compact operator as the maximal length in the graph of critical atoms.
△ Less
Submitted 24 October, 2023;
originally announced October 2023.
-
Individual based SIS models on (not so) dense large random networks
Authors:
Jean-François Delmas,
Paolo Frasca,
Federica Garin,
Viet Chi Tran,
Aurélien Velleret,
Pierre-André Zitt
Abstract:
Starting from a stochastic individual-based description of an SIS epidemic spreading on a random network, we study the dynamics when the size n of the network tends to infinity. We recover in the limit an infinite-dimensional integro-differential equation studied by Delmas, Dronnier and Zitt (2022) for an SIS epidemic propagating on a graphon. Our work covers the case of dense and sparse graphs, p…
▽ More
Starting from a stochastic individual-based description of an SIS epidemic spreading on a random network, we study the dynamics when the size n of the network tends to infinity. We recover in the limit an infinite-dimensional integro-differential equation studied by Delmas, Dronnier and Zitt (2022) for an SIS epidemic propagating on a graphon. Our work covers the case of dense and sparse graphs, provided that the number of edges grows faster than n, but not the case of very sparse graphs with O(n) edges. In order to establish our limit theorem, we have to deal with both the convergence of the random graphs to the graphon and the convergence of the stochastic process spreading on top of these random structures: in particular, we propose a coupling between the process of interest and an epidemic that spreads on the complete graph but with a modified infection rate.
Keywords: random graph, mathematical models of epidemics, measure-valued process, large network limit, limit theorem, graphon.
△ Less
Submitted 31 July, 2023; v1 submitted 26 February, 2023;
originally announced February 2023.
-
Vaccinating according to the maximal endemic equilibrium achieves herd immunity
Authors:
Jean-François Delmas,
Dylan Dronnier,
Pierre-André Zitt
Abstract:
We consider the simple epidemiological SIS model for a general heterogeneous population introduced by Lajmanovich and Yorke (1976) in finite dimension, and its infinite dimensional generalization we introduced in previous works. In this model the basic reproducing number $R_0$ is given by the spectral radius of an integral operator. If $R_0>1$, then there exists a maximal endemic equilibrium. In t…
▽ More
We consider the simple epidemiological SIS model for a general heterogeneous population introduced by Lajmanovich and Yorke (1976) in finite dimension, and its infinite dimensional generalization we introduced in previous works. In this model the basic reproducing number $R_0$ is given by the spectral radius of an integral operator. If $R_0>1$, then there exists a maximal endemic equilibrium. In this very general heterogeneous SIS model, we prove that vaccinating according to the profile of this maximal endemic equilibrium ensures herd immunity. Moreover, this vaccination strategy is critical: the resulting effective reproduction number is exactly equal to one. As an application, we estimate that if $R_0 = 2$ in an age-structured community with mixing rates fitted to social activity, applying this strategy would require approximately 29% less vaccine doses than the strategy which consists in vaccinating uniformly a proportion $1 - 1/R_0$ of the population.
From a dynamical systems point of view, we prove that the non-maximality of an equilibrium $g$ is equivalent to its linear instability in the original dynamics, and to the linear instability of the disease-free state in the modified dynamics where we vaccinate according to $g$.
△ Less
Submitted 20 November, 2022;
originally announced November 2022.
-
The effective reproduction number: convexity, concavity and invariance
Authors:
Jean-François Delmas,
Dylan Dronnier,
Pierre-André Zitt
Abstract:
Motivated by the question of optimal vaccine allocation strategies in heterogeneous population for epidemic models, we study various properties of the \emph{effective reproduction number}. In the simplest case, given a fixed, non-negative matrix $K$, this corresponds mathematically to the study of the spectral radius $R_e(η)$ of the matrix product $\mathrm{Diag}(η)K$, as a function of…
▽ More
Motivated by the question of optimal vaccine allocation strategies in heterogeneous population for epidemic models, we study various properties of the \emph{effective reproduction number}. In the simplest case, given a fixed, non-negative matrix $K$, this corresponds mathematically to the study of the spectral radius $R_e(η)$ of the matrix product $\mathrm{Diag}(η)K$, as a function of $η\in\mathbb{R}_+^n$. The matrix $K$ and the vector $η$ can be interpreted as a next-generation operator and a vaccination strategy. This can be generalized in an infinite dimensional case where the matrix $K$ is replaced by a positive integral compact operator, which is composed with a multiplication by a non-negative function $η$.
We give sufficient conditions for the function $R_e$ to be convex or a concave. Eventually, we provide equivalence properties on models which ensure that the function $R_e$ is unchanged.
△ Less
Submitted 21 November, 2022;
originally announced November 2022.
-
Optimal vaccinations: Cordons sanitaires, reducible population and optimal rays
Authors:
Jean-François Delmas,
Dylan Dronnier,
Pierre-André Zitt
Abstract:
We consider the bi-objective problem of allocating doses of a (perfect) vaccine to an infinite-dimensional metapopulation in order to minimize simultaneously the vaccination cost and the effective reproduction number $R_e$, which is defined as the spectral radius of the effective next-generation operator.
In this general framework, we prove that a cordon sanitaire, that is, a strategy that effec…
▽ More
We consider the bi-objective problem of allocating doses of a (perfect) vaccine to an infinite-dimensional metapopulation in order to minimize simultaneously the vaccination cost and the effective reproduction number $R_e$, which is defined as the spectral radius of the effective next-generation operator.
In this general framework, we prove that a cordon sanitaire, that is, a strategy that effectively disconnects the non-vaccinated population, might not be optimal, but it is still better than the "worst" vaccination strategies. Inspired by graph theory, we also compute the minimal cost which ensures that no infection occurs using independent sets. Using Frobenius decomposition of the whole population into irreducible sub-populations, we give some explicit formulae for optimal ("best" and "worst") vaccinations strategies. Eventually, we provide some sufficient conditions for a scaling of an optimal strategy to still be optimal.
△ Less
Submitted 15 December, 2022; v1 submitted 15 September, 2022;
originally announced September 2022.
-
Optimal vaccination: Various (counter) intuitive examples
Authors:
Jean-François Delmas,
Dylan Dronnier,
Pierre-André Zitt
Abstract:
In previous articles, we formalized the problem of optimal allocation strategies for a (perfect) vaccine in an infinite-dimensional metapopulation model. The aim of the current paper is to illustrate this theoretical framework with multiple examples where one can derive the analytic expression of the optimal strategies. We discuss in particular the following points: whether or not it is possible t…
▽ More
In previous articles, we formalized the problem of optimal allocation strategies for a (perfect) vaccine in an infinite-dimensional metapopulation model. The aim of the current paper is to illustrate this theoretical framework with multiple examples where one can derive the analytic expression of the optimal strategies. We discuss in particular the following points: whether or not it is possible to vaccinate optimally when the vaccine doses are given one at a time (greedy vaccination strategies); the effect of assortativity (that is, the tendency to have more contacts with similar individuals) on the shape of optimal vaccination strategies; the particular case where everybody has the same number of neighbors.
△ Less
Submitted 16 December, 2021;
originally announced December 2021.
-
Asymptotics of the Inertia Moments and the Variance Conjecture in Schatten Balls
Authors:
Benjamin Dadoun,
Matthieu Fradelizi,
Olivier Guédon,
Pierre-André Zitt
Abstract:
We study the first and second orders of the asymptotic expansion, as the dimension goes to infinity, of the moments of the Hilbert-Schmidt norm of a uniformly distributed matrix in the p-Schatten unit ball. We consider the case of matrices with real, complex or quaternionic entries, self-adjoint or not. When p > 3, this asymptotic expansion allows us to establish a generalized version of the varia…
▽ More
We study the first and second orders of the asymptotic expansion, as the dimension goes to infinity, of the moments of the Hilbert-Schmidt norm of a uniformly distributed matrix in the p-Schatten unit ball. We consider the case of matrices with real, complex or quaternionic entries, self-adjoint or not. When p > 3, this asymptotic expansion allows us to establish a generalized version of the variance conjecture for the family of p-Schatten unit balls of self-adjoint matrices.
△ Less
Submitted 16 February, 2022; v1 submitted 15 November, 2021;
originally announced November 2021.
-
Effective reproduction number: Convexity, invariance and cordons sanitaires
Authors:
Jean-François Delmas,
Dylan Dronnier,
Pierre-André Zitt
Abstract:
We consider the problem of optimal allocation strategies for a (perfect) vaccine in an infinite-metapopulation model (including SIS, SIR, SEIR,. . .), when the loss function is given by the effective reproduction number Re, which is defined as the spectral radius of the effective next generation matrix (in finite dimension) or more generally of the effective next generation operator (in infinite d…
▽ More
We consider the problem of optimal allocation strategies for a (perfect) vaccine in an infinite-metapopulation model (including SIS, SIR, SEIR,. . .), when the loss function is given by the effective reproduction number Re, which is defined as the spectral radius of the effective next generation matrix (in finite dimension) or more generally of the effective next generation operator (in infinite dimension). We give sufficient conditions for Re to be a convex or a concave function of the vaccination strategy. Then, following a previous work, we consider the bi-objective problem of minimizing simultaneously the cost and the loss of the vaccination strategies. In particular, we prove that a cordon sanitaire might not be optimal, but it is still better than the "worst" vaccination strategies. Inspired by the graph theory, we compute the minimal cost which ensures that no infection occurs using independent sets. Using Frobenius decomposition of the whole population into irreducible sub-populations, we give some explicit formulae for optimal ("best" and "worst") vaccinations strategies. Eventually, we provide equivalence properties on models which ensure that the function Re is unchanged.
△ Less
Submitted 25 October, 2021;
originally announced October 2021.
-
Targeted Vaccination Strategies for an Infinite-dimensional SIS Model
Authors:
Jean-François Delmas,
Dylan Dronnier,
Pierre-André Zitt
Abstract:
We formalize and study the problem of optimal allocation strategies for a (perfect) vaccine in the infinite-dimensional SIS model. The question may be viewed as a bi-objective minimization problem, where one tries to minimize simultaneously the cost of the vaccination, and a loss that may be either the effective reproduction number, or the proportion of the infected population in the endemic state…
▽ More
We formalize and study the problem of optimal allocation strategies for a (perfect) vaccine in the infinite-dimensional SIS model. The question may be viewed as a bi-objective minimization problem, where one tries to minimize simultaneously the cost of the vaccination, and a loss that may be either the effective reproduction number, or the proportion of the infected population in the endemic state. We prove the existence of Pareto optimal strategies, describe the corresponding Pareto frontier in both cases, and study its convexity and stability properties. We also show that vaccinating according to the profile of the endemic state is a critical allocation, in the sense that, if the initial reproduction number is larger than 1, then this vaccination strategy yields an effective reproduction number equal to 1.
△ Less
Submitted 26 August, 2021; v1 submitted 17 March, 2021;
originally announced March 2021.
-
Exact targeting of Gibbs distributions using velocity-jump processes
Authors:
Pierre Monmarché,
Mathias Rousset,
Pierre-André Zitt
Abstract:
This work introduces and studies a new family of velocity jump Markov processes directly amenable to exact simulation with the following two properties: i) trajectories converge in law when a time-step parameter vanishes towards a given Langevin or Hamil-tonian dynamics; ii) the stationary distribution of the process is always exactly given by the product of a Gaussian (for velocities) by any targ…
▽ More
This work introduces and studies a new family of velocity jump Markov processes directly amenable to exact simulation with the following two properties: i) trajectories converge in law when a time-step parameter vanishes towards a given Langevin or Hamil-tonian dynamics; ii) the stationary distribution of the process is always exactly given by the product of a Gaussian (for velocities) by any target log-density whose gradient is pointwise computabe together with some additional explicit appropriate upper bound. The process does not exhibit any velocity reflections (jump sizes can be controlled) and is suitable for the 'factorization method'. We provide a rigorous mathematical proof of: i) the small time-step convergence towards Hamiltonian/Langevin dynamics, as well as ii) the exponentially fast convergence towards the target distribution when suitable noise on velocity is present. Numerical implementation is detailed and illustrated.
△ Less
Submitted 14 September, 2020; v1 submitted 21 August, 2020;
originally announced August 2020.
-
An Infinite-Dimensional SIS Model
Authors:
Jean-François Delmas,
Dylan Dronnier,
Pierre-André Zitt
Abstract:
In this article, we introduce an infinite-dimensional deterministic SIS model which takes into account the heterogeneity of the infections and the social network among a large population. We study the long-time behavior of the dynamic. We identify the basic reproduction number $R_0$ which determines whether there exists a stable endemic steady state (super-critical case: $R_0>1$) or if the only e…
▽ More
In this article, we introduce an infinite-dimensional deterministic SIS model which takes into account the heterogeneity of the infections and the social network among a large population. We study the long-time behavior of the dynamic. We identify the basic reproduction number $R_0$ which determines whether there exists a stable endemic steady state (super-critical case: $R_0>1$) or if the only equilibrium is disease-free (critical and sub-critical case: $R_0\leq1$). As an application of this general study, we prove that the so-called ``leaky'' and ``all-or-nothing'' vaccination mechanism have the same effect on $R_0$. This framework is also very natural and intuitive to model lockdown policies and study their impact.
△ Less
Submitted 15 June, 2020;
originally announced June 2020.
-
Convergence of metadynamics: discussion of the adiabatic hypothesis
Authors:
Benjamin Jourdain,
Tony Lelièvre,
Pierre-André Zitt
Abstract:
By drawing a parallel between metadynamics and self interacting models for polymers, we study the longtime convergence of the original metadynamics algorithm in the adiabatic setting, namely when the dynamics along the collective variables decouples from the dynamics along the other degrees of freedom. We also discuss the bias which is introduced when the adiabatic assumption does not holds.
By drawing a parallel between metadynamics and self interacting models for polymers, we study the longtime convergence of the original metadynamics algorithm in the adiabatic setting, namely when the dynamics along the collective variables decouples from the dynamics along the other degrees of freedom. We also discuss the bias which is introduced when the adiabatic assumption does not holds.
△ Less
Submitted 18 April, 2019;
originally announced April 2019.
-
Ergodicity of the zigzag process
Authors:
Joris Bierkens,
Gareth Roberts,
Pierre-André Zitt
Abstract:
The zigzag process is a Piecewise Deterministic Markov Process which can be used in a MCMC framework to sample from a given target distribution. We prove the convergence of this process to its target under very weak assumptions, and establish a central limit theorem for empirical averages under stronger assumptions on the decay of the target measure. We use the classical "Meyn-Tweedie" approach. T…
▽ More
The zigzag process is a Piecewise Deterministic Markov Process which can be used in a MCMC framework to sample from a given target distribution. We prove the convergence of this process to its target under very weak assumptions, and establish a central limit theorem for empirical averages under stronger assumptions on the decay of the target measure. We use the classical "Meyn-Tweedie" approach. The main difficulty turns out to be the proof that the process can indeed reach all the points in the space, even if we consider the minimal switching rates.
△ Less
Submitted 18 January, 2019; v1 submitted 28 December, 2017;
originally announced December 2017.
-
A Weak Overdamped Limit Theorem for Langevin Processes
Authors:
Mathias Rousset,
Yushun Xu,
Pierre-André Zitt
Abstract:
In this paper, we prove convergence in distribution of Langevin processes in the overdamped asymptotics. The proof relies on the classical perturbed test function (or corrector) method, which is used both to show tightness in path space, and to identify the extracted limit with a martingale problem. The result holds assuming the continuity of the gradient of the potential energy, and a mild contro…
▽ More
In this paper, we prove convergence in distribution of Langevin processes in the overdamped asymptotics. The proof relies on the classical perturbed test function (or corrector) method, which is used both to show tightness in path space, and to identify the extracted limit with a martingale problem. The result holds assuming the continuity of the gradient of the potential energy, and a mild control of the initial kinetic energy.
△ Less
Submitted 8 March, 2019; v1 submitted 28 September, 2017;
originally announced September 2017.
-
On the persistence regime for Lotka-Volterra in randomly fluctuating environments
Authors:
Florent Malrieu,
Pierre-André Zitt
Abstract:
In this note, we study the long time behavior of Lotka-Volterra systems whose coefficients vary randomly. Bena{ï}m and Lobry (2015) recently established that randomly switching between two environments that are both favorable to the same species may lead to different regimes: extinction of one species or the other, or persistence of both species. Our purpose here is to describe more accurately the…
▽ More
In this note, we study the long time behavior of Lotka-Volterra systems whose coefficients vary randomly. Bena{ï}m and Lobry (2015) recently established that randomly switching between two environments that are both favorable to the same species may lead to different regimes: extinction of one species or the other, or persistence of both species. Our purpose here is to describe more accurately the range of parameters leading to these regimes, and the support of the invariant probability measure in case of persistence.
△ Less
Submitted 14 June, 2017; v1 submitted 29 January, 2016;
originally announced January 2016.
-
Functional inequalities for Gaussian convolutions of compactly supported measures: explicit bounds and dimension dependence
Authors:
Jean-Baptiste Bardet,
Nathaël Gozlan,
Florent Malrieu,
Pierre-André Zitt
Abstract:
The aim of this paper is to establish various functional inequalities for the convolution of a compactly supported measure and a standard Gaussian distribution on Rd. We especially focus on getting good dependence of the constants on the dimension. We prove that the Poincar{é} inequality holds with a dimension-free bound. For the logarithmic Sobolev inequality, we improve the best known results (Z…
▽ More
The aim of this paper is to establish various functional inequalities for the convolution of a compactly supported measure and a standard Gaussian distribution on Rd. We especially focus on getting good dependence of the constants on the dimension. We prove that the Poincar{é} inequality holds with a dimension-free bound. For the logarithmic Sobolev inequality, we improve the best known results (Zimmermann, JFA 2013) by getting a bound that grows linearly with the dimension. We also establish transport-entropy inequalities for various transport costs.
△ Less
Submitted 9 July, 2015;
originally announced July 2015.
-
A Wright-Fisher model with indirect selection
Authors:
Ludovic Goudenège,
Pierre-André Zitt
Abstract:
We study a generalization of the Wright--Fisher model in which some individuals adopt a behavior that is harmful to others without any direct advantage for themselves. This model is motivated by studies of spiteful behavior in nature, including several species of parasitoid hymenoptera in which sperm-depleted males continue to mate de- spite not being fertile.
We first study a single reproductive…
▽ More
We study a generalization of the Wright--Fisher model in which some individuals adopt a behavior that is harmful to others without any direct advantage for themselves. This model is motivated by studies of spiteful behavior in nature, including several species of parasitoid hymenoptera in which sperm-depleted males continue to mate de- spite not being fertile.
We first study a single reproductive season, then use it as a building block for a generalized Wright--Fisher model. In the large population limit, for male-skewed sex ratios, we rigorously derive the convergence of the renormalized process to a diffusion with a frequency-dependent selection and genetic drift. This allows a quantitative comparison of the indirect selective advantage with the direct one classically considered in the Wright--Fisher model.
From the mathematical point of view, each season is modeled by a mix between samplings with and without replacement, and analyzed by a sort of "reverse numerical analysis", viewing a key recurrence relation as a discretization scheme for a PDE. The diffusion approximation is then obtained by classical methods.
△ Less
Submitted 17 March, 2015; v1 submitted 16 January, 2014;
originally announced January 2014.
-
Piecewise deterministic Markov process - recent results
Authors:
Romain Azaïs,
Jean-Baptiste Bardet,
Alexandre Genadot,
Nathalie Krell,
Pierre-André Zitt
Abstract:
We give a short overview of recent results on a specific class of Markov process: the Piecewise Deterministic Markov Processes (PDMPs). We first recall the definition of these processes and give some general results. On more specific cases such as the TCP model or a model of switched vector fields, better results can be proved, especially as regards long time behaviour. We continue our review with…
▽ More
We give a short overview of recent results on a specific class of Markov process: the Piecewise Deterministic Markov Processes (PDMPs). We first recall the definition of these processes and give some general results. On more specific cases such as the TCP model or a model of switched vector fields, better results can be proved, especially as regards long time behaviour. We continue our review with an infinite dimensional example of neuronal activity. From the statistical point of view, these models provide specific challenges: we illustrate this point with the example of the estimation of the distribution of the inter-jumping times. We conclude with a short overview on numerical methods used for simulating PDMPs.
△ Less
Submitted 24 September, 2013;
originally announced September 2013.
-
First order global asymptotics for confined particles with singular pair repulsion
Authors:
Djalil Chafaï,
Nathael Gozlan,
Pierre-André Zitt
Abstract:
We study a physical system of $N$ interacting particles in $\mathbb{R}^d$, $d\geq1$, subject to pair repulsion and confined by an external field. We establish a large deviations principle for their empirical distribution as $N$ tends to infinity. In the case of Riesz interaction, including Coulomb interaction in arbitrary dimension $d>2$, the rate function is strictly convex and admits a unique mi…
▽ More
We study a physical system of $N$ interacting particles in $\mathbb{R}^d$, $d\geq1$, subject to pair repulsion and confined by an external field. We establish a large deviations principle for their empirical distribution as $N$ tends to infinity. In the case of Riesz interaction, including Coulomb interaction in arbitrary dimension $d>2$, the rate function is strictly convex and admits a unique minimum, the equilibrium measure, characterized via its potential. It follows that almost surely, the empirical distribution of the particles tends to this equilibrium measure as $N$ tends to infinity. In the more specific case of Coulomb interaction in dimension $d>2$, and when the external field is a convex or increasing function of the radius, then the equilibrium measure is supported in a ring. With a quadratic external field, the equilibrium measure is uniform on a ball.
△ Less
Submitted 7 September, 2014; v1 submitted 29 April, 2013;
originally announced April 2013.
-
Qualitative properties of certain piecewise deterministic Markov processes
Authors:
Michel Benaïm,
Stéphane Le Borgne,
Florent Malrieu,
Pierre-André Zitt
Abstract:
We study a class of Piecewise Deterministic Markov Processes with state space Rd x E where E is a finite set. The continuous component evolves according to a smooth vector field that is switched at the jump times of the discrete coordinate. The jump rates may depend on the whole position of the process. Working under the general assumption that the process stays in a compact set, we detail a possi…
▽ More
We study a class of Piecewise Deterministic Markov Processes with state space Rd x E where E is a finite set. The continuous component evolves according to a smooth vector field that is switched at the jump times of the discrete coordinate. The jump rates may depend on the whole position of the process. Working under the general assumption that the process stays in a compact set, we detail a possible construction of the process and characterize its support, in terms of the solutions set of a differential inclusion. We establish results on the long time behaviour of the process, in relation to a certain set of accessible points, which is shown to be strongly linked to the support of invariant measures. Under Hörmander-type bracket conditions, we prove that there exists a unique invariant measure and that the processes converges to equilibrium in total variation. Finally we give examples where the bracket condition does not hold, and where there may be one or many invariant measures, depending on the jump rates between the flows.
△ Less
Submitted 7 April, 2014; v1 submitted 18 April, 2012;
originally announced April 2012.
-
Recursive estimation of the conditional geometric median in Hilbert spaces
Authors:
Hervé Cardot,
Peggy Cénac,
Pierre-André Zitt
Abstract:
A recursive estimator of the conditional geometric median in Hilbert spaces is studied. It is based on a stochastic gradient algorithm whose aim is to minimize a weighted L1 criterion and is consequently well adapted for robust online estimation. The weights are controlled by a kernel function and an associated bandwidth. Almost sure convergence and L2 rates of convergence are proved under general…
▽ More
A recursive estimator of the conditional geometric median in Hilbert spaces is studied. It is based on a stochastic gradient algorithm whose aim is to minimize a weighted L1 criterion and is consequently well adapted for robust online estimation. The weights are controlled by a kernel function and an associated bandwidth. Almost sure convergence and L2 rates of convergence are proved under general conditions on the conditional distribution as well as the sequence of descent steps of the algorithm and the sequence of bandwidths. Asymptotic normality is also proved for the averaged version of the algorithm with an optimal rate of convergence. A simulation study confirms the interest of this new and fast algorithm when the sample sizes are large. Finally, the ability of these recursive algorithms to deal with very high-dimensional data is illustrated on the robust estimation of television audience profiles conditional on the total time spent watching television over a period of 24 hours.
△ Less
Submitted 14 April, 2012;
originally announced April 2012.
-
Quantitative ergodicity for some switched dynamical systems
Authors:
Michel Benaïm,
Stéphane Le Borgne,
Florent Malrieu,
Pierre-André Zitt
Abstract:
We provide quantitative bounds for the long time behavior of a class of Piecewise Deterministic Markov Processes with state space Rd \times E where E is a finite set. The continuous component evolves according to a smooth vector field that switches at the jump times of the discrete coordinate. The jump rates may depend on the whole position of the process. Under regularity assumptions on the jump…
▽ More
We provide quantitative bounds for the long time behavior of a class of Piecewise Deterministic Markov Processes with state space Rd \times E where E is a finite set. The continuous component evolves according to a smooth vector field that switches at the jump times of the discrete coordinate. The jump rates may depend on the whole position of the process. Under regularity assumptions on the jump rates and stability conditions for the vector fields we provide explicit exponential upper bounds for the convergence to equilibrium in terms of Wasserstein distances. As an example, we obtain convergence results for a stochastic version of the Morris-Lecar model of neurobiology.
△ Less
Submitted 6 December, 2012; v1 submitted 9 April, 2012;
originally announced April 2012.
-
On the stability of planar randomly switched systems
Authors:
Michel Benaïm,
Stéphane Le Borgne,
Florent Malrieu,
Pierre-André Zitt
Abstract:
Consider the random process (Xt) solution of dXt/dt = A(It) Xt where (It) is a Markov process on {0,1} and A0 and A1 are real Hurwitz matrices on R2. Assuming that there exists lambda in (0, 1) such that (1 - λ)A0 + \lambdaA1 has a positive eigenvalue, we establish that the norm of Xt may converge to 0 or infinity, depending on the the jump rate of the process I. An application to product of rando…
▽ More
Consider the random process (Xt) solution of dXt/dt = A(It) Xt where (It) is a Markov process on {0,1} and A0 and A1 are real Hurwitz matrices on R2. Assuming that there exists lambda in (0, 1) such that (1 - λ)A0 + \lambdaA1 has a positive eigenvalue, we establish that the norm of Xt may converge to 0 or infinity, depending on the the jump rate of the process I. An application to product of random matrices is studied. This paper can be viewed as a probabilistic counterpart of the paper "A note on stability conditions for planar switched systems" by Balde, Boscain and Mason.
△ Less
Submitted 9 April, 2012;
originally announced April 2012.
-
Total variation estimates for the TCP process
Authors:
Jean-Baptiste Bardet,
Alejandra Christen,
Arnaud Guillin,
Florent Malrieu,
Pierre-André Zitt
Abstract:
The TCP window size process appears in the modeling of the famous Transmission Control Protocol used for data transmission over the Internet. This continuous time Markov process takes its values in [0, \infty), is ergodic and irreversible. The sample paths are piecewise linear deterministic and the whole randomness of the dynamics comes from the jump mechanism. The aim of the present paper is to p…
▽ More
The TCP window size process appears in the modeling of the famous Transmission Control Protocol used for data transmission over the Internet. This continuous time Markov process takes its values in [0, \infty), is ergodic and irreversible. The sample paths are piecewise linear deterministic and the whole randomness of the dynamics comes from the jump mechanism. The aim of the present paper is to provide quantitative estimates for the exponential convergence to equilibrium, in terms of the total variation and Wasserstein distances.
△ Less
Submitted 21 September, 2012; v1 submitted 29 December, 2011;
originally announced December 2011.
-
Efficient and fast estimation of the geometric median in Hilbert spaces with an averaged stochastic gradient algorithm
Authors:
Hervé Cardot,
Peggy Cénac,
Pierre-André Zitt
Abstract:
With the progress of measurement apparatus and the development of automatic sensors it is not unusual anymore to get thousands of samples of observations taking values in high dimension spaces such as functional spaces. In such large samples of high dimensional data, outlying curves may not be uncommon and even a few individuals may corrupt simple statistical indicators such as the mean trajectory…
▽ More
With the progress of measurement apparatus and the development of automatic sensors it is not unusual anymore to get thousands of samples of observations taking values in high dimension spaces such as functional spaces. In such large samples of high dimensional data, outlying curves may not be uncommon and even a few individuals may corrupt simple statistical indicators such as the mean trajectory. We focus here on the estimation of the geometric median which is a direct generalization of the real median and has nice robustness properties. The geometric median being defined as the minimizer of a simple convex functional that is differentiable everywhere when the distribution has no atoms, it is possible to estimate it with online gradient algorithms. Such algorithms are very fast and can deal with large samples. Furthermore they also can be simply updated when the data arrive sequentially. We state the almost sure consistency and the L2 rates of convergence of the stochastic gradient estimator as well as the asymptotic normality of its averaged version. We get that the asymptotic distribution of the averaged version of the algorithm is the same as the classic estimators which are based on the minimization of the empirical loss function. The performances of our averaged sequential estimator, both in terms of computation speed and accuracy of the estimations, are evaluated with a small simulation study. Our approach is also illustrated on a sample of more 5000 individual television audiences measured every second over a period of 24 hours.
△ Less
Submitted 20 May, 2011; v1 submitted 22 January, 2011;
originally announced January 2011.
-
Poincaré inequalities and hitting times
Authors:
Patrick Cattiaux,
Arnaud Guillin,
Pierre-André Zitt
Abstract:
Equivalence of the spectral gap, exponential integrability of hitting times and Lyapunov conditions are well known. We give here the correspondance (with quantitative results) for reversible diffusion processes. As a consequence, we generalize results of Bobkov in the one dimensional case on the value of the Poincaré constant for logconcave measures to superlinear potentials. Finally, we study var…
▽ More
Equivalence of the spectral gap, exponential integrability of hitting times and Lyapunov conditions are well known. We give here the correspondance (with quantitative results) for reversible diffusion processes. As a consequence, we generalize results of Bobkov in the one dimensional case on the value of the Poincaré constant for logconcave measures to superlinear potentials. Finally, we study various functional inequalities under different hitting times integrability conditions (polynomial, ...). In particular, in the one dimensional case, ultracontractivity is equivalent to a bounded Lyapunov condition.
△ Less
Submitted 23 December, 2010;
originally announced December 2010.
-
Super Poincaré inequalities, Orlicz norms and essential spectrum
Authors:
Pierre-André Zitt
Abstract:
We prove some results about the super Poincaré inequality (SPI) and its relation to the spectrum of an operator: we show that it can be alternatively written with Orlicz norms instead of L1 norms, and we use this to give an alternative proof that a bound on the bottom of the essential spectrum implies a SPI. Finally, we apply these ideas to give a spectral proof of the log Sobolev inequality for…
▽ More
We prove some results about the super Poincaré inequality (SPI) and its relation to the spectrum of an operator: we show that it can be alternatively written with Orlicz norms instead of L1 norms, and we use this to give an alternative proof that a bound on the bottom of the essential spectrum implies a SPI. Finally, we apply these ideas to give a spectral proof of the log Sobolev inequality for the Gaussian measure.
△ Less
Submitted 25 October, 2009;
originally announced October 2009.
-
Resonances for a diffusion with small noise
Authors:
Markus Klein,
Pierre-André Zitt
Abstract:
We study resonances for the generator of a diffusion with small noise in $R^d$ :$ L_ε= -εΔ+ \nabla F \cdot \nabla$, when the potential F grows slowly at infinity (typically as a square root of the norm). The case when F grows fast is well known, and under suitable conditions one can show that there exists a family of exponentially small eigenvalues, related to the wells of F . We show that, for…
▽ More
We study resonances for the generator of a diffusion with small noise in $R^d$ :$ L_ε= -εΔ+ \nabla F \cdot \nabla$, when the potential F grows slowly at infinity (typically as a square root of the norm). The case when F grows fast is well known, and under suitable conditions one can show that there exists a family of exponentially small eigenvalues, related to the wells of F . We show that, for an F with a slow growth, the spectrum is R+, but we can find a family of resonances whose real parts behave as the eigenvalues of the "quick growth" case, and whose imaginary parts are small.
△ Less
Submitted 1 May, 2008;
originally announced May 2008.
-
Functional inequalities and uniqueness of the Gibbs measure -- from log-Sobolev to Poincaré
Authors:
Pierre-André Zitt
Abstract:
In a statistical mechanics model with unbounded spins, we prove uniqueness of the Gibbs measure under various assumptions on finite volume functional inequalities. We follow the approach of G. Royer (1999) and obtain uniqueness by showing convergence properties of a Glauber-Langevin dynamics. The result was known when the measures on the box $[-n,n]^d$ (with free boundary conditions) satisfied t…
▽ More
In a statistical mechanics model with unbounded spins, we prove uniqueness of the Gibbs measure under various assumptions on finite volume functional inequalities. We follow the approach of G. Royer (1999) and obtain uniqueness by showing convergence properties of a Glauber-Langevin dynamics. The result was known when the measures on the box $[-n,n]^d$ (with free boundary conditions) satisfied the same logarithmic Sobolev inequality. We generalize this in two directions: either the constants may be allowed to grow sub-linearly in the diameter, or we may suppose a weaker inequality than log-Sobolev, but stronger than Poincaré. We conclude by giving a heuristic argument showing that this could be the right inequalities to look at.
△ Less
Submitted 14 February, 2007;
originally announced February 2007.
-
Annealing diffusions in a slowly growing potential
Authors:
Pierre-André Zitt
Abstract:
We consider a continuous analogue of the simulated annealing algorithm in $R^d$. We prove a convergence result, under hypotheses weaker than the usual ones. In particular, we cover cases where the gradient of the potential goes to zero at infinity. The proof follows an idea of L. Miclo, but we replace the Poincaré and log-Sobolev inequalities (which do not hold in our setting) by weak Poincaré i…
▽ More
We consider a continuous analogue of the simulated annealing algorithm in $R^d$. We prove a convergence result, under hypotheses weaker than the usual ones. In particular, we cover cases where the gradient of the potential goes to zero at infinity. The proof follows an idea of L. Miclo, but we replace the Poincaré and log-Sobolev inequalities (which do not hold in our setting) by weak Poincaré inequalities. We estimate the latter with measure-capacity criteria. We show that, despite the absence of a spectral gap, the convergence still holds for the "classical" schedule t = c/ ln(t), if c is bigger than a constant related to the potential.
△ Less
Submitted 6 July, 2006;
originally announced July 2006.