fitPareto: Fit a Pareto Distribution to Binned Data

In bda: Binned Data Analysis

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

bda documentation built on Sept. 11, 2024, 9:01 p.m.