fit.Pareto {bda}

Fit a Pareto Distribution to Binned Data

Description

Fit a Pareto distribution to binned data.

Usage

  fit.Pareto(x, xm, method='mle')

Arguments

x

grouped data

xm

The location parameter: lower bound of the support of the distribution

method

fitting method: 'mle'=maximum likelihood estimate, 'percentile'=percentile matching.

Value

xm

fitted location parameter

alpha

fitted scale parameter

Examples

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'))