On optimal stopping of autoregressive sequences

Similar documents
On Optimal Stopping Problems with Power Function of Lévy Processes

Stochastic Process (ENPC) Monday, 22nd of January 2018 (2h30)

Markov processes and queueing networks

Non-homogeneous random walks on a semi-infinite strip

Ernesto Mordecki 1. Lecture III. PASI - Guanajuato - June 2010

Harmonic Functions and Brownian motion

Nonparametric regression with martingale increment errors

Statistical Models and Algorithms for Real-Time Anomaly Detection Using Multi-Modal Data

1.3 Convergence of Regular Markov Chains

Exercises with solutions (Set D)

Asymptotic properties of the maximum likelihood estimator for a ballistic random walk in a random environment

Lecture 10. Theorem 1.1 [Ergodicity and extremality] A probability measure µ on (Ω, F) is ergodic for T if and only if it is an extremal point in M.

Lecture 5. If we interpret the index n 0 as time, then a Markov chain simply requires that the future depends only on the present and not on the past.

Latent voter model on random regular graphs

Ergodic Theorems. Samy Tindel. Purdue University. Probability Theory 2 - MA 539. Taken from Probability: Theory and examples by R.

Reflected Brownian Motion

Lecturer: Olga Galinina

MS&E 321 Spring Stochastic Systems June 1, 2013 Prof. Peter W. Glynn Page 1 of 10

Analyzing Aggregated AR(1) Processes

Mixing time for a random walk on a ring

Chapter 3 sections. SKIP: 3.10 Markov Chains. SKIP: pages Chapter 3 - continued

Other properties of M M 1

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

Introduction General Framework Toy models Discrete Markov model Data Analysis Conclusion. The Micro-Price. Sasha Stoikov. Cornell University

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

Exact Simulation of the Stationary Distribution of M/G/c Queues

A Note on the Central Limit Theorem for a Class of Linear Systems 1

IEOR 6711, HMWK 5, Professor Sigman

Irreducibility. Irreducible. every state can be reached from every other state For any i,j, exist an m 0, such that. Absorbing state: p jj =1

Harmonic Functions and Brownian Motion in Several Dimensions

NEW FRONTIERS IN APPLIED PROBABILITY

Albert N. Shiryaev Steklov Mathematical Institute. On sharp maximal inequalities for stochastic processes

Exponential martingales: uniform integrability results and applications to point processes

8. Statistical Equilibrium and Classification of States: Discrete Time Markov Chains

Bayesian Quickest Change Detection Under Energy Constraints

On an Effective Solution of the Optimal Stopping Problem for Random Walks

Chapter 3 sections. SKIP: 3.10 Markov Chains. SKIP: pages Chapter 3 - continued

Lecture 6. 2 Recurrence/transience, harmonic functions and martingales

Modern Discrete Probability Branching processes

1 General problem. 2 Terminalogy. Estimation. Estimate θ. (Pick a plausible distribution from family. ) Or estimate τ = τ(θ).

Practical conditions on Markov chains for weak convergence of tail empirical processes

Lecture 7. 1 Notations. Tel Aviv University Spring 2011

Positive Harris Recurrence and Diffusion Scale Analysis of a Push Pull Queueing Network. Haifa Statistics Seminar May 5, 2008

CALCULATION OF RUIN PROBABILITIES FOR A DENSE CLASS OF HEAVY TAILED DISTRIBUTIONS

Single Index Quantile Regression for Heteroscedastic Data

The Modified Keifer-Weiss Problem, Revisited

On semilinear elliptic equations with measure data

2 Mathematical Model, Sequential Probability Ratio Test, Distortions

MULTIVARIATE DISCRETE PHASE-TYPE DISTRIBUTIONS

Problem Points S C O R E Total: 120

T 1. The value function v(x) is the expected net gain when using the optimal stopping time starting at state x:

Exercises. T 2T. e ita φ(t)dt.

A Conditional Approach to Modeling Multivariate Extremes

A.Piunovskiy. University of Liverpool Fluid Approximation to Controlled Markov. Chains with Local Transitions. A.Piunovskiy.

A Stochastic Online Sensor Scheduler for Remote State Estimation with Time-out Condition

A New Penalty-SQP Method

Lecture 17 Brownian motion as a Markov process

Generalized Fibonacci Numbers and Blackwell s Renewal Theorem

The Contact Process on the Complete Graph with Random, Vertex-dependent, Infection Rates

8.1 Concentration inequality for Gaussian random matrix (cont d)

Applied Stochastic Processes

X t = a t + r t, (7.1)

Plan Martingales cont d. 0. Questions for Exam 2. More Examples 3. Overview of Results. Reading: study Next Time: first exam

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539

Computer Intensive Methods in Mathematical Statistics

The Mabinogion Sheep Problem

FUNCTIONAL CENTRAL LIMIT (RWRE)

Stochastic Processes

Stochastic process. X, a series of random variables indexed by t

1 Math 285 Homework Problem List for S2016

process on the hierarchical group

Mixing Times and Hitting Times

An Introduction to Probability Theory and Its Applications

Uniformly Uniformly-ergodic Markov chains and BSDEs

Module 6:Random walks and related areas Lecture 24:Random woalk and other areas. The Lecture Contains: Random Walk.

Convergence Rate of Markov Chains

MATH 56A: STOCHASTIC PROCESSES CHAPTER 6

On the principle of smooth fit in optimal stopping problems

Random Walks Conditioned to Stay Positive

On detection of unit roots generalizing the classic Dickey-Fuller approach

Preliminary Exam: Probability 9:00am 2:00pm, Friday, January 6, 2012

Convex Optimization of Graph Laplacian Eigenvalues

2. Variance and Covariance: We will now derive some classic properties of variance and covariance. Assume real-valued random variables X and Y.

Maximum Likelihood Estimation

Markov Switching Models

arxiv: v2 [q-fin.cp] 9 Jun 2011

Analytical properties of scale functions for spectrally negative Lévy processes

Chapter 7. Markov chain background. 7.1 Finite state space

Probabilistic Logics and the Synthesis of Reliable Organisms from Unreliable Components

Stability of the two queue system

Lecture 12. F o s, (1.1) F t := s>t

Irregular Birth-Death process: stationarity and quasi-stationarity

Math 6810 (Probability) Fall Lecture notes

ABC methods for phase-type distributions with applications in insurance risk problems

MFM Practitioner Module: Quantitiative Risk Management. John Dodson. October 14, 2015

Perturbed Proximal Gradient Algorithm

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

The G/GI/N Queue in the Halfin-Whitt Regime I: Infinite Server Queue System Equations

Lectures on Stochastic Stability. Sergey FOSS. Heriot-Watt University. Lecture 4. Coupling and Harris Processes

for all f satisfying E[ f(x) ] <.

Transcription:

On optimal stopping of autoregressive sequences Sören Christensen (joint work with A. Irle, A. Novikov) Mathematisches Seminar, CAU Kiel September 24, 2010 Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 1 / 22

Outline 1 AR(1) processes 2 Basic results on optimal stopping 3 Threshold times as optimal stopping times 4 Phasetype distributions and overshoot 5 Continuous fit condition 6 Summary Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 1 / 22

Outline AR(1) processes 1 AR(1) processes 2 Basic results on optimal stopping 3 Threshold times as optimal stopping times 4 Phasetype distributions and overshoot 5 Continuous fit condition 6 Summary Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 2 / 22

AR(1) processes Setting AR(1) process X n = λx n 1 + Z n, X 0 = x, Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 3 / 22

AR(1) processes Setting AR(1) process X n = λx n 1 + Z n, X 0 = x, λ (0, 1), Z 1, Z 2,... iid. Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 3 / 22

AR(1) processes Setting AR(1) process X n = λx n 1 + Z n, X 0 = x, λ (0, 1), Z 1, Z 2,... iid. X n = (1 λ)x n 1 n + L n, L n = n i=1 Z i Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 3 / 22

Setting AR(1) processes AR(1) process X n = λx n 1 + Z n, X 0 = x, λ (0, 1), Z 1, Z 2,... iid. X n = (1 λ)x n 1 n + L n, X n = λ n X 0 + n 1 i=0 λi Z n i. L n = n i=1 Z i Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 3 / 22

AR(1) processes Path of an AR(1) process Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 4 / 22

AR(1) processes Threshold time, overshoot In many applications such as statistical surveillance it is of interest to find the joint distribution of τ a = inf{n N 0 : X n a} X τa a (threshold time) (overshoot). a X τa a τ a Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 5 / 22

Outline Basic results on optimal stopping 1 AR(1) processes 2 Basic results on optimal stopping 3 Threshold times as optimal stopping times 4 Phasetype distributions and overshoot 5 Continuous fit condition 6 Summary Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 6 / 22

Basic results on optimal stopping Optimal stopping of Markov processes Let (X n ) n N0 be a Markov process on E, g : E R, ρ (0, 1). Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 7 / 22

Basic results on optimal stopping Optimal stopping of Markov processes Let (X n ) n N0 be a Markov process on E, g : E R, ρ (0, 1). Find a stopping time τ that maximizes E x (ρ τ g(x τ )). Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 7 / 22

Basic results on optimal stopping Optimal stopping of Markov processes Let (X n ) n N0 be a Markov process on E, g : E R, ρ (0, 1). Find a stopping time τ that maximizes Solution : E x (ρ τ g(x τ )). v(x) := sup E x (ρ τ g(x τ )) τ S := {x E : v(x) g(x)} Optimal stopping time: value function. optimal stopping set. τ = inf{n N 0 : X n S}. Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 7 / 22

Outline Threshold times as optimal stopping times 1 AR(1) processes 2 Basic results on optimal stopping 3 Threshold times as optimal stopping times 4 Phasetype distributions and overshoot 5 Continuous fit condition 6 Summary Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 8 / 22

Threshold times as optimal stopping times Elementary approach to optimal stopping of AR(1) processes 1 Reduce the set of potential optimal stopping times to threshold time. Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 9 / 22

Threshold times as optimal stopping times Elementary approach to optimal stopping of AR(1) processes 1 Reduce the set of potential optimal stopping times to threshold time. 2 Find the joint distribution of (τ a, X τa a) to obtain an expression for E x (ρ τa g(x τa )). Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 9 / 22

Threshold times as optimal stopping times Elementary approach to optimal stopping of AR(1) processes 1 Reduce the set of potential optimal stopping times to threshold time. 2 Find the joint distribution of (τ a, X τa a) to obtain an expression for E x (ρ τa g(x τa )). 3 Find the optimal boundary a. Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 9 / 22

Threshold times as optimal stopping times Novikov-Shiryaev problem Let g(x) = (x + ) α, α > 0. Consider sup τ E x (ρ τ (X τ + ) α ) = sup E(ρ τ ((λ τ x + X τ ) + ) α ). τ Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 10 / 22

Threshold times as optimal stopping times Novikov-Shiryaev problem Let g(x) = (x + ) α, α > 0. Consider sup τ E x (ρ τ (X τ + ) α ) = sup E(ρ τ ((λ τ x + X τ ) + ) α ). τ Proposition If (X n ) n N0 is a AR(1) process, then there exists a such that the optimal stopping set S is given by S = [a, ) Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 10 / 22

Proof Threshold times as optimal stopping times For simplicity: Z 1 0. State-space: [0, ) Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 11 / 22

Proof Threshold times as optimal stopping times For simplicity: Z 1 0. State-space: [0, ) Then v(x) g(x) = sup τ E(ρ τ (λ τ x + X τ ) α ) x α = sup E(ρ τ (λ τ + X τ /x) α ) τ Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 11 / 22

Proof Threshold times as optimal stopping times For simplicity: Z 1 0. State-space: [0, ) Then v(x) g(x) = sup τ E(ρ τ (λ τ x + X τ ) α ) x α = sup E(ρ τ (λ τ + X τ /x) α ) τ Hence S = {x : v(x) g(x)} is of the form [a, ). Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 11 / 22

Outline Phasetype distributions and overshoot 1 AR(1) processes 2 Basic results on optimal stopping 3 Threshold times as optimal stopping times 4 Phasetype distributions and overshoot 5 Continuous fit condition 6 Summary Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 12 / 22

Motivation Phasetype distributions and overshoot If (X n ) n N0 is a random walk (or an AR(1) process) with Exp(β)-distributed jumps. Then τ a, X τa a are independent and X τa a d = Exp(β). Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 13 / 22

Phasetype distributions and overshoot Distributions of phasetype (J t ) t 0 : Markov chain in continuous time on {1,..., m} { }. 1,..., m: transient, absorbing. Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 14 / 22

Phasetype distributions and overshoot Distributions of phasetype (J t ) t 0 : Markov chain in continuous time on {1,..., m} { }. 1,..., m: transient, absorbing. ( ) Generator: ˆQ 1 Q q =, Q R 0 0 m m, q = Q... 1 ˆα = (α 1,..., α m, 0) initial distribution Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 14 / 22

Phasetype distributions and overshoot Distributions of phasetype (J t ) t 0 : Markov chain in continuous time on {1,..., m} { }. 1,..., m: transient, absorbing. ( ) Generator: ˆQ 1 Q q =, Q R 0 0 m m, q = Q... 1 ˆα = (α 1,..., α m, 0) initial distribution η := inf{t 0 : J t = } PH(Q, α) := P ηˆα is called phasetype distribution with parameters Q, α. Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 14 / 22

Phasetype distributions and overshoot Properties of phasetype distributions density: f (t) = αe Qt q Laplace transform: Eˆα (e sη ) = α( si Q) 1 q, s C with R(s) < δ. closed under mixture and convolution. Phasetype distributions are dense in all distributions on (0, ). Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 15 / 22

Phasetype distributions and overshoot Properties of phasetype distributions density: f (t) = αe Qt q Laplace transform: Eˆα (e sη ) = α( si Q) 1 q, s C with R(s) < δ. closed under mixture and convolution. Phasetype distributions are dense in all distributions on (0, ). Ex: m = 1, Q = ( λ) PH(Q, (1)) := Exp(λ). Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 15 / 22

Phasetype distributions and overshoot Phasetype distribution and overshoot Generalization: If (X n ) n N0 is AR(1) process with innovations Z n = S n T n, S n, T n 0 independent, S n PH(Q, α), then E x (ρ τa g(x τa )) = m E x (ρ τa 1 Ki )E(g(a + R i )), i=1 R i PH(Q, e i ), K i events. Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 16 / 22

Phasetype distributions and overshoot Phasetype distribution and overshoot Generalization: If (X n ) n N0 is AR(1) process with innovations Z n = S n T n, S n, T n 0 independent, S n PH(Q, α), then E x (ρ τa g(x τa )) = m E x (ρ τa 1 Ki )E(g(a + R i )), i=1 R i PH(Q, e i ), K i events. How to find E x (ρ τa 1 Ki )? Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 16 / 22

Phasetype distributions and overshoot Phasetype distribution and overshoot Generalization: If (X n ) n N0 is AR(1) process with innovations Z n = S n T n, S n, T n 0 independent, S n PH(Q, α), then E x (ρ τa g(x τa )) = m E x (ρ τa 1 Ki )E(g(a + R i )), i=1 R i PH(Q, e i ), K i events. How to find E x (ρ τa 1 Ki )? Answer: Use the stationary distribution, complex martingales and complex analysis. Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 16 / 22

Outline Continuous fit condition 1 AR(1) processes 2 Basic results on optimal stopping 3 Threshold times as optimal stopping times 4 Phasetype distributions and overshoot 5 Continuous fit condition 6 Summary Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 17 / 22

Continuous fit condition How to find the find the optimal threshold a? For each a R we consider v a (x) = E x (ρ τa g(x τa )). Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 18 / 22

Continuous fit condition How to find the find the optimal threshold a? For each a R we consider v a (x) = E x (ρ τa g(x τa )). Choose a such that v a (a ) = g(a ). Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 18 / 22

Continuous fit Continuous fit condition a 0.7 is the unique solution to the transcendental equation g(a ) = v a (a ). Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 19 / 22

Continuous fit condition Smooth and continuous fit condition Principle for continuous time processes Let g be differentiable, then the following can be expected: 1 If the process X enters the interior of the optimal stopping set S immediately after starting on a S, then the value function v is smooth at a. Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 20 / 22

Continuous fit condition Smooth and continuous fit condition Principle for continuous time processes Let g be differentiable, then the following can be expected: 1 If the process X enters the interior of the optimal stopping set S immediately after starting on a S, then the value function v is smooth at a. 2 If the process X doesn t enter the interior of the optimal stopping set S immediately after starting on a S, then the value function v is continuous at a Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 20 / 22

Continuous fit condition Smooth and continuous fit condition Principle for continuous time processes Let g be differentiable, then the following can be expected: 1 If the process X enters the interior of the optimal stopping set S immediately after starting on a S, then the value function v is smooth at a. 2 If the process X doesn t enter the interior of the optimal stopping set S immediately after starting on a S, then the value function v is continuous at a The second principle holds in our situation. Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 20 / 22

Outline Summary 1 AR(1) processes 2 Basic results on optimal stopping 3 Threshold times as optimal stopping times 4 Phasetype distributions and overshoot 5 Continuous fit condition 6 Summary Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 21 / 22

Summary Summary We studied optimal stopping problems driven by autoregressive processes. Elementary arguments reduces the problem to finding an optimal threshold time. The joint distribution of (τ a, X τa a) can be obtained for phasetype innovations. The optimal threshold can be found using the continuous fit principle. Sören Christensen (CAU Kiel) Optimal stopping of AR(1) sequences 20.09.10 22 / 22