This function can be used to estimate the shape
parameter using the
Maximum Likelihood Estimator method (Newman 2005). It can be used to
obtain biased and unbiased estimates of the shape and scale parameters as
well as the confidence interval for the shape parameter for the biased estimates.
alpha_mle(dat, biased = TRUE, significance = NULL)
vector of observations
TRUE/FALSE to indicate biased or unbiased estimates
level of significance
A list of the following form:
Estimate of the shape parameter of the data
Upper error bound of the estimate of shape
Lower error bound of the estimate of shape
Estimate of the scale parameter of the data (which is taken to be the minimum of the data)
Newman MEJ (2005). "Power Laws, Pareto Distributions And Zipf's Law." Contemporary Physics, 46, 323-351.
# NOT RUN {
x <- generate_pareto(10000, 5, 2)
alpha_mle(x, TRUE, 0.05)
x <- generate_pareto(10000, 5, 2)
alpha_mle(x, FALSE)
# }
Run the code above in your browser using DataLab