pareto_tailplot {tailplots} | R Documentation |
Plot the estimated t and the corresponding confidence intervals
Description
This function produces a tail plot for the estimate \hat{t}
over a range of thresholds for a given sample, including confidence intervals computed by one of three methods (unbiased, bootstrap or jackknife). The function also allows a choice between original and log scale.
Usage
pareto_tailplot(
x,
method = c("unbiased", "bootstrap", "jackknife"),
R = 1000,
conf.level = 0.95,
ci.points = 101,
xscale = "b"
)
Arguments
x |
a vector containing the sample data. |
method |
the method used for computing the confidence intervals (options include unbiased variance estimator, jackknife, and bootstrap). |
R |
the number of the bootstrap replicates. |
conf.level |
the confidence level for the interval. |
ci.points |
the number of thresholds used in the calculation of the confidence intervals. |
xscale |
the scale of the x-axis (options include "o" = original, "l" = log scale, "b" = both). |
Details
For more details about the estimator \hat{t}
and the computation of the confidence intervals see pareto_tail.
Value
A plot showing the estimated t(u)
versus threshold u
, optionally on a logarithmic x-axis and including confidence intervals. Note that on the right side of the plot, one can observe the corresponding alpha values, which indicate the shape parameter of the Pareto distribution associated with the estimated t-values.
References
Klar, B. (2024). A Pareto tail plot without moment restrictions. The American Statistician. doi:10.1080/00031305.2024.2413081
Examples
x <- actuar::rpareto1(1e3, shape=1, min=1)
pareto_tailplot(x, method="unbiased", xscale="o")