laeken (version 0.3.3)

replaceTail: Replace observations under a Pareto model

Description

Replace observations under a Pareto model for the upper tail with values drawn from the fitted distribution.

Usage

replaceTail(x, ...)

  ## S3 method for class 'paretoTail':
replaceTail(x, all = TRUE, ...)

  replaceOut(x, ...)

Arguments

x

an object of class "paretoTail" (see paretoTail).

all

a logical indicating whether all observations in the upper tail should be replaced or only those flagged as outliers.

...

additional arguments to be passed down.

Value

• A numeric vector consisting mostly of the original values, but with observations in the upper tail replaced with values from the fitted Pareto distribution.

code

replaceTail(x, all = FALSE, ...

Details

replaceOut(x, ...)

See Also

paretoTail, reweightOut, shrinkOut

Examples

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)