This paper studies the dependence structure and bounds of several basic prototypical parallel queueing systems with correlated arrival processes to different queues. The marked feature of our systems is that each queue viewed alone is a standard single-server queuing system extensively studied in the literature, but those queues are statistically dependent due to correlated arrival streams. The major difficulty in analysing those systems is that the presence of correlation makes the explicit computation of a joint performance measure either intractable or computationally intensive. In addition, it is not well understood how and in what sense arrival correlation will improve or deteriorate a system performance measure.

The objective of this paper is to provide a better understanding of the dependence structure of correlated queueing systems and to derive computable bounds for the statistics of a joint performance measure. In this paper, we obtain conditions on arrival processes under which a performance measure in two systems can be compared, in the sense of orthant and supermodular orders, among different queues and over different arrival times. Such strong comparison results enable us to study both spatial dependence (dependence among different queues) and temporal dependence (dependence over different time instances) for a joint performance measure. Further, we derive a variety of upper and lower bounds for the statistics of a stationary joint performance measure. Finally, we apply our results to synchronized queueing systems, using the ideas combined from the theory of orthant and supermodular dependence orders and majorization with respect to weighted trees (Xu and Li (2000)). Our results reveal how a performance measure can be affected, favourably or adversely, by different types of dependencies.

We study the (two-sided) exit time and position of a random walk outside boundaries which are regularly varying functions of smaller order at infinity than the square root. A natural domain of interest is those random walks which are attracted without centring to a normal law, or are relatively stable. These are shown to have ‘stable’ exit positions, in that the overshoot of the curved boundary is of smaller order of magnitude (in probability) than the boundary, as the boundary expands. Surprisingly, this remains true regardless of the shape of the boundary. Furthermore, within the same natural domain of interest, norming of the exit position by, for example, the square root of the exit time (in the finite-variance case), produces limiting distributions which are computable from corresponding functionals of Brownian motion. We give a functional limit theorem for attraction of normed sums to general infinitely divisible random variables, as a means of making this, and more general, computations. These kinds of theorems have applications in sequential analysis, for example.

We show that the stationary distribution of a two-dimensional stochastic fluid network with (possibly dependent) Lévy inputs does not have product form other than in truly obvious cases. This is in contrast to queueing networks, where product form exists for non-obvious situations in which the inputs are independent, and for Brownian networks, where it typically exists for cases where the driving processes are actually dependent.

We study a reliability system subject to shocks generated by a renewal point process. When a shock occurs, components fail independently of each other with equal probabilities that are random numbers drawn from a distribution that may differ from shock to shock. We first consider the case of a parallel system and derive closed expressions for the Laplace-Stieltjes transform and the expectation of the time to system failure and for its density in the case that the distribution function of the renewal process possesses a density. We then treat a more general system structure, which has some very important special cases, such as k-out-of-n:F systems, and derive analogous formulae.

We investigate some preservation properties of two nonparametric classes of survival distributions and their duals, under appropriate reliability operations. The aging properties defining these nonparametric classes are based on comparing the mean life of a new unit to the mean residual life function of the asymptotic remaining survival time of the unit under repeated perfect repairs. They are motivated from a point of view that realistic notions of degradation, applicable to repairable systems, should be based on contrasting some aspect of the remaining life of a repairable unit (under a given repair strategy, such as renewals) to the life of a new unit.

We consider an interacting particle system in which each site of the d-dimensional integer lattice can be in state 0, 1, or 2. Our aim is to model the spread of disease in plant populations, so think of 0 = vacant, 1 = healthy plant, 2 = infected plant. A vacant site becomes occupied by a plant at a rate which increases linearly with the number of plants within range R, up to some saturation level, F 1, above which the rate is constant. Similarly, a plant becomes infected at a rate which increases linearly with the number of infected plants within range M, up to some saturation level, F 2. An infected plant dies (and the site becomes vacant) at constant rate δ. We discuss coexistence results in one and two dimensions. These results depend on the relative dispersal ranges for plants and disease.

We consider a fully connected queueing network in which customers have one direct and many alternative routes through the network, and where customer routeing is dynamic. We obtain an asymptotically optimal routeing policy, taking the limit as the number of queues of the network increases. We observe that good policies route customers directly, unless there is a danger of servers becoming idle, in which case customers should be routed alternatively so as to avoid such idleness, and this leads to policies that perform well in moderate-sized networks.

We formulate and verify an almost-sure lattice renewal theorem for branching random walks, whose non-lattice analogue is originally due to Nerman. We also identify the limit in these renewal theorems (both lattice and non-lattice) as the limit of Kingman's well-known martingale multiplied by a deterministic factor.

In this article we introduce generalizations of several well-known reliability bounds. These bounds are based on arbitrary partitions of the family of minimal path or cut sets of the system and can be used for approximating the reliability of any coherent structure with i.i.d. components. An illustration is also given of how the general results can be applied for a specific reliability structure (two-dimensional consecutive-k 1 x k 2-out-of-n 1x n 2 system) along with extensive numerical calculations revealing that, in most cases, the generalized bounds perform better than other available bounds in the literature for this system.

Certain types of Markov jump processes x(t) with continuous state space and one or more absorbing states are studied. Cases where the transition rate in state x is of the form λ(x) = |x|δ in a neighbourhood of the origin in ℝd are considered, in particular. This type of problem arises from quantum physics in the study of laser cooling of atoms, and the present paper connects to recent work in the physics literature. The main question addressed is that of the asymptotic behaviour of x(t) near the origin for large t. The study involves solution of a renewal equation problem in continuous state space.

In this paper we investigate the application of perfect simulation, in particular Coupling from the Past (CFTP), to the simulation of random point processes. We give a general formulation of the method of dominated CFTP and apply it to the problem of perfect simulation of general locally stable point processes as equilibrium distributions of spatial birth-and-death processes. We then investigate discrete-time Metropolis-Hastings samplers for point processes, and show how a variant which samples systematically from cells can be converted into a perfect version. An application is given to the Strauss point process.

Let ε be a (single or composite) pattern defined over a sequence of Bernoulli trials. This article presents a unified approach for the study of the joint distribution of the number S n of successes (and F n of failures) and the number X n of occurrences of ε in a fixed number of trials as well as the joint distribution of the waiting time T r till the rth occurrence of the pattern and the number S T r of successes (and F T r of failures) observed at that time. General formulae are developed for the joint probability mass functions and generating functions of (X n ,S n ), (T r ,S T r ) (and (X n ,S n ,F n ),(T r ,S T r ,F T r )) when X n belongs to the family of Markov chain imbeddable variables of binomial type. Specializing to certain success runs, scans and pattern problems several well-known results are delivered as special cases of the general theory along with some new results that have not appeared in the statistical literature before.

We argue the importance both of developing simple sufficient conditions for the stability of general multiclass queueing networks and also of assessing such conditions under a range of assumptions on the weight of the traffic flowing between service stations. To achieve the former, we review a peak-rate stability condition and extend its range of application and for the latter, we introduce a generalisation of the Lu–Kumar network on which the stability condition may be tested for a range of traffic configurations. The peak-rate condition is close to exact when the between-station traffic is light, but degrades as this traffic increases.

This paper studies the connection between the dynamical and equilibrium behaviour of large uncontrolled loss networks. We consider the behaviour of the number of calls of each type in the network, and show that, under the limiting regime of Kelly (1986), all trajectories of the limiting dynamics converge to a single fixed point, which is necessarily that on which the limiting stationary distribution is concentrated. The approach uses Lyapunov techniques and involves the evolution of the transition rates of a stationary Markov process in such a way that it tends to reversibility.

Miyazawa and Taylor (1997) introduced a class of assemble-transfer batch service queueing networks which do not have tractable stationary distribution. However by assuming a certain additional arrival process at each node when it is empty, they obtain a geometric product-form stationary distribution which is a stochastic upper bound for the stationary distribution of the original network. In this paper we develop a stochastic lower bound for the original network by introducing an additional departure process at each node which tends to remove all the customers present in it. This model in combination with the aforementioned upper bound model gives a better sense for the properties of the original network.

Two related individuals are identical by descent at a genetic locus if they share the same gene copy at that locus due to inheritance from a recent common ancestor. We consider idealized continuous identity by descent (IBD) data in which IBD status is known continuously along chromosomes. IBD data contains information about the relationship between the two individuals, and about the underlying crossover processes. We present a Monte Carlo method for calculating probabilities for IBD data. The method is not restricted to Haldane's Poisson process model of crossing-over but may be used with other models including the chi-square, Kosambi renewal and Sturt models. Results of a simulation study demonstrate that IBD data can be used to distinguish between alternative models for the crossover process.

A unifying technology is introduced for finding explicit closed form expressions for joint moment generating functions of various random quantities associated with some waiting time problems. Sooner and later waiting times are covered for general discrete- and continuous-time models. The models are either Markov chains or semi-Markov processes with a finite number of states. Waiting times associated with generalized phase-type distributions, that are of interest in survival analysis and other areas, are also covered.

We derive some asymptotic results for the rate of convergence to equilibrium for the number of busy servers in an M/M/N/N queue with input rate λN and service rate 1 for N → ∞ in the ‘subcritical’ case λ ∈]0, 1[. These results improve recent contributions of Fricker, Robert and Tibi.

We investigate zero-temperature dynamics for a homogeneous ferromagnetic Ising model on the homogeneous tree of degree three (𝕋) with random (i.i.d. Bernoulli) spin configuration at time 0. Letting θ denote the probability that any particular vertex has a +1 initial spin, for infinite spin clusters do not exist at time 0 but we show that infinite ‘spin chains’ (doubly infinite paths of vertices with a common spin) exist in abundance at any time ϵ > 0. We study the structure of the subgraph of 𝕋 generated by the vertices in time-ϵ spin chains. We show the existence of a phase transition in the sense that, for some critical θc with spin chains almost surely never form for θ < θc but almost surely do form in finite time for θ > θc . We relate these results to certain quantities of physical interest, such as the t → ∞ asymptotics of the probability that any particular vertex changes spin after time t.