Large deviations for random walks under subexponentiality: the big-jump domain

Similar documents
Rare Events in Random Walks and Queueing Networks in the Presence of Heavy-Tailed Distributions

Introduction to Rare Event Simulation

Approximating the Integrated Tail Distribution

Reinsurance and ruin problem: asymptotics in the case of heavy-tailed claims

Reduced-load equivalence for queues with Gaussian input

The largest eigenvalues of the sample covariance matrix. in the heavy-tail case

On lower limits and equivalences for distribution tails of randomly stopped sums 1

Precise Large Deviations for Sums of Negatively Dependent Random Variables with Common Long-Tailed Distributions

Asymptotics and Simulation of Heavy-Tailed Processes

Rare event simulation for the ruin problem with investments via importance sampling and duality

The Fundamentals of Heavy Tails Properties, Emergence, & Identification. Jayakrishnan Nair, Adam Wierman, Bert Zwart

Asymptotic Statistics-III. Changliang Zou

Stability and Rare Events in Stochastic Models Sergey Foss Heriot-Watt University, Edinburgh and Institute of Mathematics, Novosibirsk

Corrections to the Central Limit Theorem for Heavy-tailed Probability Densities

On the inefficiency of state-independent importance sampling in the presence of heavy tails

8.1 Concentration inequality for Gaussian random matrix (cont d)

arxiv: v2 [math.pr] 22 Apr 2014 Dept. of Mathematics & Computer Science Mathematical Institute

The Convergence Rate for the Normal Approximation of Extreme Sums

Lecture 2. We now introduce some fundamental tools in martingale theory, which are useful in controlling the fluctuation of martingales.

On the estimation of the heavy tail exponent in time series using the max spectrum. Stilian A. Stoev

IEOR 6711: Stochastic Models I SOLUTIONS to the First Midterm Exam, October 7, 2008

Research Reports on Mathematical and Computing Sciences

March 1, Florida State University. Concentration Inequalities: Martingale. Approach and Entropy Method. Lizhe Sun and Boning Yang.

Limit theorems for dependent regularly varying functions of Markov chains

Asymptotics of random sums of heavy-tailed negatively dependent random variables with applications

Poisson Processes. Stochastic Processes. Feb UC3M

Other properties of M M 1

The function graphed below is continuous everywhere. The function graphed below is NOT continuous everywhere, it is discontinuous at x 2 and

arxiv: v1 [math.pr] 28 Nov 2007

The Moment Method; Convex Duality; and Large/Medium/Small Deviations

LIMITS FOR QUEUES AS THE WAITING ROOM GROWS. Bell Communications Research AT&T Bell Laboratories Red Bank, NJ Murray Hill, NJ 07974

Approximation of Heavy-tailed distributions via infinite dimensional phase type distributions

On Upper Bounds for the Tail Distribution of Geometric Sums of Subexponential Random Variables

The Subexponential Product Convolution of Two Weibull-type Distributions

The exact asymptotics for hitting probability of a remote orthant by a multivariate Lévy process: the Cramér case

On Sums of Conditionally Independent Subexponential Random Variables

Moments of the maximum of the Gaussian random walk

Rare-Event Simulation

Limiting Distributions

E cient Simulation and Conditional Functional Limit Theorems for Ruinous Heavy-tailed Random Walks

f X (x) = λe λx, , x 0, k 0, λ > 0 Γ (k) f X (u)f X (z u)du

Estimating Tail Probabilities of Heavy Tailed Distributions with Asymptotically Zero Relative Error

Asymptotic behavior for sums of non-identically distributed random variables

EFFICIENT SIMULATION FOR LARGE DEVIATION PROBABILITIES OF SUMS OF HEAVY-TAILED INCREMENTS

arxiv:math/ v2 [math.pr] 9 Oct 2007

Stability of the Defect Renewal Volterra Integral Equations

The Dynamic Analysis and Design of A Communication link with Stationary and Nonstationary Arrivals

Lecture 20: Reversible Processes and Queues

Exercises in Extreme value theory

Spring 2012 Math 541A Exam 1. X i, S 2 = 1 n. n 1. X i I(X i < c), T n =

Analysis methods of heavy-tailed data

arxiv: v1 [math.pr] 9 May 2014

Math 328 Course Notes

LARGE DEVIATION PROBABILITIES FOR SUMS OF HEAVY-TAILED DEPENDENT RANDOM VECTORS*

SOLUTIONS IEOR 3106: Second Midterm Exam, Chapters 5-6, November 8, 2012

Probability and Measure

Efficient rare-event simulation for sums of dependent random varia

Finite-time Ruin Probability of Renewal Model with Risky Investment and Subexponential Claims

Some open problems related to stability. 1 Multi-server queue with First-Come-First-Served discipline

UNIFORM LARGE DEVIATIONS FOR HEAVY-TAILED QUEUES UNDER HEAVY TRAFFIC

Free Entropy for Free Gibbs Laws Given by Convex Potentials

Stochastic-Process Limits

Analysis of an M/G/1 queue with customer impatience and an adaptive arrival process

Sub-Gaussian estimators under heavy tails

Limiting Distributions

Lecture 2: Convergence of Random Variables

STATE-DEPENDENT IMPORTANCE SAMPLING FOR REGULARLY VARYING RANDOM WALKS

Class 11 Non-Parametric Models of a Service System; GI/GI/1, GI/GI/n: Exact & Approximate Analysis.

FRACTIONAL BROWNIAN MOTION WITH H < 1/2 AS A LIMIT OF SCHEDULED TRAFFIC

Chapter 2. Poisson Processes. Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan

Control of Fork-Join Networks in Heavy-Traffic

A Diffusion Approximation for Stationary Distribution of Many-Server Queueing System In Halfin-Whitt Regime

Phenomena in high dimensions in geometric analysis, random matrices, and computational geometry Roscoff, France, June 25-29, 2012

THIELE CENTRE. Markov Dependence in Renewal Equations and Random Sums with Heavy Tails. Søren Asmussen and Julie Thøgersen

Chapter 10. Hypothesis Testing (I)

IEOR 8100: Topics in OR: Asymptotic Methods in Queueing Theory. Fall 2009, Professor Whitt. Class Lecture Notes: Wednesday, September 9.

Logarithmic scaling of planar random walk s local times

Randomly Weighted Sums of Conditionnally Dependent Random Variables

Ruin Probabilities of a Discrete-time Multi-risk Model

Large deviations of empirical processes

Uniformly Efficient Importance Sampling for the Tail Distribution of Sums of Random Variables

POISSON PROCESSES 1. THE LAW OF SMALL NUMBERS

Notes 9 : Infinitely divisible and stable laws

Lecture Notes 3 Convergence (Chapter 5)

Additive functionals of infinite-variance moving averages. Wei Biao Wu The University of Chicago TECHNICAL REPORT NO. 535

This lecture is expanded from:

University Of Calgary Department of Mathematics and Statistics

Introduction to Algorithmic Trading Strategies Lecture 10

Branching Processes II: Convergence of critical branching to Feller s CSB

Quantile-quantile plots and the method of peaksover-threshold

Self-normalized Cramér-Type Large Deviations for Independent Random Variables

18.175: Lecture 17 Poisson random variables

TAIL ASYMPTOTICS FOR A RANDOM SIGN LINDLEY RECURSION

Random matrices: Distribution of the least singular value (via Property Testing)

TOWARDS BETTER MULTI-CLASS PARAMETRIC-DECOMPOSITION APPROXIMATIONS FOR OPEN QUEUEING NETWORKS

Generalization theory

Designing load balancing and admission control policies: lessons from NDS regime

18.175: Lecture 13 Infinite divisibility and Lévy processes

BRAVO for QED Queues

Convexity Properties of Loss and Overflow Functions

Transcription:

Large deviations under subexponentiality p. Large deviations for random walks under subexponentiality: the big-jump domain Ton Dieker, IBM Watson Research Center joint work with D. Denisov (Heriot-Watt, UK) and V. Shneer (EURANDOM, The Netherlands)

Large deviations under subexponentiality p. Random walks and LD theory Let {S n = ξ 1 +... + ξ n } be a random walk, with i.i.d. step sizes ξ 1, ξ 2,... having distribution F. We often assume Eξ = 0. We are interested in asymptotics of P(S n > x), for large n and x. Light-tailed case: Bahadur-Ranga Rao asymptotics for P(S n > an) Logarithmic asymptotics on abstract spaces (Cramér s theorem), assuming the existence of exponential moments e.g. books by Varadhan, Dembo/Zeitouni, Deuschel/Stroock, etc. If Cramér s condition fails to hold, then we are in the so-called heavy-tailed framework.

Large deviations under subexponentiality p. Subexponentiality The subexponential distributions form a widely used family of heavy-tailed distributions F is subexponential if, with F(x) = 1 F(x), for all y R, lim x F(x + y) F(x) = 1 P(S 2 > x) 2F(x) as x Intuition: either ξ 1 is very large and ξ 2 moderate, or vice versa For subexponential F, P(S n > x) nf(x) as x for any n 1 single big jump principle

Large deviations under subexponentiality p. Subexponentiality: applications Internet traffic modeling; e.g., session durations Long sessions can explain long-term correlations Some models make self-similarity plausible by superposition Insurance mathematics; e.g., claim sizes Theory of ruin for insurance companies One big claim causes ruin Call durations in certain call centers

Large deviations under subexponentiality p. Motivation from queueing theory The step size ξ i is the difference between the service time of the i-th customer and its inter-arrival time (so Eξ i < 0) Let τ be the length of the busy period Baltrunas/Daley/Klüppelberg, Denisov/Shneer: For heavy-tailed service distributions, we have for some explicit constant C P(τ > n) C n P(S n > 0)

Large deviations under subexponentiality p. Subexponentiality (3) Recall: P(S n > x) nf(x) as x Main question: for which {x n } do we have P(S n > x n ) nf(x n )? We say that {x n } is a big-jump sequence. This is a LD theory for heavy tails, but we ll see an interesting connection with LD theory for light-tailed distributions at the end

Large deviations under subexponentiality p. Two important examples F with F(x) x α, α > 0, is subexponential (more generally, F(x) = l(x)x α for some slowly varying l) a finite mean/variance is not required if α > 2, known that for any ɛ > 0, x n = (1 + ɛ) (α 2)n log n is a big-jump sequence (A. Nagaev) F with F(x) e xβ, β (0, 1), is subexponential known that any {x n } with x n n 1 2(1 β) is a big-jump sequence. in particular x n = an is big-jump sequence only if β < 1/2

Large deviations under subexponentiality p. Two important examples F with F(x) x α, α > 0, is subexponential (more generally, F(x) = l(x)x α for some slowly varying l) a finite mean/variance is not required if α > 2, known that for any ɛ > 0, x n = (1 + ɛ) (α 2)n log n is a big-jump sequence (A. Nagaev) F with F(x) e xβ, β (0, 1), is subexponential known that any {x n } with x n n 1 2(1 β) is a big-jump sequence. in particular x n = an is big-jump sequence only if β < 1/2 Are these special cases of a general theory?

Large deviations under subexponentiality p. Related work A. V. and S. V. Nagaev, A. A. Borovkov, R. Doney, C. C. Heyde, I. Pinelis, L. V. Rozovskii Post-2003 contributions: Baltrunas/Daley/Klüppelberg, Borovkov/Mogulskii, Hult/Lindskog/Mikosch/Samorodnitsky, Jelenkovic/Momcilovic, Konstantinides/Mikosch, Ng/Tang/Yan/Yang, Tang

Large deviations under subexponentiality p. Four sequences {b n } is a natural-scale sequence if {S n /b n } is tight

Large deviations under subexponentiality p. Four sequences {b n } is a natural-scale sequence if {S n /b n } is tight {I n } is an insentitivity sequence if I n b n and lim sup n x I n F(x b n ) F(x) 1 = 0

Large deviations under subexponentiality p. Four sequences {b n } is a natural-scale sequence if {S n /b n } is tight {I n } is an insentitivity sequence if I n b n and lim sup n x I n {h n } is a truncation sequence if F(x b n ) F(x) 1 = 0 lim sup np(s 2 > x, ξ 1 > h n, ξ 2 > h n ) n x h n F(x) = 0

Large deviations under subexponentiality p. Four sequences {b n } is a natural-scale sequence if {S n /b n } is tight {I n } is an insentitivity sequence if I n b n and lim sup n x I n {h n } is a truncation sequence if F(x b n ) F(x) 1 = 0 lim sup np(s 2 > x, ξ 1 > h n, ξ 2 > h n ) n x h n F(x) = 0 {J n } is an h-small steps sequence if lim sup P(S n > x, ξ 1 h n,..., ξ n h n ) n x J n nf(x) = 0

Main result A lemma shows that these sequences exist if and only if F is subexponential

Main result A lemma shows that these sequences exist if and only if F is subexponential Theorem If h n = O(b n ) and h n J n, then P(S n > I n + J n ) np(s 1 > I n + J n ). It remains to find good sequences {I n } and {J n } otherwise result still valid, but weak statement

Choosing {h n } Recall: {h n } is a truncation sequence if lim sup np(s 2 > x, ξ 1 > h n, ξ 2 > h n ) n x h n F(x) = 0 If F(x) x α, then any {h n } satisfying nf(h n ) 0 is a truncation sequence If x r F(x) is subexponential for some r > 0, then any {h n } with lim sup n nh r n < is a truncation sequence

Choosing {J n }: a heuristic Recall: {J n } is an h-small steps sequence if lim sup P(S n > x, ξ 1 h n,..., ξ n h n ) n x J n nf(x) = 0 Suppose Eξ 2 <, and let J n satisfy J n 2n log[nf(j n)]. For any ɛ > 0, J n = (1 + ɛ)j n is a good choice, provided {n/j n } is a truncation sequence

The polynomial-tail example Here F(x) x α, for some α > 2. Suppose that Eξ = 0, Eξ 2 = 1. By the CLT we may set b n = n

The polynomial-tail example Here F(x) x α, for some α > 2. Suppose that Eξ = 0, Eξ 2 = 1. By the CLT we may set b n = n Any I n b n suffices: lim sup n x I n (x b n ) α x α 1 = 0

The polynomial-tail example Here F(x) x α, for some α > 2. Suppose that Eξ = 0, Eξ 2 = 1. By the CLT we may set b n = n Any I n b n suffices: lim sup n x I n (x b n ) α x α 1 = 0 Any {h n } with nf(h n ) 0 is a truncation sequence

The polynomial-tail example Here F(x) x α, for some α > 2. Suppose that Eξ = 0, Eξ 2 = 1. By the CLT we may set b n = n Any I n b n suffices: lim sup n x I n (x b n ) α x α 1 = 0 Any {h n } with nf(h n ) 0 is a truncation sequence Use the heuristic for {J n }: set (α 2)n log n and note that 2n log[nf(j n)] = 2αn log[n 1/α J n] (α 2)n log n = J n 2αn log[n 1/2 1/α ]

The Weibull example Here F(x) e xβ, for some β (0, 1). By the CLT we may set b n = n

The Weibull example Here F(x) e xβ, for some β (0, 1). By the CLT we may set b n = n To find I n we study sup x I n e (x n) β +x β 1 = 0. Since (x n) β + x β β nx β 1, we need I n n 1 2(1 β)

The Weibull example Here F(x) e xβ, for some β (0, 1). By the CLT we may set b n = n To find I n we study sup x I n e (x n) β +x β 1 = 0. Since (x n) β + x β β nx β 1, we need I n n 1 2(1 β) For any r > 0, x r F(x) is the tail of a subexponential distribution, so any {h n } with h n n 1/r is a truncation sequence

The Weibull example Here F(x) e xβ, for some β (0, 1). By the CLT we may set b n = n To find I n we study sup x I n e (x n) β +x β 1 = 0. Since (x n) β + x β β nx β 1, we need I n n 1 2(1 β) For any r > 0, x r F(x) is the tail of a subexponential distribution, so any {h n } with h n n 1/r is a truncation sequence The heuristic can be used for {J n }, yielding {J n } with J n I n

Local analogues We also study for T (0, ) P(S n (x n, x n + T]) np(s 1 (x n, x n + T]) Statements very similar, parts of the proof a lot harder!!! Use in the context of light tails: given γ > 0, subexponential F with L(γ) = e γy F(dy) <, define the RW on P through P (S 1 dx) = e γx F(dx)/L(γ). The local case allows us to conclude that P (S n > x n ) nl(γ) 1 n P (S 1 > x n )

Further examples We worked out examples with infinite variance or mean We always recovered the sharpest known big-jump sequences... but we also found big-jump sequences in new examples