Description Usage Arguments Details Value Examples
Evaluate uncertainty in the mixture of distributions of maximums assuming the underlying data generating mechanism is the 2D extremal Poisson process, but without specifying a single threshold. The mixture is across multiple thresholds.
1 2 3 | ## S3 method for class 'full_multi_fit'
fullMaxDistUncert(x, declust_obs, lt_gen, n_mc,
n_boot, n_starts, progress_tf = TRUE)
|
x |
An S3 object of class |
declust_obs |
(numeric vector) The observed data used by
|
lt_gen |
(numeric scalar) Length of each generated series. The units
(seconds, minutes, hours, etc.) should be consistent with the value of
|
n_mc |
(numeric scalar) The number of samples to draw from the distribution of the maximum |
n_boot |
(numeric scalar) The number of bootstrap replicates of the distribution of the maximum to create. |
n_starts |
(numeric scalar) The number of random starts to use in the search for the maximum |
progress_tf |
(logical scalar) Display a progress bar if TRUE, else not. |
The results of fitting many 2D extremal Poisson process are fed into
this function. The declustered data are repeatedly sampled with
replacement, and for each resampled data set the distribution of the
maximum is empirically constructed as described in
fullMaxDist.full_multi_fit
. The bootstrap replicates of the mixture
of distributions of the maximum may be used to quantify uncertainty and
construct intervals.
An S3 object of class full_max_dist_uncert_multi_thresh
,
which is a list of length n_boot
of S3 objects of class
full_max_dist_multi_thresh
.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | ## Not run:
complete_series <- -jp1tap1715wind270$value
declustered_obs <- decluster(complete_series)
full_multi_fit <- fullMultiFit(x = declustered_obs, lt = 100,
n_min = 10, n_max = 50,
weight_scale = 5, n_starts = 20)
full_multi_fit_uncert <- fullMaxDistUncert(x = full_multi_fit,
declust_obs = declustered_obs$declustered_series,
lt_gen = 200,
n_mc = 1000,
n_boot = 200)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.