xm <- 0.5
alpha <- 1.0
x <- rPareto(1000, xm, alpha)
(out <- fit.Pareto(x,method='mle'))
(out <- fit.Pareto(x,method='ls'))
xbrks <- c(0,4.5,9.5,19.5,49.5,99.5,249.5,499.5,999.5,
2499.5,4999.5,9999.5,Inf)
xhist <- binning(x, breaks=xbrks)
(out <- fit.Pareto(xhist))
(out <- fit.Pareto(xhist,method='mle'))
(out <- fit.Pareto(xhist,method='ls'))
(out <- fit.Pareto(xhist,xm=.5,method='mle'))
Run the code above in your browser using DataLab