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.

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.

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 describe methods for randomly thinning certain classes of spatial point processes. In the case of a Markov point process, the proposed method involves a dependent thinning of a spatial birth-and-death process, where clans of ancestors associated with the original points are identified, and where we simulate backwards and forwards in order to obtain the thinned process. In the case of a Cox process, a simple independent thinning technique is proposed. In both cases, the thinning results in a Poisson process if and only if the true Papangelou conditional intensity is used, and, thus, can be used as a graphical exploratory tool for inspecting the goodness-of-fit of a spatial point process model. Several examples, including clustered and inhibitive point processes, are considered.

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

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.

Consider a family of random ordered graph trees (T n )n≥1, where T n has n vertices. It has previously been established that if the associated search-depth processes converge to the normalised Brownian excursion when rescaled appropriately as n → ∞, then the simple random walks on the graph trees have the Brownian motion on the Brownian continuum random tree as their scaling limit. Here, this result is extended to demonstrate the existence of a diffusion scaling limit whenever the volume measure on the limiting real tree is nonatomic, supported on the leaves of the limiting tree, and satisfies a polynomial lower bound for the volume of balls. Furthermore, as an application of this generalisation, it is established that the simple random walks on a family of Galton-Watson trees with a critical infinite variance offspring distribution, conditioned on the total number of offspring, can be rescaled to converge to the Brownian motion on a related α-stable tree.

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.

In this paper the theory of species sampling sequences is linked to the theory of conditionally identically distributed sequences in order to enlarge the set of species sampling sequences which are mathematically tractable. The conditional identity in distribution (see Berti, Pratelli and Rigo (2004)) is a new type of dependence for random variables, which generalizes the well-known notion of exchangeability. In this paper a class of random sequences, called generalized species sampling sequences, is defined and a condition to have conditional identity in distribution is given. Moreover, two types of generalized species sampling sequence that are conditionally identically distributed are introduced and studied: the generalized Poisson-Dirichlet sequence and the generalized Ottawa sequence. Some examples are discussed.

Asymptotic independence of the components of random vectors is a concept used in many applications. The standard criteria for checking asymptotic independence are given in terms of distribution functions (DFs). DFs are rarely available in an explicit form, especially in the multivariate case. Often we are given the form of the density or, via the shape of the data clouds, we can obtain a good geometric image of the asymptotic shape of the level sets of the density. In this paper we establish a simple sufficient condition for asymptotic independence for light-tailed densities in terms of this asymptotic shape. This condition extends Sibuya's classic result on asymptotic independence for Gaussian densities.

In this paper, we consider an insurance company whose surplus (reserve) is modeled by a jump diffusion risk process. The insurance company can invest part of its surplus in n risky assets and purchase proportional reinsurance for claims. Our main goal is to find an optimal investment and proportional reinsurance policy which minimizes the ruin probability. We apply stochastic control theory to solve this problem. We obtain the closed form expression for the minimal ruin probability, optimal investment and proportional reinsurance policy. We find that the minimal ruin probability satisfies the Lundberg equality. We also investigate the effects of the diffusion volatility parameter, the market price of risk and the correlation coefficient on the minimal ruin probability, optimal investment and proportional reinsurance policy through numerical calculations.

For a time series, a plot of sample covariances is a popular way to assess its dependence properties. In this paper we give a systematic characterization of the asymptotic behavior of sample covariances of long-memory linear processes. Central and noncentral limit theorems are obtained for sample covariances with bounded as well as unbounded lags. It is shown that the limiting distribution depends in a very interesting way on the strength of dependence, the heavy-tailedness of the innovations, and the magnitude of the lags.

We present new ideas about the type of random tessellation which evolves through successive division of its cells. These ideas are developed in an intuitive way, with many pictures and only a modicum of mathematical formalism–so that the wide application of the ideas is clearly apparent to all readers. A vast number of new tessellation models, with known probability distribution for the volume of the typical cell, follow from the concepts in this paper. There are other interesting models for which results are not presented (or presented only through simulation methods), but these models have illustrative value. A large agenda of further research is opened up by the ideas in this paper.

We consider variance-optimal hedging in general continuous-time affine stochastic volatility models. The optimal hedge and the associated hedging error are determined semiexplicitly in the case that the stock price follows a martingale. The integral representation of the solution opens the door to efficient numerical computation. The setup includes models with jumps in the stock price and in the activity process. It also allows for correlation between volatility and stock price movements. Concrete parametric models will be illustrated in a forthcoming paper.

We consider a system with Poisson arrivals and independent and identically distributed service times, where requests in the system are served according to the state-dependent (Cohen's generalized) processor-sharing discipline, where each request receives a service capacity that depends on the actual number of requests in the system. For this system, we derive expressions as well as tight insensitive upper bounds for the moments of the conditional sojourn time of a request with given required service time. The bounds generalize and extend corresponding results, recently given for the single-server processor-sharing system in Cheung et al. (2006) and for the state-dependent processor-sharing system with exponential service times by the authors (2008). Analogous results hold for the waiting times. Numerical examples for the M/M/m-PS and M/D/m-PS systems illustrate the given bounds.

Consider a monoecious diploid population with nonoverlapping generations, whose size varies with time according to an irreducible, aperiodic Markov chain with states x 1 N,…,x K N, where K ≪ N. It is assumed that all matings except for selfing are possible and equally probable. At time 0 a random sample of n ≪ N genes is taken. Given two successive population sizes x j N and x i N, the numbers of gametes that individual parents contribute to offspring can be shown to be exchangeable random variables distributed as G ij . Under minimal conditions on the first three moments of G ij for all i and j, a suitable effective population size N e is derived. Then if time is recorded in a backward direction in units of 2N e generations, it can be shown that coalescent theory holds.

In this paper we propose a framework that facilitates the study of large deviations for point processes based on stationary sequences with regularly varying tails. This framework allows us to keep track both of the magnitude of the extreme values of a process and the order in which these extreme values appear. Particular emphasis is put on (infinite) linear processes with random coefficients. The proposed framework provides a fairly complete description of the joint asymptotic behavior of the large values of the stationary sequence. We apply the general result on large deviations for point processes to derive the asymptotic decay of certain probabilities related to partial sum processes as well as ruin probabilities.

Consider a sequence of exchangeable or independent binary trials ordered on a line or on a circle. The statistics denoting the number of times an F-S string of length (at least)k 1 + k 2, that is, (at least)k 1 failures followed by (at least) k 2 successes in n such trials, are studied. The associated waiting time for the rth occurrence of an F-S string of length (at least) k 1 + k 2 in linearly ordered trials is also examined. Exact formulae, lower/upper bounds and approximations are derived for their distributions. Mean values and variances of the number of occurrences of F-S strings are given in exact formulae too. Particular exchangeable and independent sequences of binary random variables, used in applied research, combined with numerical examples clarify further the theoretical results.

In this paper we obtain an approximation for the duration distribution of the excursion set generated by the minimum of two independent F random processes above a high threshold u. Moreover, we obtain a closed-form approximation for the mean duration of the conjunction of these two F processes. As an illustration, we conduct a simulation study to compare the performances of the approximated distribution and the exact distribution.