Threshold Regression for Survival Analysis: Modeling Event Times by a Stochastic Process Reaching a Boundary

Total Page:16

File Type:pdf, Size:1020Kb

Threshold Regression for Survival Analysis: Modeling Event Times by a Stochastic Process Reaching a Boundary Statistical Science 2006, Vol. 21, No. 4, 501–513 DOI: 10.1214/088342306000000330 c Institute of Mathematical Statistics, 2006 Threshold Regression for Survival Analysis: Modeling Event Times by a Stochastic Process Reaching a Boundary Mei-Ling Ting Lee and G. A. Whitmore Abstract. Many researchers have investigated first hitting times as models for survival data. First hitting times arise naturally in many types of stochastic processes, ranging from Wiener processes to Markov chains. In a survival context, the state of the underlying process repre- sents the strength of an item or the health of an individual. The item fails or the individual experiences a clinical endpoint when the process reaches an adverse threshold state for the first time. The time scale can be calendar time or some other operational measure of degradation or disease progression. In many applications, the process is latent (i.e., unobservable). Threshold regression refers to first-hitting-time models with regression structures that accommodate covariate data. The pa- rameters of the process, threshold state and time scale may depend on the covariates. This paper reviews aspects of this topic and discusses fruitful avenues for future research. Key words and phrases: Accelerated testing, calendar time, compet- ing risk, cure rate, duration, environmental studies, first hitting time, gamma process, lifetime, latent variable models, maximum likelihood, operational time, occupational exposure, Ornstein–Uhlenbeck process, Poisson process, running time, stochastic process, stopping time, sur- vival analysis, threshold regression, time-to-event, Wiener diffusion pro- cess. 1. INTRODUCTION a stochastic process, which may be latent or observ- arXiv:0708.0346v1 [stat.ME] 2 Aug 2007 able. FHT models have a long history of application Many types of lifetime, duration or time-to-event in diverse fields, including medicine, environmental data may be interpreted as first hitting times (FHT’s) science, engineering, business, economics and soci- of a boundary or threshold state by sample paths of ology. FHT models may describe the length of a hospital stay, the survival time of a transplant pa- Mei-Ling Ting Lee is Professor and Chair, Division of tient, the onset time for a cancer induced by occu- Biostatistics, School of Public Health, Ohio State University, Columbus, Ohio 43210, USA e-mail: pational exposure, the failure time of an engineering [email protected]. G. A. Whitmore is Professor, system, the depletion time of an inventory, the sur- Management Science, Desautels Faculty of vival time of a new business, the transition time for Management, McGill University, Montreal, Quebec, a stock price change and the length of a marriage. Canada H3A 165 e-mail: [email protected]. Relevant articles, covering theory and application, appear in both the lifetime data and reliability lit- This is an electronic reprint of the original article published by the Institute of Mathematical Statistics in eratures. These FHT models are gradually becom- Statistical Science, 2006, Vol. 21, No. 4, 501–513. This ing widely adopted because of their conceptual ap- reprint differs from the original in pagination and peal, realism and applicability. Recently, interest in typographic detail. them has deepened and spread, and exciting new 1 2 M.-L. T. LEE AND G. A. WHITMORE areas of application are being encountered. The po- for the process and, therefore, the FHT is usually tential applications require new conceptual view- a stopping time in the formal sense of that term in points, theoretical advances, analytical techniques stochastic process theory. Note that when the parent and methodological extensions to which the discus- process is latent, there is no direct way of observing sion returns later. the FHT event in the state space of the process. To make FHT models truly valuable in applica- In some versions of the FHT model, there is no tions, they must be capable of extension to include guarantee that process X(t) will reach the bound- regression structures. Regression structures allow the ary set , so P (S< ){< 1.} We will let S = de- effects of covariates to explain the inherent disper- note theB absence of a∞ finite hitting time with P∞(S = sion of the data, thereby taking account of vari- ) = 1 P (S< ). The later discussion will show ability and sharpening inferences. Regression struc- situations∞ − where∞ this condition is plausible and a tures also provide scientific insights into potential desirable model feature. The basic FHT model (1) causal roles of covariates in the underlying processes, assumes that is fixed in time. In some applica- boundary sets and time scales. This article reviews tions, however,B it varies with time, that is, (t). This aspects of FHT models and is concerned especially variation may be deterministic or follow aB stochastic with regression structures for FHT models, which process. will be referred to as threshold regression, or TR An exhaustive review of the first-hitting-time lit- for short. The word “threshold” refers to the fact erature is impossible within the confines of a single that the FHT is triggered by the underlying process article. Eaton and Whitmore (1977) discuss FHT’s reaching a threshold state within a boundary set, as as a general model for hospital stay. Aalen and described in more detail in the next section. Gjessing (2001) provide an excellent overview of much of this subject. Likewise, Lawless (2003) gives a com- 2. THE FIRST-HITTING-TIME (FHT) MODEL plete and compact summary of theory, models and methods (see Section 11.5, pages 518–523). Lee and A first-hitting-time (FHT) model has two basic Whitmore (2004) also provide an overview of first- components: (1) a parent stochastic process X(t),t hitting-time models for survival and time-to-event ,x with initial value X(0) = x , where{ is∈ 0 data. We will make numerous references to selected theT time∈ X } space and is the state space of theT pro- work as we proceed. There is a huge literature deal- cess and (2) a boundaryX set , where . We ing with theoretical and mathematical aspects of shall refer to the boundary setB as aB⊂X boundary, FHT models that we will not attempt to review or barrier or threshold, depending onB which is most de- incorporate. We also will not cover random growth scriptive or conventional in the context. The process curve models, such as those of Carey and Koenig X(t) may have a variety of properties, such as one (1991) and Lu and Meeker (1993), which have an or{ many} dimensions, the Markov property, continu- FHT interpretation but where the only randomness ous or discrete states, and monotonic sample paths. at the level of the individual parent process is con- Whether the sample path of the parent process is fined to a noise factor. The large literature on linear observable or latent (unobservable) is an important and nonlinear regression methods for survival data distinguishing characteristic of the FHT model. La- and reliability where the underlying models have no tent processes are the most common by far. The central FHT motivation (such as accelerated failure boundary set may also have different features as time and proportional hazards models) is also not will be illustratedB in later examples. covered. Taking the initial value X(0) = x0 of the process to lie outside the boundary set , the first hitting time of is the random variableBS defined as 3. EXAMPLES OF FIRST-HITTING-TIME B MODELS (1) S = inf t : X(t) . { ∈ B} The parent stochastic processes may take many Thus, the first hitting time is the time when the forms, from Wiener processes to Markov chains. Like- stochastic process first encounters set . We refer wise, the nature of the boundary state may vary to the state first encountered in the boundaryB set widely—for example, a fixed threshold in a Wiener by the process, that is, X(S) , as the threshold process or an absorbing state in a Markov chain. state. The boundary set defines∈B a stopping condition As the preceding description of a first-hitting-time THRESHOLD REGRESSION MODELS 3 model is quite abstract, we now list a few basic ex- (1995), Whitmore and Schenkelberg (1997) and Hor- amples to illustrate the variety encountered in ap- rocks and Thompson (2004), to name a few. Onar plications. and Padgett (2000) and Padgett and Tomlinson (2004) 1. Bernoulli process and negative binomial first extend the Wiener diffusion model to an acceler- hitting time. The number of trials S required to ated testing context. Pettit and Young (1999) set reach the mth success in a Bernoulli process Bt,t = the model in a Bayesian context. As illustrated later, 1, 2,... has a negative binomial distribution{ with the inverse Gaussian distribution has a closed-form parameters} m and p, where p is the success prob- probability density function (p.d.f.) and a compu- ability on each trial. To give this setup our stan- tationally simple cumulative distribution function dard representation, we consider a parent process (c.d.f.). Their formulas vary slightly depending on Xt,t = 0, 1, 2,... with initial value X = x = m whether the parent process is defined as rising or { } 0 0 and let Xt = x0 Bt,t = 1, 2,..., where Bt is the falling to hit the relevant boundary. preceding Bernoulli− process. The first hitting{ } time 4. Gamma process and inverse gamma first hitting time. Consider a parent process X(t),t 0 with is the first Bernoulli trial t = S for which Xt equals { ≥ } 0. The number of rocket launches required to get m initial value X(0) = x0 > 0. Let X(t)= x0 Z(t) where Z(t),t 0 is a gamma process with− scale satellites in orbit is a simple example of this FHT { ≥ } model.
Recommended publications
  • On Hitting Times for Jump-Diffusion Processes with Past Dependent Local Characteristics
    Stochastic Processes and their Applications 47 ( 1993) 131-142 131 North-Holland On hitting times for jump-diffusion processes with past dependent local characteristics Manfred Schal Institut fur Angewandte Mathematik, Universitiit Bonn, Germany Received 30 March 1992 Revised 14 September 1992 It is well known how to apply a martingale argument to obtain the Laplace transform of the hitting time of zero (say) for certain processes starting at a positive level and being skip-free downwards. These processes have stationary and independent increments. In the present paper the method is extended to a more general class of processes the increments of which may depend both on time and past history. As a result a generalized Laplace transform is obtained which can be used to derive sharp bounds for the mean and the variance of the hitting time. The bounds also solve the control problem of how to minimize or maximize the expected time to reach zero. AMS 1980 Subject Classification: Primary 60G40; Secondary 60K30. hitting times * Laplace transform * expectation and variance * martingales * diffusions * compound Poisson process* random walks* M/G/1 queue* perturbation* control 1. Introduction Consider a process X,= S,- ct + uW, where Sis a compound Poisson process with Poisson parameter A which is disturbed by an independent standard Wiener process W This is the point of view of Dufresne and Gerber (1991). One can also look on X as a diffusion disturbed by jumps. This is the point of view of Ethier and Kurtz (1986, Section 4.10). Here u and A may be zero so that the jump term or the diffusion term may vanish.
    [Show full text]
  • THE MEASURABILITY of HITTING TIMES 1 Introduction
    Elect. Comm. in Probab. 15 (2010), 99–105 ELECTRONIC COMMUNICATIONS in PROBABILITY THE MEASURABILITY OF HITTING TIMES RICHARD F. BASS1 Department of Mathematics, University of Connecticut, Storrs, CT 06269-3009 email: [email protected] Submitted January 18, 2010, accepted in final form March 8, 2010 AMS 2000 Subject classification: Primary: 60G07; Secondary: 60G05, 60G40 Keywords: Stopping time, hitting time, progressively measurable, optional, predictable, debut theorem, section theorem Abstract Under very general conditions the hitting time of a set by a stochastic process is a stopping time. We give a new simple proof of this fact. The section theorems for optional and predictable sets are easy corollaries of the proof. 1 Introduction A fundamental theorem in the foundations of stochastic processes is the one that says that, under very general conditions, the first time a stochastic process enters a set is a stopping time. The proof uses capacities, analytic sets, and Choquet’s capacibility theorem, and is considered hard. To the best of our knowledge, no more than a handful of books have an exposition that starts with the definition of capacity and proceeds to the hitting time theorem. (One that does is [1].) The purpose of this paper is to give a short and elementary proof of this theorem. The proof is simple enough that it could easily be included in a first year graduate course in probability. In Section 2 we give a proof of the debut theorem, from which the measurability theorem follows. As easy corollaries we obtain the section theorems for optional and predictable sets. This argument is given in Section 3.
    [Show full text]
  • The Ito Integral
    Notes on the Itoˆ Calculus Steven P. Lalley May 15, 2012 1 Continuous-Time Processes: Progressive Measurability 1.1 Progressive Measurability Measurability issues are a bit like plumbing – a bit ugly, and most of the time you would prefer that it remains hidden from view, but sometimes it is unavoidable that you actually have to open it up and work on it. In the theory of continuous–time stochastic processes, measurability problems are usually more subtle than in discrete time, primarily because sets of measures 0 can add up to something significant when you put together uncountably many of them. In stochastic integration there is another issue: we must be sure that any stochastic processes Xt(!) that we try to integrate are jointly measurable in (t; !), so as to ensure that the integrals are themselves random variables (that is, measurable in !). Assume that (Ω; F;P ) is a fixed probability space, and that all random variables are F−measurable. A filtration is an increasing family F := fFtgt2J of sub-σ−algebras of F indexed by t 2 J, where J is an interval, usually J = [0; 1). Recall that a stochastic process fYtgt≥0 is said to be adapted to a filtration if for each t ≥ 0 the random variable Yt is Ft−measurable, and that a filtration is admissible for a Wiener process fWtgt≥0 if for each t > 0 the post-t increment process fWt+s − Wtgs≥0 is independent of Ft. Definition 1. A stochastic process fXtgt≥0 is said to be progressively measurable if for every T ≥ 0 it is, when viewed as a function X(t; !) on the product space ([0;T ])×Ω, measurable relative to the product σ−algebra B[0;T ] × FT .
    [Show full text]
  • Tight Inequalities Among Set Hitting Times in Markov Chains
    PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000–000 S 0002-9939(XX)0000-0 TIGHT INEQUALITIES AMONG SET HITTING TIMES IN MARKOV CHAINS SIMON GRIFFITHS, ROSS J. KANG, ROBERTO IMBUZEIRO OLIVEIRA, AND VIRESH PATEL Abstract. Given an irreducible discrete-time Markov chain on a finite state space, we consider the largest expected hitting time T (α) of a set of stationary measure at least α for α ∈ (0, 1). We obtain tight inequalities among the values of T (α) for different choices of α. One consequence is that T (α) ≤ T (1/2)/α for all α < 1/2. As a corollary we have that, if the chain is lazy in a certain sense as well as reversible, then T (1/2) is equivalent to the chain’s mixing time, answering a question of Peres. We furthermore demonstrate that the inequalities we establish give an almost everywhere pointwise limiting characterisation of possible hitting time functions T (α) over the domain α ∈ (0, 1/2]. 1. Introduction Hitting times are a classical topic in the theory of finite Markov chains, with connections to mixing times, cover times and electrical network representations [5, 6]. In this paper, we consider a natural family of extremal problems for maximum expected hitting times. In contrast to most earlier work on hitting times that considered the maximum expected hitting times of individual states, we focus on hitting sets of states of at least a given stationary measure. Informally, we are interested in the following basic question: how much more difficult is it to hit a smaller set than a larger one? (We note that other, quite different extremal problems about hitting times have been considered, e.g.
    [Show full text]
  • First Passage Times of Diffusion Processes and Their Applications To
    First Passage Times of Diffusion Processes and Their Applications to Finance A thesis presented for the degree of Doctor of Philosophy Luting Li Department of Statistics The London School of Economics and Political Science United Kingdom April 2019 Declaration I certify that the thesis I have presented for examination for the Ph.D. degree of the London School of Economics and Political Science is solely my own work other than where I have clearly indicated that it is the work of others (in which case the extent of any work carried out jointly by me and any other person is clearly identified in it). The copyright of this thesis rests with the author. Quotation from it is permitted, provided that full acknowledgement is made. This thesis may not be reproduced without the prior written consent of the author. I warrant that this authorisation does not, to the best of my belief, infringe the rights of any third party. Statement of conjoint work A version of parts of Chapters 3-6 has been submitted for publication in two articles jointly authored with A. Dassios. A version of Chapter 7 has been submitted for publication jointly authored with H. Xing. i Acknowledgements Just before writing down the words appearing on this page, I sit at my desk on the 7th floor of the Columbia House, recalling those countless yet enjoyable days and nights that I have spent on pursuing knowledge and exploring the world. To me, the marvellous beauty of this journey is once in a lifetime. I would like to, first and foremost, express my sincere gratitude to my supervisors Angelos Dassios and Hao Xing, for providing me the opportunity of carrying out my research under their supervisions, for their constant support and the freedom they granted me to conduct my work according to my preferences.
    [Show full text]
  • Section 3, Simple Gaussian Processes in Continuous Time 1 Introduction
    Stochastic Calculus, Courant Institute, Fall 2014 http://www.math.nyu.edu/faculty/goodman/teaching/StochCalc2014/index.html Section 3, Simple Gaussian processes in continuous time Jonathan Goodman, October 6, 2014 1 Introduction This section is on Brownian motion and a small extension of it, the Ornstein Uhlenbeck process. Brownian motion is in some senses the simplest continuous time continuous path stochastic process. It is a reasonable yet tractable model for many stochastic processes. It is used as a mathematical way to model random noise in most continuous time continuous sample path processes. The name comes from English biologist Robert Brown. In the 1820's he was looking at pollen grains in his new more powerful microscope. Amazingly, the particles seemed to be swimming. They were all in rapid and seemingly random motion. Brown then observed fine particles of clay, which move in a similar way. Brown, and science in general, were at a loss to understand what makes small particles move about randomly in water. The first physical explanation was given by Albert Einstein in 1905. Al- though clay particles are seem small to us, they are much larger than water molecules. Einstein suggested that the particles move because they are being pushed about by random collisions with water molecules. In this picture, the fine scale random motion of a clay particle is the large scale result of many small independent random impacts from water molecules. The clay particle motions are large scale relative to water molecules. Einstein made a quantitative theory of this and used it to give one of the first physical estimates of the size of water molecules.
    [Show full text]
  • Chapter 6 Brownian Motion and Beyond
    Chapter 6 Brownian motion and beyond Brownian motion is probably the most important stochastic process in probability theory and applications. Some of its basic properties are developed in the first half of this chapter, together with a brief discussion of martingales. Based on Brownian motion, the stochastic integrals and stochastic differential equations are introduced, including an application in option pricing in a single stock market. This topic is more advanced than the rest of these notes. An effort is made to present a complete description without being bogged down in advanced theory. For a more complete treatment, the reader is referred to a standard text on stochastic differential equations such as [7]. 6.1 Brownian motion Some motivation: Brownian motion has been used to describe the erratic motion of a pollen particle suspended in a liquid that is caused by the bombardment of surrounding molecules. Viewing the motion in one dimension, one may imagine a particle moving along the real line R in a continuous random path. It is reasonable to assume that the displacements of the particle over non-overlapping time intervals are independent and the distribution of each displacement depends only on how long the interval is. For any time t, let the interval [0; t] be divided into many small intervals of equal length, then the displacement over [0; t] is the sum of the iid displacements over the small intervals. Suggested by the CLT, the position of the particle at time t should be a normal random variable. This leads to the following formal definition. Definition of Brownian motion: Let Z(t) be a real-valued process indexed by continuous time t ≥ 0.
    [Show full text]
  • Advanced Probability
    Advanced Probability Alan Sola Department of Pure Mathematics and Mathematical Statistics University of Cambridge [email protected] Michaelmas 2014 Contents 1 Conditional expectation 4 1.1 Basic objects: probability measures, σ-algebras, and random variables . .4 1.2 Conditional expectation . .8 1.3 Integration with respect to product measures and Fubini's theorem . 12 2 Martingales in discrete time 15 2.1 Basic definitions . 15 2.2 Stopping times and the optional stopping theorem . 16 2.3 The martingale convergence theorem . 20 2.4 Doob's inequalities . 22 2.5 Lp-convergence for martingales . 24 2.6 Some applications of martingale techniques . 26 3 Stochastic processes in continuous time 31 3.1 Basic definitions . 31 3.2 Canonical processes and sample paths . 36 3.3 Martingale regularization theorem . 38 3.4 Convergence and Doob's inequalities in continuous time . 40 3.5 Kolmogorov's continuity criterion . 42 3.6 The Poisson process . 44 1 4 Weak convergence 46 4.1 Basic definitions . 46 4.2 Characterizations of weak convergence . 49 4.3 Tightness . 51 4.4 Characteristic functions and L´evy'stheorem . 53 5 Brownian motion 56 5.1 Basic definitions . 56 5.2 Construction of Brownian motion . 59 5.3 Regularity and roughness of Brownian paths . 61 5.4 Blumenthal's 01-law and martingales for Brownian motion . 62 5.5 Strong Markov property and reflection principle . 67 5.6 Recurrence, transience, and connections with partial differential equations . 70 5.7 Donsker's invariance principle . 76 6 Poisson random measures and L´evyprocesses 81 6.1 Basic definitions .
    [Show full text]
  • An Introduction to Stochastic Processes in Continuous Time
    An Introduction to Stochastic Processes in Continuous Time Harry van Zanten November 8, 2004 (this version) always under construction ii Preface iv Contents 1 Stochastic processes 1 1.1 Stochastic processes 1 1.2 Finite-dimensional distributions 3 1.3 Kolmogorov's continuity criterion 4 1.4 Gaussian processes 7 1.5 Non-di®erentiability of the Brownian sample paths 10 1.6 Filtrations and stopping times 11 1.7 Exercises 17 2 Martingales 21 2.1 De¯nitions and examples 21 2.2 Discrete-time martingales 22 2.2.1 Martingale transforms 22 2.2.2 Inequalities 24 2.2.3 Doob decomposition 26 2.2.4 Convergence theorems 27 2.2.5 Optional stopping theorems 31 2.3 Continuous-time martingales 33 2.3.1 Upcrossings in continuous time 33 2.3.2 Regularization 34 2.3.3 Convergence theorems 37 2.3.4 Inequalities 38 2.3.5 Optional stopping 38 2.4 Applications to Brownian motion 40 2.4.1 Quadratic variation 40 2.4.2 Exponential inequality 42 2.4.3 The law of the iterated logarithm 43 2.4.4 Distribution of hitting times 45 2.5 Exercises 47 3 Markov processes 49 3.1 Basic de¯nitions 49 3.2 Existence of a canonical version 52 3.3 Feller processes 55 3.3.1 Feller transition functions and resolvents 55 3.3.2 Existence of a cadlag version 59 3.3.3 Existence of a good ¯ltration 61 3.4 Strong Markov property 64 vi Contents 3.4.1 Strong Markov property of a Feller process 64 3.4.2 Applications to Brownian motion 68 3.5 Generators 70 3.5.1 Generator of a Feller process 70 3.5.2 Characteristic operator 74 3.6 Killed Feller processes 76 3.6.1 Sub-Markovian processes 76 3.6.2
    [Show full text]
  • Strict Local Martingales and Bubbles
    The Annals of Applied Probability 2015, Vol. 25, No. 4, 1827–1867 DOI: 10.1214/14-AAP1037 c Institute of Mathematical Statistics, 2015 STRICT LOCAL MARTINGALES AND BUBBLES By Constantinos Kardaras, Dorte¨ Kreher1 and Ashkan Nikeghbali London School of Economics and Political Science, Humboldt-Universit¨at zu Berlin and Universit¨at Z¨urich This paper deals with asset price bubbles modeled by strict local martingales. With any strict local martingale, one can associate a new measure, which is studied in detail in the first part of the paper. In the second part, we determine the “default term” apparent in risk-neutral option prices if the underlying stock exhibits a bubble modeled by a strict local martingale. Results for certain path dependent options and last passage time formulas are given. 1. Introduction. The goal of this paper is to determine the influence of asset price bubbles on the pricing of derivatives. Asset price bubbles have been studied extensively in the economic literature looking for explanations of why they arise, but have only recently gained attention in mathemati- cal finance by Cox and Hobson [5], Pal and Protter [30], and Jarrow et al. [20–22]. When an asset price bubble exists, the market price of the asset is higher than its fundamental value. From a mathematical point of view, this is the case when the stock price process is modeled by a positive strict local martingale under the equivalent local martingale measure. Here, by a strict local martingale, we understand a local martingale, which is not a true mar- tingale. Strict local martingales were first studied in the context of financial mathematics by Delbaen and Schachermayer [6].
    [Show full text]
  • Universal Hitting Time Statistics for Integrable Flows
    UNIVERSAL HITTING TIME STATISTICS FOR INTEGRABLE FLOWS CARL P. DETTMANN, JENS MARKLOF, AND ANDREAS STROMBERGSSON¨ To Professors D. Ruelle and Ya.G. Sinai on the occasion of their 80th birthday Abstract. The perceived randomness in the time evolution of \chaotic" dynamical systems can be characterized by universal probabilistic limit laws, which do not depend on the fine features of the individual system. One important example is the Poisson law for the times at which a particle with random initial data hits a small set. This was proved in various settings for dynamical systems with strong mixing properties. The key result of the present study is that, despite the absence of mixing, the hitting times of integrable flows also satisfy universal limit laws which are, however, not Poisson. We describe the limit distributions for \generic" integrable flows and a natural class of target sets, and illustrate our findings with two examples: the dynamics in central force fields and ellipse billiards. The convergence of the hitting time process follows from a new equidistribution theorem in the space of lattices, which is of independent interest. Its proof exploits Ratner's measure classification theorem for unipotent flows, and extends earlier work of Elkies and McMullen. 1. Introduction Let (M; F ; ν) be a probability space and consider a measure-preserving dynamical system (1.1) 't : M!M: A fundamental question is how often a trajectory with random initial data x 2 M intersects a given target set D 2 F within time t. If D is fixed, this problem has led to many important developments in ergodic theory, which show that, if 't is sufficiently \chaotic" (e.g., partially hyperbolic), the number of intersections satisfies a central limit theorem and more general invariance principles.
    [Show full text]
  • Conformal Invariance of the Loop-Erased Random Walk Intensity
    Conformal invariance of the loop-erased random walk intensity Christian Beneˇs ∗1, Fredrik Johansson Viklundy2, and Gregory F. Lawlerz3 1Brooklyn College of the City University of New York 2Uppsala University and Columbia University 3University of Chicago January 13, 2014 Abstract The loop-erased random walk (LERW) edge intensity is simply the probability that the LERW passes through a given edge. We show that this normalized probability (in the bulk) of a planar LERW between two boundary points of a square grid approximation of a simply connected domain converges in the lattice size scaling limit to an explicit conformally covariant quantity which coincides with the SLE2 Green's function. Contents 1 Introduction and outline of proof 2 1.1 Introduction . .2 1.2 Acknowledgments . .3 1.3 Notation and set-up . .4 1.4 Outline of proof and a key idea . .6 2 Preliminaries 9 2.1 Green's functions . .9 2.1.1 Loops and loop measures . 10 2.1.2 Loop-erased walk . 11 2.2 Brownian loop measure . 12 2.3 Poisson kernels . 13 2.4 KMT coupling . 15 3 The combinatorial identity: proof of Theorem 1.3 15 4 Comparison of loop measures: proof of Theorem 1.4 19 ∗[email protected] [email protected] [email protected] 1 5 Asymptotics of Λ: proof of Theorem 1.5 24 5.1 Continuum functions . 24 5.2 Asymptotics of ΛA(0; a).................................. 26 5.3 Poisson kernel convergence in the slit square: Proof of Theorem 5.5 . 31 5.4 From the Slit Square to @A: Proof of Proposition 5.7 .
    [Show full text]