Replace observations under a Pareto model
Replace observations under a Pareto model for the upper tail with values drawn from the fitted distribution.
replaceTail(x, ...) ## S3 method for class 'paretoTail' replaceTail(x, all = TRUE, ...) replaceOut(x, ...)
x |
an object of class |
... |
additional arguments to be passed down. |
all |
a logical indicating whether all observations in the upper tail should be replaced or only those flagged as outliers. |
replaceOut(x, ...{}) is a simple wrapper for replaceTail(x,
all = FALSE, ...{}).
A numeric vector consisting mostly of the original values, but with observations in the upper tail replaced with values from the fitted Pareto distribution.
Andreas Alfons
A. Alfons and M. Templ (2013) Estimation of Social Exclusion Indicators from Complex Surveys: The R Package laeken. Journal of Statistical Software, 54(15), 1–25. URL http://www.jstatsoft.org/v54/i15/
A. Alfons, M. Templ, P. Filzmoser (2013) Robust estimation of economic indicators from survey samples based on Pareto tail modeling. Journal of the Royal Statistical Society, Series C, 62(2), 271–286.
data(eusilc)
## gini coefficient without Pareto tail modeling
gini("eqIncome", weights = "rb050", data = eusilc)
## gini coefficient with Pareto tail modeling
# estimate threshold
ts <- paretoScale(eusilc$eqIncome, w = eusilc$db090,
groups = eusilc$db030)
# estimate shape parameter
fit <- paretoTail(eusilc$eqIncome, k = ts$k,
w = eusilc$db090, groups = eusilc$db030)
# replacement of outliers
eqIncome <- replaceOut(fit)
gini(eqIncome, weights = eusilc$rb050)
# replacement of whole tail
eqIncome <- replaceTail(fit)
gini(eqIncome, weights = eusilc$rb050)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.