Giter Site home page Giter Site logo

saadlabyad / aslsd Goto Github PK

View Code? Open in Web Editor NEW
13.0 2.0 2.0 12.86 MB

Parametric estimation of multivariate Hawkes processes with general kernels.

License: BSD 3-Clause "New" or "Revised" License

Python 100.00%
stochastic-gradient-descent sgd hawkes-process monte-carlo

aslsd's Introduction

aslsd

Parametric estimation of multivariate Hawkes processes with general kernels.

This code is under active development, and is part of our paper Gradient-based estimation of linear Hawkes processes with general kernels.

Description

Multivariate Hawkes processes (MHP) are a fundamental class of point processes with self-excitation. When estimating parameters for these processes, a difficulty is that the two main error functionals, the log-likelihood and the least squares error (LSE), as well as the evaluation of their gradients, have a quadratic complexity in the number of observed events. In practice, this prohibits the use of exact gradient-based algorithms for parameter estimation in many settings. Furthermore, MHP models are not designed for non-stationary training data, and they cannot incorporate event information besides their timestamps: we introduce the marked time-dependent linear Hawkes (MTLH) model to overcome these limitations.

We construct an adaptive stratified sampling estimator of the gradient of the LSE of Hawkes models. This results in the ASLSD algorithm, a fast parametric estimation method for MHP and MTLH with general kernels, applicable to large datasets, which compares favourably with existing methods. We evaluate our algorithm on synthetic and real-world data.

This code implements

  • the ASLSD algorithm for the estimation of MHP and MTLH models;
  • exact cluster based simulation of MHP and MTLH models;
  • residual analysis and other evaluation metrics for MHP and MTLH models.

Dependencies

This code requires Python 3.7 or newer, as well as:

Citation

This code is part of our paper Gradient-based estimation of linear Hawkes processes with general kernels.

If you use this code as part of a scientific publication, please acknowledge our paper.

@article{cartea2021gradient,
  title={Gradient-based estimation of linear Hawkes processes with general kernels},
  author={Cartea, {\'A}lvaro and Cohen, Samuel N and Labyad, Saad},
  journal={arXiv preprint arXiv:2111.10637},
  year={2021}
}

aslsd's People

Contributors

saadlabyad avatar stochasticsam avatar

Stargazers

 avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar

Watchers

 avatar  avatar

Forkers

zshwuhan azztt

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    ๐Ÿ–– Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. ๐Ÿ“Š๐Ÿ“ˆ๐ŸŽ‰

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google โค๏ธ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.