No Arabic abstract
The stationary asymptotic properties of the diffusion limit of a multi-type branching process with neutral mutations are studied. For the critical and subcritical processes the interesting limits are those of quasi-stationary distributions conditioned on non-extinction. Limiting distributions for supercritical and critical processes are found to collapse onto rays aligned with stationary eigenvectors of the mutation rate matrix, in agreement with known results for discrete multi-type branching processes. For the sub-critical process the quasi-stationary distribution is obtained to first order in the overall mutation rate, which is assumed to be small. The sampling distribution over allele types for a sample of given finite size is found to agree to first order in mutation rates with the analogous sampling distribution for a Wright-Fisher diffusion with constant population size.
We show that the quasi-stationary distribution of the subcritical contact process on $mathbb{Z}^d$ is unique. This is in contrast with other processes which also do not come down from infinity, like stable queues and Galton-Watson, and it seems to be the first such example.
A population is considered stationary if the growth rate is zero and the age structure is constant. It thus follows that a population is considered non-stationary if either its growth rate is non-zero and/or its age structure is non-constant. We propose three properties that are related to the stationary population identity (SPI) of population biology by connecting it with stationary populations and non-stationary populations which are approaching stationarity. One of these important properties is that SPI can be applied to partition a population into stationary and non-stationary components. These properties provide deeper insights into cohort formation in real-world populations and the length of the duration for which stationary and non-stationary conditions hold. The new concepts are based on the time gap between the occurrence of stationary and non-stationary populations within the SPI framework that we refer to as Oscillatory SPI and the Amplitude of SPI. This article will appear in Bulletin of Mathematical Biology (Springer)
Given a domain G, a reflection vector field d(.) on the boundary of G, and drift and dispersion coefficients b(.) and sigma(.), let L be the usual second-order elliptic operator associated with b(.) and sigma(.). Under suitable assumptions that, in particular, ensure that the associated submartingale problem is well posed, it is shown that a probability measure $pi$ on bar{G} is a stationary distribution for the corresponding reflected diffusion if and only if $pi (partial G) = 0$ and $int_{bar{G}} L f (x) pi (dx) leq 0$ for every f in a certain class of test functions. Moreover, the assumptions are shown to be satisfied by a large class of reflected diffusions in piecewise smooth multi-dimensional domains with possibly oblique reflection.
Consider a reflecting diffusion in a domain in $R^d$ that acquires drift in proportion to the amount of local time spent on the boundary of the domain. We show that the stationary distribution for the joint law of the position of the reflecting process and the value of the drift vector has a product form. Moreover, the first component is the symmetrizing measure on the domain for the reflecting diffusion without inert drift, and the second component has a Gaussian distribution. We also consider processes where the drift is given in terms of the gradient of a potential.
We build a sequence of empirical measures on the space D(R_+,R^d) of R^d-valued c`adl`ag functions on R_+ in order to approximate the law of a stationary R^d-valued Markov and Feller process (X_t). We obtain some general results of convergence of this sequence. Then, we apply them to Brownian diffusions and solutions to Levy driven SDEs under some Lyapunov-type stability assumptions. As a numerical application of this work, we show that this procedure gives an efficient way of option pricing in stochastic volatility models.