Boundary-crossing probabilities for the Brownian motion and Poisson processes and techniques for computing the power of the Kolmogorov-Smirnov test

1971 ◽  
Vol 8 (03) ◽  
pp. 431-453 ◽  
Author(s):  
J. Durbin

Let w(t), 0 ≦ t ≦ ∞, be a Brownian motion process, i.e., a zero-mean separable normal process with Pr{w(0) = 0} = 1, E{w(t 1)w(t 2)}= min (t 1, t 2), and let a, b denote the boundaries defined by y = a(t), y = b(t), where b(0) < 0 < a(0) and b(t) < a(t), 0 ≦ t ≦ T ≦ ∞. A basic problem in many fields such as diffusion theory, gambler's ruin, collective risk, Kolmogorov-Smirnov statistics, cumulative-sum methods, sequential analysis and optional stopping is that of calculating the probability that a sample path of w(t) crosses a or b before t = T. This paper shows how this probability may be computed for sufficiently smooth boundaries by numerical solution of integral equations for the first-passage distribution functions. The technique used is to approximate the integral equations by linear recursions whose coefficients are estimated by linearising the boundaries within subintervals. The results are extended to cover the tied-down process subject to the condition w(1) = 0. Some related results for the Poisson process and the sample distribution function are given. The procedures suggested are exemplified numerically, first by computing the probability that the tied-down Brownian motion process crosses a particular curved boundary for which the true probability is known, and secondly by computing the finite-sample and asymptotic powers of the Kolmogorov-Smirnov test against a shift in mean of the exponential distribution.

1971 ◽  
Vol 8 (3) ◽  
pp. 431-453 ◽  
Author(s):  
J. Durbin

Let w(t), 0 ≦ t ≦ ∞, be a Brownian motion process, i.e., a zero-mean separable normal process with Pr{w(0) = 0} = 1, E{w(t1)w(t2)}= min (t1, t2), and let a, b denote the boundaries defined by y = a(t), y = b(t), where b(0) < 0 < a(0) and b(t) < a(t), 0 ≦ t ≦ T ≦ ∞. A basic problem in many fields such as diffusion theory, gambler's ruin, collective risk, Kolmogorov-Smirnov statistics, cumulative-sum methods, sequential analysis and optional stopping is that of calculating the probability that a sample path of w(t) crosses a or b before t = T. This paper shows how this probability may be computed for sufficiently smooth boundaries by numerical solution of integral equations for the first-passage distribution functions. The technique used is to approximate the integral equations by linear recursions whose coefficients are estimated by linearising the boundaries within subintervals. The results are extended to cover the tied-down process subject to the condition w(1) = 0. Some related results for the Poisson process and the sample distribution function are given. The procedures suggested are exemplified numerically, first by computing the probability that the tied-down Brownian motion process crosses a particular curved boundary for which the true probability is known, and secondly by computing the finite-sample and asymptotic powers of the Kolmogorov-Smirnov test against a shift in mean of the exponential distribution.


1992 ◽  
Vol 29 (02) ◽  
pp. 291-304 ◽  
Author(s):  
J. Durbin ◽  
D. Williams

An expression for the first-passage density of Brownian motion to a curved boundary is expanded as a series of multiple integrals. Bounds are given for the error due to truncation of the series when the boundary is wholly concave or wholly convex. Extensions to the Brownian bridge and to continuous Gauss–Markov processes are given. The series provides a practical method for calculating the probability that a sample path crosses the boundary in a specified time-interval to a high degree of accuracy. A numerical example is given.


1980 ◽  
Vol 85 ◽  
pp. 221-222
Author(s):  
M. Buchholz ◽  
Th. Schmidt-Kaler

The radial mass distribution (obtained by counting stars in strips) of the real cluster is compared successively to the distribution functions of a simulated cluster of 100 stars, each of which corresponds to a certain dynamical age, Tdyn, The value of Tdyn, belonging to the function most similar to the observed one is taken to be the dynamical age of the cluster. The radius is given in units of R1/2 (sphere containing half of the total mass); this unit is nearly time-independent. The difference between the distribution functions is measured by the maximum Δmax of the Kolmogorov-Smirnov test which is free from assumptions on the form of the distributions. The minimum in the plot Δmax vs Tdyn, indicates the age of the cluster. It is then converted into an absolute age, Tabs (in years), by The error due to the dynamical theory (limited number of distribution functions, etc.) is estimated at 12%, the error due to the uncertainty of diameter and mass of the cluster is about 30%. Unreliable results were obtained in case of strongly inhomogeneous reddening of the cluster. As an example, the plot of the test values for NGC 457 is given in Figure 1.


1993 ◽  
Vol 30 (01) ◽  
pp. 17-27
Author(s):  
Aimé Lachal

Let be the Brownian motion process starting at the origin, its primitive and Ut = (Xt+x + ty, Bt + y), , the associated bidimensional process starting from a point . In this paper we present an elementary procedure for re-deriving the formula of Lefebvre (1989) giving the Laplace–Fourier transform of the distribution of the couple (σ α, Uσa ), as well as Lachal's (1991) formulae giving the explicit Laplace–Fourier transform of the law of the couple (σ ab, Uσab ), where σ α and σ ab denote respectively the first hitting time of from the right and the first hitting time of the double-sided barrier by the process . This method, which unifies and considerably simplifies the proofs of these results, is in fact a ‘vectorial' extension of the classical technique of Darling and Siegert (1953). It rests on an essential observation (Lachal (1992)) of the Markovian character of the bidimensional process . Using the same procedure, we subsequently determine the Laplace–Fourier transform of the conjoint law of the quadruplet (σ α, Uσa, σb, Uσb ).


2000 ◽  
Vol 23 (11) ◽  
pp. 759-776 ◽  
Author(s):  
Seung Jun Chang ◽  
Soon Ja Kang ◽  
David Skoug

We use a generalized Brownian motion process to define a generalized Feynman integral and a conditional generalized Feynman integral. We then establish the existence of these integrals for various functionals. Finally we use the conditional generalized Feynman integral to derive a Schrödinger integral equation.


2013 ◽  
Vol 2013 ◽  
pp. 1-13 ◽  
Author(s):  
ZaiTang Huang ◽  
ChunTao Chen

We study the stability, attractors, and bifurcation of stochastic Rayleigh-van der Pol equations with jumps. We first established the stochastic stability and the large deviations results for the stochastic Rayleigh-van der Pol equations. We then examine the existence limit circle and obtain some new random attractors. We further establish stochastic bifurcation of random attractors. Interestingly, this shows the effect of the Poisson noise which can stabilize or unstabilize the system which is significantly different from the classical Brownian motion process.


1993 ◽  
Vol 30 (1) ◽  
pp. 17-27 ◽  
Author(s):  
Aimé Lachal

Let be the Brownian motion process starting at the origin, its primitive and Ut = (Xt+x + ty, Bt + y), , the associated bidimensional process starting from a point . In this paper we present an elementary procedure for re-deriving the formula of Lefebvre (1989) giving the Laplace–Fourier transform of the distribution of the couple (σ α, Uσa), as well as Lachal's (1991) formulae giving the explicit Laplace–Fourier transform of the law of the couple (σ ab, Uσab), where σ α and σ ab denote respectively the first hitting time of from the right and the first hitting time of the double-sided barrier by the process . This method, which unifies and considerably simplifies the proofs of these results, is in fact a ‘vectorial' extension of the classical technique of Darling and Siegert (1953). It rests on an essential observation (Lachal (1992)) of the Markovian character of the bidimensional process .Using the same procedure, we subsequently determine the Laplace–Fourier transform of the conjoint law of the quadruplet (σ α, Uσa, σb, Uσb).


Sign in / Sign up

Export Citation Format

Share Document