We propose a stochastic modelling of the PCR amplification process by a size-dependent branching process starting as a supercritical Bienaymé-Galton-Watson transient phase and then having a saturation near-critical size-dependent phase. This model allows us to estimate the probability of replication of a DNA molecule at each cycle of a single PCR trajectory with a very good accuracy.

We propose a new convergence criterion for the stochastic algorithm for the optimization of probabilities (SAOP) described in an earlier paper. The criterion is based on the dissection principle for irreducible finite Markov chains.

In this paper we study different types of planar random motions (performed with constant velocity) with three directions, defined by the vectors dj = (cos(2πj/3), sin(2πj/3)) for j = 0, 1, 2, changing at Poisson-paced times. We examine the cyclic motion (where the change of direction is deterministic), the completely uniform motion (where at each Poisson event each direction can be taken with probability ) and the symmetrically deviating case (where the particle can choose all directions except that taken before the Poisson event). For each of the above random motions we derive the explicit distribution of the position of the particle, by using an approach based on order statistics. We prove that the densities obtained are solutions of the partial differential equations governing the processes. We are also able to give the explicit distributions on the boundary and, for the case of the symmetrically deviating motion, we can write it as the distribution of a telegraph process. For the symmetrically deviating motion we use a generalization of the Bose-Einstein statistics in order to determine the distribution of the triple (N0, N1, N2) (conditional on N(t) = k, with N0 + N1 + N2 = N(t) + 1, where N(t) is the number of Poisson events in [0, t]), where Nj denotes the number of times the direction dj (j = 0, 1, 2) is taken. Possible extensions to four directions or more are briefly considered.

We analyse the role of Euler summation in a numerical inversion algorithm for Laplace transforms due to Abate and Whitt called the EULER algorithm. Euler summation is shown to accelerate convergence of a slowly converging truncated Fourier series; an explicit bound for the approximation error is derived that generalizes a result given by O'Cinneide. An enhanced inversion algorithm called EULER-GPS is developed using a new variant of Euler summation. The algorithm EULER-GPS makes it possible to accurately invert transforms of functions with discontinuities at arbitrary locations. The effectiveness of the algorithm is verified through numerical experiments. Besides numerical transform inversion, the enhanced algorithm is applicable to a wide range of other problems where the goal is to recover point values of a piecewise-smooth function from the Fourier series.

This paper investigates the finite- and infinite-time ruin probabilities in a discrete-time stochastic economic environment. Under the assumption that the insurance risk - the total net loss within one time period - is extended-regularly-varying or rapidly-varying tailed, various precise estimates for the ruin probabilities are derived. In particular, some estimates obtained are uniform with respect to the time horizon, and so apply in the case of infinite-time ruin.

Discrete time-series models are commonly used to represent economic and physical data. In decision making and system control, the first-passage time and level-crossing probabilities of these processes against certain threshold levels are important quantities. In this paper, we apply an integral-equation approach together with the state-space representations of time-series models to evaluate level-crossing probabilities for the AR(p) and ARMA(1,1) models and the mean first passage time for AR(p) processes. We also extend Novikov's martingale approach to ARMA(p,q) processes. Numerical schemes are used to solve the integral equations for specific examples.