In this paper we study the properties of total time on test transforms of order n and examine their applications in reliability analysis. It is shown that the successive transforms produce either distributions with increasing or bathtub-shaped failure rates or distributions with decreasing or upside bathtub-shaped failure rates. The ageing properties of the baseline distribution is compared with those of transformed distributions, and a partial order based on nth-order transforms and their implications are discussed.

The present paper generalises some results for spectrally negative Lévy processes to the setting of Markov additive processes (MAPs). A prominent role is assumed by the first passage times, which will be determined in terms of their Laplace transforms. These have the form of a phase-type distribution, with a rate matrix that can be regarded as an inverse function of the cumulant matrix. A numerically stable iteration to compute this matrix is given. The theory is first developed for MAPs without positive jumps and then extended to include positive jumps having phase-type distributions. Numerical and analytical examples show agreement with existing results in special cases.

Weak local linear approximations have played a prominent role in the construction of effective inference methods and numerical integrators for stochastic differential equations. In this note two weak local linear approximations for stochastic differential equations with jumps are introduced as a generalization of previous ones. Their respective order of convergence is obtained as well.

We consider an infectious disease spreading along the edges of a network which may have significant clustering. The individuals in the population have heterogeneous infectiousness and/or susceptibility. We define the out-transmissibility of a node to be the marginal probability that it would infect a randomly chosen neighbor given its infectiousness and the distribution of susceptibility. For a given distribution of out-transmissibility, we find the conditions which give the upper (or lower) bounds on the size and probability of an epidemic, under weak assumptions on the transmission properties, but very general assumptions on the network. We find similar bounds for a given distribution of in-transmissibility (the marginal probability of being infected by a neighbor). We also find conditions giving global upper bounds on the size and probability. The distributions leading to these bounds are network independent. In the special case of networks with high girth (locally tree-like), we are able to prove stronger results. In general, the probability and size of epidemics are maximal when the population is homogeneous and minimal when the variance of in- or out-transmissibility is maximal.

We consider a stochastic SIS model for the spread of an epidemic amongst a population of n individuals that are equally spaced upon the circumference of a circle. Whilst infectious, an individual, i say, makes both local and global infectious contacts at the points of homogeneous Poisson point processes. Global contacts are made uniformly at random with members of the entire population, whilst local contacts are made according to a contact distribution centred upon the infective. Individuals at the end of their infectious period return to the susceptible state and can be reinfected. The emphasis of the paper is on asymptotic results as the population size n → ∞. Therefore, a contact process with global infection is introduced representing the limiting behaviour as n → ∞ of the circle epidemics. A branching process approximation for the early stages of the epidemic is derived and the endemic equilibrium of a major outbreak is obtained. Furthermore, assuming exponential infectious periods, the probability of a major epidemic outbreak and the proportion of the population infectious in the endemic equilibrium are shown to satisfy the same equation which characterises the epidemic process.

We establish a functional large deviation principle and a functional moderate deviation principle for Markov-modulated risk models with reinsurance by constructing an exponential martingale approach. Lundberg's estimate of the ruin time is also presented.

In this paper we study the α-stable continuous-state branching processes (for α ∈ (1, 2]) and the α-stable continuous-state branching processes conditioned never to become extinct in the light of positive self-similarity. Understanding the interaction of the Lamperti transformation for continuous-state branching processes and the Lamperti transformation for positive, self-similar Markov processes gives access to a number of explicit results concerning the paths of α-stable continuous-state branching processes and α-stable continuous-state branching processes conditioned never to become extinct.

In this paper we consider discrete-time multidimensional Markov chains having a block transition probability matrix which is the sum of a matrix with repeating block rows and a matrix of upper-Hessenberg, quasi-Toeplitz structure. We derive sufficient conditions for the existence of the stationary distribution, and outline two algorithms for calculating the stationary distribution.

In the framework of patterns in random texts, the Markov chain embedding techniques consist of turning the occurrences of a pattern over an order-m Markov sequence into those of a subset of states into an order-1 Markov chain. In this paper we use the theory of language and automata to provide space-optimal Markov chain embedding using the new notion of pattern Markov chains (PMCs), and we give explicit constructive algorithms to build the PMC associated to any given pattern problem. The interest of PMCs is then illustrated through the exact computation of P-values whose complexity is discussed and compared to other classical asymptotic approximations. Finally, we consider two illustrative examples of highly degenerated pattern problems (structured motifs and PROSITE signatures), which further illustrate the usefulness of our approach.

It is shown that, when expressing arguments in terms of their logarithms, the Laplace transform of a function is related to the antiderivative of this function by a simple convolution. This allows efficient numerical computations of moment generating functions of positive random variables and their inversion. The application of the method is straightforward, apart from the necessity to implement it using high-precision arithmetics. In numerical examples the approach is demonstrated to be particularly useful for distributions with heavy tails, such as lognormal, Weibull, or Pareto distributions, which are otherwise difficult to handle. The computational efficiency compared to other methods is demonstrated for an M/G/1 queueing problem.

In this paper we consider a compound Poisson risk model where the insurer earns credit interest at a constant rate if the surplus is positive and pays out debit interest at another constant rate if the surplus is negative. Absolute ruin occurs at the moment when the surplus first drops below a critical value (a negative constant). We study the asymptotic properties of the absolute ruin probability of this model. First we investigate the asymptotic behavior of the absolute ruin probability when the claim size distribution is light tailed. Then we study the case where the common distribution of claim sizes are heavy tailed.

The upper tail behaviour is explored for a stopped random product ∏j=1NXj, where the factors are positive and independent and identically distributed, and N is the first time one of the factors occupies a subset of the positive reals. This structure is motivated by a heavy-tailed analogue of the factorial n!, called the factoid of n. Properties of the factoid suggested by computer explorations are shown to be valid. Two topics about the determination of the Zipf exponent in the rank-size law for city sizes are discussed.

We specify a model for a catastrophe loss index, where the initial estimate of each catastrophe loss is reestimated immediately by a positive martingale starting from the random time of loss occurrence. We consider the pricing of catastrophe insurance options written on the loss index and obtain option pricing formulae by applying Fourier transform techniques. An important advantage is that our methodology works for loss distributions with heavy tails, which is the appropriate tail behavior for catastrophe modeling. We also discuss the case when the reestimation factors are given by positive affine martingales and provide a characterization of positive affine local martingales.

We study the convergence of at-the-money implied volatilities to the spot volatility in a general model with a Brownian component and a jump component of finite variation. This result is a consequence of the robustness of the Black-Scholes formula and of the central limit theorem for martingales.

We consider a sequential rule, where an item is chosen into the group, such as a university faculty member, only if his/her score is better than the average score of those already belonging to the group. We study four variables: the average score of the members of the group after k items have been selected, the time it takes (in terms of the number of observed items) to assemble a group of k items, the average score of the group after n items have been observed, and the number of items kept after the first n items have been observed. We develop the relationships between these variables, and obtain their asymptotic behavior as k (respectively, n) tends to ∞. The assumption throughout is that the items are independent and identically distributed with a continuous distribution. Though knowledge of this distribution is not needed to implement the selection rule, the asymptotic behavior does depend on the distribution. We study in some detail the exponential, Pareto, and beta distributions. Generalizations of the ‘better than average’ rule to the β better than average rules are also considered. These are rules where an item is admitted to the group only if its score is better than β times the present average of the group, where β > 0.

A sequence of independent Bernoulli random variables with success probabilities a / (a + b + k − 1), k = 1, 2, 3, …, is embedded in a marked Poisson process with intensity 1. Using this, conditional Poisson limits follow for counts of failure strings.

We investigate the problem of using a Riemannian sum with random subintervals to approximate the iterated Itô integral ∫wdw - or, equivalently, solving the corresponding stochastic differential equation by Euler's method with variable step sizes. In the past this task has been used as a counterexample to illustrate that variable step sizes must be used with extreme caution in stochastic numerical analysis. This article establishes a class of variable step size schemes which do work.

For exponential open and closed queueing networks, we investigate the internal dependence structure, compare the internal dependence for different networks, and discuss the relation of correlation formulae to the existence of spectral gaps and comparison of asymptotic variances. A central prerequisite for the derived theorems is stochastic monotonicity of the networks. The dependence structure of network processes is described by concordance order with respect to various classes of functions. Different networks with the same first-order characteristics are compared with respect to their second-order properties. If a network is perturbed by changing the routeing in a way which holds the routeing equilibrium fixed, the resulting perturbations of the network processes are evaluated.

We study the number of collisions, Xn, of an exchangeable coalescent with multiple collisions (Λ-coalescent) which starts with n particles and is driven by rates determined by a finite characteristic measure η(dx) = x−2Λ(dx). Via a coupling technique, we derive limiting laws of Xn, using previous results on regenerative compositions derived from stick-breaking partitions of the unit interval. The possible limiting laws of Xn include normal, stable with index 1 ≤ α < 2, and Mittag-Leffler distributions. The results apply, in particular, to the case when η is a beta(a − 2, b) distribution with parameters a > 2 and b > 0. The approach taken allows us to derive asymptotics of three other functionals of the coalescent: the absorption time, the length of an external branch chosen at random from the n external branches, and the number of collision events that occur before the randomly selected external branch coalesces with one of its neighbours.

The failure rate of a mixture of even the most standard distributions used in reliability can have a complicated shape. However, failure rates of mixtures of two carefully selected distributions will have the well-known bathtub shape. Here we show that mixtures of whole families of distribtions can have a bathtub-shaped failure rate.