tstab.gpd | R Documentation |
This function computes the maximum likelihood estimate at each provided threshold and plots the estimates (pointwise), along with 95 or else from 1000 independent draws from the posterior distribution under vague independent normal prior on the log-scale and shape. The latter two methods better reflect the asymmetry of the estimates than the Wald confidence intervals.
tstab.gpd(
xdat,
thresh,
method = c("wald", "profile", "post"),
level = 0.95,
plot = TRUE,
...
)
xdat |
a vector of observations |
thresh |
a vector of candidate thresholds at which to compute the estimates. |
method |
string indicating the method for computing confidence or credible intervals.
Must be one of |
level |
confidence level of the intervals. Default to 0.95. |
plot |
logical; should parameter stability plots be displayed? Default to |
... |
additional arguments passed to |
a list with components
threshold
: vector of numerical threshold values.
mle
: matrix of modified scale and shape maximum likelihood estimates.
lower
: matrix of lower bounds for the confidence or credible intervals.
upper
: matrix of lower bounds for the confidence or credible intervals.
method
: method for the confidence or coverage intervals.
plots of the modified scale and shape parameters, with pointwise confidence/credible intervals
and an invisible data frame containing the threshold thresh
and the modified scale and shape parameters.
Leo Belzile
gpd.fitrange
dat <- abs(rnorm(10000))
u <- qnorm(seq(0.9,0.99, by= 0.01))
tstab.gpd(xdat = dat, thresh = u)
## Not run:
tstab.gpd(xdat = dat, thresh = u, method = "profile")
tstab.gpd(xdat = dat, thresh = u, method = "post")
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.