Consider the following migration process based on a closed network of N queues with $K_N$ customers. Each station is a $\cdot$/M/$\infty$ queue with service (or migration) rate $\mu$. Upon departure, a customer is routed independently and uniformly at random to another station. In addition to migration, these customers are subject to a susceptible–infected–susceptible (SIS) dynamics. That is, customers are in one of two states: I for infected, or S for susceptible. Customers can swap their state either from I to S or from S to I only in stations. More precisely, at any station, each susceptible customer becomes infected with the instantaneous rate $\alpha Y$ if there are Y infected customers in the station, whereas each infected customer recovers and becomes susceptible with rate $\beta$. We let N tend to infinity and assume that $\lim_{N\to \infty} K_N/N= \eta $, where $\eta$ is a positive constant representing the customer density. The main problem of interest concerns the set of parameters of such a system for which there exists a stationary regime where the epidemic survives in the limiting system. The latter limit will be referred to as the thermodynamic limit. We use coupling and stochastic monotonicity arguments to establish key properties of the associated Markov processes, which in turn allow us to give the structure of the phase transition diagram of this thermodynamic limit with respect to $\eta$. The analysis of the Kolmogorov equations of this SIS model reduces to that of a wave-type PDE for which we have found no explicit solution. This plain SIS model is one among several companion stochastic processes that exhibit both random migration and contagion. Two of them are discussed in the present paper as they provide variants to the plain SIS model as well as some bounds and approximations. These two variants are the departure-on-change-of-state (DOCS) model and the averaged-infection-rate (AIR) model, which both admit closed-form solutions. The AIR system is a classical mean-field model where the infection mechanism based on the actual population of infected customers is replaced by a mechanism based on some empirical average of the number of infected customers in all stations. The latter admits a product-form solution. DOCS features accelerated migration in that each change of SIS state implies an immediate departure. This model leads to another wave-type PDE that admits a closed-form solution. In this text, the main focus is on the closed stochastic networks and their limits. The open systems consisting of a single station with Poisson input are instrumental in the analysis of the thermodynamic limits and are also of independent interest. This class of SIS dynamics has incarnations in virtually all queueing networks of the literature.

We study the geometric and topological features of U-statistics of order k when the k-tuples satisfying geometric and topological constraints do not occur frequently. Using appropriate scaling, we establish the convergence of U-statistics in vague topology, while the structure of a non-degenerate limit measure is also revealed. Our general result shows various limit theorems for geometric and topological statistics, including persistent Betti numbers of Čech complexes, the volume of simplices, a functional of the Morse critical points, and values of the min-type distance function. The required vague convergence can be obtained as a result of the limit theorem for point processes induced by U-statistics. The latter convergence particularly occurs in the $\mathcal M_0$ -topology.

Across a wide variety of applications, the self-exciting Hawkes process has been used to model phenomena in which the history of events influences future occurrences. However, there may be many situations in which the past events only influence the future as long as they remain active. For example, a person spreads a contagious disease only as long as they are contagious. In this paper, we define a novel generalization of the Hawkes process that we call the ephemerally self-exciting process. In this new stochastic process, the excitement from one arrival lasts for a randomly drawn activity duration, hence the ephemerality. Our study includes exploration of the process itself as well as connections to well-known stochastic models such as branching processes, random walks, epidemics, preferential attachment, and Bayesian mixture models. Furthermore, we prove a batch scaling construction of general, marked Hawkes processes from a general ephemerally self-exciting model, and this novel limit theorem both provides insight into the Hawkes process and motivates the model contained herein as an attractive self-exciting process in its own right.

Network dynamics with point-process-based interactions are of paramount modeling interest. Unfortunately, most relevant dynamics involve complex graphs of interactions for which an exact computational treatment is impossible. To circumvent this difficulty, the replica-mean-field approach focuses on randomly interacting replicas of the networks of interest. In the limit of an infinite number of replicas, these networks become analytically tractable under the so-called ‘Poisson hypothesis’. However, in most applications this hypothesis is only conjectured. In this paper we establish the Poisson hypothesis for a general class of discrete-time, point-process-based dynamics that we propose to call fragmentation-interaction-aggregation processes, and which are introduced here. These processes feature a network of nodes, each endowed with a state governing their random activation. Each activation triggers the fragmentation of the activated node state and the transmission of interaction signals to downstream nodes. In turn, the signals received by nodes are aggregated to their state. Our main contribution is a proof of the Poisson hypothesis for the replica-mean-field version of any network in this class. The proof is obtained by establishing the propagation of asymptotic independence for state variables in the limit of an infinite number of replicas. Discrete-time Galves–Löcherbach neural networks are used as a basic instance and illustration of our analysis.

Innovations in medicine provide us longer and healthier life, leading lower mortality. Sooner rather than later, much greater longevity would be possible for us due to artificial intelligence advances in health care. Similarly, Advanced Driver Assistance Systems (ADAS) in highly automated vehicles may reduce or even eventually eliminate accidents by perceiving dangerous situations, which would minimize the number of accidents and lead to fewer loss claims for insurance companies. To model the survivor function capturing greater longevity as well as the number of claims reflecting less accidents in the long run, in this paper, we study a Cox process whose intensity process is piecewise-constant and decreasing. We derive its ultimate distributional properties, such as the Laplace transform of intensity integral process, the probability generating function of point process, their associated moments and cumulants, and the probability of no more claims for a given time point. In general, this simple model may be applicable in many other areas for modeling the evolution of gradually disappearing events, such as corporate defaults, dividend payments, trade arrivals, employment of a certain job type (e.g., typists) in the labor market, and release of particles. In particular, we discuss some potential applications to insurance.

Consider a homogeneous Poisson point process of the Euclidean plane and its Voronoi tessellation. The present note discusses the properties of two stationary point processes associated with the latter and depending on a parameter $\theta$ . The first is the set of points that belong to some one-dimensional facet of the Voronoi tessellation and such that the angle with which they see the two nuclei defining the facet is $\theta$ . The main question of interest on this first point process is its intensity. The second point process is that of the intersections of the said tessellation with a straight line having a random orientation. Its intensity is well known. The intersection points almost surely belong to one-dimensional facets. The main question here concerns the Palm distribution of the angle with which the points of this second point process see the two nuclei associated with the facet. We will give answers to these two questions and briefly discuss their practical motivations. We also discuss natural extensions to three dimensions.

We consider the problem of numerical integration when the sampling nodes form a stationary point process on the real line. In previous papers it was argued that a naïve Riemann sum approach can cause a severe variance inflation when the sampling points are not equidistant. We show that this inflation can be avoided using a higher-order Newton–Cotes quadrature rule which exploits smoothness properties of the integrand. Under mild assumptions, the resulting estimator is unbiased and its variance asymptotically obeys a power law as a function of the mean point distance. If the Newton–Cotes rule is of sufficiently high order, the exponent of this law turns out to only depend on the point process through its mean point distance. We illustrate our findings with the stereological estimation of the volume of a compact object, suggesting alternatives to the well-established Cavalieri estimator.

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.

In this paper, we introduce a new large family of Lévy-driven point processes with (and without) contagion, by generalising the classical self-exciting Hawkes process and doubly stochastic Poisson processes with non-Gaussian Lévy-driven Ornstein–Uhlenbeck-type intensities. The resulting framework may possess many desirable features such as skewness, leptokurtosis, mean-reverting dynamics, and more importantly, the ‘contagion’ or feedback effects, which could be very useful for modelling event arrivals in finance, economics, insurance, and many other fields. We characterise the distributional properties of this new class of point processes and develop an efficient sampling method for generating sample paths exactly. Our simulation scheme is mainly based on the distributional decomposition of the point process and its intensity process. Extensive numerical implementations and tests are reported to demonstrate the accuracy and effectiveness of our scheme. Moreover, we use portfolio risk management as an example to show the applicability and flexibility of our algorithms.

We consider a branching random walk on a multitype (with Q types of particles), supercritical Galton–Watson tree which satisfies the Kesten–Stigum condition. We assume that the displacements associated with the particles of type Q have regularly varying tails of index $\alpha$, while the other types of particles have lighter tails than the particles of type Q. In this paper we derive the weak limit of the sequence of point processes associated with the positions of the particles in the nth generation. We verify that the limiting point process is a randomly scaled scale-decorated Poisson point process using the tools developed by Bhattacharya, Hazra, and Roy (2018). As a consequence, we obtain the asymptotic distribution of the position of the rightmost particle in the nth generation.

We introduce and study the class of branching-stable point measures, which can be seen as an analog of stable random variables when the branching mechanism for point measures replaces the usual addition. In contrast with the classical theory of stable (Lévy) processes, there exists a rich family of branching-stable point measures with a negative scaling exponent, which can be described as certain Crump‒Mode‒Jagers branching processes. We investigate the asymptotic behavior of their cumulative distribution functions, that is, the number of atoms in (-∞, x] as x→∞, and further depict the genealogical lineage of typical atoms. For both results, we rely crucially on the work of Biggins (1977), (1992).

In this paper we propose a model for biological neural nets where the activity of the network is described by Hawkes processes having a variable length memory. The particularity in this paper is that we deal with an infinite number of components. We propose a graphical construction of the process and build, by means of a perfect simulation algorithm, a stationary version of the process. To implement this algorithm, we make use of a Kalikow-type decomposition technique. Two models are described in this paper. In the first model, we associate to each edge of the interaction graph a saturation threshold that controls the influence of a neuron on another. In the second model, we impose a structure on the interaction graph leading to a cascade of spike trains. Such structures, where neurons are divided into layers, can be found in the retina.

Consider a family of Boolean models, indexed by integers n≥1. The nth model features a Poisson point process in ℝn of intensity e{nρn}, and balls of independent and identically distributed radii distributed like X̅n√n. Assume that ρn→ρ as n→∞, and that X̅n satisfies a large deviations principle. We show that there then exist the three deterministic thresholds τd, the degree threshold, τp, the percolation probability threshold, and τv, the volume fraction threshold, such that, asymptotically as n tends to ∞, we have the following features. (i) For ρ<τd, almost every point is isolated, namely its ball intersects no other ball; (ii) for τd<ρ<τp, the mean number of balls intersected by a typical ball converges to ∞ and nevertheless there is no percolation; (iii) for τp<ρ<τv, the volume fraction is 0 and nevertheless percolation occurs; (iv) for τd<ρ<τv, the mean number of balls intersected by a typical ball converges to ∞ and nevertheless the volume fraction is 0; (v) for ρ>τv, the whole space is covered. The analysis of this asymptotic regime is motivated by problems in information theory, but it could be of independent interest in stochastic geometry. The relations between these three thresholds and the Shannon‒Poltyrev threshold are discussed.

We consider random variables observed at arrival times of a renewal process, which possibly depends on those observations and has regularly varying steps with infinite mean. Due to the dependence and heavy-tailed steps, the limiting behavior of extreme observations until a given time t tends to be rather involved. We describe the asymptotics and extend several partial results which appeared in this setting. The theory is applied to determine the asymptotic distribution of maximal excursions and sojourn times for continuous-time random walks.

In the full rectangular version of Gilbert's planar tessellation (see Gilbert (1967), Mackisack and Miles (1996), and Burridge et al. (2013)), lines extend either horizontally (with east- and west-growing rays) or vertically (north- and south-growing rays) from seed points which form a stationary Poisson point process, each ray stopping when it meets another ray that has blocked its path. In the half-Gilbert rectangular version, east- and south-growing rays, whilst having the potential to block each other, do not interact with west and north rays, and vice versa. East- and south-growing rays have a reciprocality of blocking, as do west and north. In this paper we significantly expand and simplify the half-Gilbert analytic results that we gave in Burridge et al. (2013). We also show how the idea of reciprocality of blocking between rays can be used in a much wider context, with rays not necessarily orthogonal and with seeds that produce more than two rays.

In this paper, a unique dataset of improvised explosive device attacks during “The Troubles” in Northern Ireland (NI) is analysed via a Hawkes process model. It is found that this past dependent model is a good fit to improvised explosive device attacks yielding key insights about the nature of terrorism in NI. We also present a novel approach to quantitatively investigate some of the sociological theory surrounding the Provisional Irish Republican Army which challenges previously held assumptions concerning changes seen in the organisation. Finally, we extend our use of the Hawkes process model by considering a multidimensional version which permits both self and mutual-excitations. This allows us to test how the Provisional Irish Republican Army responded to past improvised explosive device attacks on different geographical scales from which we find evidence for the autonomy of the organisation over the six counties of NI and Belfast. By incorporating a second dataset concerning British Security Force (BSF) interventions, the multidimensional model allows us to test counter-terrorism (CT) operations in NI where we find subsequent increases in violence.

In this paper, we explore some of the various issues that may occur in attempting to fit a dynamical systems (either agent- or continuum-based) model of urban crime to data on just the attack times and locations. We show how one may carry out a regression analysis for the model described by Short et al. (2008, Math. Mod. Meth. Appl. Sci.) by using simulated attack data from the agent-based model. It is discussed how one can incorporate the attack data into the partial differential equations for the expected attractiveness to burgle and the criminal density to predict crime rates between attacks. Using this predicted crime rate, we derive a likelihood function that one can maximise in order to fit parameters and/or initial conditions for the model. We focus on carrying out data assimilation for two different parameter regions, namely in the case where stationary and non-stationary crime hotspots form. It is found that the likelihood function is ‘flat’ for large ranges of parameters, and that this has major implications for crime forecasting. Hence, we look at how one might carry out a goodness-of-fit and forecasting analysis for crime rates given the range of parameter fits. We show how one can use the Kolmogorov–Smirnov statistic to assess the goodness-of-fit. The dynamical systems analysis of the partial differential equations proves invaluable to understanding how the crime rate forecasts depend on the parameters and their sensitivity. Finally, we outline several interesting directions for future research in this area where we believe that the combination of dynamical systems modelling, analysis, and data assimilation can prove effective in developing policing strategies for urban crime.

Consider a real-valued discrete-time stationary and ergodic stochastic process, called the noise process. For each dimension n, we can choose a stationary point process in ℝn and a translation invariant tessellation of ℝn. Each point is randomly displaced, with a displacement vector being a section of length n of the noise process, independent from point to point. The aim is to find a point process and a tessellation that minimizes the probability of decoding error, defined as the probability that the displaced version of the typical point does not belong to the cell of this point. We consider the Shannon regime, in which the dimension n tends to ∞, while the logarithm of the intensity of the point processes, normalized by dimension, tends to a constant. We first show that this problem exhibits a sharp threshold: if the sum of the asymptotic normalized logarithmic intensity and of the differential entropy rate of the noise process is positive, then the probability of error tends to 1 with n for all point processes and all tessellations. If it is negative then there exist point processes and tessellations for which this probability tends to 0. The error exponent function, which denotes how quickly the probability of error goes to 0 in n, is then derived using large deviations theory. If the entropy spectrum of the noise satisfies a large deviations principle, then, below the threshold, the error probability goes exponentially fast to 0 with an exponent that is given in closed form in terms of the rate function of the noise entropy spectrum. This is obtained for two classes of point processes: the Poisson process and a Matérn hard-core point process. New lower bounds on error exponents are derived from this for Shannon's additive noise channel in the high signal-to-noise ratio limit that hold for all stationary and ergodic noises with the above properties and that match the best known bounds in the white Gaussian noise case.

In this paper we consider the stationary Poisson Boolean model with spherical grains and propose a family of nonparametric estimators for the radius distribution. These estimators are based on observed distances and radii, weighted in an appropriate way. They are ratio unbiased and asymptotically consistent for a growing observation window. We show that the asymptotic variance exists and is given by a fairly explicit integral expression. Asymptotic normality is established under a suitable integrability assumption on the weight function. We also provide a short discussion of related estimators as well as a simulation study.

Consider the following random spatial network: in a large disk, construct a network using a stationary and isotropic Poisson line process of unit intensity. Connect pairs of points using the network, with initial/final segments of the connecting path formed by travelling off the network in the opposite direction to that of the destination/source. Suppose further that connections are established using ‘near geodesics’, constructed between pairs of points using the perimeter of the cell containing these two points and formed using only the Poisson lines not separating them. If each pair of points generates an infinitesimal amount of traffic divided equally between the two connecting near geodesics, and if the Poisson line pattern is conditioned to contain a line through the centre, then what can be said about the total flow through the centre? In Kendall (2011) it was shown that a scaled version of this flow has asymptotic distribution given by the 4-volume of a region in 4-space, constructed using an improper anisotropic Poisson line process in an infinite planar strip. Here we construct a more amenable representation in terms of two ‘seminal curves’ defined by the improper Poisson line process, and establish results which produce a framework for effective simulation from this distribution up to an L1 error which tends to 0 with increasing computational effort.