The no restart random walk (NRRW) is a random network growth model driven by a random walk that builds the graph while moving on it, adding and connecting a new leaf node to the current position of the walker every s steps. We show a fundamental dichotomy in NRRW with respect to the parity of s: for ${s}=1$ we prove that the random walk is transient and non-leaf nodes have degrees bounded above by an exponential distribution; for s even we prove that the random walk is recurrent and non-leaf nodes have degrees bounded below by a power law distribution. These theoretical findings highlight and confirm the diverse and rich behaviour of NRRW observed empirically.

For spectrally negative Lévy processes, we prove several fluctuation results involving a general draw-down time, which is a downward exit time from a dynamic level that depends on the running maximum of the process. In particular, we find expressions of the Laplace transforms for the two-sided exit problems involving the draw-down time. We also find the Laplace transforms for the hitting time and creeping time over the running-maximum related draw-down level, respectively, and obtain an expression for a draw-down associated potential measure. The results are expressed in terms of scale functions for the spectrally negative Lévy processes.

We extend the construction principle of phase-type (PH) distributions to allow for inhomogeneous transition rates and show that this naturally leads to direct probabilistic descriptions of certain transformations of PH distributions. In particular, the resulting matrix distributions enable the carrying over of fitting properties of PH distributions to distributions with heavy tails, providing a general modelling framework for heavy-tail phenomena. We also illustrate the versatility and parsimony of the proposed approach in modelling a real-world heavy-tailed fire insurance dataset.

As an extension of a central limit theorem established by Svante Janson, we prove a Berry–Esseen inequality for a sum of independent and identically distributed random variables conditioned by a sum of independent and identically distributed integer-valued random variables.

A large and sparse random graph with independent exponentially distributed link weights can be used to model the propagation of messages or diseases in a network with an unknown connectivity structure. In this article we study an extended setting where, in addition, the nodes of the graph are equipped with nonnegative random weights which are used to model the effect of boundary delays across paths in the network. Our main results provide approximative formulas for typical first passage times, typical flooding times, and maximum flooding times in the extended setting, over a time scale logarithmic with respect to the network size.

We introduce an idealised model for overland flow generated by rain falling on a hillslope. Our prime motivation is to show how the coalescence of runoff streams promotes the total generation of runoff. We show that, for our model, as the rate of rainfall increases in relation to the soil infiltration rate there is a distinct phase change. For low rainfall (the subcritical case) only the bottom of the hillslope contributes to the total overland runoff, while for high rainfall (the supercritical case) the whole slope contributes and the total runoff increases dramatically. We identify the critical point at which the phase change occurs, and show how it depends on the degree of coalescence. When there is no stream coalescence the critical point occurs when the rainfall rate equals the average infiltration rate, but when we allow coalescence the critical point occurs when the rainfall rate is less than the average infiltration rate, and increasing the amount of coalescence increases the total expected runoff.

We introduce variance-optimal semi-static hedging strategies for a given contingent claim. To obtain a tractable formula for the expected squared hedging error and the optimal hedging strategy we use a Fourier approach in a multidimensional factor model. We apply the theory to set up a variance-optimal semi-static hedging strategy for a variance swap in the Heston model, which is affine, in the 3/2 model, which is not, and in a market model including jumps.

We prove existence and uniqueness of a stationary distribution and absolute regularity for nonlinear GARCH and INGARCH models of order (p, q). In contrast to previous work we impose, besides a geometric drift condition, only a semi-contractive condition which allows us to include models which would be ruled out by a fully contractive condition. This results in a subgeometric rather than the more usual geometric decay rate of the mixing coefficients. The proofs are heavily based on a coupling of two versions of the processes.

For a spectrally negative Lévy process X, killed according to a rate that is a function ω of its position, we complement the recent findings of [12] by analysing (in greater generality) the exit probability of the one-sided upwards passage problem. When ω is strictly positive, this problem is related to the determination of the Laplace transform of the first passage time upwards for X that has been time-changed by the inverse of the additive functional $$\int_0^ \cdot \omega ({X_u}){\kern 1pt} {\rm{d}}u$$. In particular, our findings thus shed extra light on related results concerning first passage times downwards (resp. upwards) of continuous-state branching processes (resp. spectrally negative positive self-similar Markov processes).

We study the distribution of the negative Wiener–Hopf factor for a class of two-sided jump Lévy processes whose positive jumps have a rational Laplace transform. The positive Wiener–Hopf factor for this class of processes was studied by Lewis and Mordecki (2008). Here we obtain a formula for the Laplace transform of the negative Wiener–Hopf factor, as well as an explicit expression for its probability density in terms of sums of convolutions of known functions. Under additional regularity conditions on the Lévy measure of the studied processes, we also provide asymptotic results as $u\to-\infty$ for the distribution function F(u) of the negative Wiener–Hopf factor. We illustrate our results in some particular examples.

We consider decision problems with arbitrary action spaces, deterministic transitions, and infinite time horizon. In the usual setup when probability measures are countably additive, a general version of Kuhn’s theorem implies under fairly general conditions that for every mixed strategy of the decision maker there exists an equivalent behavior strategy. We examine to what extent this remains valid when probability measures are only assumed to be finitely additive. Under the classical approach of Dubins and Savage (2014), we prove the following statements: (1) If the action space is finite, every mixed strategy has an equivalent behavior strategy. (2) Even if the action space is infinite, at least one optimal mixed strategy has an equivalent behavior strategy. The approach by Dubins and Savage turns out to be essentially maximal: these two statements are no longer valid if we take any extension of their approach that considers all singleton plays.

We present the first algorithm that samples maxn≥0{Sn − nα}, where Sn is a mean zero random walk, and nα with $\alpha \in ({1 \over 2},1)$ defines a nonlinear boundary. We show that our algorithm has finite expected running time. We also apply this algorithm to construct the first exact simulation method for the steady-state departure process of a GI/GI/∞ queue where the service time distribution has infinite mean.

We show that several machine learning estimators, including square-root least absolute shrinkage and selection and regularized logistic regression, can be represented as solutions to distributionally robust optimization problems. The associated uncertainty regions are based on suitably defined Wasserstein distances. Hence, our representations allow us to view regularization as a result of introducing an artificial adversary that perturbs the empirical distribution to account for out-of-sample effects in loss estimation. In addition, we introduce RWPI (robust Wasserstein profile inference), a novel inference methodology which extends the use of methods inspired by empirical likelihood to the setting of optimal transport costs (of which Wasserstein distances are a particular case). We use RWPI to show how to optimally select the size of uncertainty regions, and as a consequence we are able to choose regularization parameters for these machine learning estimators without the use of cross validation. Numerical experiments are also given to validate our theoretical findings.

We study negative association for mixed sampled point processes and show that negative association holds for such processes if a random number of their points fulfils the ultra log-concave (ULC) property. We connect the negative association property of point processes with directionally convex dependence ordering, and show some consequences of this property for mixed sampled and determinantal point processes. Some applications illustrate the general theory.

We derive moment identities for the stochastic integrals of multiparameter processes in a random-connection model based on a point process admitting a Papangelou intensity. The identities are written using sums over partitions, and they reduce to sums over non-flat partition diagrams if the multiparameter processes vanish on diagonals. As an application, we obtain general identities for the moments of k-hop counts in the random-connection model, which simplify the derivations available in the literature.

We study the asymptotic behaviour of a class of small-noise diffusions driven by fractional Brownian motion, with random starting points. Different scalings allow for different asymptotic properties of the process (small-time and tail behaviours in particular). In order to do so, we extend some results on sample path large deviations for such diffusions. As an application, we show how these results characterise the small-time and tail estimates of the implied volatility for rough volatility models, recently proposed in mathematical finance.

It is well understood that a supercritical continuous-state branching process (CSBP) is equal in law to a discrete continuous-time Galton–Watson process (the skeleton of prolific individuals) whose edges are dressed in a Poissonian way with immigration which initiates subcritical CSBPs (non-prolific mass). Equally well understood in the setting of CSBPs and superprocesses is the notion of a spine or immortal particle dressed in a Poissonian way with immigration which initiates copies of the original CSBP, which emerges when conditioning the process to survive eternally. In this article we revisit these notions for CSBPs and put them in a common framework using the well-established language of (coupled) stochastic differential equations (SDEs). In this way we are able to deal simultaneously with all types of CSBPs (supercritical, critical, and subcritical) as well as understanding how the skeletal representation becomes, in the sense of weak convergence, a spinal decomposition when conditioning on survival. We have two principal motivations. The first is to prepare the way to expand the SDE approach to the spatial setting of superprocesses, where recent results have increasingly sought the use of skeletal decompositions to transfer results from the branching particle setting to the setting of measure valued processes. The second is to provide a pathwise decomposition of CSBPs in the spirit of genealogical coding of CSBPs via Lévy excursions, albeit precisely where the aforesaid coding fails to work because the underlying CSBP is supercritical.

Consider a random vector $\textbf{U}$ whose distribution function coincides in its upper tail with that of an Archimedean copula. We report the fact that the conditional distribution of $\textbf{U}$ , conditional on one of its components, has under a mild condition on the generator function independent upper tails, no matter what the unconditional tail behavior is. This finding is extended to Archimax copulas.

The signature of a coherent system has been studied extensively in the recent literature. Signatures are particularly useful in the comparison of coherent or mixed systems under a variety of stochastic orderings. Also, certain signature-based closure and preservation theorems have been established. For example, it is now well known that certain stochastic orderings are preserved from signatures to system lifetimes when components have independent and identical distributions. This applies to the likelihood ratio order, the hazard rate order, and the stochastic order. The point of departure of the present paper is the question of whether or not a similar preservation result will hold for the mean residual life order. A counterexample is provided which shows that the answer is negative. Classes of distributions for the component lifetimes for which the latter implication holds are then derived. Connections to the theory of order statistics are also considered.

This paper considers two variants of M/G/1 queues with impatient customers, which are denoted by M/G/1+Gw and M/G/1+Gs. In the M/G/1+Gw queue customers have deadlines for their waiting times, and they leave the system immediately if their services do not start before the expiration of their deadlines. On the other hand, in the M/G/1+Gs queue customers have deadlines for their sojourn times, where customers in service also immediately leave the system when their deadlines expire. In this paper we derive comparison results for performance measures of these models. In particular, we show that if the service time distribution is new better than used in expectation, then the loss probability in the M/G/1+Gs queue is greater than that in the M/G/1+Gw queue.