The multi-level Monte Carlo method proposed by Giles (2008) approximates the expectation of some functionals applied to a stochastic process with optimal order of convergence for the mean-square error. In this paper a modified multi-level Monte Carlo estimator is proposed with significantly reduced computational costs. As the main result, it is proved that the modified estimator reduces the computational costs asymptotically by a factor (p / α)2 if weak approximation methods of orders α and p are applied in the case of computational costs growing with the same order as the variances decay.

We consider the stochastic Allen-Cahn equation perturbed by smooth additive Gaussian noise in a spatial domain with smooth boundary in dimension d ≤ 3, and study the semidiscretization in time of the equation by an implicit Euler method. We show that the method converges pathwise with a rate O(Δtγ) for any γ < ½. We also prove that the scheme converges uniformly in the strong Lp-sense but with no rate given.

We consider the time behaviour associated to the sequential Monte Carlo estimate of the backward interpretation of Feynman-Kac formulae. This is particularly of interest in the context of performing smoothing for hidden Markov models. We prove a central limit theorem under weaker assumptions than adopted in the literature. We then show that the associated asymptotic variance expression for additive functionals grows at most linearly in time under hypotheses that are weaker than those currently existing in the literature. The assumptions are verified for some hidden Markov models.

A compound Poisson process whose randomized time is an independent Poisson process is called a compound Poisson process with Poisson subordinator. We provide its probability distribution, which is expressed in terms of the Bell polynomials, and investigate in detail both the special cases in which the compound Poisson process has exponential jumps and normal jumps. Then for the iterated Poisson process we discuss some properties and provide convergence results to a Poisson process. The first-crossing time problem for the iterated Poisson process is finally tackled in the cases of (i) a decreasing and constant boundary, where we provide some closed-form results, and (ii) a linearly increasing boundary, where we propose an iterative procedure to compute the first-crossing time density and survival functions.

We deal with a random graph model evolving in discrete time steps by duplicating and deleting the edges of randomly chosen vertices. We prove the existence of an almost surely asymptotic degree distribution, with stretched exponential decay; more precisely, the proportion of vertices of degree d tends to some positive number cd > 0 almost surely as the number of steps goes to ∞, and cd ~ (eπ)1/2d1/4e-2√d holds as d → ∞.

We provide explicit expressions for the constants involved in the characterisation of ergodicity of subgeometric Markov chains. The constants are determined in terms of those appearing in the assumed drift and one-step minorisation conditions. The results are fundamental for the study of some algorithms where uniform bounds for these constants are needed for a family of Markov kernels. Our results accommodate also some classes of inhomogeneous chains.

In this paper we study a generalized coupon collector problem, which consists of determining the distribution and the moments of the time needed to collect a given number of distinct coupons that are drawn from a set of coupons with an arbitrary probability distribution. We suppose that a special coupon called the null coupon can be drawn but never belongs to any collection. In this context, we obtain expressions for the distribution and the moments of this time. We also prove that the almost-uniform distribution, for which all the nonnull coupons have the same drawing probability, is the distribution which minimizes the expected time to obtain a fixed subset of distinct coupons. This optimization result is extended to the complementary distribution of the time needed to obtain the full collection, proving by the way this well-known conjecture. Finally, we propose a new conjecture which expresses the fact that the almost-uniform distribution should minimize the complementary distribution of the time needed to obtain any fixed number of distinct coupons.

This paper concerns discrete-time Markov decision chains with denumerable state and compact action sets. Besides standard continuity requirements, the main assumption on the model is that it admits a Lyapunov function ℓ. In this context the average reward criterion is analyzed from the sample-path point of view. The main conclusion is that if the expected average reward associated to ℓ2 is finite under any policy then a stationary policy obtained from the optimality equation in the standard way is sample-path average optimal in a strong sense.

This paper deals with the first passage optimality and variance minimisation problems of discrete-time Markov decision processes (MDPs) with varying discount factors and unbounded rewards/costs. First, under suitable conditions slightly weaker than those in the previous literature on the standard (infinite horizon) discounted MDPs, we establish the existence and characterisation of the first passage expected-optimal stationary policies. Second, to further distinguish the expected-optimal stationary policies, we introduce the variance minimisation problem, prove that it is equivalent to a new first passage optimality problem of MDPs, and, thus, show the existence of a variance-optimal policy that minimises the variance over the set of all first passage expected-optimal stationary policies. Finally, we use a computable example to illustrate our main results and also to show the difference between the first passage optimality here and the standard discount optimality of MDPs in the previous literature.

In this paper we consider a class of stochastic processes based on binomial observations of continuous-time, Markovian population models. We derive the conditional probability mass function of the next binomial observation given a set of binomial observations. For this purpose, we first find the conditional probability mass function of the underlying continuous-time Markovian population model, given a set of binomial observations, by exploiting a conditional Bayes' theorem from filtering, and then use the law of total probability to find the former. This result paves the way for further study of the stochastic process introduced by the binomial observations. We utilize our results to show that binomial observations of the simple birth process are non-Markovian.

We consider a production-inventory model operating in a stochastic environment that is modulated by a finite state continuous-time Markov chain. When the inventory level reaches zero, an order is placed from an external supplier. The costs (purchasing and holding costs) are modulated by the state at the order epoch time. Applying a matrix analytic approach, fluid flow techniques, and martingales, we develop methods to obtain explicit equations for these cost functionals in the discounted case and under the long-run average criterion. Finally, we extend the model to allow backlogging.

The concept of a signature is a useful tool in the analysis of semicoherent systems with continuous, and independent and identically distributed component lifetimes, especially for the comparison of different system designs and the computation of the system reliability. For such systems, we provide conversion formulae between the signature and the reliability function through the corresponding vector of dominations and we derive efficient algorithms for the computation of any of these concepts from any other. We also show how the signature can be easily computed from the reliability function via basic manipulations such as differentiation, coefficient extraction, and integration.

We show how to determine if a given vector can be the signature of a system on a finite number of components and, if so, exhibit such a system in terms of its structure function. The method employs combinatorial results from the theory of (finite) simplicial complexes, and provides a full characterization of signature vectors using a theorem of Kruskal (1963) and Katona (1968). We also show how the same approach can provide new combinatorial proofs of further results, e.g. that the signature vector of a system cannot have isolated zeroes. Finally, we prove that a signature with all nonzero entries must be a uniform distribution.

Full likelihood inference under Kingman's coalescent is a computationally challenging problem to which importance sampling (IS) and the product of approximate conditionals (PAC) methods have been applied successfully. Both methods can be expressed in terms of families of intractable conditional sampling distributions (CSDs), and rely on principled approximations for accurate inference. Recently, more general Λ- and Ξ-coalescents have been observed to provide better modelling fits to some genetic data sets. We derive families of approximate CSDs for finite sites Λ- and Ξ-coalescents, and use them to obtain ‘approximately optimal’ IS and PAC algorithms for Λ-coalescents, yielding substantial gains in efficiency over existing methods.

We consider a dynamic metapopulation involving one large population of size N surrounded by colonies of size εNN, usually called peripheral isolates in ecology, where N → ∞ and εN → 0 in such a way that εNN → ∞. The main population, as well as the colonies, independently send propagules to found new colonies (emigration), and each colony independently, eventually merges with the main population (fusion). Our aim is to study the genealogical history of a finite number of lineages sampled at stationarity in such a metapopulation. We make assumptions on model parameters ensuring that the total outer population has size of the order of N and that each colony has a lifetime of the same order. We prove that under these assumptions, the scaling limit of the genealogical process of a finite sample is a censored coalescent where each lineage can be in one of two states: an inner lineage (belonging to the main population) or an outer lineage (belonging to some peripheral isolate). Lineages change state at constant rate and (only) inner lineages coalesce at constant rate per pair. This two-state censored coalescent is also shown to converge weakly, as the landscape dynamics accelerate, to a time-changed Kingman coalescent.

In this paper a stochastic failure model for a system with stochastically dependent competing failures is analyzed. The system is subject to two types of failure: degradation failure and catastrophic failure. Both types of failure share an initial common source: an external shock process. This implies that they are stochastically dependent. In our developments of the model, the type of dependency between the two kinds of failure will be characterized. Conditional properties of the two competing risks are also investigated. These properties are the fundamental basis for the development of the maintenance strategy studied in this paper. Considering this maintenance strategy, the long-run average cost rate is derived and the optimal maintenance policy is discussed.

In this paper we consider general counting processes stopped at a random time T, independent of the process. Provided that T has the decreasing failure rate (DFR) property, we present sufficient conditions on the arrival times so that the number of events occurring before T preserves the DFR property of T. In particular, when the interarrival times are independent, we consider applications concerning the DFR property of the stationary number of customers waiting in queue for specific queueing models.

Under proportional transaction costs, a price process is said to have a consistent price system, if there is a semimartingale with an equivalent martingale measure that evolves within the bid-ask spread. We show that a continuous, multi-asset price process has a consistent price system, under arbitrarily small proportional transaction costs, if it satisfies a natural multi-dimensional generalization of the stickiness condition introduced by Guasoni (2006).

In this paper we use the Siegert formula to derive alternative expressions for the moments of the first passage time of the Ornstein-Uhlenbeck process through a constant threshold. The expression for the nth moment is recursively linked to the lower-order moments and consists of only n terms. These compact expressions can substantially facilitate (numerical) applications also for higher-order moments.

Consider the problem of drawing random variates (X1, …, Xn) from a distribution where the marginal of each Xi is specified, as well as the correlation between every pair Xi and Xj. For given marginals, the Fréchet-Hoeffding bounds put a lower and upper bound on the correlation between Xi and Xj. Any achievable correlation between Xi and Xj is a convex combination of these bounds. We call the value λ(Xi, Xj) ∈ [0, 1] of this convex combination the convexity parameter of (Xi, Xj) with λ(Xi, Xj) = 1 corresponding to the upper bound and maximal correlation. For given marginal distributions functions F1, …, Fn of (X1, …, Xn), we show that λ(Xi, Xj) = λij if and only if there exist symmetric Bernoulli random variables (B1, …, Bn) (that is {0, 1} random variables with mean ½) such that λ(Bi, Bj) = λij. In addition, we characterize completely the set of convexity parameters for symmetric Bernoulli marginals in two, three, and four dimensions.