The two-dimensional periodogram has been proposed as an estimator of the spectral density of a real, homogeneous, random field defined over a regular lattice on the plane. In the present paper, results pertaining to the asymptotic distributional properties of such a periodogram are obtained. These results generalize some of the work of Hannan (1960), Walker (1965) and more directly Olshen (1967a,b) concerning asymptotic theory for the periodogram of a stationary time series. Although extension of asymptotic theory for one-dimensional periodograms to a parallel theory for two-dimensional periodograms is not completely straightforward (one runs into problems akin to the problems encountered in extending the theory of one-dimensional trigonometric series to two dimensions), further extensions to asymptotic theory for p-dimensional periodograms (p > 2) is easily accomplished by an obvious mimicry of the definitions, theorems and proofs for the two-dimensional case. Since the notation required for the p-dimensional case is rather cumbersome, we have chosen to give results only for two-dimensional periodograms.

We consider I fluid queues in parallel. Each fluid queue has a deterministic inflow with a constant rate. At a random instant subject to a Poisson process, random amounts of fluids are simultaneously reduced. The requested amounts for the reduction are subject to a general I-dimensional distribution. The queues with inventories that are smaller than the requests are emptied. Stochastic upper bounds are considered for the stationary distribution of the joint buffer contents. Our major interest is in finding exponential product-form bounds, which turn out to have the appropriate decay rates with respect to certain linear combinations of buffer contents.

A method is provided for numerical evaluation, with any given accuracy, of the probability that at least p% of the genetic material from an individual's chromosomal segment survives to the next generation. Relevant MAPLE® V codes, for automated implementation of such evaluation, are also provided. The genomic continuum model, with Haldane's model for the crossover process, is assumed.

We derive an integral expression for the expected coverage of an n-dimensional target when the value density function of the target, the damage function, and the distribution of the damage function center around the target center are spherically symmetric. This integral expression is transformed by means of Parseval's theorem for Hankel transforms to an integral of the product of the Hankel transforms of order ½n − 1 of all underlying functions.

We consider a multicommodity flow problem on a complete graph whose edges have random, independent, and identically distributed capacities. We show that, as the number of nodes tends to infinity, the maximum utility, given by the average of a concave function of each commodity flow, has an almost-sure limit. Furthermore, the asymptotically optimal flow uses only direct and two-hop paths, and can be obtained in a distributed manner.

Under mild assumptions, we show that the exact convergence rate in total variation is also exact in weaker Wasserstein distances for the Metropolis–Hastings independence sampler. We develop a new upper and lower bound on the worst-case Wasserstein distance when initialized from points. For an arbitrary point initialization, we show that the convergence rate is the same and matches the convergence rate in total variation. We derive exact convergence expressions for more general Wasserstein distances when initialization is at a specific point. Using optimization, we construct a novel centered independent proposal to develop exact convergence rates in Bayesian quantile regression and many generalized linear model settings. We show that the exact convergence rate can be upper bounded in Bayesian binary response regression (e.g. logistic and probit) when the sample size and dimension grow together.

Let $(A_i)_{i \geq 0}$ be a finite-state irreducible aperiodic Markov chain and f a lattice score function such that the average score is negative and positive scores are possible. Define $S_0\coloneqq 0$ and $S_k\coloneqq \sum_{i=1}^k f(A_i)$ the successive partial sums, $S^+$ the maximal non-negative partial sum, $Q_1$ the maximal segmental score of the first excursion above 0, and $M_n\coloneqq \max_{0\leq k\leq\ell\leq n} (S_{\ell}-S_k)$ the local score, first defined by Karlin and Altschul (1990). We establish recursive formulae for the exact distribution of $S^+$ and derive a new approximation for the tail behaviour of $Q_1$ , together with an asymptotic equivalence for the distribution of $M_n$ . Computational methods are explicitly presented in a simple application case. The new approximations are compared with those proposed by Karlin and Dembo (1992) in order to evaluate improvements, both in the simple application case and on the real data examples considered by Karlin and Altschul (1990).

We consider a Lévy process Y(t) that is not continuously observed, but rather inspected at Poisson( $\omega$ ) moments only, over an exponentially distributed time $T_\beta$ with parameter $\beta$ . The focus lies on the analysis of the distribution of the running maximum at such inspection moments up to $T_\beta$ , denoted by $Y_{\beta,\omega}$ . Our main result is a decomposition: we derive a remarkable distributional equality that contains $Y_{\beta,\omega}$ as well as the running maximum process $\bar Y(t)$ at the exponentially distributed times $T_\beta$ and $T_{\beta+\omega}$ . Concretely, $\overline{Y}(T_\beta)$ can be written as the sum of two independent random variables that are distributed as $Y_{\beta,\omega}$ and $\overline{Y}(T_{\beta+\omega})$ . The distribution of $Y_{\beta,\omega}$ can be identified more explicitly in the two special cases of a spectrally positive and a spectrally negative Lévy process. As an illustrative example of the potential of our results, we show how to determine the asymptotic behavior of the bankruptcy probability in the Cramér–Lundberg insurance risk model.

This paper builds a mixture representation of the reliability function of the conditional residual lifetime of a coherent system in terms of the reliability functions of conditional residual lifetimes of order statistics. Some stochastic ordering properties for the conditional residual lifetime of a coherent system with independent and identically distributed components are obtained, based on the stochastically ordered coefficient vectors.

Consider an absolutely continuous distribution on [0, ∞) with finite mean μ and hazard rate function h(t) ≤ b for all t. For bμ close to 1, we would expect F to be approximately exponential. In this paper we obtain sharp bounds for the Kolmogorov distance between F and an exponential distribution with mean μ, as well as between F and an exponential distribution with failure rate b. We apply these bounds to several examples. Applications are presented to geometric convolutions, birth and death processes, first-passage times, and to decreasing mean residual life distributions.

Two classes of finite and infinite moving-average sequences of bivariate random vectors are considered. The first class has bivariate exponential marginals while the second class has bivariate geometric marginals. The theory of positive dependence is used to show that in various cases the two classes consist of associated random variables. Association is then applied to establish moment inequalities and to obtain approximations to some joint probabilities of the bivariate processes.

For the delayed Bernoulli feedback queue with first come–first served discipline under weak assumptions a relationship for the generating functions of the joint queue-length distribution at various points in time is given. A decomposition for the generating function of the stationary total queue length distribution has been proven. The Laplace-Stieltjes transform of the stationary joint workload distribution function is represented by its marginal distributions. The arrival process is Poisson, renewal or arbitrary stationary, respectively. The service times can form an i.i.d. sequence at each queue. Different kinds of product form of the generating function of the joint queue-length distribution are discussed.

We are concerned with the variation of the supercritical nearest-neighbours contact process such that first infection occurs at a lower rate; it is known that the process survives with positive probability. Regarding the rightmost infected of the process started from one site infected and conditioned to survive, we specify a sequence of space-time points at which its behaviour regenerates and, thus, obtain the corresponding strong law and central limit theorem. We also extend complete convergence in this case.

Suppose that π is a policy for resource allocation in a stochastic environment and π ∗ is an optimal policy. Two existing procedures for policy evaluation are described and compared. Both of these evaluate π by means of upper bounds on R(π ∗) – R(π), the total reward lost when making resource allocations according to π rather than π∗. The bounds developed by these two methods are called Type 1 and Type 2. We demonstrate by example that neither of these procedures dominates the other in the sense of always yielding tighter bounds. A modification to Type 2 bounds is proposed resulting in an improved procedure which always dominates the Type 1 approach.

We introduce a definition of long range dependence of random processes and fields on an (unbounded) index space $T\subseteq \mathbb{R}^d$ in terms of integrability of the covariance of indicators that a random function exceeds any given level. This definition is specifically designed to cover the case of random functions with infinite variance. We show the value of this new definition and its connection to limit theorems via some examples including subordinated Gaussian as well as random volatility fields and time series.

Let {Zn}n≥0 be a random walk with a negative drift and independent and identically distributed increments with heavy-tailed distribution, and let M = supn≥0Zn be its supremum. Asmussen and Klüppelberg (1996) considered the behavior of the random walk given that M > x for large x, and obtained a limit theorem, as x → ∞, for the distribution of the quadruple that includes the time τ = τ(x) to exceed level x, position Zτ at this time, position Zτ-1 at the prior time, and the trajectory up to it (similar results were obtained for the Cramér-Lundberg insurance risk process). We obtain here several extensions of this result to various regenerative-type models and, in particular, to the case of a random walk with dependent increments. Particular attention is given to describing the limiting conditional behavior of τ. The class of models includes Markov-modulated models as particular cases. We also study fluid models, the Björk-Grandell risk process, give examples where the order of τ is genuinely different from the random walk case, and discuss which growth rates are possible. Our proofs are purely probabilistic and are based on results and ideas from Asmussen, Schmidli and Schmidt (1999), Foss and Zachary (2002), and Foss, Konstantopoulos and Zachary (2007).

The standard uniformization technique for continuous-time Markov chains is generalized to non-exponential stochastic networks.

An anticipation process is a continuous-time Markov process, having an element of deterministic drift in its transition structure. The theory of such processes, developed here from an elementary (Chapman-Kolmogorov) point of view, draws together a number of threads of theoretical and applied probability.

Suppose that the state variables x = (x1,…,xn)′ where the dot refers to derivatives with respect to time t, and u ∊ U is a vector of controls. The object is to transfer x to x1 by choosing the controls so that the functional takes on its minimum value J(x) called the Bellman function (although we shall define it in a different way). The Dynamic Programming Principle leads to the maximisation with respect to u of and equality is obtained upon maximisation.

We consider the two-dimensional process {X(t), V(t)} where {V(t)} is Brownian motion with drift, and {X(t)} is its integral. In this note we derive the joint density function of T and V(T) where T is the time at which the process {X(t)} first returns to its initial value. A series expansion of the marginal density of T is given in the zero-drift case. When V(0) and the drift are both positive there is a positive probability that {Χ (t)} never returns to its initial value. We show how this probability grows for small drift. Finally, using the Kontorovich–Lebedev transform pair we obtain the escape probability explicitly for arbitrary values of the drift parameter.