Sequentially performs the entropy difference (ED) test or the multiplier or parametric bootstrap score tests for the GEVr model.
gevrSeqTests(
data,
bootnum = NULL,
method = c("ed", "pbscore", "multscore"),
information = c("expected", "observed"),
allowParallel = FALSE,
numCores = 1
)
Data should be contain n rows, each a GEVr observation.
If method equals 'pbscore' or 'multscore', the number of bootstrap simulations to use.
Which test to run: ED test (ed), multiplier (multscore) or parametric bootstrap (pbscore) score test.
To use expected (default) or observed information in the score tests.
If method equals 'pbscore', should the parametric boostrap procedure be run in parallel or not. Defaults to false.
If allowParallel is true, specify the number of cores to use.
Function returns a dataframe containing the test statistics, estimates, and p-value results of the sequential tests.
Value of r to be tested.
Raw p-values from the individual tests at each value of r.
Transformed p-values according to the ForwardStop stopping rule.
Transformed p-values according to the StrongStop stopping rule.
Returned test statistics of each individual test.
Estimated location parameter for the given r.
Estimated scale parameter for the given r.
Estimated shape parameter for the given r.
GEVr data (in matrix x) should be of the form \(x[i,1] > x[i, 2] > \cdots > x[i, r]\) for each observation \(i = 1, \ldots, n\). See function `pSeqStop' for details on transformed p-values.
# NOT RUN {
x <- rgevr(200, 5, loc = 0.5, scale = 1, shape = 0.25)
gevrSeqTests(x, method = "ed")
# }
Run the code above in your browser using DataLab