Past Stochastic Analysis Seminar

Abstract: First order asymptotic of scalar renewal sequences with infinite mean characterized by regular variation has been classified in the 60's (Garsia and Lamperti). In the recent years, the question of higher order asymptotic for renewal sequences with infinite mean was motivated by obtaining 'mixing rates' for dynamical systems with infinite measure. In this talk I will present the recent results we have obtained on higher order asymptotic for renewal sequences with infinite mean and their consequences for error rates in certain limit theorems (such as arcsine law for null recurrent Markov processes).

Abstract: We consider different notions of rough paths on manifolds and study some of the relations between these definitions. Furthermore, we explore extensions to manifolds modelled along infinite dimensional Banach spaces.

Abstract: We consider random walks and Lévy processes in the free nilpotent Lie group as rough paths. For any p > 1, we completely characterise (almost) all Lévy processes whose sample paths have finite p-variation, provide a Lévy-Khintchine formula for the characteristic function of the signature of a Lévy process treated as a rough path, and give sufficient conditions under which a sequence of random walks converges weakly to a Lévy process in rough path topologies. At the heart of our analysis is a criterion for tightness of p-variation for a collection of càdlàg strong Markov processes. We demonstrate applications of our results to weak convergence of stochastic flows.

Abstract: A diffusion process on a Riemannian manifold whose generator is one half of the Laplacian is called a Brownian motion. The mean local time of Brownian motion on a hypersurface will be considered, as will the situation in which a Brownian motion is conditioned to arrive in a fixed submanifold at a fixed positive time. Doing so provides motivation for the remainder of the talk, in which a probabilistic formula for the integral of the heat kernel over a submanifold is proved and used to deduce lower bounds, an asymptotic relation and derivative estimates applicable to the conditioned process.

Abstract: The martingale optimal transport aims to optimally transfer a probability measure to another along the class of martingales. This problem is mainly motivated by the robust superhedging of exotic derivatives in financial mathematics, which turns out to be the corresponding Kantorovich dual. In this paper we consider the continuous-time martingale transport on the Skorokhod space of cadlag paths. Similar to the classical setting of optimal transport, we introduce different dual problems and establish the corresponding dualities by a crucial use of the S-topology and the dynamic programming principle. This is a joint work with Gaoyue Guo and Nizar Touzi.

Abstract: Joint work with Oleg Zaboronsky (Warwick).

Some one dimensional nearest neighbour particle systems are examples of Pfaffian point processes - where all intensities are determined by a single kernel.In some cases these kernels have appeared in the random matrix literature (where the points are the positions of eigenvalues). We are attempting to use random matrix tools on the particle sytems, and particle tools on the random matrices.

Abstract:   Consider the square $[0,n]^2$ with points from a Poisson point process of intensity 1 distributed within it. In a seminal work, Baik, Deift and Johansson proved that the number of points $L_n$ (length) on a maximal increasing path (an increasing path that contains the most number of points), when properly centered and scaled, converges to the Tracy-Widom distribution. Later Johansson showed that all maximal paths lie within the strip of width $n^{\frac{2}{3} +\epsilon}$ around the diagonal with probability tending to 1 as $n \to \infty$. We shall discuss recent work on the Gaussian behaviour of the length $L_n^{(\gamma)}$ of a maximal increasing path restricted to lie within a strip of width $n^{\gamma}, \gamma< \frac{2}{3}$.

 

TBC

There is a simple way to “glue together” a coupled pair of continuum random trees to produce a topological sphere. The sphere comes equipped with a measure and a space-filling curve (which describes the “interface” between the trees). We present an explicit and canonical way to embed the sphere into the Riemann sphere. In this embedding, the measure is Liouville quantum gravity with parameter gamma in (0,2), and the curve is space-filling version of SLE with kappa=16/gamma^2. Based on joint work with Bertrand Duplantier and Scott Sheffield

We consider a controlled system, in which an input $X: [0, T] \rightarrow E:= \mathbb{R}^{d}$ is a continuous but potentially highly oscillatory path and the corresponding output $Y$ is the line integral along $X$, for some unknown function $f: E \rightarrow E$. The rough paths theory provides a general framework to answer the question on which mild condition of $X$ and $f$, the integral $I(X)$ is well defined. It is robust enough to allow to treat stochastic integrals in a deterministic way. In this paper we are interested in identification of controlled systems of this type. The difficulty comes from the high dimensionality caused by the input of a function type. We propose novel adaptive and non-parametric algorithms to learn the functional relationship between the  input and the output from the data by carefully choosing the feature set of paths based on the rough paths theory and applying linear regression techniques. The algorithms is demonstrated on a financial application where the task is to predict the P$\&$L of the unknown trading strategy.

Abstract:In the first part of the talk, using classical hypergeometric identities, I will compute the harmonic measure of finite intervals and their complementaries for the Lévy stable process on the line. This gives a simple and unified proof of several results by Blumenthal-Getoor-Ray, Rogozin, and Kyprianou-Pardo-Watson. In the second part of the talk, I will consider the two-dimensional Markov process based on the stable Lévy process and its area process. I will give two explicit formulae for the harmonic measure of the split complex plane. These formulae allow to compute the persistence exponent of the stable area process, solving a problem raised by Zhan Shi. This is based on two joint works with Christophe Profeta.

Abstract: It has been recently shown that rough volatility models reproduce very well the statistical properties of low frequency financial data. In such models, the volatility process is driven by a fractional Brownian motion with Hurst parameter of order 0.1. The goal of this talk is to explain how such fractional dynamics can be obtained from the behaviour of market participants at the microstructural scales.

Using limit theorems for Hawkes processes, we show that a rough volatility naturally arises in the presence of high frequency trading combined with metaorders splitting. This is joint work with Thibault Jaisson.

: Estimating volatility from recent high frequency data, we revisit the question of the smoothness of the volatility process. Our main result is that log-volatility behaves essentially as a fractional Brownian motion with Hurst exponent H of order 0.1, at any reasonable time scale.

This leads us to adopt the fractional stochastic volatility (FSV) model of Comte and Renault.

We call our model Rough FSV (RFSV) to underline that, in contrast to FSV, H<1/2.

We demonstrate that our RFSV model is remarkably consistent with financial time series data; one application is that it enables us to obtain improved forecasts of realized volatility.

Furthermore, we find that although volatility is not long memory in the RFSV model, classical statistical procedures aiming at detecting volatility persistence tend to conclude the presence of long memory in data generated from it.

This sheds light on why long memory of volatility has been widely accepted as a stylized fact.

Finally, we provide a quantitative market microstructure-based foundation for our findings, relating the roughness of volatility to high frequency trading and order splitting.

This is joint work with Jim Gatheral and Thibault Jaisson.

We work in the context of Markovian rough paths associated to a class of uniformly subelliptic Dirichlet forms and prove an almost-Gaussian tail-estimate for the accumulated local p-variation functional, which has been introduced and studied by Cass, Litterer and Lyons. We comment on the significance of these estimates to a range of currently-studied problems, including the recent results of Ni Hao, and Chevyrev and Lyons.

The main goal of the talk is to set up a framework for constructing the likelihood for discretely observed RDEs. The main idea is to contract a function mapping the discretely observed data to the corresponding increments of the driving noise. Once this is known, the likelihood of the observations can be written as the likelihood of the increments of the corresponding noise times the Jacobian correction.

Constructing a function mapping data to noise is equivalent to solving the inverse problem of looking for the input given the output of the Ito map corresponding to the RDE. First, I simplify the problem by assuming that the driving noise is linear between observations. Then, I will introduce an iterative process and show that it converges in p-variation to the piecewise linear path X corresponding to the observations. Finally, I will show that the total error in the likelihood construction is bounded in p-variation.

Multiplicative chaos theory originated from the study of turbulence by Kolmogorov in the 1940s and it was mathematically founded by Kahane in the 1980s. Recently the theory has drawn much of attention due to its connection to SLEs and statistical physics.  In this talk I shall present some recent development of multiplicative chaos theory, as well as its applications to Liouville quantum gravity.

We consider the problem of large-scale regression where both the number of predictors, p, and the number of observations, n, may be in the order of millions or more. Computing a simple OLS or ridge regression estimator for such data, though potentially sensible from a purely statistical perspective (if n is large enough), can be a real computational challenge. One recent approach to tackling this problem in the common situation where the matrix of predictors is sparse, is to first compress the data by mapping it to an n by L matrix with L << p, using a scheme called b-bit min-wise hashing (Li and König, 2011). We study this technique from a theoretical perspective and obtain finite-sample bounds on the prediction error of regression following such data compression, showing how it exploits the sparsity of the data matrix to achieve good statistical performance. Surprisingly, we also find that a main effects model in the compressed data is able to approximate an interaction model in the original data. Fitting interactions requires no modification of the compression scheme, but only a higher-dimensional mapping with a larger L.
This is joint work with Nicolai Meinshausen (ETH Zürich).

The model of random interlacements is a one-parameter family of random subsets of $\Z^d$, which locally describes the trace of a simple random walk on a $d$-dimensional torus running up to time $u$ times its volume. Here, $u$ serves as an intensity parameter.

Its complement, the so-called vacant set, has been show to undergo a non-trivial percolation phase transition in $u$, i.e., there is $u_*(d)\in (0,\infty)$ such that for all $u<u_*(d)$ the vacant set has a unique infinite connected component (supercritical phase), while for $u>u_*(d)$ all connected components are finite.

So far all results regarding geometric properties of this infinite connected component have been proven under the assumption that $u$ is close to zero. 

I will discuss a recent result, which states that throughout most of the supercritical phase simple random walk on the infinite connected component is transient, provided that the dimension is high enough.

This is joint work with Alexander Drewitz

Abstract: L\'evy processes are increasingly popular for modelling stochastic process data with jump behaviour. In practice statisticians only observe discretely sampled increments of the process, leading to a statistical inverse problem. To understand the jump behaviour of the process one needs to make inference on the infinite-dimensional parameter given by the L\'evy measure. We discuss recent developments in the analysis of this problem, including in particular functional limit theorems for commonly used estimators of the generalised distribution function of the L\'evy measure, and their application to statistical uncertainty quantification methodology (confidence bands and tests). 

We consider the problem of minimising the commute or shuttle time for a diffusion between the endpoints of an interval. The control is the scale function for the diffusion. We show that the dynamic version of the problem has the same solution as the static version if we start at an end point and consider the much harder case where the starting point is in the interior.

tba

tba

 We describe how cross-kernel matrices, that is, kernel matrices between the data and a custom chosen set of `feature spanning points' can be used for learning. The main potential of cross-kernel matrices is that (a) they provide Nyström-type speed-ups for kernel learning without relying on subsampling, thus avoiding potential problems with sampling degeneracy, while preserving the usual approximation guarantees and the attractive linear scaling of standard Nyström methods and (b) the use of non-square matrices for kernel learning provides a non-linear generalization of the singular value decomposition and singular features. We present a novel algorithm, Ideal PCA (IPCA), which is a cross-kernel matrix variant of PCA, showcasing both advantages: we demonstrate on real and synthetic data that IPCA allows to (a) obtain kernel PCA-like features faster and (b) to extract novel features of empirical advantage in non-linear manifold learning and classification.

tba

The Renormalization Group (RG) was pioneered by the physicist Kenneth Wilson in the early 70's and since then it has become a fundamental tool in physics. RG remains the most general philosophy for understanding how many models in statistical mechanics behave near their critical point but implementing RG analysis in a mathematically rigorous way remains quite challenging.

I will describe how analysis of RG flows translate into statements about continuum limits, universality, and cross-over phenomena - as a concrete example I will speak about some joint work with Abdelmalek Abdesselam and Gianluca Guadagni.

Spectral volume and surface measures via the Dixmier trace for local symmetric Dirichlet spaces with Weyl type eigenvalue asymptotics

The purpose of this talk is to present the author's recent results of on an

operator theoretic way of looking atWeyl type Laplacian eigenvalue asymptotics

for local symmetric Dirichlet spaces.

For the Laplacian on a d-dimensional Riemannian manifoldM, Connes' trace

theorem implies that the linear functional  coincides with

(a constant multiple of) the integral with respect to the Riemannian volume

measure of M, which could be considered as an operator theoretic paraphrase

of Weyl's Laplacian eigenvalue asymptotics. Here  denotes a Dixmier trace,

which is a trace functional de_ned on a certain ideal of compact operators on

a Hilbert space and is meaningful e.g. for compact non-negative self-adjoint

operators whose n-th largest eigenvalue is comparable to 1/n.

The first main result of this talk is an extension of this fact in the framework

of a general regular symmetric Dirichlet space satisfying Weyl type asymptotics

for the trace of its associated heat semigroup, which was proved for Laplacians

on p.-c.f. self-simiar sets by Kigami and Lapidus in 2001 under a rather strong

assumption.

Moreover, as the second main result of this talk it is also shown that, given a

local regular symmetric Dirichlet space with a sub-Gaussian heat kernel upper

bound and a (sufficiently regular) closed subset S, a “spectral surface measure"

on S can be obtained through a similar linear functional involving the Lapla-

cian with Dirichlet boundary condition on S. In principle, corresponds to the

second order term for the eigenvalue asymptotics of this Dirichlet Laplacian, and

when the second order term is explicitly known it is possible to identify  For

example, in the case of the usual Laplacian on Rd and a Lipschitz hypersurface S,is a constant multiple of the usual surface measure on S.

Maximal couplings are couplings of Markov processes where the tail probabilities of the coupling time attain the total variation lower bound (Aldous bound) uniformly for all time. Markovian couplings are coupling strategies where neither process is allowed to look into the future of the other before making the next transition. These are easier to describe and play a fundamental role in many branches of probability and analysis. Hsu and Sturm proved that the reflection coupling of Brownian motion is the unique Markovian maximal coupling (MMC) of Brownian motions starting from two different points. Later, Kuwada proved that to have a MMC for Brownian motions on a Riemannian manifold, the manifold should have a reflection structure, and thus proved the first result connecting this purely probabilistic phenomenon (MMC) to the geometry of the underlying space.

Motivated by stochastic models for order books in stock exchanges we consider stochastic partial differential equations with a free boundary condition. Such equations can be considered generalizations of the classic (deterministic) Stefan problem of heat condition in a two-phase medium. 

Extending results by Kim, Zheng & Sowers we allow for non-linear boundary interaction, general Robin-type boundary conditions and fairly general drift and diffusion coefficients. Existence of maximal local and global solutions is established by transforming the equation to a fixed-boundary problem and solving a stochastic evolution equation in suitable interpolation spaces. Based on joint work with Marvin Mueller.

'We study an optimal switching problem with a state constraint: the controller is only allowed to choose strategies that keep the controlled diffusion in a closed domain. We prove that the value function associated to the weak formulation of this problem is the limit of the value function associated to an unconstrained switching problem with penalized coefficients, as the penalization parameter goes to infinity. This convergence allows to set a dynamic programming principle for the constrained switching problem. We then prove that the value function is a constrained viscosity solution to a system of variational inequalities (SVI for short). We finally prove that the value function is the maximal solution to this SVI. All our results are obtained without any regularity assumption on the constraint domain.’

Lawler, Schramm and Werner studied 2-point chordal restriction measures and gave several constructions using SLE tools.

It is possible to characterize 3-point chordal restriction measures in a similar manner. Their boundaries are SLE(8/3)-like curves with a slightly different drift term.

@email

Cubature is the business of describing a probability measure in terms of an empirical measure sharing its support with the original measure, of small support, and with identical integrals for a class of functions (eg polynomials with degree less than k). 

Applying cubature to already discrete sets of scenarios provides a powerful tool for scenario management and summarising data.  We refer to this process as recombination. It is a feasible operation in real time and has lead to high accuracy pde solvers.

The practical complexity of this operation has changed! By a factor corresponding to the dimension of the space of polynomials. 

We discuss the algorithm and give home computed examples of nested sparse grids with only positive weights in moderate dimensions (eg degree 1-8 in dimension 7).  Positive weights have significant advantage over signed ones when available.
 

   Stéphane Mallat

   Ecole Normale Superieure

Learning functionals in high dimension requires to find sources of regularity and invariants, to reduce dimensionality. Stability to actions of diffeomorphisms is a strong property satisfied by many physical functionals and most signal classification problems. We introduce a scattering operator in a path space, calculated with iterated multiscale wavelet transforms, which is invariant to rigid movements and stable to diffeomorphism actions. It provides a Euclidean embedding of geometric distances and a representation of stationary random processes. Applications will be shown for image classification and to learn quantum chemistry energy functionals.

In this talk we study the large time behaviour of some semilinear parabolic PDEs by a purely probabilistic approach. For that purpose, we show that the solution of a backward stochastic differential equation (BSDE) in finite horizon $T$ taken at initial time behaves like a linear term in $T$ shifted with a solution of the associated ergodic BSDE taken at inital time. Moreover we give an explicit rate of convergence: we show that the following term in the asymptotic expansion has an exponential decay. This is a Joint work with Ying Hu and Pierre-Yves Meyer from Rennes (IRMAR - France).

In this talk, we develop rough integration with jumps, offering a pathwise view on stochastic integration against cadlag processes.  A class of Marcus-like rough paths is introduced,which contains D. Williams’ construction of stochastic area for Lévy processes. We then established a Lévy–Khintchine type formula for the expected signature, based on“Marcus(canonical)"stochastic calculus. This calculus fails for non-Marcus-like Lévy rough paths and we treat the general case with Hunt’ theory of Lie group valued Lévy processes is made.

In this talk we will present some recent developments in the theory of ambit fields with a particular focuson limit theorems.
Ambit fields is a tempo-spatial class of models, which has been originally introduced by Barndorff-Nielsen and Schmiegel in the context of turbulence,
but found applications also in biology and finance. Its purely temporal analogue, Levy semi-stationary processes, has a continuous moving average structure
with an additional multiplicative random input (volatility or intermittency). We will briefly describe the main challenges of ambit stochastics, which
include questions from stochastic analysis, statistics and numerics. We will then focus on certain type of high frequency functionals typically called power variations.
We show some surprising non-standard limit theorems, which strongly depend on the driving Levy process. The talk is based on joint work with O.E. Barndorff-Nielsen, A. Basse-O'Connor,
J.M. Corcuera and R. Lachieze-Rey. 

Motivated by stochastic models for order books in stock exchanges we consider stochastic partial differential equations with a free boundary condition. Such equations can be considered generalizations of the classic (deterministic) Stefan problem of heat condition in a two-phase medium. 

Extending results by Kim, Zheng & Sowers we allow for non-linear boundary interaction, general Robin-type boundary conditions and fairly general drift and diffusion coefficients. Existence of maximal local and global solutions is established by transforming the equation to a fixed-boundary problem and solving a stochastic evolution equation in suitable interpolation spaces. Based on joint work with Marvin Mueller.

@email 

I will describe the Spatial Lambda-Fleming-Viot process, which is a model of evolution in a spatial continuum, and discuss the time and spatial scales on which selectively advantageous genes propagate through space. The appropriate scaling depends on the dimension of space, resulting in three distinct cases; d=1, d=2 and d>=3. In d=1 the limiting genealogy is the Brownian net whereas, by contrast, in d=2 local interactions give rise to a delicate damping mechanism and result in a finite limiting branching rate. This is joint work with Alison Etheridge and Daniel Straulino.

The theory of regularity structures allows one to give a meaning to several stochastic PDEs, including the Parabolic Anderson Model. So far, these equations have been considered on a torus. The goal of this talk is to explain how one can define the PAM on the whole space R^3. This is a joint work with Martin Hairer.

In the Erdös-Rényi random graph process, starting from an empty graph, in each step a new random edge is added to the evolving graph. One of its most interesting features is the `percolation phase transition': as the ratio of the number of edges to vertices increases past a certain critical density, the global structure changes radically, from only small components to a single giant component plus small ones.

In this talk we consider Achlioptas processes, which have become a key example for random graph processes with dependencies between the edges.

Starting from an empty graph these proceed as follows: in each step two potential edges are chosen uniformly at random, and using some rule one of them is selected and added to the evolving graph. We discuss why, for a large class of rules, the percolation phase transition is qualitatively comparable to the classical Erdös-Rényi process.

Based on joint work with Oliver Riordan.

We give a necessary and sufficient condition for a map defined on a compact, quasiconvex and simply-connected space to factor through a tree. This condition can be checked using currents. In particular if the target is some Euclidean space and the map is H\"older continuous with exponent bigger than 1/2, such maps can be characterized by the vanishing of some integrals over the winding number. Moreover, this shows that if the target is the Heisenberg group equipped with the Carnot-Carath\'eodory metric and the H\"older exponent of the map is bigger than 2/3, the map factors through a tree.

The idea of A-free group, where A is a discrete ordered abelian group, has been introduced by Myasnikov, Remeslennikov and Serbin. It generalises the construction of free groups. A proof will be outlined that a group is A-free for some A if and only if it acts freely and without inversions on a \lambda-tree, where \lambda is an arbitrary ordered abelian group.

The Taylor expansion of a controlled differential equation suggests that the solution at time 1 depends on the driving path only through the latter's iterated integrals up to time 1, if the vector field is infinitely differentiable. Hambly and Lyons proved that this remains true for Lipschitz vector fields if the driving path has bounded total variation. We extend the Hambly-Lyons result for weakly geometric rough paths in finite dimension. Joint work with X. Geng, T. Lyons and D. Yang.    

We consider the numerical approximation of Kolmogorov equations arising in the context of option pricing under L\'evy models and beyond in a multivariate setting. The existence and uniqueness of variational solutions of the partial integro-differential equations (PIDEs) is established in Sobolev spaces of fractional or variable order.

Most discretization methods for the considered multivariate models suffer from the curse of dimension which impedes an efficient solution of the arising systems. We tackle this problem by the use of sparse discretization methods such as classical sparse grids or tensor train techniques. Numerical examples in multiple space dimensions confirm the efficiency of the described methods.

This work is joint with Peter Sarnak.

It is shown that the topologies and nestings of the zero and nodal sets of random (Gaussian) band limited functions have universal laws of distribution. Qualitative features of the supports of these distributions are determined. In particular the results apply to random monochromatic waves and to random real algebraic hyper-surfaces in projective space.

We study infinite (random) systems of interacting particles living in a Euclidean space X and possessing internal parameter (spin) in R¹. Such systems are described by Gibbs measures on the space Γ(X,R¹) of marked configurations in X (with marks in R¹). For a class of pair interactions, we show the occurrence of phase transition, i.e. non-uniqueness of the corresponding Gibbs measure, in both 'quenched' and 'annealed' counterparts of the model.

A model of a financial market is complete if any payoff can be obtained as the terminal value of a self-financing trading strategy. It is well known that numerous models, for example stochastic volatility models, are however incomplete. We present conditions, which, in a general diffusion framework, guarantee that in such cases the market of primitive assets enlarged with an appropriate number of traded derivative contracts is complete. From a purely mathematical point of view we prove an integral representation theorem which guarantees that every local Q-martingale can be represented as a stochastic integral with respect to the vector of primitive assets and derivative contracts.

Abstract: We introduce a new framework for defining integration against rough path. This framework generalizes rough integral, and gives a natural explanation of some of the regularity requirements in rough path theory.

When $\partial D$ is spatially homogeneous, we show that we can recover the lower and upper Minkowski dimensions of $\partial D$ from the sort time behaviour of $E(s)$. Furthermore, when the Minkowski dimension exists, finer geometric fluctuations can be recovered and $E(s)$ is controlled by $s^\alpha e^{f(\log(1/s))}$ for small $s$, for some $\alpha \in (0, \infty)$ and some regularly varying function $f$. The function $f$ is not constant is general and carries some geometric information.

When $\partial D$ is statistically self-similar, the Minkowski dimension and content of $\partial D$ typically exist and can be recovered from $E(s)$. Furthermore, $E(s)$ has an almost sure expansion $E(s) = c s^{\alpha} N_\infty + o(s^\alpha)$ for small $s$, for some $c$ and $\alpha \in (0, \infty)$ and some positive random variable $N_\infty$ with unit expectation arising as the limit of some martingale. In some cases, we can show that the fluctuations around this almost sure behaviour are governed by a central limit theorem, and conjecture that this is true more generally.

This is based on joint work with David Croydon and Ben Hambly.