We consider upper‒lower (UL) (and lower‒upper (LU)) factorizations of the one-step transition probability matrix of a random walk with the state space of nonnegative integers, with the condition that both upper and lower triangular matrices in the factorization are also stochastic matrices. We provide conditions on the free parameter of the UL factorization in terms of certain continued fractions such that this stochastic factorization is possible. By inverting the order of the factors (also known as a Darboux transformation) we obtain a new family of random walks where it is possible to state the spectral measures in terms of a Geronimus transformation. We repeat this for the LU factorization but without a free parameter. Finally, we apply our results in two examples; the random walk with constant transition probabilities, and the random walk generated by the Jacobi orthogonal polynomials. In both situations we obtain urn models associated with all the random walks in question.

We present a new composite quadrature rule which is exact for polynomials of degree 2N + K – 1 with N abscissas at each subinterval and K boundary conditions. The corresponding orthogonal polynomials are introduced and the analytic formulae for abscissas and weight functions are presented. Numerical results show that the new quadrature rule is more efficient, compared with classical ones.

We consider a birth–death process {X(t),t≥0} on the positive integers for which the origin is an absorbing state with birth coefficients λn,n≥0, and death coefficients μn,n≥0. If we define A=∑n=1∞ 1/λnπn and S=∑n=1∞ (1/λnπn)∑i=n+1∞ πi, where {πn,n≥1} are the potential coefficients, it is a well-known fact (see van Doorn (1991)) that if A=∞ and S<∞, then λC>0 and there is precisely one quasistationary distribution, namely, {aj(λC)}, where λC is the decay parameter of {X(t),t≥0} in C={1,2,...} and aj(x)≡μ1-1πjxQj(x), j=1,2,.... In this paper we prove that there is a unique quasistationary distribution that attracts all initial distributions supported in C, if and only if the birth–death process {X(t),t≥0} satisfies bothA=∞ and S<∞. That is, for any probability measure M={mi, i=1,2,...}, we have limt→∞ℙM(X(t)=j∣ T>t)= aj(λC), j=1,2,..., where T=inf{t≥0 : X(t)=0} is the extinction time of {X(t),t≥0} if and only if the birth–death process {X(t),t≥0} satisfies both A=∞ and S<∞.

Suppose that n identical particles evolve according to the same marginal Markov chain. In this setting we study chains such as the Ehrenfest chain that move a prescribed number of randomly chosen particles at each epoch. The product chain constructed by this device inherits its eigenstructure from the marginal chain. There is a further chain derived from the product chain called the composition chain that ignores particle labels and tracks the numbers of particles in the various states. The composition chain in turn inherits its eigenstructure and various properties such as reversibility from the product chain. The equilibrium distribution of the composition chain is multinomial. The current paper proves these facts in the well-known framework of state lumping and identifies the column eigenvectors of the composition chain with the multivariate Krawtchouk polynomials of Griffiths. The advantages of knowing the full spectral decomposition of the composition chain include (a) detailed estimates of the rate of convergence to equilibrium, (b) construction of martingales that allow calculation of the moments of the particle counts, and (c) explicit expressions for mean coalescence times in multi-person random walks. These possibilities are illustrated by applications to Ehrenfest chains, the Hoare and Rahman chain, Kimura's continuous-time chain for DNA evolution, a light bulb chain, and random walks on some specific graphs.

In this paper we develop a constructive approach to studying continuously and discretely sampled functionals of Lévy processes. Estimates for the rate of convergence of the discretely sampled functionals to the continuously sampled functionals are derived, reducing the study of the latter to that of the former. Laguerre reduction series for the discretely sampled functionals are developed, reducing their study to that of the moment generating function of the pertinent Lévy processes and to that of the moments of these processes in particular. The results are applied to questions of contingent claim valuation, such as the explicit valuation of Asian options, and illustrated in the case of generalized inverse Gaussian Lévy processes.

On the basis of the Karlin and McGregor result, which states that the transition probability functions of a birth and death process can be expressed via the introduction of an orthogonal polynomial system and a spectral measure, we investigate in this paper how the Laplace transforms and the distributions of different transient characteristics related to excursions of a birth and death process can be expressed by means of the basic orthogonal polynomial system and the spectral measure. This allows us in particular to give a probabilistic interpretation of the series introduced by Stieltjes to study the convergence of the fundamental continued fraction associated with the system. Throughout the paper, we pay special attention to the case when the birth and death process is ergodic. Under the assumption that the spectrum of the spectral measure is discrete, we show how the distributions of different random variables associated with excursions depend on the fundamental continued fraction, the orthogonal polynomial system and the spectral measure.

Let α and β be any two real numbers and let be the Jacobi polynomial sequences. For any non-zero real number a, is an orthogonal polynomial sequence with respect to convolution if and only if either (i) b = 1, α = 0 and β + 1 is not equal to a negative integer or (ii) b = — 1, β = 0 and α + 1 is not equal to a negative integer.