This function uses the Hill's Estimator to estimate the shape parameter
of a given set of data. (Nair et al. 2019; Pokorna 2016; Hill 1975)
It is especially useful when the data is known not to follow an exact
Pareto distribution but the tail of the data does. Thus, the specification
of k
, the k
th largest observation, allows to specify the
point from where Pareto-like behavior may be seen. It is also possible to
specify the value at which the tail begins.
When k=n
, the Hill's Estimator returns the same estimate as
alpha_mle
with a warning notifying the user.
alpha_hills(dat, k, value = FALSE)
vector of observations
number of observations / value equal to or greater than to consider for tail
(TRUE/FALSE) indicating if the value which is specified in "k" (TRUE)
A list of the following form:
Estimate of the shape parameter of the data
Estimate of the scale parameter of the data (which is taken to be the minimum of the data)
Nair J, Wierman A, Zwart B (2019). "The Fundamentals Of Heavy Tails: Properties, Emergence, And Identification." http://users.cms.caltech.edu/ adamw/heavytails.html.
Pokorna M (2016). Estimation and Application of the Tail Index. Bachelor's thesis, Charles University in Prague, Faculty of Social Sciences, Institute of Economic Studies.
Hill B (1975). "A Simple General Approach To Inference About The Tail Of A Distribution."The Annals of Statistics, 3(5), 1163-1174.
# NOT RUN {
x <- generate_pareto(10000, 5, 2)
alpha_hills(x, 400)
# }
Run the code above in your browser using DataLab