Skip to contents

Package

exuber-package exuber
exuber: Econometric Analysis of Explosive Time Series

Recursive Augmented Dickey Fuller

Estimation & Critical Values

radf()
Recursive Augmented Dickey-Fuller Test
radf_mc_cv() radf_mc_distr()
Monte Carlo Critical Values
radf_wb_cv() radf_wb_distr()
Wild Bootstrap Critical Values
radf_wb_cv2() radf_wb_distr2()
Wild Bootstrap Critical Values
radf_sb_cv() radf_sb_distr()
Panel Sieve Bootstrap Critical Values
radf_crit
Stored Monte Carlo Critical Values

Heteroskedasticity-robust (time-transformed)

Bootstrap-free alternative to radf_wb_cv() under time-varying volatility

radf_tt()
Time-Transformed Test for Explosive Bubbles under Non-stationary Volatility
radf_tt_cv()
Monte Carlo critical values for the time-transformed test (STADF/GSTADF)

Volatility-robust (other routes)

Further tests robust to time-varying innovation variance; see docs/enhancements/volatility-robustness.md

radf_sbz_cv()
SBZ Weighted Least Squares Bubble Test with Union-of-Rejections
radf_kp()
Kernel-Purged Heteroskedasticity-Robust PSY Test
radf_sign()
Sign-Based Bubble Test (sPWY / sPSY)
radf_sign_cv()
Monte Carlo Critical Values for the Sign-Based Test
radf_sign_dm()
Recursively Demeaned Sign-Based Bubble Test (s-bar-PWY / s-bar-PSY)
radf_sign_dm_cv()
Monte Carlo Critical Values for the Recursively Demeaned Sign-Based Test
ssu_test()
Stochastic Unit Root Bubble Test (Kurozumi & Nishi 2025)
radf_svadf()
SV-ADF Asymmetric-Threshold Bubble Dating (Sarkar & Wells 2026)

Dating and root inference

Start/end/recovery dates and confidence intervals on the explosive root; see docs/enhancements/dating-and-root-inference.md

explosive_root()
Estimate the Explosive Autoregressive Root over a Sub-Sample
root_ci()
Confidence Interval and Doubling Time for an Explosive Root
root_ci_datestamp()
Root Confidence Intervals for Every Datestamped Episode
dating_pdc()
Sequential Sample-Splitting Bubble Dating (PDC/KS)
radf_recovery()
Reverse-Regression Dating of Crisis Origination and Market Recovery
radf_recovery_cv()
Monte Carlo Critical Values for Reverse-Regression Recovery Dating
dating_hls()
SSR/BIC Bubble Dating (Harvey, Leybourne & Sollis 2017)
dating_hlw()
Multi-Bubble SSR/BIC Dating (Harvey, Leybourne & Whitehouse 2020)
dating_knp()
Bias-Corrected Single-Bubble Dating (Kejriwal, Nguyen & Perron 2025)

Real-time monitoring

Sequential/real-time bubble detection; see docs/enhancements/monitoring.md

radf_monitor()
Real-Time Monitoring for Explosive Bubbles
monitor_cusum()
CUSUM Real-Time Monitoring for Explosive Bubbles
lbi_test()
Locally Best Invariant Test for a Bubble (Breitung & Diegel 2025)
monitor_lbi()
Sequential LBI Monitoring for an Unknown Bubble Start Date (Breitung & Diegel 2025)

Multivariate / panel bubble tests

Panel and cross-series tests; see docs/enhancements/multivariate.md

radf_common()
Common-Bubble Detection via PCA + PSY
radf_common_cv()
Critical Values for the Common-Bubble (PCA + PSY) Test
cobubble_test()
Test for Co-explosive Behaviour Between Two Series
contagion_reg()
Bubble Contagion Regression (Greenaway-McGrevy & Phillips 2016)

Alternative paradigms

Non-ADF-family approaches; see docs/enhancements/alternative-paradigms.md

quantile_test()
Quantile Unit Root Test for Bubble Detection (Global Test)
monitor_quantile()
QPWY Recursive Quantile Monitoring (Wu, Shi & Wu 2025)

Analysis

summary(<radf_obj>)
Summarizing radf models
diagnostics()
Diagnostics on hypothesis testing
datestamp()
Date-stamping periods of mildly explosive behavior

Tidying

tidy(<radf_obj>) augment(<radf_obj>)
Tidy a radf_obj object
tidy(<radf_cv>) augment(<radf_cv>)
Tidy a radf_cv object
tidy(<radf_distr>)
Tidy a radf_distr object
tidy(<ds_radf>)
Tidy a ds_radf object

Plotting

autoplot(<radf_obj>) autoplot2(<radf_obj>) shade()
Plotting radf models
autoplot(<radf_distr>)
Plotting a radf_distr object
autoplot(<ds_radf>)
Plotting a ds_radf object
scale_exuber_manual() theme_exuber()
Exuber scale and theme functions

Simulation

sim_psy1()
Simulation of a single-bubble process
sim_psy2()
Simulation of a two-bubble process
sim_ps1()
Simulation of a single-bubble process with multiple forms of collapse regime
sim_blan()
Simulation of a Blanchard (1979) / Rotermann-Wilfling (2018) bubble process
sim_evans()
Simulation of an Evans (1991) bubble process
sim_div()
Simulation of dividends

Simulation (additional DGPs)

Data-generating processes for axes exuber’s original sim_*() functions don’t cover – time-varying volatility, non-Gaussian innovations, long memory, and distinct multi-series/branching mechanisms; see docs/enhancements/simulation-dgps.md

sim_innov()
Simulate innovations with heavy-tailed/skewed marginal distributions
sim_vol_garch()
Simulate GARCH(1,1)/TGARCH(1,1) innovations
sim_vol_cir()
Simulate CIR-type stochastic-volatility innovations
sim_vol_sv()
Simulate AR(1) lognormal stochastic-volatility innovations
sim_fi()
Simulate fractionally-integrated (long-memory) innovations
sim_common()
Simulation of a latent common-factor bubble across multiple series
sim_coexplosive()
Simulation of a bivariate co-explosive pair
sim_tree()
Simulation of a stochastic branching-tree bubble
sim_mar()
Simulation of a mixed causal-noncausal AR(1,1) bubble
sim_msbubble()
Simulation of a Markov-switching present-value bubble
sim_falsebubble()
Simulation of a deterministic technology-adoption "false bubble" null

Helpers

psy_minw() psy_ds()
Helper functions in accordance to PSY(2015)
index() `index<-`()
Retrieve/Replace the index
series_names() `series_names<-`()
Retrieve/Replace series names
ps_tb()
Helper function to find tb from the Phillips and Shi (2020)

Extra tidying methods

tidy_join() augment_join()
Tidy into a joint model
tidy_join(<radf_obj>) augment_join(<radf_obj>)
Tidy into a joint model