TOPICS
Search

Hawkes Process


The term "Hawkes process" encompasses a number of point process notions; while many are similar, some are rather different. There are also different formulations for univariate and multivariate point processes.

In some literature, a univariate Hawkes process is defined to be a self-exciting temporal point process N whose conditional intensity function lambda=lambda(t) is defined to be

 lambda(t)=mu(t)+sum_(i:tau_i<t)nu(t-tau_i),
(1)

where mu(t) is the background rate of the process N, where tau_i are the points in time occurring prior to time t, and where nu is a function which governs the clustering density of N. The function nu is sometimes called the exciting function or the excitation function of N. Similarly, some authors (Mehrdad and Zhu 2014) denote the conditional intensity function lambda by lambda_t and rewrite the addend in (1) as

 sum_(i:tau_i<t)nu(t-tau_i)=int_(-infty)^tnu(t-s)N(ds).
(2)

The processes upon which Hawkes himself made the most progress were univariate self-exciting temporal point processes N whose conditional intensity function lambda depends linearly on the past arrivals (Hawkes 1971). As a result, some authors refer to such processes as Hawkes processes. In general, however, such lambda behavior is typically specified, i.e., processes for which lambda is linear are referred to as linear Hawkes processes and are differentiated from their nonlinear counterparts whose conditional intensity functions lambda are nonlinear (Mehrdad and Zhu 2014).

Still other authors consider two alternative brands of univariate Hawkes processes, one the so-called intensity-based Hawkes process and the other the so-called cluster-based version, which are equivalent though are studied in different contexts (Dassios and Zhao 2013). In this case, the intensity-based process is a temporal point process N_t={tau_i}_(i=1,2,...) on R^+ which has a nonnegative exponentially-decaying F_t-stochastic intensity lambda_t of the form

 lambda_t=a+(lambda_0-a)e^(-deltat)+sum_(0<=tau_i<t)Y_ie^(-delta(t-tau_i))
(3)

for t>=0, where {F_t}_(t>=0) is a history of the process N_t with respect to which {lambda_t}_(t>=0) is adapted, a>=0 is the constant reversion level, lambda_0>0 is the initial intensity at time t=0, delta>0 is the constant rate of exponential decay, {Y_i}_(i=1,2,...) are the sizes of the self-excited jumps viewed as independent random variables distributed according to some distribution function G(y), y>0, and each new mark Y_i is independent of its arrival time tau_i and the preceding history. The mark affects subsequent arrival times, so the full sequences {tau_i} and {Y_i} are generally dependent. This is equivalent to the cluster-based version in which (3) is viewed as a marked Poisson cluster process C={(tau_i,Y_i)}_(i=1,2,...), the only difference being that from the cluster-based perspective:

1. The set I subset {tau_i} subset R^+ consists of elements known as immigrants which are distributed as an inhomogeneous Poisson process with rate

 a+(lambda_0-a)e^(-deltat)
(4)

for t>=0.

2. The set {Y_i:tau_i in I} of marks associated to the immigrants I are independent of the immigrants and are distributed as independent random variables according to some distribution G.

3. Each immigrant tau_m generates a single cluster C_m independent of other clusters where here, each C_m is viewed as a random set subject to a certain branching structure (Dassios and Zhao 2013) which satisfies the property that C= union _mC_m.

In addition to these ambiguities, several authors (e.g., Mehrdad and Zhu 2014), have made generalized versions of the univariate process described in (1) which are still referred to as Hawkes processes. One example includes adapting (1) so that the process has different exciting functions, the result of which is a collection (N^n)_(n in N) of non-explosive simple point processes for which:

1. N^0 is an inhomogeneous Poisson process with intensity gamma_0(t) at time t;

2. For every n in N, N^n is a simple point process with intensity

 lambda_t^n=int_0^tgamma_n(t-s)N^(n-1)(ds)=sum_(tau in N^(n-1),0<tau<t)gamma_n(t-tau);
(5)

3. For every n in {0,1,2,...}, N^(n+1) is an inhomogeneous Poisson process with intensity lambda^(n+1) conditional on N^0,N^1,...,N^n.

In this context, provided that the superposition N=sum_(n=0)^(infty)N^n is non-explosive, meaning that it has finitely many points in every bounded time interval almost surely, N is said to be a univariate Hawkes point process with excitation functions (gamma_n)_(n in {0,1,2,...}) while N^0 is called the immigrant process and N^n the nth generation offspring process (Mehrdad and Zhu 2014). Note that when gamma_0=mu and when gamma_n=nu for any n in N, this extended model reduces to the classical linear model (2).

Due to the vast usage of the term Hawkes process among univariate point processes, one expects there to be room for an equally large number of definitions of multivariate Hawkes processes. Surprisingly, however, the most common use of the term is assigned to a relatively straightforward extension of equations (1) and (2), whereby one says that a multivariate d-dimensional counting process N=(N_1,...,N_d) taking values in N^d is a multivariate Hawkes process whenever the associated conditional intensity function (lambda_1,...,lambda_d) defined by

 P(N_i has a jump in [t,t+dt]|F_t)=lambda_(i,t)dt+o(dt)
(6)

for i=1, ..., d has the form

 lambda_(i,t)=mu_i+int_0^tsum_(j=1)^dphi_(ij)(t-s)N_j(ds)
(7)

(Bacry et al. 2012). Here, P stands for probability, F_t is the sigma-algebra generated by N up to present time t, mu_i in R^+, and phi_(ij):R^+->R^+ for i=1,...,d.

It is worth noting, however, that precisely as with the univariate case, some authors distinguish between different "types" of multivariate Hawkes processes (Liniger 2009) while other authors define completely separate types of multivariate functions to be multivariate Hawkes processes (Carlsson et al. 2007).


See also

Conditional Intensity Function, Multivariate Function, Point Process, Poisson Process, Probability, Self-Exciting Point Process, Stochastic, Temporal Point Process, Univariate Function

Portions of this entry contributed by Christopher Stover

Explore with Wolfram|Alpha

References

Bacry, E.; Delattre, S.; Hoffmann, M.; and Muzy, J. F. "Scaling Limits for Hawkes Processes and Application to Financial Statistics." 3 Feb 2012. https://arxiv.org/abs/1202.0842v1.Carlsson, J.; Foo, M. C.; Lee, H. H.; and Shek, H. "High Frequency Trade Prediction with Bivariate Hawkes Process." 2007. https://users.iems.northwestern.edu/~armbruster/2007msande444/report1b.pdf.Dassios, A. and Zhao, H. "Exact Simulation of Hawkes Process with Exponentially Decaying Intensity." Electron. Commun. Probab. 18, 1-13, 2013.Hawkes, A. G. "Spectra of Some Self-Exciting and Mutually Exciting Point Processes." Biometrika 58, 83-90, 1971.Hawkes, A. G. and Adamopoulos, L. "Cluster Models for Earthquakes: Regional Comparisons." Bull. Int. Statist. Inst. 45, 454-461, 1973.Hawkes, A. G. and Oakes, D. "A Cluster Process Representation of a Self-Exciting Process." J. Appl. Prob. 11, 493-503, 1974.Liniger, T. J. "Multivariate Hawkes Processes." Doctor of Sciences dissertation. Zürich, Switzerland: Swiss Federal Institute of Technology Zurich, 2009. https://doi.org/10.3929/ethz-a-006037599.Mehrdad, B. and Zhu, L. "On the Hawkes Process with Different Exciting Functions." 2014. https://arxiv.org/abs/1403.0994v1.Numberphile. Featuring H. Fry. "The Mathematics of Crime and Terrorism." Feb. 3, 2016. https://www.youtube.com/watch?v=lCjspXB5F4A.

Referenced on Wolfram|Alpha

Hawkes Process

Cite this as:

Weisstein, Eric W., with contributions by Christopher Stover. "Hawkes Process." From MathWorld--A Wolfram Resource. https://mathworld.wolfram.com/HawkesProcess.html

Subject classifications