We investigate the probabilistic evolution of a near-critical bisexual branching process with mating depending on the number of couples in the population. We determine sufficient conditions which guarantee either the almost sure extinction of such a process or its survival with positive probability. We also establish some limiting results concerning the sequences of couples, females, and males, suitably normalized. In particular, gamma, normal, and degenerate distributions are proved to be limit laws. The results also hold for bisexual Bienaymé–Galton–Watson processes, and can be adapted to other classes of near-critical bisexual branching processes.

We give a general construction of sequential games among multiple players, as well as a construction of the composition of sequential games. We obtain new properties of the optimal class of win-by-k games, including closure under composition and independence between the winner of the game and the number of points played. We obtain new results on the asymptotic efficiency of the n-point, win-by-k games.

Ebrahimi (2001) proposes an interesting model for assessing a system's reliability, expressing the failure time of the system in terms of a deterioration process and covariates. He also provides illustrative examples and gives some properties of the model. In this note, we give conditions for negative ageing of the system's lifetime, and we correct one of his statements. Moreover, for the lifetimes of two systems of the same kind, some stochastic comparisons are presented.

In a defaultable market, an investor trades having only partial information about the behavior of the market. Taking into account the intraday stock movements, the risky asset prices are modelled by marked point processes. Their dynamics depend on an unobservable process, representing the amount of news reaching the market. This is a marked point process, which may have common jump times with the risky asset price processes. The problem of hedging a defaultable claim is studied. In order to discuss all these topics, in this paper we examine stochastic control problems using backward stochastic differential equations (BSDEs) and filtering techniques. The goal of this paper is to construct a sequence of functions converging to the value function, each of these is the unique solution of a suitable BSDE.

In this paper we study efficient simulation algorithms for estimating P(X›x), where X is the total time of a job with ideal time T that needs to be restarted after a failure. The main tool is importance sampling, where a good importance distribution is identified via an asymptotic description of the conditional distribution of T given X›x. If T≡t is constant, the problem reduces to the efficient simulation of geometric sums, and a standard algorithm involving a Cramér-type root, γ(t), is available. However, we also discuss an algorithm that avoids finding the root. If T is random, particular attention is given to T having either a gamma-like tail or a regularly varying tail, and to failures at Poisson times. Different types of conditional limit occur, in particular exponentially tilted Gumbel distributions and Pareto distributions. The algorithms based upon importance distributions for T using these asymptotic descriptions have bounded relative error as x→∞ when combined with the ideas used for a fixed t. Nevertheless, we give examples of algorithms carefully designed to enjoy bounded relative error that may provide little or no asymptotic improvement over crude Monte Carlo simulation when the computational effort is taken into account. To resolve this problem, an alternative algorithm using two-sided Lundberg bounds is suggested.

We study the weak convergence of the extremes of supercritical branching Lévy processes $\{\mathbb{X}_t, t \ge0\}$ whose spatial motions are Lévy processes with regularly varying tails. The result is drastically different from the case of branching Brownian motions. We prove that, when properly renormalized, $\mathbb{X}_t$ converges weakly. As a consequence, we obtain a limit theorem for the order statistics of $\mathbb{X}_t$.

The strong law of large numbers of the Marcinkiewicz–Zygmund type is established for the total population in a subcritical branching process with immigration. The moment convergence of the total population is obtained under appropriate moment conditions on the offspring distribution and the immigration distribution.

Consider a machine, which may or may not have a defect, and the probability q that this machine is defective is unknown. In order to determine whether the machine is defective, it is tested. On each test, the defect is found with probability p, if it has not been found yet. Performing n tests costs cn dollars and there is a fine of 1 dollar if there is a defect and it is not found on the tests. When should we stop testing, in order to minimize the cost?

This problem is treated in a minimax setting: we try to find a strategy that works well, even for ‘bad’ q's. It turns out that the minimax optimal stop rule can be unexpectedly complicated. For example, if p = 1/2 and cn = cn = 0.25n, then the optimal rule is to start by performing one test. If a defect is found we stop, otherwise we perform a second test. If a defect is found, then again we stop, else toss a coin and stop if this shows heads. If we still have not stopped, a third and last test is performed.

We consider the sum Sd of record values in a sequence of independent random variables that are uniformly distributed on 1,…,d. This sum can be interpreted as the total amount of time spent in record lifetimes in the standard renewal theoretic setup. We investigate the distributional limit of Sd and some related quantities as d→∞. Some explicit values are given for d=6, a case that can be interpreted as a simple game of chance.

Technology developed in a predecessor paper (Chen and Seneta (1996)) is applied to provide, in a unified manner, a sharpening of bivariate Bonferroni-type bounds on P(v1≥r, v2≥u) obtained by Galambos and Lee (1992; upper bound) and Chen and Seneta (1986; lower bound).

In a planar percolation model, faces of the underlying graph, as well as the sites and bonds, may be viewed as random elements. With this viewpoint, Whitney duality allows construction of a planar dual percolation model for each planar percolation model, which applies to mixed models with sites, bonds, and faces open or closed at random. Using self-duality for percolation models on the square lattice, information is obtained about the percolative region in the mixed model.

We consider the question of an optimal transaction between two investors to minimize their risks. We define a dynamic entropic risk measure using backward stochastic differential equations related to a continuous-time single jump process. The inf-convolution of dynamic entropic risk measures is a key transformation in solving the optimization problem.

This paper illustrates a simple method of deriving distributions for samples from cell populations in exponential growth. The distribution of the elapsed proportion of the DNA-synthetic phase of the cell cycle is derived and used to find the distribution of DNA content in a random sample of cells. Allowing for a normally-distributed measurement error, the resulting distribution agrees well with DNA distributions observed empirically. The statistical analysis of DNA distributions, and the interpretation of the distributional parameters in terms of cell population kinetics, are discussed and illustrated with an example.

In first-passage percolation (FPP) models, the passage time Tℓ from the origin to the point ℓeℓ satisfies f(ℓ) := ETℓ = μℓ + o(ℓ½+ε), where μ ∊ (0,∞) denotes the time constant. Yet, for lattice FPP, it is not known rigorously that f(ℓ) is eventually monotonically increasing. Here, for the Poisson-based Euclidean FPP of Howard and Newman (Prob. Theory Relat. Fields108 (1997), 153–170), we prove an explicit formula for f′(ℓ). In all dimensions, for certain values of the model's only parameter we have f′(ℓ) ≥ C > 0 for large ℓ.

Under appropriate long-range dependence conditions, it is well known that the joint distribution of the number of exceedances of several high levels is asymptotically compound Poisson. Here we investigate the structure of a cluster of exceedances for stationary sequences satisfying a suitable local dependence condition, under which it is only necessary to get certain limiting probabilities, easy to compute, in order to obtain limiting results for the highest order statistics, exceedance counts and upcrossing counts.

Most generalizations of the classical occupancy model involve non-homogeneous shot assignment probabilities, but retain the independence of the individual shot assignments. Hence, these models are associated with non-homogeneous Poisson processes. The present article discusses a generalization in which the shot assignments are not independent, but which result in clustering of the shots. Conditions are given under which this clustered occupancy model converges to a Poisson cluster process. Limiting distributions for the number of empty cells are obtained for various allocation intensities when the total number of shots is deterministic as well as random. In particular, it is shown that when the allocation is sparse, then the limiting distribution of the number of empty cells is compound Poisson.

We consider a percolation model on the d-dimensional Euclidean space which consists of spheres centred at the points of a Poisson point process of intensity ?. The radii of the spheres are random and are chosen independently and identically according to a distribution of a positive random variable. We show that the percolation function is continuous everywhere except perhaps at the critical point. Further, we show that the percolation functions converge to the appropriate percolation function except at the critical point when the radius random variables are uniformly bounded and converge weakly to another bounded random variable.

As an extension of a central limit theorem established by Svante Janson, we prove a Berry–Esseen inequality for a sum of independent and identically distributed random variables conditioned by a sum of independent and identically distributed integer-valued random variables.

We consider the problem of efficient simulation estimation of the density function at the tails, and the probability of large deviations for a sum of independent, identically distributed (i.i.d.), light-tailed, and nonlattice random vectors. The latter problem besides being of independent interest, also forms a building block for more complex rare event problems that arise, for instance, in queueing and financial credit risk modeling. It has been extensively studied in the literature where state-independent, exponential-twisting-based importance sampling has been shown to be asymptotically efficient and a more nuanced state-dependent exponential twisting has been shown to have a stronger bounded relative error property. We exploit the saddle-point-based representations that exist for these rare quantities, which rely on inverting the characteristic functions of the underlying random vectors. These representations reduce the rare event estimation problem to evaluating certain integrals, which may via importance sampling be represented as expectations. Furthermore, it is easy to identify and approximate the zero-variance importance sampling distribution to estimate these integrals. We identify such importance sampling measures and show that they possess the asymptotically vanishing relative error property that is stronger than the bounded relative error property. To illustrate the broader applicability of the proposed methodology, we extend it to develop an asymptotically vanishing relative error estimator for the practically important expected overshoot of sums of i.i.d. random variables.