bootconfr {extremogram}

Confidence bands for the sample return time extremogram

Description

The function estimates confidence bands for the sample return time extremogram using the stationary bootstrap.

Usage

bootconfr(x, R, l, maxlag, uplevel = 1, lowlevel = 0, type, par,
  start = 1, cutoff = 1, alpha = 0.05)

Arguments

x

Univariate time series (a vector).

R

Number of bootstrap replications (an integer).

l

Mean block size for stationary bootstrap or mean of the geometric distribution used to generate resampling blocks (an integer that is not longer than the length of the time series).

maxlag

Number of lags to include in the extremogram (an integer)

uplevel

Quantile of the time series to indicate a upper tail extreme event (a number between 0 and 1, default is 1).

lowlevel

Quantile of the time series to indicate a lower tail extreme event (a number between 0 and 1, default is 0).

type

Extremogram type (see function extremogramr).

par

If par = 1, the bootstrap replication procedure will be parallelized. If par = 0, no parallelization will be used.

start

The lag that the extremogram plots starts at (an integer not greater than maxlag, default is 1).

cutoff

The cutoff of the y-axis on the plot (a number between 0 and 1, default is 1).

alpha

Significance level for the confidence bands (a number between 0 and 1, default is 0.05).

Value

Returns a plot of the confidence bands for the sample return time extremogram.

References

1. Davis, R. A., Mikosch, T., & Cribben, I. (2012). Towards estimating extremal serial dependence via the bootstrapped extremogram. Journal of Econometrics,170(1), 142-152.

2. Davis, R. A., Mikosch, T., & Cribben, I. (2011). Estimating extremal dependence in univariate and multivariate time series via the extremogram.arXiv preprint arXiv:1107.5592.

Examples

# generate a GARCH(1,1) process
omega    = 1
alpha    = 0.1
beta     = 0.6
n        = 1000
uplevel  = 0.95
lowlevel = 0.05
type     = 3
maxlag   = 70
df       = 3
R        = 10
l        = 30
par      = 0
G = extremogram:::garchsim(omega,alpha,beta,n,df)

extremogramr(G, type, maxlag, uplevel, lowlevel, 1, 1)
bootconfr(G, R, l, maxlag, uplevel, lowlevel, type, par, 1, 1, 0.05)