powered by
The maximum likelihood estimate of b is the minimum of x and the maximum likelihood estimate of a is 1/(mean(log(x)) - log(b)).
b
x
a
1/(mean(log(x)) - log(b))
mlpareto(x, na.rm = FALSE, ...)
mlpareto returns an object of class
mlpareto
univariateML. This is a named numeric vector with maximum likelihood estimates for a and b and the following attributes:
univariateML
model
The name of the model.
density
The density associated with the estimates.
logLik
The loglikelihood at the maximum.
support
The support of the density.
n
The number of observations.
call
The call as captured my match.call
match.call
a (non-empty) numeric vector of data values.
logical. Should missing values be removed?
currently affects nothing.
For the density function of the Pareto distribution see Pareto.
Johnson, N. L., Kotz, S. and Balakrishnan, N. (1995) Continuous Univariate Distributions, Volume 1, Chapter 20. Wiley, New York.
Pareto for the Pareto density.
mlpareto(precip)
Run the code above in your browser using DataLab