Search

In this paper, new conditions for the stability of V-geometrically ergodic Markov chains are introduced. The results are based on an extension of the standard perturbation theory formulated by Keller and Liverani. The continuity and higher regularity properties are investigated. As an illustration, an asymptotic expansion of the invariant probability measure for an autoregressive model with independent and identically distributed noises (with a nonstandard probability density function) is obtained.

In this paper we quantify the efficiency of parallelism in systems that are prone to failures and exhibit power law processing delays. We characterize the performance of two prototype schemes of parallelism, redundant and split, in terms of both the power law exponent and exact asymptotics of the delay distribution tail. We also develop the optimal splitting scheme which ensures that split always outperforms redundant.

Sharp tail bounds for the sum of d random variables with given marginal distributions and arbitrary dependence structure have been known since Makarov (1981) and Rüschendorf (1982) for d=2 and, in some examples, for d≥3. Based on a duality result, dual bounds have been introduced in Embrechts and Puccetti (2006b). In the homogeneous case, F 1=···=F n , with monotone density, sharp tail bounds were recently found in Wang and Wang (2011). In this paper we establish the sharpness of the dual bounds in the homogeneous case under general conditions which include, in particular, the case of monotone densities and concave densities. We derive the corresponding optimal couplings and also give an effective method to calculate the sharp bounds.

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

In this paper we introduce a variation on Glaser's method for determining the shape of the failure rate function of a mixture. It has often been seen that the shape of the failure rate depends on the mixing parameter q. Our method provides an explanation for this phenomenon. We then illustrate our technique with the mixture of an exponential and a gamma density for all possible cases.

In this paper we explore the first passage times of constant-elasticity-of-variance (CEV) processes with two-sided reflecting barriers. The explicit Laplace transforms of the first passage times are derived. Our results can include analytic formulae concerning Laplace transforms of first passage times of reflected Ornstein–Uhlenbeck processes, reflected geometric Brownian motions, and reflected square-root processes.

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 X i with E(X i )=µ, and where there is an edge between i and j with probability 1-e-X i X j /(µ 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.

In this paper we study the functional central limit theorem (CLT) for stationary Markov chains with a self-adjoint operator and general state space. We investigate the case when the variance of the partial sum is not asymptotically linear in n, and establish that conditional convergence in distribution of partial sums implies the functional CLT. The main tools are maximal inequalities that are further exploited to derive conditions for tightness and convergence to the Brownian motion.

In this paper we deal with an M/G/1 vacation system with the sojourn time (wait plus service) limit and two typical vacation rules, i.e. multiple and single vacation rules. Using the level crossing approach, we derive recursive equations for the steady-state distributions of the virtual waiting times in M/G/1 vacation systems with a general vacation time and two vacation rules.

The risk processes considered in this paper are generated by an underlying Markov process with a regenerative structure and an independent sequence of independent and identically distributed claims. Between the arrivals of claims the process increases at a rate which is a nonnegative function of the present value of the Markov process. The intensity for a claim to occur is another nonnegative function of the value of the Markov process. The claim arrival times are the regeneration times for the Markov process. Two-sided claims are allowed, but the distribution of the positive claims is assumed to have a Laplace transform that is a rational function. The main results describe the joint Laplace transform of the time at ruin and the deficit at ruin. The method used consists in finding partial eigenfunctions for the generator of the joint process consisting of the Markov process and the accumulated claims process, a joint process which is also Markov. These partial eigenfunctions are then used to find a martingale that directly leads to an expression for the desired Laplace transform. In the final section, three examples are given involving different types of the underlying Markov process.

Nearly fifty years after the introduction of skew Brownian motion by Itô and McKean (1963), the first passage time distribution remains unknown. In this paper we first generalize results of Pitman and Yor (2011) and Csáki and Hu (2004) to derive formulae for the distribution of ranked excursion heights of skew Brownian motion, and then use these results to derive the first passage time distribution.

Let (X, J) denote a Markov-modulated Brownian motion (MMBM) and denote its supremum process by S. For some a > 0, let σ(a) denote the time when the reflected process Y := S - X first surpasses the level a. Furthermore, let σ−(a) denote the last time before σ(a) when X attains its current supremum. In this paper we shall derive the joint distribution of S σ(a), σ−(a), and σ(a), where the latter two will be given in terms of their Laplace transforms. We also provide some remarks on scale matrices for MMBMs with strictly positive variation parameters. This extends recent results for spectrally negative Lévy processes to MMBMs. Due to well-known fluid embedding and state-dependent killing techniques, the analysis applies to Markov additive processes with phase-type jumps as well. The result is of interest to applications such as the dividend problem in insurance mathematics and the buffer overflow problem in queueing theory. Examples will be given for the former.

The notion of the signature is a basic concept and a powerful tool in the analysis of networks and reliability systems of binary type. An appropriate definition of this concept has recently been introduced for systems that have ν possible states (with ν ≥ 3). In this paper we analyze in detail several properties and the most relevant aspects of such a general definition. For simplicity's sake, we focus our attention on the case ν = 3. Our analysis will however provide a number of hints for understanding the basic aspects of the general case.

A discrete-time SIS model is presented that allows individuals in the population to vary in terms of their susceptibility to infection and their rate of recovery. This model is a generalisation of the metapopulation model presented in McVinish and Pollett (2010). The main result of the paper is a central limit theorem showing that fluctuations in the proportion of infected individuals around the limiting proportion converges to a Gaussian random variable when appropriately rescaled. In contrast to the case where there is no variation amongst individuals, the limiting Gaussian distribution has a nonzero mean.

This paper is devoted to the perfect simulation of a stationary process with an at most countable state space. The process is specified through a kernel, prescribing the probability of the next state conditional to the whole past history. We follow the seminal work of Comets, Fernández and Ferrari (2002), who gave sufficient conditions for the construction of a perfect simulation algorithm. We define backward coalescence times for these kind of processes, which allow us to construct perfect simulation algorithms under weaker conditions than in Comets, Fernández and Ferrari (2002). We discuss how to construct backward coalescence times (i) by means of information depths, taking into account some a priori knowledge about the histories that occur; and (ii) by identifying suitable coalescing events.

In this paper we determine the distributions of occupation times of a Markov-modulated Brownian motion (MMBM) in separate intervals before a first passage time or an exit from an interval. We derive the distributions in terms of their Laplace transforms, and we also distinguish between occupation times in different phases. For MMBMs with strictly positive variation parameters, we further propose scale functions.

International audience; Modeling insurance risks is a task that received an increasing attention because of Solvency Capital Requirements. The ruin probability has become a standard risk measure to assess regulatory capital. In this paper we focus on discrete time models for nite time horizon. Several results are available in the literature allowing to calibrate the ruin probability by means of the sum of the tail probabilities of individual claim amounts. The aim of this work is to obtain asymptotics for such probabilities under multivariate regularly variation and, more precisely, to derive them from Breiman's Theorem extensions. We thus exhibit new situations where the ruin probability admits computable equivalents. Consequences are also derived in terms of the Value-at-Risk.

The theory of L 2-spectral gaps for reversible Markov chains has been studied by many authors. In this paper we consider positive recurrent general state space Markov chains with stationary transition probabilities. Replacing the assumption of reversibility with a weaker assumption, we still obtain a simple necessary and sufficient condition for the spectral gap property of the associated Markov operator in terms of the isoperimetric constant. We show that this result can be applied to a large class of Markov chains, including those that are related to positive recurrent finite-range random walks on Z.

In this paper we analyze the so-called Parisian ruin probability, which arises when the surplus process stays below 0 longer than a fixed amount of time ζ > 0. We focus on a general spectrally negative Lévy insurance risk process. For this class of processes, we derive an expression for the ruin probability in terms of quantities that can be calculated explicitly in many models. We find its Cramér-type and convolution-equivalent asymptotics when reserves tend to ∞. Finally, we analyze some explicit examples.

Sequential decision problems can often be modeled as Markov decision processes. Classical solution approaches assume that the parameters of the model are known. However, model parameters are usually estimated and uncertain in practice. As a result, managers are often interested in how estimation errors affect the optimal solution. In this paper we illustrate how sensitivity analysis can be performed directly for a Markov decision process with uncertain reward parameters using the Bellman equations. In particular, we consider problems involving (i) a single stationary parameter, (ii) multiple stationary parameters, and (iii) multiple nonstationary parameters. We illustrate the applicability of this work through a capacitated stochastic lot-sizing problem.

We investigate the Zagreb index, one of the topological indices, of random recursive trees in this paper. Through a recurrence equation, the first two moments of Zn, the Zagreb index of a random recursive tree of size n, are obtained. We also show that the random process {Zn − E[Zn], n ≥ 1} is a martingale. Then the asymptotic normality of the Zagreb index of a random recursive tree is given by an application of the martingale central limit theorem. Finally, two other topological indices are also discussed in passing.

In this paper we provide some sufficient conditions to stochastically compare linear combinations of independent random variables. The main results extend those given in Proschan (1965), Ma (1998), Zhao et al. (2011), and Yu (2011). In particular, we propose a new sufficient condition to compare the peakedness of linear combinations of independent random variables which may have heavy-tailed properties.

In this paper we consider an extension of the two-dimensional risk model introduced in Avram, Palmowski and Pistorius (2008a). To this end, we assume that there are two insurers. The first insurer is subject to claims arising from two independent compound Poisson processes. The second insurer, which can be viewed as a different line of business of the same insurer or as a reinsurer, covers a proportion of the claims arising from one of these two compound Poisson processes. We derive the Laplace transform of the time until ruin of at least one insurer when the claim sizes follow a general distribution. The surplus level of the first insurer when the second insurer is ruined first is discussed at the end in connection with some open problems.

We construct an integrated probabilistic model to capture interactions between atoms of a nanocomponent. We then use this model to assess reliabilities of nanocomponents with different structures. Several properties of our proposed model are also described under a sparseness condition. The model is an extension of our previous model based on Markovian random field theory. The proposed integrated model is flexible in that pairwise relationship information among atoms as well as features of individual atoms can be easily incorporated. An important feature that distinguishes the integrated probabilistic model from our previous model is that the integrated approach uses all available sources of information with different weights for different types of interaction. In this paper we consider the nanocomponent at a fixed moment of time, say the present moment, and we assume that the present state of the nanocomponent depends only on the present states of its atoms.

In this paper we consider a general Lévy process X reflected at a downward periodic barrier At and a constant upper barrier K, giving a process VKt=Xt +LAt−LKt. We find the expression for a loss rate defined by lK=ELK1 and identify its asymptotics as K→∞ when X has light-tailed jumps and EX1

Recently, there has been considerable interest in the calculation of decay rates for models that can be viewed as quasi-birth-and-death (QBD) processes with infinitely many phases. In this paper we make a contribution to this endeavour by considering some classes of models in which the transition function is not homogeneous in the phase direction. We characterize the range of decay rates that are compatible with the dynamics of the process away from the boundary. In many cases, these rates can be attained by changing the transition structure of the QBD process at level 0. Our approach, which relies on the use of orthogonal polynomials, is an extension of that in Motyer and Taylor (2006) for the case where the generator has homogeneous blocks. © 2011 Applied Probability Trust.

Nonspatial stochastic models of populations subject to catastrophic events result in the common conclusion that the survival probability of the population is nondecreasing with respect to the random number of individuals removed at each catastrophe. The purpose of this paper is to prove that such a monotonic relationship is not true for simple spatial models based on Harris' contact processes, whose dynamics are described by hypergraph structures rather than traditional graph structures. More precisely, it is proved that, for a wide range of parameters, the destruction of (infinite) hyperplanes does not affect the existence of a nontrivial invariant measure, whereas the destruction of large (finite) cubes drives the population to extinction, a result that we depict by using the biological picture: forest fires are more devastating than tornadoes. This indicates that the geometry of the subsets struck by catastrophes is somewhat more important than their area, thus the need to consider spatial rather than nonspatial models in this context.

To enhance the performance of a system, a common practice employed by reliability engineers is to use redundant components in the system. In this paper we compare lifetimes of series (parallel) systems arising out of different allocations of one or two standby redundancies. These comparisons are made with respect to the increasing concave (convex) order, the hazard rate order, and the stochastic precedence order. The main results extend some related conclusions in the literature.

A nanocomponent is a collection of atoms arranged to a specific design in order to achieve a desired function with an acceptable performance and reliability. The type of atoms, the manner in which they are arranged within the nanocomponent, and their interrelationship have a direct effect on the nanocomponent's reliability (survival) function. In this paper we propose models based on the notion of a copula that are used to describe the relationship between the atoms of a nanocomponent. Having defined these models, we go on to construct a ‘nanocomponent’ model in order to obtain the reliability function of a nanocomponent.

Quasi-stationary distributions, as discussed in Darroch and Seneta (1965), have been used in biology to describe the steady state behaviour of population models which, while eventually certain to become extinct, nevertheless maintain an apparent stochastic equilibrium for long periods. These distributions have some drawbacks: they need not exist, nor be unique, and their calculation can present problems. In this paper, we give biologically plausible conditions under which the quasi-stationary distribution is unique, and can be closely approximated by distributions that are simple to compute.

Let {X i } i=1 n be a sequence of random variables with two possible outcomes, denoted 0 and 1. Define a random variable S n,m to be the maximum number of 1s within any m consecutive trials in {X i } i=1 n . The random variable S n,m is called a discrete scan statistic and has applications in many areas. In this paper we evaluate the distribution of discrete scan statistics when {X i } i=1 n consists of exchangeable binary trials. We provide simple closed-form expressions for both conditional and unconditional distributions of S n,m for 2m ≥ n. These results are also new for independent, identically distributed Bernoulli trials, which are a special case of exchangeable trials.

In this paper we prove that the stationary distribution of populations in genetic algorithms focuses on the uniform population with the highest fitness value as the selective pressure goes to ∞ and the mutation probability goes to 0. The obtained sufficient condition is based on the work of Albuquerque and Mazza (2000), who, following Cerf (1998), applied the large deviation principle approach (Freidlin-Wentzell theory) to the Markov chain of genetic algorithms. The sufficient condition is more general than that of Albuquerque and Mazza, and covers a set of parameters which were not found by Cerf.

It can be seen that a mixture of an exponential distribution and a gamma distribution with increasing failure rate for the right choice of parameters can yield a distribution with a bathtub-shaped failure rate. In this paper we consider a continuous mixture of exponentials and a continuous mixture of gammas with increasing failure rates and show that the resulting mixture has a bathtub-shaped failure rate.

In this paper we study the Wiener-Hopf factorization for a class of L\'evy processes with double-sided jumps, characterized by the fact that the density of the L\'evy measure is given by an infinite series of exponential functions with positive coefficients. We express the Wiener-Hopf factors as infinite products over roots of a certain transcendental equation, and provide a series representation for the distribution of the supremum/infimum process evaluated at an independent exponential time. We also introduce five eight-parameter families of L\'evy processes, defined by the fact that the density of the L\'evy measure is a (fractional) derivative of the theta-function, and we show that these processes can have a wide range of behavior of small jumps. These families of processes are of particular interest for applications, since the characteristic exponent has a simple expression, which allows efficient numerical computation of the Wiener-Hopf factors and distributions of various functionals of the process., Comment: 12 pages, published online at http://projecteuclid.org/euclid.jap/1294170516

The purpose of this paper is to study an optimal stopping problem with constraints for a Markov chain with general state space by using the convex analytic approach. The costs are assumed to be nonnegative. Our model is not assumed to be transient or absorbing and the stopping time does not necessarily have a finite expectation. As a consequence, the occupation measure is not necessarily finite, which poses some difficulties in the analysis of the associated linear program. Under a very weak hypothesis, it is shown that the linear problem admits an optimal solution, guaranteeing the existence of an optimal stopping strategy for the optimal stopping problem with constraints.

In this paper we study the moment generating function order and the new better than used in the moment generating function order (NBUMG) life distributions. A closure property of this order under an independent random sum is deduced, and stochastic comparisons among the block replacement policy, the age replacement policy, the complete repair policy, and the minimal repair policy of an NBUMG component are investigated.

A sequence of random variables is said to be extended negatively dependent (END) if the tails of its finite-dimensional distributions in the lower-left and upper-right corners are dominated by a multiple of the tails of the corresponding finite-dimensional distributions of a sequence of independent random variables with the same marginal distributions. The goal of this paper is to establish the strong law of large numbers for a sequence of END and identically distributed random variables. In doing so we derive some new inequalities of large deviation type for the sums of END and identically distributed random variables being suitably truncated. We also show applications of our main result to risk theory and renewal theory.

In this paper we study the asymptotic normality of discrete-time Markov control processes in Borel spaces, with possibly unbounded cost. Under suitable hypotheses, we show that the cost sequence is asymptotically normal. As a special case, we obtain a central limit theorem for (noncontrolled) Markov chains.

This paper is concerned with properties of the algebraic degree of the Laplace-Stieltjes transform of phase-type (PH) distributions. The main problem of interest is: given a PH generator, how do we find the maximum and the minimum algebraic degrees of all irreducible PH representations with that PH generator? Based on the matrix exponential (ME) order of ME distributions and the spectral polynomial algorithm, a method for computing the algebraic degree of a PH distribution is developed. The maximum algebraic degree is identified explicitly. Using Perron-Frobenius theory of nonnegative matrices, a lower bound and an upper bound on the minimum algebraic degree are found, subject to some conditions. Explicit results are obtained for special cases.

In this paper we obtain several mixture representations of the reliability function of the inactivity time of a coherent system under the condition that the system has failed at time t (> 0) in terms of the reliability functions of inactivity times of order statistics. Some ordering properties of the inactivity times of coherent systems with independent and identically distributed components are obtained, based on the stochastically ordered coefficient vectors between systems.

The topic of the present paper is a generalized St Petersburg game in which the distribution of the payoff X is given by P(X = sr (k-1)/α) = pq k-1, k = 1, 2,…, where p + q = 1, s = 1 / p, r = 1 / q, and 0 < α ≤ 1. For the case in which α = 1, we extend Feller's classical weak law and Martin-Löf's theorem on convergence in distribution along the 2 n -subsequence. The analog for 0 < α < 1 turns out to converge in distribution to an asymmetric stable law with index α. Finally, some limit theorems for polynomial and geometric size total gains, as well as for extremes, are given.

In the present paper an almost-sure renewal theorem for branching random walks (BRWs) on the real line is formulated and established. The theorem constitutes a generalization of Nerman's theorem on the almost-sure convergence of Malthus normed supercritical Crump-Mode-Jagers branching processes counted with general characteristic and Gatouras' almost-sure renewal theorem for BRWs on a lattice.

In this paper we consider the problem of identifiability for the two-state Markovian arrival process (MAP2). In particular, we show that the MAP2 is not identifiable, providing the conditions under which two different sets of parameters induce identical stationary laws for the observable process.

Li and Shaked (2007) introduced the family of generalized total time on test transform (TTT) stochastic orders, which is parameterized by a real function h that can be used to capture the preferences of a decision maker. It is natural to look for properties of these orders when there is an uncertainty in determining the appropriate function h. In this paper we study these orders when h is nondecreasing. We note that all these orders are location independent, and we characterize the dispersive order, and the location-independent riskier order, by means of the generalized TTT orders with nondecreasing h. Further properties, which strengthen known properties of the dispersive order, are given. A useful nontrivial closure property of the generalized TTT orders with nondecreasing h is obtained. Applications in poverty comparisons, risk management, and reliability theory are described.

In this paper we consider the probabilities of finite- and infinite-time absolute ruins in the renewal risk model with constant premium rate and constant force of interest. In the particular case of the compound Poisson model, explicit asymptotic expressions for the finite- and infinite-time absolute ruin probabilities are given. For the general renewal risk model, we present an asymptotic expression for the infinite-time absolute ruin probability. Conditional distributions of Poisson processes and probabilistic techniques regarding randomly weighted sums are employed in the course of this study.

In this paper we study the stability of queueing systems with impatient customers and a single server operating under a FIFO (first-in-first-out) discipline. We first give a sufficient condition for the existence of a stationary workload in the case of impatience until the beginning of service. We then provide a weaker condition of existence on an enriched probability space using the theory of Anantharam et al. (1997), (1999). The case of impatience until the end of service is also investigated.

In this paper we propose a general bivariate random effect model with special emphasis on frailty models and environmental effect models, and present some stochastic comparisons. The relationship between the conditional and the unconditional hazard gradients are derived and some examples are provided. We investigate how the well-known stochastic orderings between the distributions of two frailties translate into the orderings between the corresponding survival functions. These results are used to obtain the properties of the bivariate multiplicative model and the shared frailty model.

The aim of this paper is to study the stochastic monotonicity and continuity properties of the extinction time of Bellman-Harris branching processes depending on their reproduction laws. Moreover, we show their applications in an epidemiological context, obtaining an optimal criterion to establish the proportion of susceptible individuals in a given population that must be vaccinated in order to eliminate an infectious disease. First the spread of infection is modelled by a Bellman-Harris branching process. Finally, we provide a simulation-based method to determine the optimal vaccination policies.

System signatures are useful tools in the study and comparison of coherent systems. In this paper, we define and study a similar concept, called the joint signature, for two coherent systems which share some components. Under an independent and identically distributed assumption on component lifetimes, a pseudo-mixture representation based on this joint signature is obtained for the joint distribution of the lifetimes of both systems. Sufficient conditions are given based on the respective joint signatures of two pairs of systems, each with shared components, to ensure various forms of bivariate stochastic orderings between the joint lifetimes of the two pairs of systems. The present investigation is focused on the joint behavior of pairs of systems which have one or more shared components. Our main goal is to characterize the joint (bivariate) distribution of the lifetimes T\ and T2 of the two systems. Our approach involves the expansion of the scope of the theory of system signatures. We will develop representation results for the joint distribution G(t\, t2) = P(T\ < t\, T2 < t2) in terms of the newly defined joint signature of the two systems. Before reviewing the characteristics of system signatures in the univariate case that will be relevant to the present study, and before we proceed with the development of joint signatures and their applications, we will begin with a brief discussion of the type of systems and scenarios from which the present work is motivated. A frequently encountered example of systems with shared components occurs in networked computing in which a server (say a file server or Web server) is used in tandem with several individual computers. It is typical that departments within a company or university will store almost all of the files for the department's individual PCs on one central server. If the central server goes down, the PCs with local disks may retain certain limited capabilities, while other PCs may not work at all. The performance of any given pair of PCs will depend on the

In a two-person red-and-black game, each player holds an integral amount of chips. At each stage of the game, each player can bet any integral amount in his possession, winning the chips of his opponent with a probability which is a function of the ratio of his bet to the sum of both players' bets and is called a win probability function. Both players seek to maximize the probability of winning the entire fortune of his opponent. In this paper we propose two new models. In the first model, at each stage, there is a positive probability that two players exchange their bets. In the second model, the win probability functions are stage dependent. In both models, we obtain suitable conditions on the win probability functions such that it is a Nash equilibrium for the subfair player to play boldly and for the superfair player to play timidly.