We consider the likelihood ratio tests to detect an epidemic alternative in the following two cases of normal observations: (1) the alternative specifies a square wave drift in the mean value of an i.i.d. sequence; (2) the alternative permits a square wave drift in the intercept of a simple linear regression model. To develop the approximations for the significance levels leads us to consider boundary-crossing problems of some two-dimensional discrete-time Gaussian fields. By the method which was proposed originally by Woodroofe (1976) and adapted to study maxima of some random fields by Siegmund (1988), some large deviations for the conditional non-linear boundary-crossing probabilities are developed. Some results of Monte Carlo experiments confirm the accuracy of these approximations.

In this paper the class of cyclostationary Gaussian random processes is studied. Basic asymptotics are given for the class of Gaussian processes that are centered and differentiable in mean square. Then, under certain conditions on the non-degeneration of the centered cyclostationary Gaussian process with integrable covariance functions, the Gnedenko-type limit formula is established for and all x > 0.

The well-known Cameron–Martin formula allows us to calculate the mathematical expectation where Ws is a Wiener process. This paper extends this result to the case of piecewise continuous martingales. As a particular case the mathematical expectations of a functional of generalized Ornstein– Uhlenbeck processes and pure jump processes are calculated.

We give explicit expressions for the Slepian model process of non-stationary Gaussian processes following level crossings and local maxima. We also include a detailed analysis of the high-level case.

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.

Let γ t and δ t denote the residual life at t and current life at t, respectively, of a renewal process , with the sequence of interarrival times. We prove that, given a function G, under mild conditions, as long as holds for a single positive integer n, then is a Poisson process. On the other hand, for a delayed renewal process with the residual life at t, we find that for some fixed positive integer n, if is independent of t, then is an arbitrarily delayed Poisson process. We also give some corresponding results about characterizing the common distribution function F of the interarrival times to be geometric when F is discrete. Finally, we obtain some characterization results based on the total life or independence of γ t and δ t.

This article provides the theoretical basis of the virtual customer method or positive rare perturbation (RPA) method of sensitivity analysis, and in particular gives a short proof of the light traffic derivative result of Reiman and Simon [5] based on Campbell's formula. As a by-product, we obtain the archetypal H = λG formula associated with a stationary quantity of a queueing system.

By using the alternating projection theorem of J. von Neumann, we obtain explicit formulae for the best linear interpolator and interpolation error of missing values of a stationary process. These are expressed in terms of multistep predictors and autoregressive parameters of the process. The key idea is to approximate the future by a finite-dimensional space.

A unified exposition of random Johnson–Mehl tessellations in d-dimensional Euclidean space is presented. In particular, Johnson-Mehl tessellations generated by time-inhomogeneous Poisson processes and nucleation-exclusion models are studied. The ‘practical' cases d = 2 and d = 3 are discussed in detail. Several new results are established, including first- and second-order moments of various characteristics for both Johnson–Mehl tesselations and sectional Johnson–Mehl tessellations.

This paper considers sets of points from a Poisson process in the plane, chosen to be close together, and their properties. In particular, the perimeter of the convex hull of such a point set is investigated. A number of different models for the selection of such points are considered, including a simple nearest-neighbour model. Extensions to marked processes and applications to modelling animal territories are discussed.

It is well known that most commonly used discrete distributions fail to belong to the domain of maximal attraction for any extreme value distribution. Despite this negative finding, C. W. Anderson showed that for a class of discrete distributions including the negative binomial class, it is possible to asymptotically bound the distribution of the maximum. In this paper we extend Anderson's result to discrete-valued processes satisfying the usual mixing conditions for extreme value results for dependent stationary processes. We apply our result to obtain bounds for the distribution of the maximum based on negative binomial autoregressive processes introduced by E. McKenzie and Al-Osh and Alzaid. A simulation study illustrates the bounds for small sample sizes.

We investigate the asymptotic behaviour of empirical processes truncated outside an interval about the (1 – s(n)/n)-quantile where s(n) → ∞ and s(n)/n → 0 as the sample size n tends to ∞. It is shown that extreme value (Poisson) processes and, alternatively, the homogeneous Poisson process may serve as approximations if certain von Mises conditions hold.

We show that using the FIFO service discipline at single server stations with ILR (increasing likelihood ratio) service time distributions in networks of monotone queues results in stochastically earlier departures throughout the network. The converse is true at stations with DLR (decreasing likelihood ratio) service time distributions. We use these results to establish the validity of the following comparisons:

(i) The throughput of a closed network of FIFO single-server queues will be larger (smaller) when the service times are ILR (DLR) rather than exponential with the same means.

(ii) The total stationary number of customers in an open network of FIFO single-server queues with Poisson external arrivals will be stochastically smaller (larger) when the service times are ILR (DLR) rather than exponential with the same means.

We also give a surprising counterexample to show that although FIFO stochastically maximizes the number of departures by any time t from an isolated single-server queue with IHR (increasing hazard rate, which is weaker than ILR) service times, this is no longer true for networks of more than one queue. Thus the ILR assumption cannot be relaxed to IHR.

Finally, we consider multiclass networks of exponential single-server queues, where the class of a customer at a particular station determines its service rate at that station, and show that serving the customer with the highest service rate (which is SEPT — shortest expected processing time first) results in stochastically earlier departures throughout the network, among all preemptive work-conserving policies. We also show that a cµ rule stochastically maximizes the number of non-defective service completions by any time t when there are random, agreeable, yields.

We derive an asymptotic relation for the expected value of the stopping time , where Sk is a random walk with negative drift.

The paper gives a spectral representation for a class of random fields which are bounded in mean square almost surely. A characterisation of the corresponding spectral measure in the representation is obtained based on Beurling's duality theory and generalised Fourier transforms. A representation for the covariance function of asymptotically stationary random fields is also derived.

The Hammersley–Clifford theorem gives the form that the joint probability density (or mass) function of a Markov random field must take. Its exponent must be a sum of functions of variables, where each function in the summand involves only those variables whose sites form a clique. From a statistical modeling point of view, it is important to establish the converse result, namely, to give the conditional probability specifications that yield a Markov random field. Besag (1974) addressed this question by developing a one-parameter exponential family of conditional probability models. In this article, we develop new models for Markov random fields by establishing sufficient conditions for the conditional probability specifications to yield a Markov random field.

This paper considers the undershoot of a general continuous-time risk process with dependent increments under a certain initial level. The increments are given by the locations and amounts of claims which are described by a stationary marked point process. Under a certain balance condition, it is shown that the distribution of the undershoot depends only on the mark distribution and on the intensity of the underlying point process, but not on the form of its distribution. In this way an insensitivity property is extended which has been proved in Björk and Grandell [3] for the ruin probability, i.e. for the probability that after a finite time interval the initial level will be crossed from above.

Let ψ be a functional of the sample path of a stochastic system driven by a Poisson process with rate λ . It is shown in a very general setting that the expectation of ψ,E λ [ψ], is an analytic function of λ under certain moment conditions. Instead of following the straightforward approach of proving that derivatives of arbitrary order exist and that the Taylor series converges to the correct value, a novel approach consisting in a change of measure argument in conjunction with absolute monotonicity is used. Functionals of non-homogeneous Poisson processes and Wiener processes are also considered and applications to light traffic derivatives are briefly discussed.

The dependence of coincidence of the global, local and pairwise Markov properties on the underlying undirected graph is examined. The pairs of these properties are found to be equivalent for graphs with some small excluded subgraphs. Probabilistic representations of the corresponding conditional independence structures are discussed.

The use of taboo probabilities in Markov chains simplifies the task of calculating the queue-length distribution from data recording customer departure times and service commencement times such as might be available from automatic bank-teller machine transaction records or the output of telecommunication network nodes. For the case of Poisson arrivals, this permits the construction of a new simple exact O(n 3) algorithm for busy periods with n customers and an O(n 2 log n) algorithm which is empirically verified to be within any prespecified accuracy of the exact algorithm. The algorithm is extended to the case of Erlang-k interarrival times, and can also cope with finite buffers and the real-time estimates problem when the arrival rate is known.