Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!
Get Started for Free

fun.RMFMKL.qs

Fit FMKL generalised lambda distribution to data set using quantile matching


Description

This function fits FMKL generalised lambda distribution to data set using quantile matching

Usage

fun.RMFMKL.qs(data, fmkl.init = c(-0.25, 1.5), leap = 3, FUN = "runif.sobol", 
trial.n = 100, len = 1000, type = 7, no = 10000)

Arguments

data

Dataset to be fitted

fmkl.init

Initial values for FMKL distribution optimization, c(-0.25,1.5) tends to work well.

leap

Scrambling (0,1,2,3) for the Sobol sequence for the distribution fit. See scrambling/leap argument for runif.sobol, runif.halton or QUnif.

FUN

A character string of either "runif.sobol" (default), "runif.halton" or "QUnif".

trial.n

Number of evenly spaced quantile ranging from 0 to 1 to be used in the checking phase, to find the best set of initial values for optimisation, this is intended to be lower than len to speed up the fitting algorithm. Default is 100.

len

Number of evenly spaced quantile ranging from 0 to 1 to be used, default is 1000

type

Type of quantile to be used, default is 7, see quantile

no

Number of initial random values to find the best initial values for optimisation.

Details

This function provides quantile matching fitting scheme for FMKL 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 FMKL generalised lambda distribution.

Value

A vector representing four parametefmkl of the FMKL generalised lambda distribution.

Author(s)

Steve Su

References

Su (2008). Fitting GLD to data via quantile matching method. (Book chapter to appear)

See Also

Examples

## Fitting the normal distribution
# fun.RMFMKL.qs(data=rnorm(1000,2,3),fmkl.init=c(-0.25,1.5),leap=3)

GLDEX

Fitting Single and Mixture of Generalised Lambda Distributions (RS and FMKL) using Various Methods

v2.0.0.7
GPL (>= 3)
Authors
Steve Su, with contributions from: Diethelm Wuertz, Martin Maechler and Rmetrics core team members for low discrepancy algorithm, Juha Karvanen for L moments codes, Robert King for gld C codes and starship codes, Benjamin Dean for corrections and input in ks.gof code and R core team for histsu function.
Initial release
2020-02-04

We don't support your browser anymore

Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.