MLE.pareto1: Maximum Likelihood estimation of the Pareto 1 distribution,...

Description Usage Arguments Value References Examples

Description

Maximum Likelihood estimation of the Pareto 1 distribution, with weights

Usage

1
MLE.pareto1(data, weights = rep(1, length(x)), threshold = min(x))

Arguments

data

a vector of observations

weights

a vector of weights (default = 1)

threhold

the threshold parameter of the Pareto 1 distribution (mu)

Value

a list with the index alpha and k, the number of observations above threshold

References

Charpentier & Flachaire (2019) Pareto Models for Top Incomes hal-02145024

Examples

1
2
3
4
5
set.seed(123)
x <- rpareto1(100, 1, 1.5)
w <- rgamma(100, 10, 10)
estim <- MLE.pareto1(data=x, weights=w, threshold=1)
estim

freakonometrics/TopIncomes documentation built on July 7, 2019, 8:06 a.m.