An approximate expression is found for the probability density of the recurrence time between major outbreaks of a host-vector disease such as yellow fever, in which the infectious period of a host is short and the lifetime of a vector is short compared with that of a host. The typical interval between epidemics based on the modal value of this distribution is obtained. These results are the analogues of those obtained by Bartlett (1956), (1966) for the periodicity of measles epidemics in small communities.

Let A1, A2, ···, An be events on a given probability space and let Br, n be the event that exactly r of the A's occur. Let further Sk (n) be the kth binomial moment of the number of the A's which occur. A sufficient condition is given for the existence of lim P (Br,n), as n→ + ∞, in terms of limits of the Sk(n)'s and a formula is given for the limit above. This formula for the limit is similar to the sieve theorem of Takács (1967) for infinite sequences of events and in the proof we make use of Takács's analytic method. The result is immediately applicable to the limit distribution of the maximum of (dependent) random variables X1, X2, ···, Xn by choosing Aj = {Xj ≧ x}. Our main theorem is reformulated for this special case and an example is given for illustration.

A solution is obtained for the signal-noise problem X = M + Z in which M and Z are independent Gaussian processes. Conditions on the processes are given which insure that the best estimate under generalized square — error loss is the conditional mean of M given X = x. A sequence of approximators of the best estimate is also given.

The substability and the ergodicity of various queueing models are discussed. This paper considers the vector-valued queueing process Xrn = (xrn,1,xrn,2, · ··) with non-negative components and a constant initial value Xrr= a. For this, the substability is derived under simple conditions by showing the finiteness of With respect to the ergodicity, in order to make use of Borovkov's theorem, we additionally assume that the distribution of the interarrival of customers has non-bounded tail for any given past sequence.

A general branching process is extended to allow life length and reproduction probabilities to depend on randomly changing environmental states. First and second moments of the population size with respect to time are derived assuming that the environmental process is semi-Markovian. The results are similar to Markov-renewal type equations which allow, under discrete time, iterative computation of both moments of the population through time.

Stochastic dynamical system models for non-linear random vibrations are considered and their discrete-time version, non-linear time series models are introduced using the local Gaussian modelling method. Some computational problems and implications of the present method in non-linear time series analysis are discussed.

A doubly stochastic switching model previously used by the author to check against doubly stochastic ‘white-noise' models is generalized in the case of Poisson processes to any distribution of fluctuating intensity rates. The resulting Cox process is used in place of a simple Poisson process as the basis of a mixed model, which is then fitted to some data on automobile accidents previously analysed by Seal (1980) in terms of a mixed Poisson process. It is shown that there is just significant evidence in favour of some additional heterogeneity ‘within' each year for individual drivers; the nature of the data analysed (yearly time intervals) does not, however, allow discrimination between heterogeneity which is definitely attributable to a common fluctuating environment, and that which is in effect independent for each driver.

We prove a ratio identity giving a transform of the distribution of the position of a two barrier random walk as a ratio of the equivalent transforms for the one barrier random walk. This ratio identity holds for skip free random walks where the increments are dependent on both the previous increment and also the position of the random walk.

Consider a random graph where the mean degree is given and fixed. In this paper we derive the maximal size of the largest connected component in the graph. We also study the related question of the largest possible outbreak size of an epidemic occurring ‘on’ the random graph (the graph describing the social structure in the community). More precisely, we look at two different classes of random graphs. First, the Poissonian random graph in which each node i is given an independent and identically distributed (i.i.d.) random weight Xi with E(Xi)=µ, and where there is an edge between i and j with probability 1-e-XiXj/(µ n), independently of other edges. The second model is the thinned configuration model in which the n vertices of the ground graph have i.i.d. ground degrees, distributed as D, with E(D) = µ. The graph of interest is obtained by deleting edges independently with probability 1-p. In both models the fraction of vertices in the largest connected component converges in probability to a constant 1-q, where q depends on X or D and p. We investigate for which distributions X and D with given µ and p, 1-q is maximized. We show that in the class of Poissonian random graphs, X should have all its mass at 0 and one other real, which can be explicitly determined. For the thinned configuration model, D should have all its mass at 0 and two subsequent positive integers.

We study a collector's problem with K renewal arrival processes for different type items, where the objective is to collect complete sets. In particular, we derive the asymptotic distribution of the sequence of interarrival times between set completions.

In this paper we study the tail behavior of the maximum exceedance of a sequence of independent and identically distributed random variables over a random walk. For both light-tailed and heavy-tailed cases, we derive a precise asymptotic formula, which extends and unifies some existing results in the recent literature of applied probability.

We consider backward stochastic differential equations in a setting where noise is generated by a countable state, continuous time Markov chain, and the terminal value is prescribed at a stopping time. We show that, given sufficient integrability of the stopping time and a growth bound on the terminal value and BSDE driver, these equations admit unique solutions satisfying the same growth bound (up to multiplication by a constant). This holds without assuming that the driver is monotone in y, that is, our results do not require that the terminal value be discounted at some uniform rate. We show that the conditions are satisfied for hitting times of states of the chain, and hence present some novel applications of the theory of these BSDEs.

We consider a class of optimal stopping problems involving both the running maximum as well as the prevailing state of a linear diffusion. Instead of tackling the problem directly via the standard free boundary approach, we take an alternative route and present a parameterized family of standard stopping problems of the underlying diffusion. We apply this family to delineate circumstances under which the original problem admits a unique, well-defined solution. We then develop a discretized approach resulting in a numerical algorithm for solving the considered class of stopping problems. We illustrate the use of the algorithm in both a geometric Brownian motion and a mean reverting diffusion setting.

We propose a stochastic model for the failure times of items subject to two external random shocks occurring as events in an underlying bivariate counting process. This is a special formulation of the competing risks model, which is of interest in reliability theory and survival analysis. Specifically, we assume that a system, or an item, fails when the sum of the two types of shock reaches a critical random threshold. In detail, the two kinds of shock occur according to a bivariate space-fractional Poisson process, which is a two-dimensional vector of independent homogeneous Poisson processes time-changed by an independent stable subordinator. Various results are given, such as analytic hazard rates, failure densities, the probability that the failure occurs due to a specific type of shock, and the survival function. Some special cases and ageing notions related to the NBU characterization are also considered. In this way we generalize certain results in the literature, which can be recovered when the underlying process reduces to the homogeneous Poisson process.

We consider two random sequential packing processes in which spheres of unit radius are randomly attached to the surface of a fixed unit sphere. Independent random spheres are generated and added successively, provided there is no overlap with previous spheres. In model 1, the process stops when a trial sphere intersects one of the previously-accepted spheres. In model 2, random sequential packing, any such overlapping trial sphere is discarded and the next random sphere is tried, until it is impossible to add any further spheres.

Previous workers have conjectured convincingly that no exact analytical solution is possible for this type of problem. We use Monte Carlo simulation methods to estimate transition probabilities for the two models. Because some probabilities are extremely small, a simulation using independent repetitions of the model would be inefficient. We designed a branching process of conditionally binomial trials, and performed over 108 trials on a supercomputer.

Excess life distributions for discrete renewal processes may be computed by using elementary discrete Markov chain concepts involving absorption probabilities. Excess life distributions in general may then be obtained by approximating the renewal process under study by a suitably chosen sequence of discrete renewal processes. The technique is illustrated in the cases of renewal processes with interarrival distributions which are a linear combination of two exponentials and uniform [0, 1]. A related algorithm is described for computer generated approximations of excess life distributions corresponding to continuous interarrival time distributions with an increasing c.d.f. Conditions for convergence of this algorithm are examined.

We introduce the concept of the absorbing process for analysing a state process. Our aim is to show the existence of the absorbing process with probability one. This process is shown to be stationary, asymptotically stationary, periodic or a.m.s., if the input distribution has such properties. The real process is absorbed into this process so that its stability and some other properties are easily derived.

The concept of θ-stationarity for a simple second-order point process in R2 is introduced. This concept is closely related to that of isotropy. Some θ-stationary processes are defined. Techniques are given for simulating realisations of these processes. The second-order analysis of these processes which have an obvious point of reference or origin is considered. Methods are suggested for modelling spatial patterns which are realisations of such processes. These methods are illustrated using simulated data. The ideas are extended to multitype point processes.

In [14] a necessary and sufficient condition was obtained for there to exist uniquely a Q-process with a specified invariant measure, under the assumption that Q is a stable, conservative, single-exit matrix. The purpose of this note is to demonstrate that, for an arbitrary stable and conservative q-matrix, the same condition suffices for the existence of a suitable Q-process, but that this process might not be unique. A range of examples is considered, including pure-birth processes, a birth process with catastrophes, birth-death processes and the Markov branching process with immigration.

Modeling of manufacturing lines and data communications makes frequent use of tandem queues with blocking. Here we present and study such a system with a k-stage blocking scheme in which processing in each station requires attendance of the server of that station together with the servers of the next k – 1 stations. This scheme describes the conventional manufacturing and communications blocking schemes but is also representative of a wider range of applications. Explicit expressions for residence times, throughput, queueing times, waiting times and other measures of performance are obtained for the case of communications type flow and just-in-time input. Use of the results in modeling and analyzing a highway traffic flow situation is presented.