It is well known and readily seen that the maximum of n independent and uniformly on [0, 1] distributed random variables, suitably standardised, converges in total variation distance, as n increases, to the standard negative exponential distribution. We extend this result to higher dimensions by considering copulas. We show that the strong convergence result holds for copulas that are in a differential neighbourhood of a multivariate generalised Pareto copula. Sklar’s theorem then implies convergence in variational distance of the maximum of n independent and identically distributed random vectors with arbitrary common distribution function and (under conditions on the marginals) of its appropriately normalised version. We illustrate how these convergence results can be exploited to establish the almost-sure consistency of some estimation procedures for max-stable models, using sample maxima.

Hawkes processes have been widely used in many areas, but their probability properties can be quite difficult. In this paper an elementary approach is presented to obtain moments of Hawkes processes and/or the intensity of a number of marked Hawkes processes, in which the detailed outline is given step by step; it works not only for all Markovian Hawkes processes but also for some non-Markovian Hawkes processes. The approach is simpler and more convenient than usual methods such as the Dynkin formula and martingale methods. The method is applied to one-dimensional Hawkes processes and other related processes such as Cox processes, dynamic contagion processes, inhomogeneous Poisson processes, and non-Markovian cases. Several results are obtained which may be useful in studying Hawkes processes and other counting processes. Our proposed method is an extension of the Dynkin formula, which is simple and easy to use.

By a random process with immigration at random times we mean a shot noise process with a random response function (response process) in which shots occur at arbitrary random times. Such random processes generalize random processes with immigration at the epochs of a renewal process which were introduced in Iksanov et al. (2017) and bear a strong resemblance to a random characteristic in general branching processes and the counting process in a fixed generation of a branching random walk generated by a general point process. We provide sufficient conditions which ensure weak convergence of finite-dimensional distributions of these processes to certain Gaussian processes. Our main result is specialised to several particular instances of random times and response processes.

We study the problem of stopping a Brownian bridge X in order to maximise the expected value of an exponential gain function. The problem was posed by Ernst and Shepp (2015), and was motivated by bond selling with non-negative prices.

Due to the non-linear structure of the exponential gain, we cannot rely on methods used in the literature to find closed-form solutions to other problems involving the Brownian bridge. Instead, we must deal directly with a stopping problem for a time-inhomogeneous diffusion. We develop techniques based on pathwise properties of the Brownian bridge and martingale methods of optimal stopping theory, which allow us to find the optimal stopping rule and to show the regularity of the value function.

By a general shot noise process we mean a shot noise process in which the counting process of shots is arbitrary locally finite. Assuming that the counting process of shots satisfies a functional limit theorem in the Skorokhod space with a locally Hölder continuous Gaussian limit process, and that the response function is regularly varying at infinity, we prove that the corresponding general shot noise process satisfies a similar functional limit theorem with a different limit process and different normalization and centering functions. For instance, if the limit process for the counting process of shots is a Brownian motion, then the limit process for the general shot noise process is a Riemann–Liouville process. We specialize our result for five particular counting processes. Also, we investigate Hölder continuity of the limit processes for general shot noise processes.

A famous result in renewal theory is the central limit theorem for renewal processes. Since, in applications, usually only observations from a finite time interval are available, a bound on the Kolmogorov distance to the normal distribution is desirable. We provide an explicit non-uniform bound for the renewal central limit theorem based on Stein’s method and track the explicit values of the constants. For this bound the inter-arrival time distribution is required to have only a second moment. As an intermediate result of independent interest we obtain explicit bounds in a non-central Berry–Esseen theorem under second moment conditions.

We explore the first passage problem for sticky reflecting diffusion processes with double exponential jumps. The joint Laplace transform of the first passage time to an upper level and the corresponding overshoot is studied. In particular, explicit solutions are presented when the diffusion part is driven by a drifted Brownian motion and by an Ornstein–Uhlenbeck process.

In the first part of this paper we study approximations of trajectories of piecewise deterministic processes (PDPs) when the flow is not given explicitly by the thinning method. We also establish a strong error estimate for PDPs as well as a weak error expansion for piecewise deterministic Markov processes (PDMPs). These estimates are the building blocks of the multilevel Monte Carlo (MLMC) method, which we study in the second part. The coupling required by the MLMC is based on the thinning procedure. In the third part we apply these results to a two-dimensional Morris–Lecar model with stochastic ion channels. In the range of our simulations the MLMC estimator outperforms classical Monte Carlo.

We consider point process convergence for sequences of independent and identically distributed random walks. The objective is to derive asymptotic theory for the largest extremes of these random walks. We show convergence of the maximum random walk to the Gumbel or the Fréchet distributions. The proofs depend heavily on precise large deviation results for sums of independent random variables with a finite moment generating function or with a subexponential distribution.

Residential burglary is a social problem in every major urban area. As such, progress has been to develop quantitative, informative and applicable models for this type of crime: (1) the Deterministic-time-step (DTS) model [Short, D’Orsogna, Pasour, Tita, Brantingham, Bertozzi & Chayes (2008) Math. Models Methods Appl. Sci.18, 1249–1267], a pioneering agent-based statistical model of residential burglary criminal behaviour, with deterministic time steps assumed for arrivals of events in which the residential burglary aggregate pattern formation is quantitatively studied for the first time; (2) the SSRB model (agent-based stochastic-statistical model of residential burglary crime) [Wang, Zhang, Bertozzi & Short (2019) Active Particles, Vol. 2, Springer Nature Switzerland AG, in press], in which the stochastic component of the model is theoretically analysed by introduction of a Poisson clock with time steps turned into exponentially distributed random variables. To incorporate independence of agents, in this work, five types of Poisson clocks are taken into consideration. Poisson clocks (I), (II) and (III) govern independent agent actions of burglary behaviour, and Poisson clocks (IV) and (V) govern interactions of agents with the environment. All the Poisson clocks are independent. The time increments are independently exponentially distributed, which are more suitable to model individual actions of agents. Applying the method of merging and splitting of Poisson processes, the independent Poisson clocks can be treated as one, making the analysis and simulation similar to the SSRB model. A Martingale formula is derived, which consists of a deterministic and a stochastic component. A scaling property of the Martingale formulation with varying burglar population is found, which provides a theory to the finite size effects. The theory is supported by quantitative numerical simulations using the pattern-formation quantifying statistics. Results presented here will be transformative for both elements of application and analysis of agent-based models for residential burglary or in other domains.

Let X = (Xt)t≥0 be a stochastic process issued from $x \in \mathbb{R}$ that admits a marginal stationary measure v, i.e. vPtf = vf for all t ≥ 0, where $\textbf{P}_t\,f(x)= \mathbb{E}_x[f(\textbf{X}_t)]$ . In this paper, we introduce the (resp. biorthogonal) spectral projections correlation functions which are expressed in terms of projections.” Also, update first published online date, if available. into the eigenspaces of Pt (resp. and of its adjoint in the weighted Hilbert space L2 (v)). We obtain closed-form expressions involving eigenvalues, the condition number and/or the angle between the projections in the following different situations: when X = X with X = (Xt)t ≥ 0 being a Markov process, X is the subordination of X in the sense of Bochner, and X is a non-Markovian process which is obtained by time-changing X with an inverse of a subordinator. It turns out that these spectral projections correlation functions have different expressions with respect to these classes of processes which enables to identify substantial and deep properties about their dynamics. This interesting fact can be used to design original statistical tests to make inferences, for example, about the path properties of the process (presence of jumps), distance from symmetry (self-adjoint or non-self-adjoint) and short-to-long-range dependence. To reveal the usefulness of our results, we apply them to a class of non-self-adjoint Markov semigroups studied in Patie and Savov (to appear, Mem. Amer. Math. Soc., 179p), and then time-change by subordinators and their inverses.

We determine the order of magnitude of $\mathbb{E}|\sum _{n\leqslant x}f(n)|^{2q}$, where $f(n)$ is a Steinhaus or Rademacher random multiplicative function, and $0\leqslant q\leqslant 1$. In the Steinhaus case, this is equivalent to determining the order of $\lim _{T\rightarrow \infty }\frac{1}{T}\int _{0}^{T}|\sum _{n\leqslant x}n^{-it}|^{2q}\,dt$.

In particular, we find that $\mathbb{E}|\sum _{n\leqslant x}f(n)|\asymp \sqrt{x}/(\log \log x)^{1/4}$. This proves a conjecture of Helson that one should have better than squareroot cancellation in the first moment and disproves counter-conjectures of various other authors. We deduce some consequences for the distribution and large deviations of $\sum _{n\leqslant x}f(n)$.

The proofs develop a connection between $\mathbb{E}|\sum _{n\leqslant x}f(n)|^{2q}$ and the $q$th moment of a critical, approximately Gaussian, multiplicative chaos and then establish the required estimates for that. We include some general introductory discussion about critical multiplicative chaos to help readers unfamiliar with that area.

In this paper, we study the nonlinear diffusion equation associated with a particle system where the common drift depends on the rate of absorption of particles at a boundary. We provide an interpretation of this equation, which is also related to the supercooled Stefan problem, as a structural credit risk model with default contagion in a large interconnected banking system. Using the method of heat potentials, we derive a coupled system of Volterra integral equations for the transition density and for the loss through absorption. An approximation by expansion is given for a small interaction parameter. We also present a numerical solution algorithm and conduct computational tests.

We consider the behaviour of minimax recursions defined on random trees. Such recursions give the value of a general class of two-player combinatorial games. We examine in particular the case where the tree is given by a Galton–Watson branching process, truncated at some depth 2n, and the terminal values of the level 2n nodes are drawn independently from some common distribution. The case of a regular tree was previously considered by Pearl, who showed that as n → ∞ the value of the game converges to a constant, and by Ali Khan, Devroye and Neininger, who obtained a distributional limit under a suitable rescaling.

For a general offspring distribution, there is a surprisingly rich variety of behaviour: the (unrescaled) value of the game may converge to a constant, or to a discrete limit with several atoms, or to a continuous distribution. We also give distributional limits under suitable rescalings in various cases.

We also address questions of endogeny. Suppose the game is played on a tree with many levels, so that the terminal values are far from the root. To be confident of playing a good first move, do we need to see the whole tree and its terminal values, or can we play close to optimally by inspecting just the first few levels of the tree? The answers again depend in an interesting way on the offspring distribution.

We also mention several open questions.

We present a discrete-type approximation scheme to solve continuous-time optimal stopping problems based on fully non-Markovian continuous processes adapted to the Brownian motion filtration. The approximations satisfy suitable variational inequalities which allow us to construct $\varepsilon$ -optimal stopping times and optimal values in full generality. Explicit rates of convergence are presented for optimal values based on reward functionals of path-dependent stochastic differential equations driven by fractional Brownian motion. In particular, the methodology allows us to design concrete Monte Carlo schemes for non-Markovian optimal stopping time problems as demonstrated in the companion paper by Bezerra et al.

This paper generalizes the Kunita–Watanabe decomposition of an $L^2$ space. The generalization comes from using nonlinear stochastic integrals where the integrator is a family of continuous martingales bounded in $L^2$ . This result is also the solution of an optimization problem in $L^2$ . First, martingales are assumed to be stochastic integrals. Then, to get the general result, it is shown that the regularity of the family of martingales with respect to its spatial parameter is inherited by the integrands in the integral representation of the martingales. Finally, an example showing how the results of this paper, with the Clark–Ocone formula, can be applied to polynomial functions of Brownian integrals.

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.