We consider the problem of efficient estimation via simulation of first passage time probabilities for a multidimensional random walk with heavy-tailed increments. In addition to being a natural generalization to the problem of computing ruin probabilities in insurance - in which the focus is the maximum of a one-dimensional random walk with negative drift - this problem captures important features of large deviations for multidimensional heavy-tailed processes (such as the role played by the mean of the process in connection to the location of the target set). We develop a state-dependent importance sampling estimator for this class of multidimensional problems. Then, using techniques based on Lyapunov inequalities, we argue that our estimator is strongly efficient in the sense that the relative mean squared error of our estimator can be made arbitrarily small by increasing the number of replications, uniformly as the probability of interest approaches 0.

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.

We study a stochastic differential game between two insurance companies who employ reinsurance to reduce the risk of exposure. Under the assumption that the companies have large insurance portfolios compared to any individual claim size, their surplus processes can be approximated by stochastic differential equations. We formulate competition between the two companies as a game with a single payoff function which depends on the surplus processes. One company chooses a dynamic reinsurance strategy in order to maximize this expected payoff, while the other company simultaneously chooses a dynamic reinsurance strategy so as to minimize the same quantity. We describe the Nash equilibrium of this stochastic differential game and solve it explicitly for the case of maximizing/minimizing the exit probability.

We consider a Poisson cluster model, motivated by insurance applications. At each claim arrival time, modeled by the point of a homogeneous Poisson process, we start a cluster process which represents the number or amount of payments triggered by the arrival of a claim in a portfolio. The cluster process is a Lévy or truncated compound Poisson process. Given the observations of the process over a finite interval, we consider the expected value of the number and amount of payments in a future time interval. We also give bounds for the error encountered in this prediction procedure.

We introduce and motivate the study of (n + 1) × r arrays X with Bernoulli entries Xk,j and independently distributed rows. We study the distribution of which denotes the number of consecutive pairs of successes (or runs of length 2) when reading the array down the columns and across the rows. With the case r = 1 having been studied by several authors, and permitting some initial inferences for the general case r > 1, we examine various distributional properties and representations of Sn for the case r = 2, and, using a more explicit analysis, the case of multinomial and identically distributed rows. Applications are also given in cases where the array X arises from a Pólya sampling scheme.

In Reinert and Röllin (2009) a new approach - called the ‘embedding method’ - was introduced, which allows us to make use of exchangeable pairs for normal and multivariate normal approximations with Stein's method in cases where the corresponding couplings do not satisfy a certain linearity condition. The key idea is to embed the problem into a higher-dimensional space in such a way that the linearity condition is then satisfied. Here we apply the embedding to U-statistics as well as to subgraph counts in random graphs.

We analyze the notion of ‘reliability prediction’ by studying in detail a key property that is tacitly assumed to make reliability prediction possible. The analysis leads in turn to a special type of point process for which the connection of future to past can be explicitly displayed. In this type of process, the semi-renewal process, all finite-dimensional distributions are completely determined by the distribution of the time to the first event in the process. The theory provides a heretofore unappreciated unification of the two most commonly used reliability prediction models for maintained systems, namely, the renewal and revival processes. We show that familiar results from renewal theory extend and generalize to semi-renewal processes.

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 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.

We consider a Markov additive process (MAP) with phase-type jumps, starting at 0. Given a positive level u, we determine the joint distribution of the undershoot and overshoot of the first jump over the level u, the maximal level before this jump, the time of attaining this maximum, and the time between the maximum and the jump. The analysis is based on first passage times and time reversion of MAPs. A marginal of the derived distribution is the Gerber-Shiu function, which is of interest to insurance risk. Several examples serve to compare the present result with the literature.

We take a new look at transient, or time-dependent Little laws for queueing systems. Through the use of Palm measures, we show that previous laws (see Bertsimas and Mourtzinou (1997)) can be generalized. Furthermore, within this framework, a new law can be derived as well, which gives higher-moment expressions for very general types of queueing system; in particular, the laws hold for systems that allow customers to overtake one another. What is especially novel about our approach is the use of Palm measures that are induced by nonstationary point processes, as these measures are not commonly found in the queueing literature. This new higher-moment law is then used to provide expressions for all moments of the number of customers in the system in an M/G/1 preemptive last-come-first-served queue at a time t > 0, for any initial condition and any of the more famous preemptive disciplines (i.e. preemptive-resume, and preemptive-repeat with and without resampling) that are analogous to the special cases found in Abate and Whitt (1987c), (1988). These expressions are then used to derive a nice structural form for all of the time-dependent moments of a regulated Brownian motion (see Abate and Whitt (1987a), (1987b)).

A semi-numerical method is derived to compute the Laplace transform of the equilibrium busy period probability density function in a fluid queue with constant output rate when the buffer is nonempty. The input process is controlled by a continuous-time semi-Markov chain (CTSMC) with n states such that in each state the input rate is constant. The holding time in states with net positive output rate - so-called emptying states - is assumed to be an exponentially distributed random variable, whereas in states with net positive input rate - so-called filling states - it may have an arbitrary probability distribution. The result is demonstrated by applying it to various systems, including fluid queues with two on-off input sources. The latter exercise in part shows consistency with prior results but also solves the problem in the case where there are two emptying states. Numerical results are presented for selected examples which expose discontinuities in the busy period distribution when the number of emptying states changes, e.g. as a result of increasing the fluid arrival rate in one or more states of the controlling CTSMC.

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.

We provide sufficient conditions which ensure that the intrinsic martingale in the supercritical branching random walk converges exponentially fast to its limit. We include in particular the case of Galton-Watson processes so that our results can be seen as a generalization of a result given in the classical treatise by Asmussen and Hering (1983). As an auxiliary tool, we prove ultimate versions of two results concerning the exponential renewal measures which may be of interest in themselves and which correct, generalize, and simplify some earlier works.

The conditional least-squares estimators of the variances are studied for a critical branching process with immigration that allows the offspring distributions to have infinite fourth moments. We derive different forms of limiting distributions for these estimators when the offspring distributions have regularly varying tails with index α. In particular, in the case in which 2 < α < 8/3, the normalizing factor of the estimator for the offspring variance is smaller than √n, which is different from that of Winnicki (1991).

Motivated by a problem arising in the mining industry, we present a first study of the energy required to reduce a unit mass fragment by consecutively using several devices. Two devices are considered, which we represent as different stochastic fragmentation processes. Following the self-similar energy model introduced in Bertoin and Martínez (2005), we compute the average energy required to attain a size η0 with this two-device procedure. We then asymptotically compare, as η0 goes to 0 or 1, its energy requirement with that of individual fragmentation processes. In particular, we show that, for a certain range of parameters of the fragmentation processes and of their energy cost functions, the consecutive use of two devices can be asymptotically more efficient than using each of them separately, or vice versa.

We consider a random walks system on ℤ in which each active particle performs a nearest-neighbor random walk and activates all inactive particles it encounters. The movement of an active particle stops when it reaches a certain number of jumps without activating any particle. We prove that if the process relies on efficient particles (i.e. those particles with a small probability of jumping to the left) being placed strategically on ℤ, then it might survive, having active particles at any time with positive probability. On the other hand, we may construct a process that dies out eventually almost surely, even if it relies on efficient particles. That is, we discuss what happens if particles are initially placed very far away from each other or if their probability of jumping to the right tends to 1 but not fast enough.

We consider a preferential duplication model for growing random graphs, extending previous models of duplication graphs by selecting the vertex to be duplicated with probability proportional to its degree. We show that a special case of this model can be analysed using the same stochastic approximation as for vertex-reinforced random walks, and show that ‘trapping’ behaviour can occur, such that the descendants of a particular group of initial vertices come to dominate the graph.

We study the optimal liquidation strategy for a call spread in the case when an investor, who does not hedge, believes in a volatility that differs from the implied volatility. The liquidation problem is formulated as an optimal stopping problem, which we solve explicitly. We also provide a sensitivity analysis with respect to the model parameters.

Expectations of unbounded functions of dependent nonnegative integer-valued random variables are approximated by the expectations of the functions of independent copies of these random variables. The Lindeberg method is used.