SCI (version 1.0-2)

dist.start: Rough estimates for parameters of selected distributions

Description

Produces rough parameter estimates for specific distributions (distr) that are useful as starting values for maximum likelihood estimation.

Usage

dist.start(x, distr,...)
lmom.start(x, distr = c("gamma","genlog","gev","gumbel", "lnorm","norm","pe3","weibull"),...)
mom.start(x, distr = c("gamma","gumbel","logis","lnorm","norm", "weibull"),...)

Arguments

x

numeric vector

distr

A character string "name" naming a distribution for which the corresponding density function (dname), the corresponding distribution function (pname) and the quantile function (qname) must be defined (see for example GammaDist)

...

arguments passed to other functions, currently not used.

Value

named list, names correspond to distribution parameters. In case of failure, the same list with NA values is returned.

Details

lmom.start uses L-moments for parameter estimation. In most cases it relies on functionality of the lmomco package. Currently available distributions are: "gamma", "genlog", "gev", "gumbel", "logis", "lnorm", "norm", "pe3", "weibull".

mom.start uses moments (e.g. mean, standard deviation) for parameter estimation. Some estimates are precise, others only approximations that provide reasonable starting values. Currently available distributions are: "gamma", "gumbel", "logis", "lnorm", "norm", "weibull". dist.start calls first lmom.start to estimate parameters. In case of failure mom.start is called, hopefully producing reasonable parameter estimates.

Examples

lmom.start(rgamma(100,shape=0.5,rate=1),"gamma")

mom.start(rgamma(100,shape=0.5,rate=1),"gamma")

dist.start(rgamma(100,shape=0.5,rate=1),"gamma")