Showing papers in "arXiv: Probability in 2007"

G -Expectation, G -Brownian Motion and Related Stochastic Calculus of Itô Type

Abstract: We introduce a notion of nonlinear expectation --G--expectation-- generated by a nonlinear heat equation with infinitesimal generator G. We first discuss the notion of G-standard normal distribution. With this nonlinear distribution we can introduce our G-expectation under which the canonical process is a G--Brownian motion. We then establish the related stochastic calculus, especially stochastic integrals of Ito's type with respect to our G--Brownian motion and derive the related Ito's formula. We have also give the existence and uniqueness of stochastic differential equation under our G-expectation. As compared with our previous framework of g-expectations, the theory of G-expectation is intrinsic in the sense that it is not based on a given (linear) probability space.

Compressed Sensing and Redundant Dictionaries

Abstract: This article extends the concept of compressed sensing to signals that are not sparse in an orthonormal basis but rather in a redundant dictionary. It is shown that a matrix, which is a composition of a random matrix of certain type and a deterministic dictionary, has small restricted isometry constants. Thus, signals that are sparse with respect to the dictionary can be recovered via Basis Pursuit from a small number of random measurements. Further, thresholding is investigated as recovery algorithm for compressed sensing and conditions are provided that guarantee reconstruction with high probability. The different schemes are compared by numerical experiments.

The Littlewood-Offord Problem and invertibility of random matrices

Abstract: We prove two basic conjectures on the distribution of the smallest singular value of random n times n matrices with independent entries. Under minimal moment assumptions, we show that the smallest singular value is of order n^{-1/2}, which is optimal for Gaussian matrices. Moreover, we give a optimal estimate on the tail probability. This comes as a consequence of a new and essentially sharp estimate in the Littlewood-Offord problem: for i.i.d. random variables X_k and real numbers a_k, determine the probability P that the sum of a_k X_k lies near some number v. For arbitrary coefficients a_k of the same order of magnitude, we show that they essentially lie in an arithmetic progression of length 1/p.

Differential equation approximations for Markov chains

Abstract: We formulate some simple conditions under which a Markov chain may be approximated by the solution to a differential equation, with quantifiable error probabilities. The role of a choice of coordinate functions for the Markov chain is emphasised. The general theory is illustrated in three examples: the classical stochastic epidemic, a population process model with fast and slow variables, and core-finding algorithms for large random hypergraphs.

G-Brownian Motion and Dynamic Risk Measure under Volatility Uncertainty

Abstract: We introduce a new notion of G-normal distributions. This will bring us to a new framework of stochastic calculus of Ito's type (Ito's integral, Ito's formula, Ito's equation) through the corresponding G-Brownian motion. We will also present analytical calculations and some new statistical methods with application to risk analysis in finance under volatility uncertainty. Our basic point of view is: sublinear expectation theory is very like its special situation of linear expectation in the classical probability theory. Under a sublinear expectation space we still can introduce the notion of distributions, of random variables, as well as the notions of joint distributions, marginal distributions, etc. A particularly interesting phenomenon in sublinear situations is that a random variable Y is independent to X does not automatically implies that X is independent to Y. Two important theorems have been proved: The law of large number and the central limit theorem.

Rate of Converrgence for ergodic continuous Markov processes : Lyapunov versus Poincare

Abstract: We study the relationship between two classical approaches for quantitative ergodic properties : the first one based on Lyapunov type controls and popularized by Meyn and Tweedie, the second one based on functional inequalities (of Poincar\'e type). We show that they can be linked through new inequalities (Lyapunov-Poincar\'e inequalities). Explicit examples for diffusion processes are studied, improving some results in the literature. The example of the kinetic Fokker-Planck equation recently studied by H\'erau-Nier, Helffer-Nier and Villani is in particular discussed in the final section.

Martingale proofs of many-server heavy-traffic limits for Markovian queues

Abstract: This is an expository review paper illustrating the ``martingale method'' for proving many-server heavy-traffic stochastic-process limits for queueing models, supporting diffusion-process approximations. Careful treatment is given to an elementary model -- the classical infinite-server model $M/M/\infty$, but models with finitely many servers and customer abandonment are also treated. The Markovian stochastic process representing the number of customers in the system is constructed in terms of rate-1 Poisson processes in two ways: (i) through random time changes and (ii) through random thinnings. Associated martingale representations are obtained for these constructions by applying, respectively: (i) optional stopping theorems where the random time changes are the stopping times and (ii) the integration theorem associated with random thinning of a counting process. Convergence to the diffusion process limit for the appropriate sequence of scaled queueing processes is obtained by applying the continuous mapping theorem. A key FCLT and a key FWLLN in this framework are established both with and without applying martingales.

Strong invariance principles for dependent random variables

Abstract: We establish strong invariance principles for sums of stationary and ergodic processes with nearly optimal bounds. Applications to linear and some nonlinear processes are discussed. Strong laws of large numbers and laws of the iterated logarithm are also obtained under easily verifiable conditions.

Anticipated backward stochastic differential equations

Abstract: In this paper we discuss new types of differential equations which we call anticipated backward stochastic differential equations (anticipated BSDEs). In these equations the generator includes not only the values of solutions of the present but also the future. We show that these anticipated BSDEs have unique solutions, a comparison theorem for their solutions, and a duality between them and stochastic differential delay equations.

The largest eigenvalues of finite rank deformation of large Wigner matrices: convergence and nonuniversality of the fluctuations

Abstract: In this paper, we investigate the asymptotic spectrum of complex or real Deformed Wigner matrices $(M_N)_N$ defined by $M_N=W_N/\sqrt{N}+A_N$ where $W_N$ is an $N\times N$ Hermitian (resp., symmetric) Wigner matrix whose entries have a symmetric law satisfying a Poincare inequality. The matrix $A_N$ is Hermitian (resp., symmetric) and deterministic with all but finitely many eigenvalues equal to zero. We first show that, as soon as the first largest or last smallest eigenvalues of $A_N$ are sufficiently far from zero, the corresponding eigenvalues of $M_N$ almost surely exit the limiting semicircle compact support as the size $N$ becomes large. The corresponding limits are universal in the sense that they only involve the variance of the entries of $W_N$. On the other hand, when $A_N$ is diagonal with a sole simple nonnull eigenvalue large enough, we prove that the fluctuations of the largest eigenvalue are not universal and vary with the particular distribution of the entries of $W_N$.

On the optimal dividend problem for a spectrally negative L\'{e}vy process

Abstract: In this paper we consider the optimal dividend problem for an insurance company whose risk process evolves as a spectrally negative L\'{e}vy process in the absence of dividend payments. The classical dividend problem for an insurance company consists in finding a dividend payment policy that maximizes the total expected discounted dividends. Related is the problem where we impose the restriction that ruin be prevented: the beneficiaries of the dividends must then keep the insurance company solvent by bail-out loans. Drawing on the fluctuation theory of spectrally negative L\'{e}vy processes we give an explicit analytical description of the optimal strategy in the set of barrier strategies and the corresponding value function, for either of the problems. Subsequently we investigate when the dividend policy that is optimal among all admissible ones takes the form of a barrier strategy.

Vacant Set of Random Interlacements and Percolation

Abstract: We introduce a model of random interlacements made of a countable collection of doubly infinite paths on Z^d, d bigger or equal to 3. A non-negative parameter u measures how many trajectories enter the picture. This model describes in the large N limit the microscopic structure in the bulk, which arises when considering the disconnection time of a discrete cylinder with base a d-1-dimensional discrete torus of side-length N, or the set of points visited by simple random walk on the d-dimensional discrete torus of side-length N by times of order uN^d. We study the percolative properties of the vacant set left by the interlacement at level u, which is an infinite, connected, translation invariant random subset of Z^d. We introduce a critical value such that the vacant set percolates for u below the critical value, and does not percolate for u above the critical value. Our main results show that the critical value is finite when d is bigger or equal to 3, and strictly positive when d is bigger or equal to 7.

Differential equations driven by rough paths: an approach via discrete approximation

Abstract: A theory of differential equations driven by a non-differentiable path has recently been developed by Lyons. We develop an alternative approach to this theory, using (modified Euler approximations), and investigate its applicability to stochastic differential equations driven by Brownian motion. We also give some other examples showing that the main results are reasonably sharp.

Negative dependence and the geometry of polynomials

Abstract: We introduce the class of {\em strongly Rayleigh} probability measures by means of geometric properties of their generating polynomials that amount to the stability of the latter. This class covers important models such as determinantal measures (e.g. product measures, uniform random spanning tree measures) and distributions for symmetric exclusion processes. We show that strongly Rayleigh measures enjoy all virtues of negative dependence and we also prove a series of conjectures due to Liggett, Pemantle, and Wagner, respectively. Moreover, we extend Lyons' recent results on determinantal measures and we construct counterexamples to several conjectures of Pemantle and Wagner on negative dependence and ultra log-concave rank sequences.

Integral Formulas for the Asymmetric Simple Exclusion Process

Abstract: In this paper we obtain general integral formulas for probabilities in the asymmetric simple exclusion process (ASEP) on the integer lattice with nearest neighbor hopping rates p to the right and q=1-p to the left. For the most part we consider an N-particle system but for certain of these formulas we can take the limit as N goes to infinity. First we obtain, for the N-particle system, a formula for the probability of a configuration at time t, given the initial configuration. For this we use Bethe Ansatz ideas to solve the master equation, extending a result of Schuetz for the case N=2. The main results of the paper, derived from this, are integral formulas for the probability, for given initial configuration, that the m'th left-most particle is at x at time t. In one of these formulas we can take the limit as N goes to infinity, and it gives the probability for an infinite system where the initial configuration is bounded on one side. For the special case of the totally asymmetric simple exclusion process (TASEP) our formulas reduce to the known ones.

Multivariate normal approximation with Stein's method of exchangeable pairs under a general linearity condition

Abstract: In this paper we establish a multivariate exchangeable pairs approach within the framework of Stein's method to assess distributional distances to potentially singular multivariate normal distributions. By extending the statistics into a higher-dimensional space, we also propose an embedding method which allows for a normal approximation even when the corresponding statistics of interest do not lend themselves easily to Stein's exchangeable pairs approach. To illustrate the method, we provide the examples of runs on the line as well as double-indexed permutation statistics.

Fluctuations of eigenvalues and second order Poincar\'e inequalities

Abstract: Linear statistics of eigenvalues in many familiar classes of random matrices are known to obey gaussian central limit theorems. The proofs of such results are usually rather difficult, involving hard computations specific to the model in question. In this article we attempt to formulate a unified technique for deriving such results via relatively soft arguments. In the process, we introduce a notion of `second order Poincar\'e inequalities': just as ordinary Poincar\'e inequalities give variance bounds, second order Poincar\'e inequalities give central limit theorems. The proof of the main result employs Stein's method of normal approximation. A number of examples are worked out, some of which are new. One of the new results is a CLT for the spectrum of gaussian Toeplitz matrices.

Lectures on two-dimensional critical percolation

Abstract: These are the notes corresponding to the course given at the IAS-Park City graduate summer school in July 2007.

A fractional generalization of the Poisson processes

Abstract: It is our intention to provide via fractional calculus a generalization of the pure and compound Poisson processes, which are known to play a fundamental role in renewal theory, without and with reward, respectively. We first recall the basic renewal theory including its fundamental concepts like waiting time between events, the survival probability, the counting function. If the waiting time is exponentially distributed we have a Poisson process, which is Markovian. However, other waiting time distributions are also relevant in applications, in particular such ones with a fat tail caused by a power law decay of its density. In this context we analyze a non-Markovian renewal process with a waiting time distribution described by the Mittag-Leffler function. This distribution, containing the exponential as particular case, is shown to play a fundamental role in the infinite thinning procedure of a generic renewal process governed by a power asymptotic waiting time. We then consider the renewal theory with reward that implies a random walk subordinated to a renewal process.

SLE and the free field: Partition functions and couplings

Abstract: Schramm-Loewner Evolutions ($\SLE$) are random curves in planar simply connected domains; the massless (Euclidean) free field in such a domain is a random distribution. Both have conformal invariance properties in law. In the present article, some relations between the two objects are studied. We establish identities of partition functions between different versions of $\SLE$ and the free field with appropriate boundary conditions; this involves $\zeta$-regularization and the Polyakov-Alvarez conformal anomaly formula. We proceed with a construction of couplings of $\SLE$ with the free field, showing that, in a precise sense, chordal $\SLE$ is the solution of a stochastic "differential" equation driven by the free field. Existence, uniqueness in law, and pathwise uniqueness for these SDEs are proved for general $\kappa>0$.

Differential Equations Driven by Hölder Continuous Functions of Order Greater than 1/2

Abstract: We derive estimates for the solutions to differential equations driven by a H\"older continuous function of order $\beta>1/2$. As an application we deduce the existence of moments for the solutions to stochastic partial differential equations driven by a fractional Brownian motion with Hurst parameter $H>{1/2}$.

Stochastic Heat Equation Driven by Fractional Noise and Local Time

Abstract: The aim of this paper is to study the $d$-dimensional stochastic heat equation with a multiplicative Gaussian noise which is white in space and it has the covariance of a fractional Brownian motion with Hurst parameter $% H\in (0,1)$ in time. Two types of equations are considered. First we consider the equation in the It\^{o}-Skorohod sense, and later in the Stratonovich sense. An explicit chaos development for the solution is obtained. On the other hand, the moments of the solution are expressed in terms of the exponential moments of some weighted intersection local time of the Brownian motion.

An explicit formula for the Skorokhod map on $[0,a]$

Abstract: The Skorokhod map is a convenient tool for constructing solutions to stochastic differential equations with reflecting boundary conditions In this work, an explicit formula for the Skorokhod map $\Gamma_{0,a}$ on $[0,a]$ for any $a>0$ is derived Specifically, it is shown that on the space $\mathcal{D}[0,\infty)$ of right-continuous functions with left limits taking values in $\mathbb{R}$, $\Gamma_{0,a}=\Lambda_a\circ \Gamma_0$, where $\Lambda_a:\mathcal{D}[0,\infty)\to\mathcal{D}[0,\infty)$ is defined by \[\Lambda_a(\phi)(t)=\phi(t)-\sup_{s\in[0,t]}\biggl[\bigl(\ phi(s)-a\bigr)^+\wedge\inf_{u\in[s,t]}\phi(u)\biggr]\] and $\Gamma_0:\mathcal{D}[0,\infty)\to\mathcal{D}[0,\infty)$ is the Skorokhod map on $[0,\infty)$, which is given explicitly by \[\Gamma_0(\psi)(t)=\psi(t)+\sup_{s\in[0,t]}[-\psi(s)]^+\] In addition, properties of $\Lambda_a$ are developed and comparison properties of $\Gamma_{0,a}$ are established

Law of Large Numbers and Central Limit Theorem under Nonlinear Expectations

Abstract: The law of large numbers (LLN) and central limit theorem (CLT) are long and widely been known as two fundamental results in probability theory. Recently problems of model uncertainties in statistics, measures of risk and superhedging in finance motivated us to introduce, in [4] and [5] (see also [2], [3] and references herein), a new notion of sublinear expectation, called \textquotedblleft% $G$-expectation\textquotedblright, and the related \textquotedblleft$G$-normal distribution\textquotedblright from which we were able to define G-Brownian motion as well as the corresponding stochastic calculus. The notion of G-normal distribution plays the same important rule in the theory of sublinear expectation as that of normal distribution in the classic probability theory. It is then natural and interesting to ask if we have the corresponding LLN and CLT under a sublinear expectation and, in particular, if the corresponding limit distribution of the CLT is a G-normal distribution. This paper gives an affirmative answer. The proof of our CLT is short since we borrow a deep interior estimate of fully nonlinear PDE in [6] which extended a profound result of [1] (see also [7]) to parabolic PDEs. The assumptions of our LLN and CLT can be still improved. But the discovered phenomenon plays the same important rule in the theory of nonlinear expectation as that of the classical LLN and CLT in classic probability theory.

On asymptotics of eigenvectors of large sample covariance matrix

Abstract: Let \{$X_{ij}$\}, $i,j=...,$ be a double array of i.i.d. complex random variables with $EX_{11}=0,E|X_{11}|^2=1$ and $E|X_{11}|^4<\infty$, and let $A_n=\frac{1}{N}T_n^{{1}/{2}}X_nX_n^*T_n^{{1}/{2}}$, where $T_n^{{1}/{2}}$ is the square root of a nonnegative definite matrix $T_n$ and $X_n$ is the $n\times N$ matrix of the upper-left corner of the double array. The matrix $A_n$ can be considered as a sample covariance matrix of an i.i.d. sample from a population with mean zero and covariance matrix $T_n$, or as a multivariate $F$ matrix if $T_n$ is the inverse of another sample covariance matrix. To investigate the limiting behavior of the eigenvectors of $A_n$, a new form of empirical spectral distribution is defined with weights defined by eigenvectors and it is then shown that this has the same limiting spectral distribution as the empirical spectral distribution defined by equal weights. Moreover, if \{$X_{ij}$\} and $T_n$ are either real or complex and some additional moment assumptions are made then linear spectral statistics defined by the eigenvectors of $A_n$ are proved to have Gaussian limits, which suggests that the eigenvector matrix of $A_n$ is nearly Haar distributed when $T_n$ is a multiple of the identity matrix, an easy consequence for a Wishart matrix.

Brownian excursion area, Wright's constants in graph enumeration, and other Brownian areas

Abstract: This survey is a collection of various results and formulas by different authors on the areas (integrals) of five related processes, viz.\spacefactor =1000 Brownian motion, bridge, excursion, meander and double meander; for the Brownian motion and bridge, which take both positive and negative values, we consider both the integral of the absolute value and the integral of the positive (or negative) part. This gives us seven related positive random variables, for which we study, in particular, formulas for moments and Laplace transforms; we also give (in many cases) series representations and asymptotics for density functions and distribution functions. We further study Wright's constants arising in the asymptotic enumeration of connected graphs; these are known to be closely connected to the moments of the Brownian excursion area. The main purpose is to compare the results for these seven Brownian areas by stating the results in parallel forms; thus emphasizing both the similarities and the differences. A recurring theme is the Airy function which appears in slightly different ways in formulas for all seven random variables. We further want to give explicit relations between the many different similar notations and definitions that have been used by various authors. There are also some new results, mainly to fill in gaps left in the literature. Some short proofs are given, but most proofs are omitted and the reader is instead referred to the original sources.

Some recent advances in theory and simulation of fractional diffusion processes

Abstract: To offer a view into the rapidly developing theory of fractional diffusion processes we describe in some detail three topics of present interest: (i) the well-scaled passage to the limit from continuous time random walk under power law assumptions to space-time fractional diffusion, (ii) the asymptotic universality of the Mittag-Leffler waiting time law in time-fractional processes, (iii) our method of parametric subordination for generating particle trajectories.

Strong Spherical Asymptotics for Rotor-Router Aggregation and the Divisible Sandpile

Abstract: The rotor-router model is a deterministic analogue of random walk. It can be used to define a deterministic growth model analogous to internal DLA. We prove that the asymptotic shape of this model is a Euclidean ball, in a sense which is stronger than our earlier work. For the shape consisting of $n=\omega_d r^d$ sites, where $\omega_d$ is the volume of the unit ball in $\R^d$, we show that the inradius of the set of occupied sites is at least $r-O(\log r)$, while the outradius is at most $r+O(r^\alpha)$ for any $\alpha > 1-1/d$. For a related model, the divisible sandpile, we show that the domain of occupied sites is a Euclidean ball with error in the radius a constant independent of the total mass. For the classical abelian sandpile model in two dimensions, with $n=\pi r^2$ particles, we show that the inradius is at least $r/\sqrt{3}$, and the outradius is at most $(r+o(r))/\sqrt{2}$. This improves on bounds of Le Borgne and Rossin. Similar bounds apply in higher dimensions.

Central and non-central limit theorems for weighted power variations of fractional Brownian motion

Abstract: In this paper, we prove some central and non-central limit theorems for renormalized weighted power variations of order q>=2 of the fractional Brownian motion with Hurst parameter H in (0,1), where q is an integer. The central limit holds for 1/(2q) 1-1/(2q), we show the convergence in L^2 to a stochastic integral with respect to the Hermite process of order q.

Quasi-stationary distributions and diffusion models in population dynamics

Abstract: In this paper, we study quasi-stationarity for a large class of Kolmogorov diffusions. The main novelty here is that we allow the drift to go to $- \infty$ at the origin, and the diffusion to have an entrance boundary at $+\infty$. These diffusions arise as images, by a deterministic map, of generalized Feller diffusions, which themselves are obtained as limits of rescaled birth--death processes. Generalized Feller diffusions take nonnegative values and are absorbed at zero in finite time with probability 1. An important example is the logistic Feller diffusion. We give sufficient conditions on the drift near 0 and near $+ \infty$ for the existence of quasi-stationary distributions, as well as rate of convergence in the Yaglom limit and existence of the $Q$-process. We also show that under these conditions, there is exactly one quasi-stationary distribution, and that this distribution attracts all initial distributions under the conditional evolution, if and only if $+\infty$ is an entrance boundary. In particular this gives a sufficient condition for the uniqueness of quasi-stationary distributions. In the proofs spectral theory plays an important role on $L^2$ of the reference measure for the killed process.