The geometric distribution is characterized. The lack of memory property is replaced by the constancy of the conditional variance. Then the characterization is obtained.

Caching is an effective mechanism for reducing bandwidth usage and alleviating server load. However, the use of caching entails a compromise between content freshness and refresh cost. Constant refreshing allows a high degree of content freshness at a greater cost of system resource. Conversely, too little refreshing reduces content freshness but saves the cost of resource usage. To address the freshness-cost problem, we formulate the refresh scheduling problem with a generic cost model and use this cost model to determine an optimal refresh frequency that gives the best trade-off between refresh cost and content freshness. We prove the existence and uniqueness of an optimal refresh frequency under the assumptions that the arrival of updated content is Poisson and the age-related cost monotonically increases with decreasing freshness. In addition, we provide an analytic comparison of system performance under fixed refresh scheduling and random refresh scheduling, showing that two refresh schedulings with the same average refresh frequency are mathematically equivalent in terms of the long-run average cost.

We study the durations of the avalanches in the maximal avalanche decomposition of the Bak-Sneppen evolution model. We show that all the avalanches in this maximal decomposition have infinite expectation, but only ‘barely’, in the sense that if we made the appropriate threshold a tiny amount smaller (in a certain sense), then the avalanches would have finite expectation. The first of these results is somewhat surprising, since simulations suggest finite expectations.

In this paper various notions of positive and negative dependence for bivariate stochastic processes are introduced and their interrelationship is studied. Examples are given to illustrate these concepts.

The purpose of this paper is to determine the exact distribution of the final size of an epidemic for a wide class of models of susceptible–infective–removed type. First, a nonstationary version of the classical Reed–Frost model is constructed that allows us to incorporate, in particular, random levels of resistance to infection in the susceptibles. Then, a randomized version of this nonstationary model is considered in order to take into account random levels of infectiousness in the infectives. It is shown that, in both cases, the distribution of the final number of infected individuals can be obtained in terms of Abel–Gontcharoff polynomials. The new methodology followed also provides a unified approach to a number of recent works in the literature.

For simple random walk on the integers, consider the chance that the walk has traveled distance k from its start given that its first return is at time 2n. We derive a limiting approximation accurate to order 1/n. We give a combinatorial explanation for a functional equation satisfied by the limit and show this yields the functional equation of Riemann's zeta function.

Let T denote a subset of the possible transitions between the states of a finite Markov chain and let Yk denote the time of the kth occurrence of a T-transition. Formulae are derived for the generating functions of Yk, of Yj + k — Yj and of Yj + k — Yj in the limit as j → ∞, for both discrete-time and continuoustime chains. Several particular cases are briefly discussed.

We consider a large class of exponential random graph models and prove the existence of a region of parameter space corresponding to the emergent multipartite structure, separated by a phase transition from a region of disordered graphs. An essential feature is the formalism of graph limits as developed by Lovász et al. for dense random graphs.

We consider the problem of conditioning a non-explosive birth and death process to remain positive until time T, and consider weak convergence of this conditional process as T → ∞. By a suitable almost sure construction we prove weak convergence. The almost sure construction used is of independent interest but relies heavily on the strong monotonic properties of birth and death processes.

For correlated random walks a method of transition probability matrices as an alternative to the much-used methods of probability generating functions and difference equations has been investigated in this paper. To illustrate the use of transition probability matrices for computing the various probabilities for correlated random walks, the transition probability matrices for restricted/unrestricted one-dimensional correlated random walk have been defined and used to obtain some of the probabilities.

A discrete random variable describing the number of comparisons made in a sequence of comparisons between two opponents which terminates as soon as one opponent wins m comparisons is studied. By equating two different expressions for the mean of the variable, a closed form for the incomplete beta function with equal arguments is obtained. This expression is used in deriving asymptotic (m-large) expressions for the mean and variance. The standardized variate is shown to converge to the Gaussian distribution as m→ ∞. A result corresponding to the DeMoivre-Laplace limit theorem is proved. Finally applications are made to the genetic code problem, to Banach's Match Box Problem, and to the World Series of baseball.

This work is concerned with the existence of an optimal control strategy for the long-run average continuous control problem of piecewise-deterministic Markov processes (PDMPs). In Costa and Dufour (2008), sufficient conditions were derived to ensure the existence of an optimal control by using the vanishing discount approach. These conditions were mainly expressed in terms of the relative difference of the α-discount value functions. The main goal of this paper is to derive tractable conditions directly related to the primitive data of the PDMP to ensure the existence of an optimal control. The present work can be seen as a continuation of the results derived in Costa and Dufour (2008). Our main assumptions are written in terms of some integro-differential inequalities related to the so-called expected growth condition, and geometric convergence of the post-jump location kernel associated to the PDMP. An example based on the capacity expansion problem is presented, illustrating the possible applications of the results developed in the paper.

In this paper we propose a new theory and methodology to tackle the problem of unifying Monte Carlo samples from distributed densities into a single Monte Carlo draw from the target density. This surprisingly challenging problem arises in many settings (for instance, expert elicitation, multiview learning, distributed ‘big data’ problems, etc.), but to date the framework and methodology proposed in this paper (Monte Carlo fusion) is the first general approach which avoids any form of approximation error in obtaining the unified inference. In this paper we focus on the key theoretical underpinnings of this new methodology, and simple (direct) Monte Carlo interpretations of the theory. There is considerable scope to tailor the theory introduced in this paper to particular application settings (such as the big data setting), construct efficient parallelised schemes, understand the approximation and computational efficiencies of other such unification paradigms, and explore new theoretical and methodological directions.

The Asmussen–Kroese Monte Carlo estimators of P(Sn > u) and P(SN > u) are known to work well in rare event settings, where SN is the sum of independent, identically distributed heavy-tailed random variables X1,…,XN and N is a nonnegative, integer-valued random variable independent of the Xi. In this paper we show how to improve the Asmussen–Kroese estimators of both probabilities when the Xi are nonnegative. We also apply our ideas to estimate the quantity E[(SN-u)+].

In this paper we apply the recently established Wiener-Hopf Monte Carlo simulation technique for Lévy processes from Kuznetsov et al. (2011) to path functionals; in particular, first passage times, overshoots, undershoots, and the last maximum before the passage time. Such functionals have many applications, for instance, in finance (the pricing of exotic options in a Lévy model) and insurance (ruin time, debt at ruin, and related quantities for a Lévy insurance risk process). The technique works for any Lévy process whose running infimum and supremum evaluated at an independent exponential time can be sampled from. This includes classic examples such as stable processes, subclasses of spectrally one-sided Lévy processes, and large new families such as meromorphic Lévy processes. Finally, we present some examples. A particular aspect that is illustrated is that the Wiener-Hopf Monte Carlo simulation technique (provided that it applies) performs much better at approximating first passage times than a ‘plain’ Monte Carlo simulation technique based on sampling increments of the Lévy process.

We consider a drift parameter estimation problem when the state process is a reflected fractional Brownian motion (RFBM) with a nonzero drift parameter and the observation is the associated local time process. The RFBM process arises as the key approximating process for queueing systems with long-range dependent and self-similar input processes, where the drift parameter carries the physical meaning of the surplus service rate and plays a central role in the heavy-traffic approximation theory for queueing systems. We study a statistical estimator based on the cumulative local time process and establish its strong consistency and asymptotic normality.

The basic jump-telegraph process with exponentially distributed interarrival times deserves interest in various applied fields such as financial modelling and queueing theory. Aiming to propose a more general setting, we analyse such a stochastic process when the interarrival times separating consecutive velocity changes (and jumps) have generalized Mittag-Leffler distributions, and constitute the random times of a fractional alternating Poisson process. By means of renewal theory-based issues we obtain the forward and backward transition densities of the motion in series form, and prove their uniform convergence. Specific attention is then given to the case of jumps with constant size, for which we also obtain the mean of the process. Finally, we investigate the first-passage time of the process through a constant positive boundary, providing its formal distribution and suitable lower bounds.

We obtain an expression for the correlation of the maxima of two correlated Brownian motions.