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

In freakonometrics/TopIncome: Estimation of Top Incomes Using Pareto Models

Description

Maximum Likelihood estimation of the Pareto 1 distribution, with weights

Usage

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

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

Examples

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/TopIncome documentation built on Oct. 16, 2021, 5:58 p.m.