A starship estimation method for the generalized λ distributions

A development of the 'starship' method (Owen, 1988), a computer intensive estimation method, is presented for two forms of generalized lambda distributions (g lambda d). The method can be used for the full parameter space and is flexible, allowing choice of both the form of the generalized lambda distribution and of the nature of fit required. Some examples of its use in fitting data and approximating distributions are given. Some simulation studies explore the sampling distribution of the parameter estimates produced by this method for selected values of the parameters and consider comparisons with two other methods, for one of the g lambda d distributional forms, not previously so investigated. In the forms and parameter regions available to the other methods, it is demonstrated that the starship compares favourably. Although the differences between the methods, where available, tend to disappear with larger samples, the parameter coverage, flexibility and adaptability of the starship method make it attractive. However, the paper also demonstrates that care is needed when fitting and using such quantile-defined distributional families that are rich in shape, but have complex properties.