We consider the model of Deijfen, Häggström and Bagley (2004) for competing growth of two infection types in R^d, based on the Richardson model on Z^d. Stochastic ball-shaped infection outbursts transmit the infection type of the center to all points of the ball that are not yet infected. Relevant parameters of the model are the initial infection configuration, the (type-dependent) growth rates, and the radius distribution of the infection outbursts. The main question is that of coexistence: Which values of the parameters allow the unbounded growth of both types with positive probability? Deijfen, Häggström and Bagley (2004) conjectured that the initial configuration is basically irrelevant for this question, and gave a proof for this under strong assumptions on the radius distribution, which, e.g. do not include the case of a deterministic radius. Here we give a proof that does not rely on these assumptions. One of the tools to be used is a slight generalization of the model with immune regions and delayed initial infection configurations.

Consider the region L = {(x ,y) : 0 ≤ y ≤ Clog(1 + x), x > 0} for a constant C > 0. We study the percolation and coverage properties of this region. For the coverage properties, we place a Poisson point process of intensity λ on the entire half space R₊ x R and associated with each Poisson point we place a box of a random side length ρ. Depending on the tail behaviour of the random variable ρ we exhibit a phase transition in the intensity for the eventual coverage of the region L. For the percolation properties, we place a Poisson point process of intensity λ on the region R². At each point of the process we centre a box of a random side length ρ. In the case ρ ≤ R for some fixed R > 0 we study the critical intensity λ_c of the percolation on L.

We consider statistical inference for a parametric cooperative sequential adsorption model for spatial time series data, based on maximum likelihood. We establish asymptotic normality of the maximum likelihood estimator in the thermodynamic limit. We also perform and discuss some numerical simulations of the model, which illustrate the procedure for creating confidence intervals for large samples.

A forwards induction policy is a type of greedy algorithm for Markov decision processes. It is straightforward to implement and is optimal for a large class of models, especially in stochastic resource allocation. In this paper we consider a model for the optimal allocation of resources in pre-clinical pharmaceutical research. We show that although they are not always strictly optimal, forwards induction policies perform well.

Let η = (η₁,...,η_n) be a positive random vector. If its coordinates η_i and η_j are exchangeable, i.e. the distribution of η is invariant with respect to the swap π_ij of its ith and jth coordinates, then Ef(η) = Ef(π_ijη) for all integrable functions f. In this paper we study integrable random vectors that satisfy this identity for a particular family of functions f, namely those which can be written as the positive part of the scalar product 〈u, η〉 with varying weights u. In finance such functions represent payoffs from exchange options with η being the random part of price changes, while from the geometric point of view they determine the support function of the so-called zonoid of η. If the expected values of such payoffs are π_ij-invariant, we say that η is ij-swap-invariant. A full characterisation of the swap-invariance property and its relationship to the symmetries of expected payoffs of basket options are obtained. The first of these results relies on a characterisation theorem for integrable positive random vectors with equal zonoids. Particular attention is devoted to the case of asset prices driven by Lévy processes. Based on this, concrete semi-static hedging techniques for multi-asset barrier options, such as weighted barrier swap options, weighted barrier quanto-swap options, or certain weighted barrier spread options, are suggested.

For a given bivariate Lévy process (U_t, L_t)_t≥0, distributional properties of the stationary solutions of the stochastic differential equation dV_t = V_t-dU_t + dL_t are analysed. In particular, the expectation and autocorrelation function are obtained in terms of the process (U, L) and in several cases of interest the tail behavior is described. In the case where U has jumps of size -1, necessary and sufficient conditions for the law of the solutions to be (absolutely) continuous are given.

This paper is concerned with the behaviour of a Lévy process when it crosses over a positive level, u, starting from 0, both as u becomes large and as u becomes small. Our main focus is on the time, τ_u, it takes the process to transit above the level, and in particular, on the stability of this passage time; thus, essentially, whether or not τ_u behaves linearly as u ↓ 0 or u → ∞. We also consider the conditional stability of τ_u when the process drifts to -∞ almost surely. This provides information relevant to quantities associated with the ruin of an insurance risk process, which we analyse under a Cramér condition.

We present a new perfect simulation algorithm for stationary chains having unbounded variable length memory. This is the class of infinite memory chains for which the family of transition probabilities is represented by a probabilistic context tree. We do not assume any continuity condition: our condition is expressed in terms of the structure of the context tree. More precisely, the length of the contexts is a deterministic function of the distance to the last occurrence of some determined string of symbols. It turns out that the resulting class of chains can be seen as a natural extension of the class of chains having a renewal string. In particular, our chains exhibit a visible regeneration scheme.

We consider the problem of maximizing the probability of stopping on any of the last m successes in independent Bernoulli trials with random horizon of length N, where m is a predetermined integer. A prior is given for N. It is known that, when N is degenerate, i.e. P{N = n} = 1 for a given n > m, the sum-the-multiplicative-odds theorem gives the solution and shows that the optimal rule is a threshold rule, i.e. it stops on the first success appearing after a given stage. However, when N is nondegenerate, the optimal rule is not necessarily a threshold rule. So our main concern in Section 2 is to give a sufficient condition for the optimal rule to be a threshold rule when N is a bounded random variable such that P{N ≤ n} = 1. Application will be made to the usual (discrete arrival time) secretary problem with a random number N of applicants in Section 3. When N is uniform or curtailed geometric, the optimal rules are shown to be threshold rules and their asymptotic results are obtained. We also examine, as a nonhomogeneous Poisson process model, an intermediate prior that allows N to be uniform or degenerate. In Section 4 we consider a continuous arrival time version of the secretary problem with a random number M of applicants. It is shown that, whatever the distribution of M, we can win with probability greater than or equal to u_m^*, where u_m^* is, as given in (1.4), the asymptotic win probability of the usual secretary problem when N degenerates to n and n → ∞.

In this paper we consider the scaled limit of a continuous-time random walk (CTRW) based on a Markov chain {X_n, n ≥ 0} and two observables, τ(∙) and V(∙), corresponding to the renewal times and jump sizes. Assuming that these observables belong to the domains of attraction of some stable laws, we give sufficient conditions on the chain that guarantee the existence of the scaled limits for CTRWs. An application of the results to a process that arises in quantum transport theory is provided. The results obtained in this paper generalize earlier results contained in Becker-Kern, Meerschaert and Scheffler (2004) and Meerschaert and Scheffler (2008), and the recent results of Henry and Straka (2011) and Jurlewicz, Kern, Meerschaert and Scheffler (2010), where {X_n, n ≥ 0} is a sequence of independent and identically distributed random variables.

We introduce a new point process, the dynamic contagion process, by generalising the Hawkes process and the Cox process with shot noise intensity. Our process includes both self-excited and externally excited jumps, which could be used to model the dynamic contagion impact from endogenous and exogenous factors of the underlying system. We have systematically analysed the theoretical distributional properties of this new process, based on the piecewise-deterministic Markov process theory developed in Davis (1984), and the extension of the martingale methodology used in Dassios and Jang (2003). The analytic expressions of the Laplace transform of the intensity process and the probability generating function of the point process have been derived. An explicit example of specified jumps with exponential distributions is also given. The object of this study is to produce a general mathematical framework for modelling the dependence structure of arriving events with dynamic contagion, which has the potential to be applicable to a variety of problems in economics, finance, and insurance. We provide an application of this process to credit risk, and a simulation algorithm for further industrial implementation and statistical analysis.

We consider a stochastic control model for a queueing system driven by a two-dimensional fractional Brownian motion with Hurst parameter 0 < H < 1. In particular, when H > ½, this model serves to approximate a controlled two-station tandem queueing model with heavy-tailed ON/OFF sources in heavy traffic. We establish the weak convergence results for the distribution of the state process and construct an explicit stationary state process associated with given controls. Based on suitable coupling arguments, we show that each state process couples with its stationary counterpart and we use it to represent the long-run average cost functional in terms of the stationary process. Finally, we establish the existence result of an optimal control, which turns out to be independent of the initial data.

In the so-called lightbulb process, on days r = 1,...,n, out of n lightbulbs, all initially off, exactly r bulbs, selected uniformly and independent of the past, have their status changed from off to on, or vice versa. With X the number of bulbs on at the terminal time n, an even integer, and μ = n/2, σ² = var(X), we have sup_{z ∈ R} | P((X - μ)/σ ≤ z) - P(Z ≤ z) | ≤ nΔ̅₀/2σ² + 1.64n/σ³ + 2/σ, where Z is a standard normal random variable and Δ̅₀ = 1/2√n + 1/2n + e^-n/2/3 for n ≥ 6, yielding a bound of order O(n^-1/2) as n → ∞. A similar, though slightly larger bound, holds for odd n. The results are shown using a version of Stein's method for bounded, monotone size bias couplings. The argument for even n depends on the construction of a variable X^s on the same space as X that has the X-size bias distribution, that is, which satisfies E[Xg(X)] = μE[g(X^s)] for all bounded continuous g, and for which there exists a B ≥ 0, in this case B = 2, such that X ≤ X^s ≤ X + B almost surely. The argument for odd n is similar to that for even n, but one first couples X closely to V, a symmetrized version of X, for which a size bias coupling of V to V^s can proceed as in the even case. In both the even and odd cases, the crucial calculation of the variance of a conditional expectation requires detailed information on the spectral decomposition of the lightbulb chain.