We derive a new theoretical lower bound for the expected supremum of drifted fractional Brownian motion with Hurst index $H\in(0,1)$ over a (in)finite time horizon. Extensive simulation experiments indicate that our lower bound outperforms the Monte Carlo estimates based on very dense grids for $H\in(0,\tfrac{1}{2})$. Additionally, we derive the Paley–Wiener–Zygmund representation of a linear fractional Brownian motion in the general case and give an explicit expression for the derivative of the expected supremum at $H=\tfrac{1}{2}$ in the sense of Bisewski, Dȩbicki and Rolski (2021).

Pseudo cross-variograms appear naturally in the context of multivariate Brown–Resnick processes, and are a useful tool for analysis and prediction of multivariate random fields. We give a necessary and sufficient criterion for a matrix-valued function to be a pseudo cross-variogram, and further provide a Schoenberg-type result connecting pseudo cross-variograms and multivariate correlation functions. By means of these characterizations, we provide extensions of the popular univariate space–time covariance model of Gneiting to the multivariate case.

In this paper we study the asymptotic behaviour of a random uniform parking function $\pi_n$ of size n. We show that the first $k_n$ places $\pi_n(1),\ldots,\pi_n(k_n)$ of $\pi_n$ are asymptotically independent and identically distributed (i.i.d.) and uniform on $\{1,2,\ldots,n\}$, for the total variation distance when $k_n = {\rm{o}}(\sqrt{n})$, and for the Kolmogorov distance when $k_n={\rm{o}}(n)$, improving results of Diaconis and Hicks. Moreover, we give bounds for the rate of convergence, as well as limit theorems for certain statistics such as the sum or the maximum of the first $k_n$ parking places. The main tool is a reformulation using conditioned random walks.

We consider a spatial model of cancer in which cells are points on the d-dimensional torus $\mathcal{T}=[0,L]^d$, and each cell with $k-1$ mutations acquires a kth mutation at rate $\mu_k$. We assume that the mutation rates $\mu_k$ are increasing, and we find the asymptotic waiting time for the first cell to acquire k mutations as the torus volume tends to infinity. This paper generalizes results on waiting for $k\geq 3$ mutations in Foo et al. (2020), which considered the case in which all of the mutation rates $\mu_k$ are the same. In addition, we find the limiting distribution of the spatial distances between mutations for certain values of the mutation rates.

In this paper, we analyze a two-queue random time-limited Markov-modulated polling model. In the first part of the paper, we investigate the fluid version: fluid arrives at the two queues as two independent flows with deterministic rate. There is a single server that serves both queues at constant speeds. The server spends an exponentially distributed amount of time in each queue. After the completion of such a visit time to one queue, the server instantly switches to the other queue, i.e., there is no switch-over time.

For this model, we first derive the Laplace–Stieltjes transform (LST) of the stationary marginal fluid content/workload at each queue. Subsequently, we derive a functional equation for the LST of the two-dimensional workload distribution that leads to a Riemann–Hilbert boundary value problem (BVP). After taking a heavy-traffic limit, and restricting ourselves to the symmetric case, the BVP simplifies and can be solved explicitly.

In the second part of the paper, allowing for more general (Lévy) input processes and server switching policies, we investigate the transient process limit of the joint workload in heavy traffic. Again solving a BVP, we determine the stationary distribution of the limiting process. We show that, in the symmetric case, this distribution coincides with our earlier solution of the BVP, implying that in this case the two limits (stationarity and heavy traffic) commute.

We study a stochastic differential equation with an unbounded drift and general Hölder continuous noise of order $\lambda \in (0,1)$. The corresponding equation turns out to have a unique solution that, depending on a particular shape of the drift, either stays above some continuous function or has continuous upper and lower bounds. Under some mild assumptions on the noise, we prove that the solution has moments of all orders. In addition, we provide its connection to the solution of some Skorokhod reflection problem. As an illustration of our results and motivation for applications, we also suggest two stochastic volatility models which we regard as generalizations of the CIR and CEV processes. We complete the study by providing a numerical scheme for the solution.

We consider the asset price as the weak solution to a stochastic differential equation driven by both a Brownian motion and the counting process martingale whose predictable compensator follows shot-noise and Hawkes processes. In this framework, we discuss the Esscher martingale measure where the conditions for its existence are detailed. This generalizes certain relationships not yet encountered in the literature.

We study an optimal reinsurance problem for a diffusion model, in which the drift of the claim follows an Ornstein–Uhlenbeck process. The aim of the insurer is to maximize the expected exponential utility of its terminal wealth. We consider two cases: full information and partial information. Full information occurs when the insurer directly observes the drift; partial information occurs when the insurer observes only its claims. By applying stochastic control and by solving the corresponding Hamilton–Jacobi–Bellman equations, we find the value function and the optimal reinsurance strategy under both full and partial information. We determine a relationship between the value function and reinsurance strategy under full information with the value function and reinsurance strategy under partial information.

This paper addresses the asymptotic analysis of sojourn functionals of spatiotemporal Gaussian random fields with long-range dependence (LRD) in time, also known as long memory. Specifically, reduction theorems are derived for local functionals of nonlinear transformation of such fields, with Hermite rank $m\geq 1,$ under general covariance structures. These results are proven to hold, in particular, for a family of nonseparable covariance structures belonging to the Gneiting class. For $m=2,$ under separability of the spatiotemporal covariance function in space and time, the properly normalized Minkowski functional, involving the modulus of a Gaussian random field, converges in distribution to the Rosenblatt-type limiting distribution for a suitable range of values of the long-memory parameter.

We consider a dual risk model with constant expense rate and i.i.d. exponentially distributed gains $C_i$ ($i=1,2,\dots$) that arrive according to a renewal process with general interarrival times. We add to this classical dual risk model the proportional gain feature; that is, if the surplus process just before the ith arrival is at level u, then for $a>0$ the capital jumps up to the level $(1+a)u+C_i$. The ruin probability and the distribution of the time to ruin are determined. We furthermore identify the value of discounted cumulative dividend payments, for the case of a Poisson arrival process of proportional gains. In the dividend calculations, we also consider a random perturbation of our basic risk process modeled by an independent Brownian motion with drift.

This article examines the impact of the largest claims reinsurance treaties on loss reserve of the ceding company. The largest claims reinsurance, known as LCR, and ECOMOR reinsurance treaties are considered to be the two most appropriate reinsurance treaties for large or catastrophe claims. Then, it studies the impact of such treaties on loss reserves. Through a simulation study, it shown that, under a more general situation, the LCR treaty can be a more efficient (in some sense, see below) treaty than the ECOMOR treaty for the ceding company.

We clarify and refine the definition of a reciprocal random field on an undirected graph, with the reciprocal chain as a special case, by introducing four new properties: the factorizing, global, local, and pairwise reciprocal properties, in decreasing order of strength, with respect to a set of nodes $\delta$. They reduce to the better-known Markov properties if $\delta$ is the empty set, or, with the exception of the local property, if $\delta$ is a complete set. Conditions for each reciprocal property to imply the next stronger property are derived, and it is shown that, conditionally on the values at a set of nodes $\delta_0$, all four properties are preserved for the subgraph induced by the remaining nodes, with respect to the node set $\delta\setminus\delta_0$. We note that many of the above results are new even for reciprocal chains.

In this study, we consider a class of multiple-drawing opposite-reinforcing urns with time-dependent replacement rules. The class has the symmetric property of a Friedman-type urn. We divide the class into a small-increment regime and a large-increment regime. For small-increment schemes, we prove almost-sure convergence and a central limit theorem for the proportion of white balls by stochastic approximation. For large-increment schemes, by assuming the affinity condition, we show almost-sure convergence of the proportion of white balls by martingale theory and present a way to identify the limit distribution of the proportion of white balls.

Given a spectrally negative Lévy process, we predict, in an $L_1$ sense, the last passage time of the process below zero before an independent exponential time. This optimal prediction problem generalises [2], where the infinite-horizon problem is solved. Using a similar argument as that in [24], we show that this optimal prediction problem is equivalent to solving an optimal prediction problem in a finite-horizon setting. Surprisingly (unlike the infinite-horizon problem), an optimal stopping time is based on a curve that is killed at the moment the mean of the exponential time is reached. That is, an optimal stopping time is the first time the process crosses above a non-negative, continuous, and non-increasing curve depending on time. This curve and the value function are characterised as a solution of a system of nonlinear integral equations which can be understood as a generalisation of the free boundary equations (see e.g. [21, Chapter IV.14.1]) in the presence of jumps. As an example, we numerically calculate this curve in the Brownian motion case and for a compound Poisson process with exponential-sized jumps perturbed by a Brownian motion.

We consider solutions of Lévy-driven stochastic differential equations of the form $\textrm{d} X_t=\sigma(X_{t-})\textrm{d} L_t$, $X_0=x$, where the function $\sigma$ is twice continuously differentiable and the driving Lévy process $L=(L_t)_{t\geq0}$ is either vector or matrix valued. While the almost sure short-time behavior of Lévy processes is well known and can be characterized in terms of the characteristic triplet, there is no complete characterization of the behavior of the solution X. Using methods from stochastic calculus, we derive limiting results for stochastic integrals of the form $t^{-p}\int_{0+}^t\sigma(X_{t-})\,\textrm{d} L_t$ to show that the behavior of the quantity $t^{-p}(X_t-X_0)$ for $t\downarrow0$ almost surely reflects the behavior of $t^{-p}L_t$. Generalizing $t^{{\kern1pt}p}$ to a suitable function $f\colon[0,\infty)\rightarrow\mathbb{R}$ then yields a tool to derive explicit law of the iterated logarithm type results for the solution from the behavior of the driving Lévy process.

The invariant Galton–Watson (IGW) tree measures are a one-parameter family of critical Galton–Watson measures invariant with respect to a large class of tree reduction operations. Such operations include the generalized dynamical pruning (also known as hereditary reduction in a real tree setting) that eliminates descendant subtrees according to the value of an arbitrary subtree function that is monotone nondecreasing with respect to an isometry-induced partial tree order. We show that, under a mild regularity condition, the IGW measures are attractors of arbitrary critical Galton–Watson measures with respect to the generalized dynamical pruning. We also derive the distributions of height, length, and size of the IGW trees.

We introduce a variant of Shepp’s classical urn problem in which the optimal stopper does not know whether sampling from the urn is done with or without replacement. By considering the problem’s continuous-time analog, we provide bounds on the value function and, in the case of a balanced urn (with an equal number of each ball type), an explicit solution is found. Surprisingly, the optimal strategy for the balanced urn is the same as in the classical urn problem. However, the expected value upon stopping is lower due to the additional uncertainty present.

We give an example of an FIID vertex-labeling of ${\mathbb T}_3$ whose marginals are uniform on $[0,1]$, and if we delete the edges between those vertices whose labels are different, then some of the remaining clusters are infinite. We also show that no such process can be finitary.

We present a reflection principle for a wide class of symmetric random motions with finite velocities. We propose a deterministic argument which is then applied to trajectories of stochastic processes. In the case of symmetric correlated random walks and the symmetric telegraph process, we provide a probabilistic result recalling the classical reflection principle for Brownian motion, but where the initial velocity has a crucial role. In the case of the telegraph process we also present some consequences which lead to further reflection-type characteristics of the motion.