Thresholding Bursty Time Series

Models extremes of 'bursty' time series via Continuous Time Random Exceedances (CTRE). See , K. Hees, S. Nayak, P.Straka, 2018.

Travis-CI BuildStatus CRAN_Status_Badge Downloads DL_Total

This R package provides tools to analyse extremes of ‘bursty’ time series. Burstiness is characterized by heavy-tailed inter-arrival times and scale-free event dynamics. The CTRE model captures burstiness by generalizing the Poisson process to a fractional Poisson process, with Mittag-Leffler inter-arrival times. Parameter estimates are read off from stability plots, and goodness of fit is assessed via diagnostic plots; see the Shiny app below.


“Peaks Over Threshold for Bursty Time Series”, Katharina Hees, Smarak Nayak, Peter Straka (2018).

Shiny App

The package comes with two examples of bursty time series: solar flare magnitudes and bitcoin trading volumes. For parameter estimates of the Mittag-Leffler distribution, see the tab “Exceedance Times”. CTRE model assumptions are checked via a QQ plot of the Mittag-Leffler distribution; an empirical copula plot checking for dependence between inter-arrival times and magnitudes; and a plot of the autocorrelation function for the two series (interarrival times and magnitudes). For the standard POT model plots, see the “Exceedances” tab.

Install from GitHub


Run shiny app

You can run the above Shiny app from within RStudio:


Package usage

You can


Reference manual

It appears you don't have a PDF plugin for this browser. You can click here to download the reference manual.


0.1.0 by Peter Straka, 9 months ago

Report a bug at

Browse source code at

Authors: Katharina Hees [aut] , Peter Straka [aut, cre]

Documentation:   PDF Manual  

GPL-3 license

Imports assertthat, graphics, MittagLeffleR, plyr, stats, tea, utils, zoo

Suggests evmix, magrittr, shiny

See at CRAN