fun.RPRS.mm: Fit RS generalised lambda distribution to data set using moment
matching
Description
This function fits RS generalised lambda distribution to data set
using moment matching
Usage
fun.RPRS.mm(data, rs.init = c(-1.5, 1.5), leap = 3, FUN = "runif.sobol",
no = 10000)
Value
A vector representing four parameters of the RS generalised lambda distribution.
Arguments
data
Dataset to be fitted
rs.init
Initial values for RS distribution optimization,
c(-1.5,1.5) tends to work well.
leap
See scrambling argument in fun.gen.qrn.
FUN
A character string of either "runif.sobol" (default), "runif.sobol.owen",
"runif.halton" or "QUnif".
no
Number of initial random values to find the best initial values for
optimisation.
Author
Steve Su
Details
This function provides method of moment fitting scheme for RS GLD. Note this
function can fail if there are no defined percentiles from the data set or if
the initial values do not lead to a valid RS generalised lambda distribution.
This function is based on scheme detailed in the literature below but it has
been modified by the author (Steve Su).
References
Karian, Z. and E. Dudewicz (2000). Fitting Statistical
Distributions: The Generalized Lambda Distribution and Generalised Bootstrap
Methods. New York, Chapman and Hall.