It is known [8] that a certain class of bond-decorated graphs exhibits multiple AB percolation phase transitions. Sufficient conditions are given under which the corresponding AB percolation critical probabilities may be identified as points of intersection of the graph of a certain polynomial with the boundary of the percolative region of an associated two-parameter bond-site percolation model on the underlying undecorated graph. The main result of the article is used to prove that the graphs in [8] exhibit multiple AB percolation critical probabilities. The possibility of identifying AB percolation critical exponents with corresponding limits for the bond-site model is discussed.

A large deviation expansion is used to evaluate the impact of errors in a geographical database on the computation of travel times. We work in the framework of discrete random variables and improve a theorem by Book to solve this problem. Simulations are provided to illustrate the methodology.

In Gadag and Rajarshi (1987), we studied a bivariate (multitype) branching process based on infinite and finite lines of descent, of particles of a supercritical one-dimensional (multitype) Galton-Watson branching process (GWBP). In this paper, we discuss a few more meaningful and interesting univariate and multitype branching processes, based on exact progeny lengths of particles in a GWBP. Our constructions relax the assumption of supercriticality made in Gadag and Rajarshi (1987). We investigate some finite-time and asymptotic results of these processes in some details and relate them to the original process. These results are then used to propose new and better estimates of the offspring mean. An illustration based on the branching process of the white male population of the USA is also given. We believe that our work offers a rather finer understanding of the branching property.

Let {Z(t), t ≧ 0} be an age-dependent subcritical branching process. In this paper we show that if the lifetime distribution is subexponential, EZα (t) ~ EZ(t) (t →∞) for every α ≧ 1. If furthermore the lifetime distribution has a subexponential density, a rate of convergence result in the above relation is established.

Examples are given of point processes that are non-stationary but have stationary forward recurrence time distributions. They are obtained by modification of stationary Poisson and renewal processes.

Consider a spectrally negative risk process where, on ruin, the deficit is immediately paid, and the process restarts from 0. When the process reaches a threshold b, all the surplus above b is paid as dividend. Applying the theory of exit times for a spectrally negative Lévy process and its reflection at the maximum and at the minimum, we obtain recursive formulae for the following moments. (i) The moments of the discounted loss until the process reaches b. This is equivalent to the moments of the discounted dividends in the dual model under the barrier strategy. (ii) The moments of the discounted loss for models with and without a dividend barrier for the infinite horizon. (iii) The moments of the discounted dividends for the infinite horizon.

We have found a mistake in the proofs of Navarro (2008, Theorem 2.3(b) and 2.3(c)) due to misapplication of properties of hazard rate and likelihood ratio orders. In this paper we show with an example that the stated results do not hold. This example is interesting since it proves some unexpected properties for these orderings under the formation of coherent systems. The result stated in Navarro (2008, Theorem 2.3(a)) for the usual stochastic order is correct.

A transformation is introduced to stabilize the variance of the largest value in a sample from a normal distribution. The transformed distribution is found to approximate closely to the limiting Gumbel form for all sample sizes. The practical implications of the result are discussed.

Vasicek (1977) proved that among all queueing disciplines that do not change the departure process of the queue, FIFO and LIFO yield, respectively, the smallest and the largest expectation of any given convex function of the service delay. In this note we further show that, if arriving customers join the queue stochastically ‘closer' to the server(s), then the expected value of any convex function of service delay is larger. As a more interesting result, we also show that if the function under consideration is concave, then the conclusion will be exactly the opposite. This result indicates that LIFO will be the best discipline if the delay cost is an increasing function but at a diminishing rate.

Comparison results for Markov processes with respect to function-class-induced (integral) stochastic orders have a long history. The most general results so far for this problem have been obtained based on the theory of evolution systems on Banach spaces. In this paper we transfer the martingale comparison method, known for the comparison of semimartingales to Markovian semimartingales, to general Markov processes. The basic step of this martingale approach is the derivation of the supermartingale property of the linking process, giving a link between the processes to be compared. This property is achieved using the characterization of Markov processes by the associated martingale problem in an essential way. As a result, the martingale comparison method gives a comparison result for Markov processes under a general alternative but related set of regularity conditions compared to the evolution system approach.

Let X1, X2,… and Y1, Y2,… be two sequences of absolutely continuous, independent and identically distributed (i.i.d.) random variables with equal means E(Xi)=E(Yi), i=1,2,… In this work we provide upper bounds for the total variation and Kolmogorov distances between the distributions of the partial sums ∑i=1nXi and ∑i=1nYi. In the case where the distributions of the Xis and the Yis are compared with respect to the convex order, the proposed upper bounds are further refined. Finally, in order to illustrate the applicability of the results presented, we consider specific examples concerning gamma and normal approximations.

In the present paper, results given in Natvig (1990) are generalized to a multistate, strongly coherent, nonrepairable system of independent components by considering the reduction in remaining system time above a certain state due to a jump downwards of a component. This reduction also equals the increase in remaining system time above a certain state due to a minimal repair of the component at its time of jump downwards. The expected value of the sum of such reductions/increases for the different possible jumps downwards of the component is the building block of the Natvig measure of the importance of the component in the multistate case. Hence, the whole distributions of these reductions/increases are arrived at, not only the expectations, throwing more light on the consequences for the system of the deterioration of the components.

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.

Consider a device which begins operating at time t = 0 and then operates continuously in time. The operation of the device is satisfactory so long as a physical parameter X(t) of the device remains less than a critical value c > 0 (X(0) = 0). The device does not operate satisfactorily (i.e., it is “failed”) at a time t if X(t) ≧ c. Suppose that X(t) is a random variable with mean EX(t) = ta (a> 0), i.e., the device parameter X drifts with a linear rate a.

We study a production-inventory system in which input is deterministic, and its rate is the controlled parameter. The output is a compound Poisson process with exponential jump-size distribution. The parameters of the output process are unknown. The following adaptive policy is used: At each time t the unknown parameters are estimated by maximum likelihood method and the input rate which would be optimal if these were the true values of the parameter is used. The issues investigated are (1) the convergence of maximum likelihood estimates to the true value, and (2) the asymptotic properties of the cost under the adaptive policy.

In this paper we consider the factor analysis for Lévy-driven multivariate price models with stochastic volatility. Our main aim is to provide conditions on the volatility process under which we can possibly reduce the dimension of the driving Lévy motion. We find that these conditions depend on a particular form of the multivariate Lévy process. In some settings we concentrate on nondegenerate symmetric α-stable Lévy motions.

In this paper we present a set of results relating to the occupation time α(t) of a process X(·). The first set of results concerns exact characterizations of α(t), e.g. in terms of its transform up to an exponentially distributed epoch. In addition, we establish a central limit theorem (entailing that a centered and normalized version of α(t)∕t converges to a zero-mean normal random variable as t→∞) and the tail asymptotics of ℙ(α(t)∕t≥q). We apply our findings to spectrally positive Lévy processes reflected at the infimum and establish various new occupation time results for the corresponding model.

Distributions of some random variables have been characterized by independence of certain functions of these random variables. For example, let X and Y be two independent and identically distributed random variables having the gamma distribution. Laha showed that U = X + Y and V = X | Y are also independent random variables. Lukacs showed that U and V are independently distributed if, and only if, X and Y have the gamma distribution. Ferguson characterized the exponential distribution in terms of the independence of X – Y and min (X, Y). The best-known of these characterizations is that first proved by Kac which states that if random variables X and Y are independent, then X + Y and X – Y are independent if, and only if, X and Y are jointly Gaussian with the same variance. In this paper, Kac's hypotheses have been somewhat modified. In so doing, we obtain a larger class of distributions which we shall call class λ1. A subclass λ0 of λ1 enjoys many nice properties of the Gaussian distribution, in particular, in non-linear filtering.

We analyze here the optimal premium rates for an insurance company for the case where the claims process is a stochastic process with unknown parameters. The beliefs of the insurance company about these parameters are represented by some prior distribution function, and the premium rates are determined on the basis of this distribution. The insurer learns about the unknown parameters from the past behavior of the claims process, and revises the prior distribution whenever new information is obtained (i.e., a claims stream is observed).

The learning process introduces multi-period dynamic considerations to the problem of the rate determination, since the rates determined at any period affect the claims process, and thus the information available for future decisions. We show how the optimal rates can be obtained by using some dynamic programming techniques. We also prove that the optimal rates obtained in our case are lower than those determined by an insurance company that does not revise its distribution.