Some Examples. Uniform motion. Poisson processes on the real line

Similar documents
Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

PDE Notes, Lecture #11

Separation of Variables

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k

Convergence of Random Walks

Introduction to the Vlasov-Poisson system

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

Levy Process and Infinitely Divisible Law

DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10

Tractability results for weighted Banach spaces of smooth functions

Linear First-Order Equations

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

Differentiation Rules Derivatives of Polynomials and Exponential Functions

Davar Khoshnevisan. Topics in Probability: Lévy Processes Math ; Spring 2011

CLARK-OCONE FORMULA BY THE S-TRANSFORM ON THE POISSON WHITE NOISE SPACE

6 General properties of an autonomous system of two first order ODE

Table of Common Derivatives By David Abraham

The Ehrenfest Theorems

Self-normalized Martingale Tail Inequality

PROBABILITY THEORY EXPONENTIAL FUNCTIONAL OF LÉVY PROCESSES: GENERALIZED WEIERSTRASS PRODUCTS AND WIENER-HOPF FACTORIZATION P. PATIE AND M.

Least-Squares Regression on Sparse Spaces

The total derivative. Chapter Lagrangian and Eulerian approaches

3.6. Let s write out the sample space for this random experiment:

Lecture 10: October 30, 2017

arxiv: v4 [math.pr] 27 Jul 2016

Additive Lévy Processes

N may be the number of photon interactions in a given volume or the number of radioactive disintegrations in a given time. Expectation value:

Symbolic integration with respect to the Haar measure on the unitary groups

On conditional moments of high-dimensional random vectors given lower-dimensional projections

Topic 7: Convergence of Random Variables

Euler equations for multiple integrals

Sturm-Liouville Theory

Implicit Differentiation

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods

1. Aufgabenblatt zur Vorlesung Probability Theory

Lower Bounds for the Smoothed Number of Pareto optimal Solutions

Function Spaces. 1 Hilbert Spaces

Some functions and their derivatives

Introduction to Markov Processes

arxiv: v1 [math.pr] 4 Feb 2016

Monotonicity for excited random walk in high dimensions

SYMPLECTIC GEOMETRY: LECTURE 3

Generalization of the persistent random walk to dimensions greater than 1

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Notes on Lie Groups, Lie algebras, and the Exponentiation Map Mitchell Faulk

x 2 2x 8 (x 4)(x + 2)

Analysis IV, Assignment 4

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

7.1 Support Vector Machine

Asymptotic estimates on the time derivative of entropy on a Riemannian manifold

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

An Optimal Algorithm for Bandit and Zero-Order Convex Optimization with Two-Point Feedback

WUCHEN LI AND STANLEY OSHER

Noether s theorem applied to classical electrodynamics

On a limit theorem for non-stationary branching processes.

Convergence of random variables, and the Borel-Cantelli lemmas

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold

STABLE DISTRIBUTIONS

Relative Entropy and Score Function: New Information Estimation Relationships through Arbitrary Additive Perturbation

12.11 Laplace s Equation in Cylindrical and

Lower bounds on Locality Sensitive Hashing

Tutorial on Maximum Likelyhood Estimation: Parametric Density Estimation

Quantum Mechanics in Three Dimensions

1 Math 285 Homework Problem List for S2016

Nonlinear Schrödinger equation with a white-noise potential: Phase-space approach to spread and singularity

Agmon Kolmogorov Inequalities on l 2 (Z d )

1 Heisenberg Representation

High-Dimensional p-norms

Ramsey numbers of some bipartite graphs versus complete graphs

Sharp Thresholds. Zachary Hamaker. March 15, 2010

LECTURE NOTES ON DVORETZKY S THEOREM

Chapter 4. Electrostatics of Macroscopic Media

Problem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs

STUDENT S COMPANIONS IN BASIC MATH: THE FOURTH. Trigonometric Functions

Global Solutions to the Coupled Chemotaxis-Fluid Equations

u!i = a T u = 0. Then S satisfies

The Exact Form and General Integrating Factors

Math 342 Partial Differential Equations «Viktor Grigoryan

SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES

Generalized Tractability for Multivariate Problems

II. First variation of functionals

Chapter 6: Energy-Momentum Tensors

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.

cosh x sinh x So writing t = tan(x/2) we have 6.4 Integration using tan(x/2) = 2 2t 1 + t 2 cos x = 1 t2 We will revisit the double angle identities:

The Generalized Incompressible Navier-Stokes Equations in Besov Spaces

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS

Jointly continuous distributions and the multivariate Normal

Solution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010

Equilibrium in Queues Under Unknown Service Times and Service Value

arxiv: v1 [cs.it] 21 Aug 2017

cosh x sinh x So writing t = tan(x/2) we have 6.4 Integration using tan(x/2) 2t 1 + t 2 cos x = 1 t2 sin x =

Asymptotic determination of edge-bandwidth of multidimensional grids and Hamming graphs

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Lecture XVI: Symmetrical spacetimes

Lecture 6: Calculus. In Song Kim. September 7, 2011

Inverse Problems = Quest for Information

The Three-dimensional Schödinger Equation

Transcription:

Some Examples Our immeiate goal is to see some examples of Lévy processes, an/or infinitely-ivisible laws on. Uniform motion Choose an fix a nonranom an efine X := for all (1) Then, {X } is a [nonranom] Lévy process with Lévy triple ( ). The process {X } enotes uniform motion in the irection of. Poisson processes on the real line If N = Poiss(λ) for some λ>, then Ee ξn = = e ξ e λ λ! = exp λ 1 e ξ (2) That is, E exp(ξn) = exp( Ψ(ξ)), where Ψ(ξ) = 1 e ξ ξ1l (1) ( ) () an () := λδ {1} (). The corresponing Lévy process is calle the Poisson process with intensity parameter λ. 9

1 2. Some Examples Nonstanar Brownian motion with rift The Lévy triple ( I), where I enotes the ( ) ientity matrix, belongs to a vector of i.i.. stanar-normal ranom variables, an the corresponing Lévy process is [stanar] -imensional Brownian motion. We can generalize this example easily: Choose an fix a vector, an a ( ) matrix σ, an consier the Lévy triple (σ). The corresponing Lévy exponent is Ψ(ξ) =( ξ)+ 1 2 σξ2 Therefore, Ψ is the Lévy exponent of ranom vector X in if an only if X = + σz where Z is a vector of i.i.. stanar-normal ranom variables. The corresponing Lévy process is escribe by W := + σb, where B is stanar Brownian motion [check!]. The th coorinate of W is a Brownian motion with mean an variance 2 := (σ σ), an the coorinates of W are not in general inepenent. Since lim (W /) = a.s. by the law of large numbers for Brownian motion, is calle the rift of the nonstanar Brownian motion W. Isotropic stable laws Choose an fix a number α. An isotropic stable law of inex α is the infinitely-ivisible istribution with Lévy exponent Ψ(ξ) =ξ α, where is a fixe constant. The corresponing Lévy process is calle the isotropic stable process with inex α. We consier only ranom vectors with Lévy exponent exp( ξ α ) in this iscussion. Of course, =leas to Ψ, which is the exponent of the infinitely ivisible, but trivial, ranom variable X. Therefore, we stuy only =. Also, we nee exp( Ψ(ξ)) = exp{ ξ α } 1, an this means that cannot be negative. Finally, Ψ() =, an hence we have α>. Lemma 1 (Lévy). Ψ(ξ) =ξ α is a Lévy exponent iff α ( 2]. The case α =2is the Gaussian case we just saw. An α =1is also noteworthy; the resulting istribution is the isotropic [or symmetric, in imension one] Cauchy istribution whose ensity is () = Γ +1 (+1)/2 2 1+ 2 π (+1)/2 /2 2 for all (3) Proof of Lemma 1. Exercise 8 below shows that Ψ(ξ) = O(ξ 2 ) as ξ ; therefore, α 2. Since α =2is Gaussian, we limit ourselves to α<2.

Isotropic stable laws 11 Next let us consier α ( 2) [an of course >]. In that case, a change of variables shows that for all ξ, 1 e ξ (ξ )1l (1) () +α = (1 cos(ξ )) +α (4) ξ α [The first ientity is justifie because the left-most integral is a raial function of ξ, an hence real.] Therefore, we can choose C ( ) so that () =C (+α) is the Lévy measure with exponent Ψ(ξ) = exp( ξ α ) iff (3) on page 3 hols. Sometimes, a reasonable knowlege of the Lévy exponent of an infinitelyivisible law yiels insight into its structure. Here is a first example; Skoroho (1961) contains much more precise [an very useful] estimates of the tails of stable laws. Proposition 2. If X has an isotropic stable istribution with α ( 2), then for all β>, E(X β ) < iff β<α. In other wors, except in the case that α =2, the ecay of the tail of an isotropic stable law is slow. It is also possible to say more about the tails of the istrbution (Theorem 14, page 5). But that requires a more sophisticate analysis. Proof. Because E exp( X) = exp( α )=Ecos( X) for some >, we may apply (4) in the following form: 1 e ξ +α = (1 cos(ξ )) +α ξα (5) eplace ξ by X an take expectations to obtain E X β E (1 cos( X)) +β = Now 1 1 e α β+1 1 β+1 < 1 e α 1 e α β+1 for every β> +β Therefore, E(X β ) < iff 1 (1 e α )/ β+1 <. The result follows from this an the fact that (θ/2) 1 e θ θ for all θ ( 1).

12 2. Some Examples The asymmetric Cauchy istribution on the line We can specialize the preceing example to = α =1an obtain the symmetric Cauchy law µ on the line. Of course the ensity of µ is known as well. 1 But more significantly, we have learne that the Lévy triple of µ is ( ), where () 2. This suggests that perhaps we can create an asymmetric variation of the Cauchy law by consiering a Lévy triple of the form ( ) for a Lévy measure of the form () = 1 2 1l ()()+ 2 2 1l ( )() where 1 = 2 are both positive an is selecte carefully. This can in fact be one, as we will see next. Theorem 3. For every 1 2 > an θ [ 2/π 2/π] there exists an an infinitely-ivisible Borel probability measure µ on such that: (i) The Lévy triple of µ has the form ( ) for as above; an (ii) ˆµ(ξ) = exp( ξ θξ log ξ ) for ξ \{}. 2 The measure µ given above is calle the Cauchy istribution on with asymmetry parameter θ [an scale parameter ]. When θ =, µ is the symmetric Cauchy law. When θ =2/π, µ is calle the completely asymmetric Cauchy law on. The Gamma istribution on the half line It is easy to see that the Gamma (αλ) istribution on ( ) is infinitely ivisible for every α λ >. Next we ientify its Lévy triple. Proposition 4. If µ is a Gamma (αλ) istribution on + for some α λ >, then µ is infinitely ivisible with :=, σ :=, an Lévy measure an Lévy exponent Ψ respectively given by () = αe λ 1l () () Ψ(ξ) =α log 1 ξ λ where log enotes the principle branch of the logarithm. 1 It is () =1/{π(1 + 2 )} for <<. 2 I am making a little fuss about ξ \{}, since the function ξ log ξ is efine only for ξ \{}. But of course ˆµ() = 1 because µ is a probability measure. Alternatively, we can efine the function ξ log ξ cotinuously on all of [ ) by letting log be lim ξ ξ= ξ log ξ =. If so, then the state formula for ˆµ(ξ) is vali for all ξ.

Problems for Lecture 2 13 Aing inepenent Lévy processes Finally, let me mention the following evice which can be use to generate new Lévy processes from ol ones. Lemma 5. If {X } an { X } are inepenent Lévy processes on with respective triples ( σ ) an ( σ ), then {X + Y } is a Lévy process on with Lévy triple (AΣ M), where A := + an M := +, an Σ can be chosen in any fashion as long as Σ Σ=σ σ + σ σ. This is elementary, an can be checke by irectly verifying the efining properties of Lévy processes. However, I emphasize that: (i) There is always such a Σ; 3 an (ii) even though Σ might not be efine uniquely, Σ Σ is Problems for Lecture 2 1. Prove that if := 1 2 ( sin ) 1 2 sin, then for all ξ>, 1 e ξ + ξ1l (1) () 2 = πξ + ξ ln ξ + ξ 2 Deuce Theorem 3 from this ientity. 2. This problem outlines the proof of Proposition 4. (1) Suppose : ( ) has a continuous erivative L 1 (), an (+) := lim () an ( ) := lim () exist an are finite. Prove that for all λ ρ >, (λ) ((λ + ρ)) = (+) ( ) ln 1+ ρ λ This is the Frullani integral ientity. (2) Prove that for all λ> an ξ, 1 ξ λ = exp e λ 1 e ξ an euce Proposition 4 from this ientity. 3 (Stable scaling). Prove that if X is an isotropic stable process in with inex α ( 2], then Y := 1/α X efines a Lévy process for every >fixe. Explicitly compute the Lévy exponent of Y. Is the same result true when epens on? What happens if you consier instea an asymmetric Cauchy process on the line? 3 This follows reaily from the fact that σ σ + σ σ is a nonnegative-efinite matrix.