Has anyone had problems with the Weibull distribution using the ExtDist Package?
From the documentation:
Parameter Estimation for a distribution with unknown shape parameters
Example from: Rinne (2009) Dataset p.338 and example pp.418-419
Parameter estimates are given as shape = 99.2079 and scale = 2.5957.
The log-likelihood for this data and Rinne's parameter estimates is
-1163.278.
data <- c(35,38,42,56,58,61,63,76,81,83,86,90,99,104,113,114,117,119,141,183)
est.par <- eWeibull(X=data, method="numerical.MLE"); est.par
plot(est.par)
However when I run this I get the following output:
Parameters for the Weibull distribution.
(found using the numerical.MLE method.)
Parameter Type Estimate S.E.
shape shape 5.82976007 1.79326460
scale scale 0.06628166 0.02129258
This is clearly wrong but I am not sure if I have made a mistake or if there is a bug in the package?