We study collision probabilities concerning the simple balls-and-bins problem developed by Wendl (2003). In this article we give the factorial moment of the number of collisions. Moreover, we obtain a Poisson approximation for the number of collisions using the Chen-Stein method.

Let X and Y be two simple symmetric continuous-time random walks on the vertices of the n-dimensional hypercube, Z2n. We consider the class of co-adapted couplings of these processes, and describe an intuitive coupling which is shown to be the fastest in this class.

We study a model arising in chemistry where n elements numbered 1, 2, …, n are randomly permuted and if i is immediately to the left of i + 1 then they become stuck together to form a cluster. The resulting clusters are then numbered and considered as elements, and this process keeps repeating until only a single cluster is remaining. In this article we study properties of the distribution of the number of permutations required.

In this paper we present a robust numerical method to compute the scale function W(q)(x) of a general spectrally negative Lévy process (X, P). The method is based on the Esscher transform of measure Pν under which X is taken and the scale function is determined. This change of measure makes it possible for the scale function to be bounded and, hence, makes numerical computation easy, fast, and stable. Working under the new measure Pν and using the method of Abate and Whitt (1992) and Choudhury, Lucantoni, and Whitt (1994), we give a fast stable numerical algorithm for the computation of W(q)(x).

We describe the genealogy of individuals with infinite descent in a supercritical continuous-state branching process.

Detection of repeated sequences within complete genomes is a powerful tool to help understanding genome dynamics and species evolutionary history. To distinguish significant repeats from those that can be obtained just by chance, statistical methods have to be developed. In this paper we show that the distribution of the number of long repeats in long sequences generated by stationary Markov chains can be approximated by a Poisson distribution with explicit parameter. Thanks to the Chen-Stein method we provide a bound for the approximation error; this bound converges to 0 as soon as the length n of the sequence tends to ∞ and the length t of the repeats satisfies n2ρt = O(1) for some 0 < ρ < 1. Using this Poisson approximation, p-values can then be easily calculated to determine if a given genome is significantly enriched in repeats of length t.

We present two new stochastic volatility models in which option prices for European plain-vanilla options have closed-form expressions. The models are motivated by the well-known SABR model, but use modified dynamics of the underlying asset. The asset process is modelled as a product of functions of two independent stochastic processes: a Cox-Ingersoll-Ross process and a geometric Brownian motion. An application of the models to options written on foreign currencies is studied.

The class of processes formed as the aggregation of Ornstein-Uhlenbeck processes has proved useful in modeling time series from a number of areas and includes several interesting special cases. This paper examines the second-order properties of this class. Bounds on the one-step prediction error variance are proved and consistency of the minimum contrast estimation is demonstrated.

We compute the posterior distributions of the initial population and parameter of binary branching processes in the limit of a large number of generations. We compare this Bayesian procedure with a more naïve one, based on hitting times of some random walks. In both cases, central limit theorems are available, with explicit variances.

The TCP window size process can be modeled as a piecewise-deterministic Markov process that increases linearly and experiences downward jumps at Poisson times. We present a transient analysis of this window size process. Our main result is the Laplace transform of the transient moments. Formulae for the integer and fractional moments are derived, as well as an explicit characterization of the speed of convergence to steady state. Central to our approach are the infinitesimal generator and Dynkin's martingale.

In this paper, a fourth moment bound for partial sums of functionals of strongly ergodic Markov chains is established. This type of inequality plays an important role in the study of the empirical process invariance principle. This inequality is specially adapted to the technique of Dehling, Durieu, and Volný (2008). The same moment bound can be proved for dynamical systems whose transfer operator has some spectral properties. Examples of applications are given.

The distributions of the run occurrences for a sequence of independent and identically distributed (i.i.d.) experiments are usually obtained by combinatorial methods (see Balakrishnan and Koutras (2002, Chapter 5)) and the resulting formulae are often very tedious, while the distributions for non i.i.d. experiments are generally intractable. It is therefore of practical interest to find a suitable approximate model with reasonable approximation accuracy. In this paper we demonstrate that the negative binomial distribution is the most suitable approximate model for the number of k-runs: it outperforms the Poisson approximation, the general compound Poisson approximation as observed in Eichelsbacher and Roos (1999), and the translated Poisson approximation in Rollin (2005). In particular, its accuracy of approximation in terms of the total variation distance improves when the number of experiments increases, in the same way as the normal approximation improves in the Berry-Esseen theorem.

The representation of the reliability function of the lifetime of a coherent system as a mixture of the reliability function of order statistics associated with the lifetimes of its components is a very useful tool to study the ordering and the limiting behaviour of coherent systems. In this paper, we obtain several representations of the reliability functions of residual lifetimes of used coherent systems under two particular conditions on the status of the components or the system in terms of the reliability functions of residual lifetimes of order statistics.

A dual Markov branching process (DMBP) is by definition a Siegmund's predual of some Markov branching process (MBP). Such a process does exist and is uniquely determined by the so-called dual-branching property. Its q-matrix Q is derived and proved to be regular and monotone. Several equivalent definitions for a DMBP are given. The criteria for transience, positive recurrence, strong ergodicity, and the Feller property are established. The invariant distributions are given by a clear formulation with a geometric limit law.

The mean time to extinction of a critical Galton-Watson process with initial population size k is shown to be asymptotically equivalent to two integral transforms: one involving the kth iterate of the probability generating function and one involving the generating function itself. Relating the growth of these transforms to the regular variation of their arguments, immediately connects statements involving the regular variation of the probability generating function, its iterates at 0, the quasistationary measures, their partial sums, and the limiting distribution of the time to extinction. In the critical case of finite variance we also give the growth of the mean time to extinction, conditioned on extinction occurring by time n.

In this paper a branching process approximation for the spread of a Reed-Frost epidemic on a network with tunable clustering is derived. The approximation gives rise to expressions for the epidemic threshold and the probability of a large outbreak in the epidemic. We investigate how these quantities vary with the clustering in the graph and find that, as the clustering increases, the epidemic threshold decreases. The network is modeled by a random intersection graph, in which individuals are independently members of a number of groups and two individuals are linked to each other if and only if there is at least one group that they are both members of.

By considering k-out-of-n systems with independent and nonidentically distributed components, we discuss stochastic monotone properties of the residual life and the inactivity time. We then present some stochastic comparisons of two systems based on the residual life and inactivity time.

The reduced Markov branching process is a stochastic model for the genealogy of an unstructured biological population. Its limit behavior in the critical case is well studied for the Zolotarev-Slack regularity parameter α ∈ (0, 1]. We turn to the case of very heavy-tailed reproduction distribution α = 0 assuming that Zubkov's regularity condition holds with parameter β ∈ (0, ∞). Our main result gives a new asymptotic pattern for the reduced branching process conditioned on nonextinction during a long time interval.

We consider weak and strong survival for branching random walks on multigraphs with bounded degree. We prove that, at the strong critical value, the process dies out locally almost surely. We relate the weak critical value to a geometric parameter of the multigraph. For a large class of multigraphs (which enlarges the class of quasi-transitive or regular graphs), we prove that, at the weak critical value, the process dies out globally almost surely. Moreover, for the same class, we prove that the existence of a pure weak phase is equivalent to nonamenability. The results are extended to branching random walks on weighted graphs.

We study the susceptible-infected-susceptible model of the spread of an endemic infection. We calculate an exact expression for the mean number of transmissions for all values of the population size (N) and the infectivity. We derive the large-N asymptotic behavior for the infectivitiy below, above, and in the critical region. We obtain an analytical expression for the probability distribution of the number of transmissions, n, in the critical region. We show that this distribution has an n-3/2 singularity for small n and decays exponentially for large n. The exponent decreases with the distance from the threshold, diverging to ∞ far below and approaching 0 far above.