Source author record

Fabien Panloup

Fabien Panloup 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

14works
7topics
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

14 published item(s)

preprint2022arXiv

On the cost of Bayesian posterior mean strategy for log-concave models

In this paper, we investigate the problem of computing Bayesian estimators using Langevin Monte-Carlo type approximation. The novelty of this paper is to consider together the statistical and numerical counterparts (in a general log-concave setting). More precisely, we address the following question: given $n$ observations in $\mathbb{R}^q$ distributed under an unknown probability $\mathbb{P}_{θ^\star}$ with $θ^\star \in \mathbb{R}^d$ , what is the optimal numerical strategy and its cost for the approximation of $θ^\star$ with the Bayesian posterior mean? To answer this question, we establish some quantitative statistical bounds related to the underlying Poincaré constant of the model and establish new results about the numerical approximation of Gibbs measures by Cesaro averages of Euler schemes of (over-damped) Langevin diffusions. These last results include in particular some quantitative controls in the weakly convex case based on new bounds on the solution of the related Poisson equation of the diffusion.

preprint2020arXiv

A general drift estimation procedure for stochastic differential equations with additive fractional noise

In this paper we consider the drift estimation problem for a general differential equation driven by an additive multidimensional fractional Brownian motion, under ergodic assumptions on the drift coefficient. Our estimation procedure is based on the identification of the invariant measure, and we provide consistency results as well as some information about the convergence rate. We also give some examples of coefficients for which the identifiability assumption for the invariant measure is satisfied.

preprint2020arXiv

Sub-exponential convergence to equilibrium for Gaussian driven Stochastic Differential Equations with semi-contractive drift

The convergence to the stationary regime is studied for Stochastic Differential Equations driven by an additive Gaussian noise and evolving in a semi-contractive environment, i.e. when the drift is only contractive out of a compact set but does not have repulsive regions. In this setting, we develop a synchronous coupling strategy to obtain sub-exponential bounds on the rate of convergence to equilibrium in Wasserstein distance. Then by a coalescent coupling close to terminal time, we derive a similar bound in total variation distance.

preprint2016arXiv

Rate of convergence to equilibrium of fractional driven stochastic differential equations with some multiplicative noise

We investigate the problem of the rate of convergence to equilibrium for ergodic stochastic differential equations driven by fractional Brownian motion with Hurst parameter $H\textgreater{}1/2$ and multiplicative noise component $σ$. When $σ$ is constant and for every $H\in(0,1)$, it was proved in \cite{hairer} that, under some mean-reverting assumptions, such a process converges to its equilibrium at a rate of order $t^{-α}$ where $α\in(0,1)$ (depending on $H$). The aim of this paper is to extend such types of results to some multiplicative noise setting. More precisely, we show that we can recover such convergence rates when $H\textgreater{}1/2$ and the inverse of the diffusion coefficient $σ$ is a Jacobian matrix. The main novelty of this work is a type of extension of Foster-Lyapunov like techniques to this non-Markovian setting, which allows us to put in place an asymptotic coupling scheme such as in \cite{hairer} without resorting to deterministic contracting properties.

preprint2016arXiv

Regret bounds for Narendra-Shapiro bandit algorithms

Narendra-Shapiro (NS) algorithms are bandit-type algorithms that have been introduced in the sixties (with a view to applications in Psychology or learning automata), whose convergence has been intensively studied in the stochastic algorithm literature. In this paper, we adress the following question: are the Narendra-Shapiro (NS) bandit algorithms competitive from a \textit{regret} point of view? In our main result, we show that some competitive bounds can be obtained for such algorithms in their penalized version (introduced in \cite{Lamberton_Pages}). More precisely, up to an over-penalization modification, the pseudo-regret $\bar{R}_n$ related to the penalized two-armed bandit algorithm is uniformly bounded by $C \sqrt{n}$ (where $C$ is made explicit in the paper). \noindent We also generalize existing convergence and rates of convergence results to the multi-armed case of the over-penalized bandit algorithm, including the convergence toward the invariant measure of a Piecewise Deterministic Markov Process (PDMP) after a suitable renormalization. Finally, ergodic properties of this PDMP are given in the multi-armed case.

preprint2016arXiv

Stochastic Heavy Ball

This paper deals with a natural stochastic optimization procedure derived from the so-called Heavy-ball method differential equation, which was introduced by Polyak in the 1960s with his seminal contribution [Pol64]. The Heavy-ball method is a second-order dynamics that was investigated to minimize convex functions f . The family of second-order methods recently received a large amount of attention, until the famous contribution of Nesterov [Nes83], leading to the explosion of large-scale optimization problems. This work provides an in-depth description of the stochastic heavy-ball method, which is an adaptation of the deterministic one when only unbiased evalutions of the gradient are available and used throughout the iterations of the algorithm. We first describe some almost sure convergence results in the case of general non-convex coercive functions f . We then examine the situation of convex and strongly convex potentials and derive some non-asymptotic results about the stochastic heavy-ball method. We end our study with limit theorems on several rescaled algorithms.

preprint2016arXiv

Weighted Multilevel Langevin Simulation of Invariant Measures

We investigate a weighted Multilevel Richardson-Romberg extrapolation for the ergodic approximation of invariant distributions of diffusions adapted from the one introduced in~[Lemaire-Pagès, 2013] for regular Monte Carlo simulation. In a first result, we prove under weak confluence assumptions on the diffusion, that for any integer $R\ge2$, the procedure allows us to attain a rate $n^{\frac{R}{2R+1}}$ whereas the original algorithm convergence is at a weak rate $n^{1/3}$. Furthermore, this is achieved without any explosion of the asymptotic variance. In a second part, under stronger confluence assumptions and with the help of some second order expansions of the asymptotic error, we go deeper in the study by optimizing the choice of the parameters involved by the method. In particular, for a given $\varepsilon\textgreater{}0$, we exhibit some semi-explicit parameters for which the number of iterations of the Euler scheme required to attain a Mean-Squared Error lower than $\varepsilon^2$ is about $\varepsilon^{-2}\log(\varepsilon^{-1})$. Finally, we numerically this Multilevel Langevin estimator on several examples including the simple one-dimensional Ornstein-Uhlenbeck process but also on a high dimensional diffusion motivated by a statistical problem. These examples confirm the theoretical efficiency of the method.

preprint2015arXiv

Long time behavior of Markov processes and beyond

This note provides several recent progresses in the study of long time behavior of Markov processes. The examples presented below are related to other scientific fields as PDE's, physics or biology. The involved mathematical tools as propagation of chaos, coupling, functional inequalities, provide a good picture of the classical methods that furnish quantitative rates of convergence to equilibrium.

preprint2013arXiv

A stochastic model for speculative bubbles

This paper aims to provide a simple modelling of speculative bubbles and derive some quantitative properties of its dynamical evolution. Starting from a description of individual speculative behaviours, we build and study a second order Markov process, which after simple transformations can be viewed as a turning two-dimensional Gaussian process. Then, our main problem is to ob- tain some bounds for the persistence rate relative to the return time to a given price. In our main results, we prove with both spectral and probabilistic methods that this rate is almost proportional to the turning frequency ω of the model and provide some explicit bounds. In the continuity of this result, we build some estimators of ω and of the pseudo-period of the prices. At last, we end the paper by a proof of the quasi-stationary distribution of the process, as well as the existence of its persistence rate.

preprint2013arXiv

Approximation of stationary solutions to SDEs driven by multiplicative fractional noise

In a previous paper, we studied the ergodic properties of an Euler scheme of a stochastic differential equation with a Gaussian additive noise in order to approximate the stationary regime of such equation. We now consider the case of multiplicative noise when the Gaussian process is a fractional Brownian Motion with Hurst parameter H>1/2 and obtain some (functional) convergences properties of some empirical measures of the Euler scheme to the stationary solutions of such SDEs.

preprint2010arXiv

Approximation of stationary solutions of Gaussian driven Stochastic Differential Equations

We study sequences of empirical measures of Euler schemes associated to some non-Markovian SDEs: SDEs driven by Gaussian processes with stationary increments. We obtain the functional convergence of this sequence to a stationary solution to the SDE. Then, we end the paper by some specific properties of this stationary solution. We show that, in contrast to Markovian SDEs, its initial random value and the driving Gaussian process are always dependent. However, under an integral representation assumption, we also obtain that the past of the solution is independent to the future of the underlying innovation process of the Gaussian driving process.

preprint2009arXiv

Approximation of the distribution of a stationary Markov process with application to option pricing

We build a sequence of empirical measures on the space D(R_+,R^d) of R^d-valued càdlàg functions on R_+ in order to approximate the law of a stationary R^d-valued Markov and Feller process (X_t). We obtain some general results of convergence of this sequence. Then, we apply them to Brownian diffusions and solutions to Lévy driven SDE's under some Lyapunov-type stability assumptions. As a numerical application of this work, we show that this procedure gives an efficient way of option pricing in stochastic volatility models.

preprint2008arXiv

A connection between extreme value theory and long time approximation of SDE's

We consider a sequence $(ξ_n)_{n\ge1}$ of $i.i.d.$ random values living in the domain of attraction of an extreme value distribution. For such sequence, there exists $(a_n)$ and $(b_n)$, with $a_n>0$ and $b_n\in\ER$ for every $n\ge 1$, such that the sequence $(X_n)$ defined by $X_n=(\max(ξ_1,...,ξ_n)-b_n)/a_n$ converges in distribution to a non degenerated distribution. In this paper, we show that $(X_n)$ can be viewed as an Euler scheme with decreasing step of an ergodic Markov process solution to a SDE with jumps and we derive a functional limit theorem for the sequence $(X_n)$ from some methods used in the long time numerical approximation of ergodic SDE's.