starship
for more details.This function is the objective funciton minimised in the methods. It is a goodness of fit measure carried out on the depths of the data.
starship.obj(par, data, inverse.eps, param = "fmkl")
GeneralisedLambdaDistribution
for details on the
definitions of these parametersfmkl
uses Freimer, Mudholkar, Kollia and Lin (1988) (default).
rs
uses Ramberg and Schmeiser (1974)gld
)
is a transformation of the uniform distribution. Thus the inverse of this
transformation is the distribution function for the gld. The starship method
applies different values of the parameters of the distribution to the
distribution function, calculates the depths q corresponding to the data
and chooses the parameters that make the depths closest to a uniform
distribution.The closeness to the uniform is assessed by calculating the Anderson-Darling
goodness-of-fit test on the transformed data against the uniform, for a
sample of size length(data)
.
This function returns that objective function. It is provided as a seperate
function to allow users to carry out minimisations using optim
or other methods. The recommended method is to use the starship
function.
Ramberg, J. S. & Schmeiser, B. W. (1974), An approximate method for generating asymmetric random variables, Communications of the ACM 17, 78--82. King, R.A.R. & MacGillivray, H. L. (1999), A starship method for fitting the generalised $lambda$ distributions, Australian and New Zealand Journal of Statistics 41, 353--374
Owen, D. B. (1988), The starship, Communications in Statistics - Computation and Simulation 17, 315--323.
starship
,
starship.adaptivegrid
data <- rgl(100,0,1,.2,.2)
starship.obj(c(0,1,.2,.2),data,inverse.eps=1e-10,"fmkl")
Run the code above in your browser using DataLab