Maximum likelihood fit for Siler distribution

Do a maximum likelihood Siler fit for the input data vector x. To allow unconstrained optimization to be used, optimization is done on the transformed variable abar = log(a/a0).

fitSilerMaxLik(x, a0 = c(0.175, 1.4, 0.368 * 0.01, 0.075 * 0.001, 0.917 *
  0.1), x0 = 0, calcHessian = FALSE)

Arguments

x

Locations at which to evaluate probability density function
a0

Initial value for optimization. Default from Gage and Dyke 1986, Table 2, Level 15
x0

x-offset. The default is 0 (no offset).
calcHessian

Whether to calculate the Hessian

Value

A list consisting of the fit (on the transformed variable abar) and maximum likelihood estimate of a. Optionally, the Hessian of the log-likelihood is returned, which allows estimation of the standard errors of the maximum likelihood estimate via the observed Fisher information matrix.

See also

dsiler