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). https://arxiv.org/abs/1802.05218
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.
library("devtools")
install_github("UNSW-MATH/CTRE")
library(CTRE)
You can run the above Shiny app from within RStudio:
runCTREshiny()
You can
ctre
object from a time
series, a data frame, or two vectors.plot
thin
interarrival
,
time
and
magnitudes
MLestimates
mlqqplot
,
acf
,
empcopula
)Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.