xtibber: Test-Inversion Bootstrap for Extreme-Value Analysis

View source: R/xtibber.R

xtibberR Documentation

Test-Inversion Bootstrap for Extreme-Value Analysis

Description

Test-inversion bootstrap (TIB) for fevd class objects.

Usage

xtibber(x, type = c("return.level", "parameter"), which.one,
	tib.method = c("interp", "rm"), nuisance = "shape", B,
	test.pars, rsize, block.length = 1, shuffle = NULL,
	replace = TRUE, alpha = 0.05, qcov = NULL,
	qcov.base = NULL, stud = FALSE, step.size, tol = 1e-04,
	max.iter = 1000, keep.iters = TRUE, verbose = FALSE, ...)

Arguments

x

List object of class “fevd”.

type

character string stating whether to calculate TIB intervals for a return level or a parameter as this funciton will only calculate an interval for a single parameter/return level at a time.

which.one

number or character stating which return level or which parameter to find CIs for.

tib.method

character stating whether to estimate the TIB interval by interpolating from a series of pre-determined values of the nuisance parameter or to use the Robbins-Monroe (RM) method. See the help file for tibber from the distillery package for more information.

nuisance

character naming the nuisance parameter.

B, rsize, block.length, shuffle, replace

See the help file for booter from the distillery package for more information on these arguments.

test.pars

numeric vector giving the sequence of nuisance parameter values for the interpolation method, or a numeric vector of length two giving the starting values for the RM method.

alpha

numeric between zero and one giving the desired confidence level.

qcov

numeric matrix with rows the same length as q and columns equal to the number of parameters (+ 1 for the threshold, if a POT model). This gives any covariate values for a nonstationary model. If NULL, and model is non-stationary, only the intercept terms for modeled parameters are used, and if a non-constant threshold, only the first threshold value is used. Not used if model is stationary.

qcov.base

numeric matrix analogous to qcov. When provided, the function returns the difference in return levels between the level for the covariates in qcov and the level for covariates in qcov.base.

stud

logical if TRUE will calculate Studentized intervals (generally not profitable with the TIB method).

step.size

Used with the RM method only. Numeric giving the size of increments to use in the root-finding algorithm.

tol

Used with the RM method only. Numeric stating how close to the desired level of confidence is satisfactory.

max.iter

numeric giving the maximum number of iterations for the root-finding algorithm before giving up.

keep.iters

logical, should all of the values in the root-finding search be kept? Needed if a plot will be made.

verbose

logical, if TRUE will print progress information to the screen.

...

optional arguments to nlminb.

Details

This function provides a wrapper to the tibber function from distillery for “fevd” objects.

Value

See the help file for tibber for more information on the value

Author(s)

Eric Gilleland

References

Gilleland, E. (2020) Bootstrap methods for statistical inference. Part I: Comparative forecast verification for continuous variables. Journal of Atmospheric and Oceanic Technology, 37 (11), 2117 - 2134, doi: 10.1175/JTECH-D-20-0069.1.

Gilleland, E. (2020) Bootstrap methods for statistical inference. Part II: Extreme-value analysis. Journal of Atmospheric and Oceanic Technology, 37 (11), 2135 - 2144, doi: 10.1175/JTECH-D-20-0070.1.

See Also

fevd, tibber, booter

Examples

## Not run: 
data("ftcanmax")
fit <- fevd( Prec, data = ftcanmax )

tbfit <- xtibber( fit, which.one = 100, B = 500,
		 test.pars = seq(-0.01,0.2,,100), verbose = TRUE )

tbfit

plot( tbfit )


## End(Not run)

extRemes documentation built on May 29, 2024, 5:27 a.m.