First-hitting-time model

From Infogalactic: the planetary knowledge core
Jump to: navigation, search

Lua error in package.lua at line 80: module 'strict' not found.

In statistics, first-hitting-time models are a sub-class of survival models. The first hitting time, also called first passage time, of a set A with respect to an instance of a stochastic process is the time until the stochastic process first enters A.

More colloquially, a first passage time in a stochastic system, is the time taken for a state variable to reach a certain value. Understanding this metric allows one to further understand the physical system under observation, and as such has been the topic of research in very diverse fields, from Economics to Ecology.[1]

Examples

A common example of a first-hitting-time model is a ruin problem, such as Gambler's ruin. In this example, an entity (often described as a gambler or an insurance company) has an amount of money which varies randomly with time, possibly with some drift. The model considers the event that the amount of money reaches 0, representing bankruptcy. The model can answer questions such as the probability that this occurs within finite time, or the mean time until which it occurs.

First-hitting-time models can be applied to expected lifetimes, of patients or mechanical devices. When the process reaches an adverse threshold state for the first time, the patient dies, or the device breaks down.

First passage time of a 1D Brownian Particle

One of the simplest and omnipresent stochastic systems is that of the Brownian particle in one dimension. This system describes the motion of a particle which moves stochastically in one dimensional space, with equal probability of moving to the left or to the right. Given that Brownian motion is used often as a tool to understand more complex phenomena, it is important to understand the probability of a first passage time of the Brownian particle of reaching some position distant from its start location. This is done through the following means.

The probability density function (PDF) for a particle in one dimension is found by solving the one-dimensional diffusion equation. (This equation states that the position probability density diffuses out over time. It is analogous to say, cream in a cup of coffee if the cream was all contained within some small location initially. In the long time limit the cream has diffused throughout the entire drink evenly.) Namely,


 \frac{\partial p(x,t \mid x_{0})}{\partial t}=D\frac{\partial^2p(x,t \mid
x_{0})}{\partial x^2},

given the initial condition p(x,t={0} \mid x_{0})=\delta(x-x_{0}); where x(t) is the position of the particle at some given time, x_0 is the tagged particle's initial position, and D is the diffusion constant with the S.I. units m^2s^{-1} (an indirect measure of the particle's speed). The bar in the argument of the instantaneous probability refers to the conditional probability. The diffusion equation states that the speed at which the probability for finding the particle at x(t) is position dependent.

It can be shown that the one-dimensional PDF is


 p(x,t; x_0)=\frac{1}{\sqrt{4\pi Dt}}\exp\left(-\frac{(x-x_0)^2}{4Dt}\right).

This states that the probability of finding the particle at x(t) is Gaussian, and the width of the Gaussian is time dependent. More specifically the Full Width at Half Maximum (FWHM) – technically, this is actually the Full Duration at Half Maximum as the independent variable is time – scales like


 \rm{FWHM}\sim\sqrt{t}.

Using the PDF one is able to derive the average of a given function, L, at time t:


\langle  L(t) \rangle\equiv \int^{\infty}_{-\infty} L(x,t) p(x,t) dx,

where the average is taken over all space (or any applicable variable).

The First Passage Time Density (FPTD) is the probability that a particle has first reached a point x_c at time t. This probability density is calculable from the Survival probability (a more common probability measure in statistics). Consider the absorbing boundary condition p(x_c,t)=0 (The subscript c for the absorption point x_c is an abbreviation for cliff used in many texts as an analogy to an absorption point). The PDF satisfying this boundary condition is given by


p(x,t; x_0, x_c) = \frac{1}{\sqrt{4\pi Dt}} \left( \exp\left(-\frac{(x-x_0)^2}{4Dt}\right) - \exp\left(-\frac{(x-(2x_c-x_0))^2}{4Dt}\right) \right),

for x<x_c. The survival probability, the probability that the particle has remained at a position x < x_c for all times up to t, is given by


 S(t)\equiv\int_{-\infty}^{x_c} p(x,t; x_{0}, x_c) dx = \operatorname{erf}\left(\frac{x_c-x_{0}}{2\sqrt{D t}}\right),

where \operatorname{erf} is the error function. The relation between the Survival probability and the FPTD is as follows(the probability that a particle has reached the absorption point between times t and t+dt is f(t)dt=S(t)-S(t+dt). If one uses the first-order Taylor approximation, the definition of the FPTD follows):

 f(t)=-\frac{\partial S(t)}{\partial t}.

By using the diffusion equation and integrating by parts, the explicit FPTD is


 f(t)\equiv\frac{|x_c-x_{0}|}{\sqrt{4\pi Dt^3}} \exp\left(- \frac{(x_c-x_{0})^2}{4Dt}\right).

The first-passage time for a Brownian particle therefore follows a Lévy distribution.

For t\gg\frac{(x_c-x_{0})^2}{4D}, it follows from above that


 f(t)=\frac{\Delta x}{\sqrt{4\pi Dt^3}}\sim t^{-3/2},

where \Delta x\equiv |x_c-x_{0}|. This equation states that the probability for a Brownian particle achieving a first passage at some long time (defined in the paragraph above) becomes increasingly small, but always finite.

The first moment of the FPTD diverges (as it is a so-called heavy-tailed distribution), therefore one cannot calculate the average FPT, so instead, one can calculate the typical time, the time when the FPTD is at a maximum (\partial f/\partial t=0), i.e.,


\tau_{\rm{ty}}=\frac{\Delta x^2}{6D}.

Latent vs observable

In many real world applications, the process is latent, or unobservable. When first hitting time models are equipped with regression structures, accommodating covariate data, we call such regression structure Threshold regression. The threshold state, parameters of the process, and even time scale may depend on corresponding covariates.

A first-hitting-time (FHT) model has two underlying components: (1) a parent stochastic process \{X(t)\}\,\,, and (2) a threshold. The first hitting time is defined as the time when the stochastic process first reaches the threshold. It is very important to distinguish whether the sample path of the parent process is latent (i.e., unobservable) or observable, and such distinction is a characteristic of the FHT model. By far, latent processes are most common. To give an example, we can use a Wiener process \{X(t), t\ge0\,\}\, as the parent stochastic process. Such Wiener process can be defined with the mean parameter {\mu}\,\,, the variance parameter {\sigma^2}\,\,, and the initial value X(0)=x_0>0\,.

See also

References

  1. Redner 2001
  • Lua error in package.lua at line 80: module 'strict' not found.
  • Lua error in package.lua at line 80: module 'strict' not found.
  • Lua error in package.lua at line 80: module 'strict' not found.
  • Lua error in package.lua at line 80: module 'strict' not found.