Search

We introduce a variant of the Hawkes-fed birth-death process, in which the conditional intensity does not increase at arrivals, but at departures from the system. Since arrivals cause excitation after a delay equal to their lifetimes, we call this a delayed Hawkes process. We introduce a general family of models admitting a cluster representation containing the Hawkes, delayed Hawkes and ephemerally self-exciting processes as special cases. For this family of models, as well as their nonlinear extensions, we prove existence, uniqueness and stability. Our family of models satisfies the same FCLT as the classical Hawkes process; however, we describe a scaling limit for the delayed Hawkes process in which sojourn times are stretched out by a factor $\sqrt T$, after which time gets contracted by a factor $T$. This scaling limit highlights the effect of sojourn-time dependence. The cluster representation renders our family of models tractable, allowing for transform characterisation by a fixed-point equation and for an analysis of heavy-tailed asymptotics. In the Markovian case, for a multivariate network of delayed Hawkes birth-death processes, an explicit recursive procedure is presented to calculate the $d$th-order moments analytically. Finally, we compare the delayed Hawkes process to the regular Hawkes process in the stochastic ordering, which enables us to describe stationary distributions and heavy-traffic behaviour., Comment: 38 pages, 1 figure

Point processes model the occurrence of a countable number of random points over some support. They can model diverse phenomena, such as chemical reactions, stock market transactions and social interactions. We show that JumpProcesses.jl is a fast, general-purpose library for simulating point processes. JumpProcesses.jl was first developed for simulating jump processes via stochastic simulation algorithms (SSAs) (including Doob's method, Gillespie's methods, and Kinetic Monte Carlo methods). Historically, jump processes have been developed in the context of dynamical systems to describe dynamics with discrete jumps. In contrast, the development of point processes has been more focused on describing the occurrence of random events. In this paper, we bridge the gap between the treatment of point and jump process simulation. The algorithms previously included in JumpProcesses.jl can be mapped to three general methods developed in statistics for simulating evolutionary point processes. Our comparative exercise revealed that the library initially lacked an efficient algorithm for simulating processes with variable intensity rates. We, therefore, extended JumpProcesses.jl with a new simulation algorithm, Coevolve, that enables the rapid simulation of processes with locally-bounded variable intensity rates. It is now possible to efficiently simulate any point process on the real line with a non-negative, left-continuous, history-adapted and locally bounded intensity rate coupled or not with differential equations. This extension significantly improves the computational performance of JumpProcesses.jl when simulating such processes, enabling it to become one of the few readily available, fast, general-purpose libraries for simulating evolutionary point processes.

We study a sequence of $N-$particle mean-field systems, each driven by $N$ simple point processes $Z^{N,i}$ in a random environment. Each $Z^{N,i}$ has the same intensity $(f(X^N_{t-}))_t$ and at every jump time of $Z^{N,i},$ the process $X^N$ does a jump of height $U_i/\sqrt{N}$ where the $U_i$ are disordered centered random variables attached to each particle. We prove the convergence in distribution of $X^N$ to some limit process $\bar X$ that is solution to an SDE with a random environment given by a Gaussian variable, with a convergence speed for the finite-dimensional distributions. This Gaussian variable is created by a CLT as the limit of the patial sums of the $U_i.$ To prove this result, we use a coupling for the classical CLT relying on the result of [Koml\'os, Major and Tusn\'ady (1976)], that allows to compare the conditional distributions of $X^N$ and $\bar X$ given the random environment, with the same Markovian technics as the ones used in [Erny, L\"ocherbach and Loukianova (2022)]., Nous étudions une suite de systèmes de $N$ particules en champ-moyen, où chaque système est dirigé par $N$ processus ponctuels simples $Z^{N,i}$ dans un environnement aléatoire. Chaque Z^{N,i} a la même intensité $f(X^N_{t-})_t$ et à chaque instant de saut de $Z^{N,i,},$ le processus $X^N$ fait un saut d'amplitude $U_i/\sqrt{N}$ où $U_i$ sont des variables aléatoires de désordre centrées attachées à chaque particule. Nous prouvons la convergence en loi de $X^N$ vers un processus limite $\bar X$ qui est solution d'une EDS dans un environnement aléatoire donné par une variable gaussienne, avec une vitesse de convergence pour les lois fini-dimensionelles. Cette variable gaussienne est créé par un TCL comme la limite des sommes partielles des $U_i.$ Pour montrer ce résultat, nous utilisons un couplage pour le TCL classique qui repose sur le résultat de [Koml\'os, Major et Tusn\'ady (1976)], qui permet de comparer les lois conditionnelles de $X^N$ et de $\bar X$ étant donné les variables d'environnement, avec les mêmes techniques markoviennes que celles utilisées dans [Erny, L\"ocherbach et Loukianova (2022)].

We define a general class of random systems of horizontal and vertical weighted broken lines on the quarter plane whose distribution are proved to be translation invariant. This invariance stems from a reversibility property of the model. This class of systems generalizes several classical processes of the same kind, such as Hammersley's broken line processes involved in Last Passage Percolation theory or such as the six-vertex model for some special sets of parameters. The novelty comes here from the introduction of a weight associated with each line. The lines are initially generated by spatially homogeneous weighted Poisson Point Process and their evolution (turn, split, crossing) are ruled by a Markovian dynamics which preserves Kirchhoff's node law for the line weights at each intersection. Among others, we derive some new explicit invariant measures for some bullet models as well as new reversible properties for some six-vertex models with an external electromagnetic field., Comment: 13 figures

We present general existence and uniqueness results for marked models with pair interactions, exemplified through Gibbs point processes on path space. More precisely, we study a class of infinite-dimensional diffusions under Gibbsian interactions, in the context of marked point configurations: the starting points belong to $\mathbb{R}^d$, and the marks are the paths of Langevin diffusions. We use the entropy method to prove existence of an infinite-volume Gibbs point process and use cluster expansion tools to provide an explicit activity domain in which uniqueness holds., Comment: Added a missing term when proving the form of the interaction range. 34 pages, 6 figures

We present a general statistical methodology for analysing a Laguerre tessellation data set viewed as a realization of a marked point process model. In the first step, for the points, we use a nested sequence of multiscale processes which constitute a flexible parametric class of pairwise interaction point process models. In the second step, for the marks/radii conditioned on the points, we consider various exponential family models where the canonical sufficient statistic is based on tessellation characteristics. For each step, parameter estimation based on maximum pseudolikelihood methods is tractable. For model selection, we consider maximized log pseudolikelihood functions for models of the radii conditioned on the points. Model checking is performed using global envelopes and corresponding tests in both steps and moreover by comparing observed and simulated tessellation characteristics in the second step. We apply our methodology for a 3D Laguerre tessellation data set representing the microstructure of a polycrystalline metallic material, where simulations under a fitted model may substitute expensive laboratory experiments.

Analysis of images of sets of fibers such as myelin sheaths or skeletal muscles must account for both the spatial distribution of fibers and differences in fiber shape. This necessitates a combination of point process and shape analysis methodology. In this paper, we develop a K-function for shape-valued point processes by embedding shapes as currents, thus equipping the point process domain with metric structure inherited from a reproducing kernel Hilbert space. We extend Ripley's K-function which measures deviations from spatial homogeneity of point processes to fiber data. The paper provides a theoretical account of the statistical foundation of the K-function and its extension to fiber data, and we test the developed K-function on simulated as well as real data sets. This includes a fiber data set consisting of myelin sheaths, visualizing the spatial and fiber shape behavior of myelin configurations at different debts., Comment: 12 pages, 2 figures

We prove that, in the coupon collector's problem, the point processes given by the times of $r$-th arrivals for coupons of each type, centered and normalized in a proper way, converge toward a non-homogeneous Poisson point process. This result is then used to derive some generalizations and infinite-dimensional extensions of classical limit theorems on the topic., Comment: 11 pages

We prove limit theorems for functionals of a Poisson point process using the Malliavin calculus on the Poisson space. The target distribution is conditionally either a Gaussian vector or a Poisson random variable. The convergence is stable and our conditions are expressed in terms of the Malliavin operators. For conditionally Gaussian limits, we also obtain quantitative bounds, given for the Monge-Kantorovich transport distance in the univariate case; and for another probabilistic variational distance in higher dimension. Our work generalizes several limit theorems on the Poisson space, including the seminal works by Peccati, Sol��, Taqqu & Utzet for Gaussian approximations; and by Peccati for Poisson approximations; as well as the recently established fourth-moment theorem on the Poisson space of D��bler & Peccati. We give an application to stochastic processes., 30 pages, comments are welcome

We study weak convergence of a sequence of point processes to a scale-invariant simple point process. For a deterministic sequence $(z_n)_{n\in\mathbb{N}}$ of positive real numbers increasing to infinity as $n \to \infty$ and a sequence $(X_k)_{k\in\mathbb{N}}$ of independent non-negative integer-valued random variables, we consider the sequence of point processes \begin{equation*} \nu_n=\sum_{k=1}^\infty X_k \delta_{z_k/z_n}, \quad n\in \mathbb{N}, \end{equation*} and prove that, under some general conditions, it converges vaguely in distribution to a scale-invariant Poisson process $\eta_c$ on $(0,\infty)$ with the intensity measure having the density $ct^{-1}$, $t\in(0,\infty)$. An important motivating example from probabilistic number theory relies on choosing $X_k \sim {\rm Geom}(1-1/p_k)$ and $z_k=\log p_k$, $k\in \mathbb{N}$, where $(p_k)_{k \in \mathbb{N}}$ is an enumeration of the primes in increasing order. We derive a general result on convergence of the integrals $\int_0^1 t \nu_n(dt)$ to the integral $\int_0^1 t \eta_c(dt)$, the latter having a generalized Dickman distribution, thus providing a new way of proving Dickman convergence results. We extend our results to the multivariate setting and provide sufficient conditions for vague convergence in distribution for a broad class of sequences of point processes obtained by mapping the points from $(0,\infty)$ to $\mathbb{R}^d$ via multiplication by i.i.d. random vectors. In addition, we introduce a new class of multivariate Dickman distributions which naturally extends the univariate setting., Comment: Final version, to appear in Electronic Journal of Probability

We present an approximate Bayesian inference approach for estimating the intensity of an inhomogeneous Poisson process, where the intensity function is modelled using a Gaussian process (GP) prior via a sigmoid link function. Augmenting the model using a latent marked Poisson process and P\'olya--Gamma random variables we obtain a representation of the likelihood which is conjugate to the GP prior. We estimate the posterior using a variational free--form mean field optimisation together with the framework of sparse GPs. Furthermore, as alternative approximation we suggest a sparse Laplace's method for the posterior, for which an efficient expectation--maximisation algorithm is derived to find the posterior's mode. Both algorithms compare well against exact inference obtained by a Markov Chain Monte Carlo sampler and standard variational Gauss approach solving the same model, while being one order of magnitude faster. Furthermore, the performance and speed of our method is competitive with that of another recently proposed Poisson process model based on a quadratic link function, while not being limited to GPs with squared exponential kernels and rectangular domains., Comment: 34 pages; 6 figures

We test three common information criteria (IC) for selecting the order of a Hawkes process with an intensity kernel that can be expressed as a mixture of exponential terms. These processes find application in high-frequency financial data modelling. The information criteria are Akaike's information criterion (AIC), the Bayesian information criterion (BIC) and the Hannan-Quinn criterion (HQ). Since we work with simulated data, we are able to measure the performance of model selection by the success rate of the IC in selecting the model that was used to generate the data. In particular, we are interested in the relation between correct model selection and underlying sample size. The analysis includes realistic sample sizes and parameter sets from recent literature where parameters were estimated using empirical financial intra-day data. We compare our results to theoretical predictions and similar empirical findings on the asymptotic distribution of model selection for consistent and inconsistent IC., Comment: 28 pages, 1 figure, submitted to a special issue of Quantitative Finance. Some typos in the previous version have been now corrected

Given a collection $\mathcal L$ of $n$ points on a sphere $\mathbf{S}^2_n$ of surface area $n$, a fair allocation is a partition of the sphere into $n$ parts each of area $1$, and each associated with a distinct point of $\mathcal L$. We show that if the $n$ points are chosen uniformly at random and the partition is defined by considering the gravitational field defined by the $n$ points, then the expected distance between a point on the sphere and the associated point of $\mathcal L$ is $O(\sqrt{\log n})$. We use our result to define a matching between two collections of $n$ independent and uniform points on the sphere, and prove that the expected distance between a pair of matched points is $O(\sqrt{\log n})$, which is optimal by a result of Ajtai, Koml\'os, and Tusn\'ady., Comment: 26 pages, 5 figures

For general thinning procedures, its inverse operation, the condensing, is studied and a link to integration-by-parts formulas is established. This extends the recent results on that link for independent thinnings of point processes to general thinnings of finite point processes. In particular, the classical integration-by-parts formulas appear as the example of independent thinnings. Moreover, the representation of the splitting kernel of finite point processes in terms of its reduced Palm kernels is extended to general thinnings. This link is studied in the context of discrete random variables and yields analogue characterizations of their distributions. Results on independent thinnings are complemented by a discrete stick breaking characterization of distributions., Comment: 30p

We consider the Wiener process with drift $$ dX_t=\mu dt +\sigma d W_t $$ with initial value problem $X_0=x_0$, where $x_0 \in R$, $ \mu \in R$ and $\sigma > 0$ are parameters. By use values $(z_k)_{k \in N}$ of corresponding trajectories at a fixed positive moment $t$, the infinite-sample consistent estimates of each unknown parameter of the Wiener process with drift are constructed under assumption that all another parameters are known. Further, we propose a certain approach for estimation of unknown parameters $x_0,\mu,\sigma$ of the Wiener process with drift by use the values $(z^{(1)}_k)_{k \in N}$ and $(z^{(2)}_k)_{k \in N}$ being the results of observations on the $2k$-th and $2k+1$-th trajectories of the Wiener process with drift at moments $t_1$ and $t_2$ , respectively., Comment: 12 pages, 1 figures, 4 tables. arXiv admin note: substantial text overlap with arXiv:1608.04507

We decompose the genealogy of a general superprocess with spatially dependent branching mechanism with respect to the last individual alive (Williams decomposition). This is a generalization of the main result of Delmas and H\'{e}nard [Electron. J. Probab.,18,1-43,2013] where only superprocesses with spatially dependent quadratic branching mechanism were considered. As an application of the Williams decomposition, we prove that, for some superprocesses, the normalized total mass will converge to a point mass at its extinction time. This generalizes a result of Tribe [Ann. Probab.,20,286-311,1992] in the sense that our branching mechanism is more general.

Consider a random trigonometric polynomial $X_n: \mathbb R \to \mathbb R$ of the form $$ X_n(t) = \sum_{k=1}^n \left( \xi_k \sin (kt) + \eta_k \cos (kt)\right), $$ where $(\xi_1,\eta_1),(\xi_2,\eta_2),\ldots$ are independent identically distributed bivariate real random vectors with zero mean and unit covariance matrix. Let $(s_n)_{n\in\mathbb N}$ be any sequence of real numbers. We prove that as $n\to\infty$, the number of real zeros of $X_n$ in the interval $[s_n+a/n, s_n+ b/n]$ converges in distribution to the number of zeros in the interval $[a,b]$ of a stationary, zero-mean Gaussian process with correlation function $(\sin t)/t$. We also establish similar local universality results for the centered random vectors $(\xi_k,\eta_k)$ having an arbitrary covariance matrix or belonging to the domain of attraction of a two-dimensional $\alpha$-stable law., Comment: 20 pages, extended version. New results (including the stable case) were added

We develop a canonical framework for the study of the problem of registration of multiple point processes subjected to warping, known as the problem of separation of amplitude and phase variation. The amplitude variation of a real random function $\{Y(x):x\in[0,1]\}$ corresponds to its random oscillations in the $y$-axis, typically encapsulated by its (co)variation around a mean level. In contrast, its phase variation refers to fluctuations in the $x$-axis, often caused by random time changes. We formalise similar notions for a point process, and nonparametrically separate them based on realisations of i.i.d. copies $\{\Pi_i\}$ of the phase-varying point process. A key element in our approach is to demonstrate that when the classical phase variation assumptions of Functional Data Analysis (FDA) are applied to the point process case, they become equivalent to conditions interpretable through the prism of the theory of optimal transportation of measure. We demonstrate that these induce a natural Wasserstein geometry tailored to the warping problem, including a formal notion of bias expressing over-registration. Within this framework, we construct nonparametric estimators that tend to avoid over-registration in finite samples. We show that they consistently estimate the warp maps, consistently estimate the structural mean, and consistently register the warped point processes, even in a sparse sampling regime. We also establish convergence rates, and derive $\sqrt{n}$-consistency and a central limit theorem in the Cox process case under dense sampling, showing rate optimality of our structural mean estimator in that case., Comment: Published at http://dx.doi.org/10.1214/15-AOS1387 in the Annals of Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Hawkes (1971) introduced a powerful multivariate point process model of mutually exciting processes to explain causal structure in data. In this paper it is shown that the Granger causality structure of such processes is fully encoded in the corresponding link functions of the model. A new nonparametric estimator of the link functions based on a time-discretized version of the point process is introduced by using an infinite order autoregression. Consistency of the new estimator is derived. The estimator is applied to simulated data and to neural spike train data from the spinal dorsal horn of a rat., Comment: 20 pages, 4 figures

We study the long-range dependence (LRD) of the increments of the fractional Poisson process (FPP), the fractional negative binomial process (FNBP) and the increments of the FNBP. We first point out an error in the proof of Theorem 1 of Biard and Saussereau (2014) and prove that the increments of the FPP has indeed the short-range dependence (SRD) property, when the fractional index $\beta$ satisfies $0, Comment: 17 pages

We study the mutual regularity properties of Palm measures of point processes, and establish that a key determining factor for these properties is the rigidity-tolerance behaviour of the point process in question (for those processes that exhibit such behaviour). Thereby, we extend the results of Osada-Shirai, Bufetov and Olshanski to new ensembles, particularly those that are devoid of any determinantal structure. These include the zeroes of the standard planar Gaussian analytic function and several others., Comment: To appear in the Electronic Communications in Probability

The notion of stability can be generalised to point processes by defining the scaling operation in a randomised way: scaling a configuration by $t$ corresponds to letting such a configuration evolve according to a Markov branching particle system for -$\log t$ time. We prove that these are the only stochastic operations satisfying basic associativity and distributivity properties and we thus introduce the notion of branching-stable point processes. We characterise stable distributions with respect to local branching as thinning-stable point processes with multiplicities given by the quasi-stationary (or Yaglom) distribution of the branching process under consideration. Finally we extend branching-stability to random variables with the help of continuous branching (CB) processes, and we show that, at least in some frameworks, $\mathcal{F}$-stable integer random variables are exactly Cox (doubly stochastic Poisson) random variables driven by corresponding CB-stable continuous random variables., Comment: 31 pages. To appear in Electronic Journal of Probability

The variational principle for Gibbs point processes with general finite range interaction is proved. Namely, the Gibbs point processes are identified as the minimizers of the free excess energy equals to the sum of the specific entropy and the mean energy. The interaction is very general and includes superstable pairwise potential, finite or infinite multibody potential, geometrical interaction, hardcore interaction. The only restrictive assumption involves the finite range property.

We present some correlated fractional counting processes on a finite-time interval. This will be done by considering a slight generalization of the processes in Borges et al. (2012). The main case concerns a class of space-time fractional Poisson processes and, when the correlation parameter is equal to 0, the univariate distributions coincide with those of the space-time fractional Poisson process in Orsingher and Polito (2012). On the one hand, when we consider the time fractional Poisson process, the multivariate finite-dimensional distributions are different from those presented for the renewal process in Politi et al. (2011). We also consider a case concerning a class of fractional negative binomial processes.

In this work we consider a one-dimensional Brownian motion with constant drift moving among a Poissonian cloud of obstacles. Our main result proves convergence of the law of processes conditional on survival up to time $t$ as $t$ converges to infinity in the critical case where the drift coincides with the intensity of the Poisson process. This complements a previous result of T. Povel, who considered the same question in the case where the drift is strictly smaller than the intensity. We also show that the end point of the process conditioned on survival up to time $t$ rescaled by $\sqrt{t} $ converges in distribution to a non-trivial random variable, as $t$ tends to infinity, which is in fact invariant with respect to the drift $h>0$. We thus prove that it is sub-ballistic and estimate the speed of escape. The latter is in a sharp contrast with discrete models of dimension larger or equal to $2$ when the behaviour at criticality is ballistic, see [7], and even to many one dimensional models which exhibit ballistic behaviour at criticality, see [8].

Shot-noise processes are used in applied probability to model a variety of physical systems in, for example, teletraffic theory, insurance and risk theory, and in the engineering sciences. In this paper we prove a large deviation principle for the sample-paths of a general class of multidimensional state-dependent Poisson shot-noise processes. The result covers previously known large deviation results for one-dimensional state-independent shot-noise processes with light tails. We use the weak convergence approach to large deviations, which reduces the proof to establishing the appropriate convergence of certain controlled versions of the original processes together with relevant results on existence and uniqueness.

The purpose of this article is to develop a theory behind the occurrence of “path-integral” kernels in the study of extended determinantal point processes and non-intersecting line ensembles. Our first result shows how determinants involving such kernels arise naturally in studying ratios of partition functions and expectations of multiplicative functionals for ensembles of non-intersecting paths on weighted graphs. Our second result shows how Fredholm determinants with extended kernels (as arise in the study of extended determinantal point processes such as the Airy[subscript 2] process) are equal to Fredholm determinants with path-integral kernels. We also show how the second result applies to a number of examples including the stationary (GUE) Dyson Brownian motion, the Airy[subscript 2] process, the Pearcey process, the Airy[subscript 1] and Airy[subscript 2→1] processes, and Markov processes on partitions related to the zz-measures., National Science Foundation (U.S.) (Grant DMS- 1056390), National Science Foundation (U.S.) (Grant DMS-1208998), Clay Mathematics Institute (Clay Research Fellowship), Microsoft Research (Schramm Memorial Fellowship)

International audience; We consider the population model associated to continuous state branching processes and we are interested in the so-called Eve property that asserts the existence of an ancestor with an overwhelming progeny at large times, and more generally, in the possible behaviours of the frequencies among the population at large times. In this paper, we classify all the possible behaviours according to the branching mechanism of the continuous state branching process.

In this paper, we propose a stochastic process, which is a Cox-Ingersoll-Ross process with Hawkes jumps. It can be seen as a generalization of the classical Cox-Ingersoll-Ross process and the classical Hawkes process with exponential exciting function. Our model is a special case of the affine point processes. Laplace transforms and limit theorems have been obtained, including law of large numbers, central limit theorems and large deviations., Comment: 14 pages

In a recent development in the literature, a new temporal rainfall model, based on the Bartlett-Lewis clustering mechanism and intended for sub-hourly application, was introduced. That model replaced the rectangular rain cells of the original model with finite Poisson processes of instantaneous pulses, allowing greater variability in rainfall intensity over short intervals. In the present paper, the basic instantaneous pulse model is first extended to allow for randomly varying storm types. A systematic comparison of a number of key model variants, fitted to 5-minute rainfall data from Germany, then generates further new insights into the models, leading to the development of an additional model extension, which introduces dependence between rainfall intensity and duration in a simple way. The new model retains the original rectangular cells, previously assumed inappropriate for fine-scale data, obviating the need for the computationally more intensive instantaneous pulse model., Comment: To appear in Hydrological Sciences Journal. 36 pages, 6 tables, 10 figures

We adapt Johnson's sufficiency postulate, Carnap's prediction invariance postulate and B��ge's learn-merge invariance to the context of Papangelou processes and discuss equivalence of their generalizations, in particular their weak and strong generalizations. This discussion identifies a condition which occurs in the construction of Papangelou processes. In particular, we show that these generalizations characterize classes of Poisson and P��lya point processes.

Consider an urn model where at each step one of $q$ colors is sampled according to some probability distribution and a ball of that color is placed in an urn. The distribution of assigning balls to urns may depend on the color of the ball. Collisions occur when a ball is placed in an urn which already contains a ball of different color. Equivalently, this can be viewed as sequentially coloring a complete $q$-partite graph wherein a collision corresponds to the appearance of a monochromatic edge. Using a Poisson embedding technique, the limiting distribution of the first collision time is determined and the possible limits are explicitly described. Joint distribution of successive collision times and multi-fold collision times are also derived. The results can be used to obtain the limiting distributions of running times in various birthday problem based algorithms for solving the discrete logarithm problem, generalizing previous results which only consider expected running times. Asymptotic distributions of the time of appearance of a monochromatic edge are also obtained for other graphs., Comment: Minor revision. 35 pages, 2 figures. To appear in Annals of Applied Probability

Let X be a translation invariant point process on the complex plane and let D be a bounded open set whose boundary has zero Lebesgue measure. We ask what does the point configuration obtained by taking the points of X outside D tell us about the point configuration inside D? We show that for the Ginibre ensemble, it determines the number of points in D. For the translation-invariant zero process of a planar Gaussian Analytic Function, we show that it determines the number as well as the centre of mass of the points in D. Further, in both models we prove that the outside says "nothing more" about the inside, in the sense that the conditional distribution of the inside points, given the outside, is mutually absolutely continuous with respect to the Lebesgue measure on its supporting submanifold., Comment: To appear in Duke Mathematical Journal

Let $B_s$ be a $d$-dimensional Brownian motion and $\omega(dx)$ be an independent Poisson field on $\mathbb{R}^d$. The almost sure asymptotics for the logarithmic moment generating function [\log\math bb{E}_0\exp\biggl{\pm\theta\int_0^t\bar{V}(B_s) ds\biggr}\qquad (t\to\infty)] are investigated in connection with the renormalized Poisson potential of the form [\bar{V}(x)=\int_{\mathbb{R}^d}{\frac{1}{|y-x|^p}}[\omega(dy)-dy],\qquad x\in\mathbb{R}^d.] The investigation is motivated by some practical problems arising from the models of Brownian motion in random media and from the parabolic Anderson models., Comment: Published in at http://dx.doi.org/10.1214/11-AOP655 the Annals of Probability (http://www.imstat.org/aop/) by the Institute of Mathematical Statistics (http://www.imstat.org)

We consider non-colliding Brownian bridges starting from two points and returning to the same position. These positions are chosen such that, in the limit of large number of bridges, the two families of bridges just touch each other forming a tacnode. We obtain the limiting process at the tacnode, the "asymmetric tacnode process". It is a determinantal point process with correlation kernel given by two parameters: (1) the curvature's ratio \lambda>0 of the limit shapes of the two families of bridges, (2) a parameter \sigma controlling the interaction on the fluctuation scale. This generalizes the result for the symmetric tacnode process (\lambda=1 case)., Comment: 21 pages, 1 figure, LaTeX; Includes a further representation of the kernel

Hawkes process is a class of simple point processes that is self-exciting and has clustering effect. The intensity of this point process depends on its entire past history. It has wide applications in finance, neuroscience and many other fields. In this paper, we study the large deviations for nonlinear Hawkes processes. The large deviations for linear Hawkes processes has been studied by Bordenave and Torrisi. In this paper, we prove first a large deviation principle for a special class of nonlinear Hawkes processes, that is, a Markovian Hawkes process with nonlinear rate and exponential exciting function, and then generalize it to get the result for sum of exponentials exciting functions. We then provide an alternative proof for the large deviation principle for a linear Hawkes process. Finally, we use an approximation approach to prove the large deviation principle for a special class of nonlinear Hawkes processes with general exciting functions., Comment: Published in at http://dx.doi.org/10.1214/14-AAP1003 the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org)

We consider a queue where the server is the Euclidean space, and the customers are random closed sets (RACS) of the Euclidean space. These RACS arrive according to a Poisson rain and each of them has a random service time (in the case of hail falling on the Euclidean plane, this is the height of the hailstone, whereas the RACS is its footprint). The Euclidean space serves customers at speed 1. The service discipline is a hard exclusion rule: no two intersecting RACS can be served simultaneously and service is in the First In First Out order: only the hailstones in contact with the ground melt at speed 1, whereas the other ones are queued; a tagged RACS waits until all RACS arrived before it and intersecting it have fully melted before starting its own melting. We give the evolution equations for this queue. We prove that it is stable for a sufficiently small arrival intensity, provided the typical diameter of the RACS and the typical service time have finite exponential moments. We also discuss the percolation properties of the stationary regime of the RACS in the queue., Comment: 26 pages

For a Borel set $A$ and a homogeneous Poisson point process $\eta$ in $\R^d$ of intensity $\lambda >0$, define the Poisson--Voronoi approximation $ A_\eta$ of $A$ as a union of all Voronoi cells with nuclei from $\eta$ lying in $A$. If $A$ has a finite volume and perimeter we find an exact asymptotic of $\E\Vol(A\Delta A_\eta)$ as $\lambda\to\infty$ where $\Vol$ is the Lebesgue measure. Estimates for all moments of $\Vol(A_\eta)$ and $\Vol(A\Delta A_\eta)$ together with their asymptotics for large $\lambda$ are obtained as well., Comment: 19 pages, minor revisions

We consider an ensemble of $N$ discrete nonintersecting paths starting from equidistant points and ending at consecutive integers. Our first result is an explicit formula for the correlation kernel that allows us to analyze the process as $N\to \infty$. In that limit we obtain a new general class of kernels describing the local correlations close to the equidistant starting points. As the distance between the starting points goes to infinity, the correlation kernel converges to that of a single random walker. As the distance to the starting line increases, however, the local correlations converge to the sine kernel. Thus, this class interpolates between the sine kernel and an ensemble of independent particles. We also compute the scaled simultaneous limit, with both the distance between particles and the distance to the starting line going to infinity, and obtain a process with number variance saturation, previously studied by Johansson., Comment: 34 pages, 9 figures; reference added, Theorem 2.1 extended, typos corrected

We consider the beta-Laguerre ensemble, a family of distributions generalizing the joint eigenvalue distribution of the Wishart random matrices. We show that the bulk scaling limit of these ensembles exists for all beta>0 for a general family of parameters and it is the same as the bulk scaling limit of the corresponding beta-Hermite ensemble., 36 pages

We study a family of continuous time Markov jump processes on strict partitions (partitions with distinct parts) preserving the distributions introduced by Borodin (1997) in connection with projective representations of the infinite symmetric group. The one-dimensional distributions of the processes (i.e., the Borodin's measures) have determinantal structure. We express the dynamical correlation functions of the processes in terms of certain Pfaffians and give explicit formulas for both the static and dynamical correlation kernels using the Gauss hypergeometric function. Moreover, we are able to express our correlation kernels (both static and dynamical) through those of the z-measures on partitions obtained previously by Borodin and Olshanski in a series of papers. The results about the fixed time case were announced in the author's note arXiv:1002.2714. A part of the present paper contains proofs of those results., Comment: AMS-LaTeX, 59 pages. v2: Added new results about connections with the z-measures and orthogonal spectral projection operators. Investigated asymptotic behaviour of the dynamical Pfaffian correlation kernel. Removed double contour integral expressions for correlation kernels to shorten the text

Suppose that red and blue points occur as independent Poisson processes of equal intensity in R^d, and that the red points are matched to the blue points via straight edges in a translation-invariant way. We address several closely related properties of such matchings. We prove that there exist matchings that locally minimize total edge length in d=1 and d>=3, but not in the strip R x [0,1]. We prove that there exist matchings in which every bounded set intersects only finitely many edges in d>=2, but not in d=1 or in the strip. It is unknown whether there exists a matching with no crossings in d=2, but we prove positive answers to various relaxations of this question. Several open problems are presented., 21 pages

If we follow an asexually reproducing population through time, then the amount of time that has passed since the most recent common ancestor (MRCA) of all current individuals lived will change as time progresses. The resulting "MRCA age" process has been studied previously when the population has a constant large size and evolves via the diffusion limit of standard Wright--Fisher dynamics. For any population model, the sample paths of the MRCA age process are made up of periods of linear upward drift with slope +1 punctuated by downward jumps. We build other Markov processes that have such paths from Poisson point processes on $\mathbb{R}_{++}\times\mathbb{R}_{++}$ with intensity measures of the form $\lambda\otimes\mu$ where $\lambda$ is Lebesgue measure, and $\mu$ (the "family lifetime measure") is an arbitrary, absolutely continuous measure satisfying $\mu((0,\infty))=\infty$ and $\mu((x,\infty))0$. Special cases of this construction describe the time evolution of the MRCA age in $(1+\beta)$-stable continuous state branching processes conditioned on nonextinction--a particular case of which, $\beta=1$, is Feller's continuous state branching process conditioned on nonextinction. As well as the continuous time process, we also consider the discrete time Markov chain that records the value of the continuous process just before and after its successive jumps. We find transition probabilities for both the continuous and discrete time processes, determine when these processes are transient and recurrent and compute stationary distributions when they exist., Comment: Published in at http://dx.doi.org/10.1214/09-AAP616 the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org)

As well as arising naturally in the study of non-intersecting random paths, random spanning trees, and eigenvalues of random matrices, determinantal point processes (sometimes also called fermionic point processes) are relatively easy to simulate and provide a quite broad class of models that exhibit repulsion between points. The fundamental ingredient used to construct a determinantal point process is a kernel giving the pairwise interactions between points: the joint distribution of any number of points then has a simple expression in terms of determinants of certain matrices defined from this kernel. In this paper we initiate the study of an analogous class of point processes that are defined in terms of a kernel giving the interaction between $2M$ points for some integer $M$. The role of matrices is now played by $2M$-dimensional "hypercubic" arrays, and the determinant is replaced by a suitable generalization of it to such arrays -- Cayley's first hyperdeterminant. We show that some of the desirable features of determinantal point processes continue to be exhibited by this generalization., Comment: 12 pages

We study properties of a sequence $��$ obtained by a randomselection of integers $n$, where $n\in��$ with probability $\varpi_{n}$, independently of the other choices. We distinguish two cases : if $\limsup_{n\to\infty}n\varpi_{n}

We study interacting systems of linear Brownian motions whose drift vector at every time point is determined by the relative ranks of the coordinate processes at that time. Our main objective has been to study the long range behavior of the spacings between the Brownian motions arranged in increasing order. For finitely many Brownian motions interacting in this manner, we characterize drifts for which the family of laws of the vector of spacings is tight, and show its convergence to a unique stationary joint distribution given by independent exponential distributions with varying means. We also study one particular countably infinite system, where only the minimum Brownian particle gets a constant upward drift, and prove that independent and identically distributed exponential spacings remain stationary under the dynamics of such a process. Some related conjectures in this direction have also been discussed., Comment: 25 pages; to appear in the Annals of Applied Probability

Consider a time-varying collection of n points on the positive real axis, modeled as exponentials of n Brownian motions whose drift vector at every time point is determined by the relative ranks of the coordinate processes at that time. If at each time point we divide the points by their sum, under suitable assumptions the rescaled point process converges to a stationary distribution (depending on n and the vector of drifts) as time goes to infinity. This stationary distribution can be exactly computed using a recent result of Pal and Pitman. The model and the rescaled point process are both central objects of study in models of equity markets introduced by Banner, Fernholz, and Karatzas. In this paper, we look at the behavior of this point process under the stationary measure as $n$ tends to infinity. Under a certain `continuity at the edge' condition on the drifts, we show that one of the following must happen: either (i) all points converge to zero, or (ii) the maximum goes to one and the rest go to zero, or (iii) the processes converge in law to a non-trivial Poisson-Dirichlet distribution. The proof employs, among other things, techniques from Talagrand's analysis of the low temperature phase of Derrida's Random Energy Model of spin glasses. The main result establishes a universality property for the BFK models and aids in explicit asymptotic computations using known results about the Poisson-Dirichlet law., Comment: 30 pages, 1 figures, to appear in PTRF

Consider an infinite random matrix $H=(h_{ij})_{0, Comment: 27 pages, 3 figures. This is a new version that incorporates the changes detailed in the Erratum that was published in the same journal

We introduce a new class of spatial-temporal point processes based on Voronoi tessellations. At each step of such a process, a point is chosen at random according to a distribution determined by the associated Voronoi cells. The point is then removed, and a new random point is added to the configuration. The dynamics are simple and intuitive and could be applied to modeling natural phenomena. We prove ergodicity of these processes under wide conditions., Comment: 17 pages, 4 figures

A random dense countable set is characterized (in distribution) by independence and stationarity. Two examples are `Brownian local minima' and `unordered infinite sample'. They are identically distributed. A framework for such concepts, proposed here, includes a wide class of random equivalence classes., Comment: 40 pages. Supersedes math.PR/0508414 and math.PR/0511011