tick.hawkes.SimuHawkesExpKernels¶

class tick.hawkes.SimuHawkesExpKernels(adjacency, decays, baseline=None, end_time=None, period_length=None, max_jumps=None, seed=None, verbose=True, force_simulation=False)[source]¶

Hawkes process with exponential kernels simulation

They are defined by the intensity:

$\forall i \in [1 \dots D], \quad \lambda_i(t) = \mu_i + \sum_{j=1}^D \int \phi_{ij}(t - s) dN_j(s)$

where

$D$ is the number of nodes

$\mu_i(t)$ are the baseline intensities

$\phi_{ij}$ are the kernels

$dN_j$ are the processes differentiates

and with an exponential parametrisation of the kernels

$\phi_{ij}(t) = \alpha_{ij} \beta_{ij} \exp (- \beta_{ij} t) 1_{t > 0}$

where $\alpha_{ij}$ is the intensity of the kernel and $\beta_{ij}$ its decay. The matrix of all $\alpha$ is called adjacency matrix.

Parameters

baseline : np.ndarray or list

The baseline of all intensities, also noted $\mu(t)$ . It might be three different types:

np.ndarray, shape=(n_nodes,) : One baseline per node is given. Hence baseline is assumed to be constant, ie. $\mu_i(t) = \mu_i$

np.ndarray, shape=(n_nodes, n_intervals) : n_intervals baselines are given per node. This assumes parameter period_length is also given. In this case baseline is piecewise constant on intervals of size period_length / n_intervals and periodic.

list of tick.base.TimeFunction, shape=(n_nodes,) : One function is given per node, ie. $\mu_i(t)$ is explicitely given.

adjacency : np.ndarray, shape=(n_nodes, n_nodes)

Intensities of exponential kernels, also named $\alpha_{ij}$

decays : float or np.ndarray, shape=(n_nodes, n_nodes)

Decays of exponential kernels, also named $\beta_{ij}$ If a float is given, all decays are equal to this float

end_time : float, default=None

Time until which this point process will be simulated

max_jumps : int, default=None

Simulation will stop if this number of jumps in reached

seed : int, default = None

The seed of the random sampling. If it is None then a random seed (different at each run) will be chosen.

force_simulation : bool, default = False

If force is not set to True, simulation won’t be run if the matrix of the L1 norm of kernels has a spectral radius greater or equal to 1 as it would be unstable

Attributes

timestamps : list of np.ndarray, size=n_nodes

A list of n_nodes timestamps arrays, each array containing the timestamps of all the jumps for this node

simulation_time : float

Time until which this point process has been simulated

n_total_jumps : int

Total number of jumps simulated

tracked_intensity : list[np.ndarray], size=n_nodes

A record of the intensity with which this point process has been simulated. Note: you must call track_intensity before simulation to record it

intensity_tracked_times : np.ndarray

The times at which intensity has been recorded. Note: you must call track_intensity before simulation to record it

intensity_track_step : float

Step with which the intensity has been recorded

__init__(adjacency, decays, baseline=None, end_time=None, period_length=None, max_jumps=None, seed=None, verbose=True, force_simulation=False)[source]¶: Initialize self. See help(type(self)) for accurate signature.