
replaceTail(x, ...) ## S3 method for class 'paretoTail':
replaceTail(x, all = TRUE, ...)
replaceOut(x, ...)
"paretoTail"
(see
paretoTail
).replaceTail(x, all = FALSE, ...
replaceOut(x, ...
)paretoTail
, reweightOut
,
shrinkOut
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)
Run the code above in your browser using DataLab