plot.starship: Plots to compare a fitted generalised lambda distribution to data

Description

Plots to compare a Generalised Lambda Distribution fitted via the starship to data

Usage

"plot"(x, data, ask = FALSE, one.page = TRUE,  breaks = "Sturges", plot.title = "default",...)

Arguments

An object of class starship.

data

Data to which the gld was fitted. Leave this as NULL if the return.data argument was TRUE in the starship call that created x.

ask

Ask for user input before next plot --- passed to par(ask). Does not permanently change this setting. Ignored if one.page is TRUE

one.page

Put the two plots on one page using par(mfrow=c(2,1)). Does not permanently change this setting.

breaks

Control the number of histogram bins --- passed to hist.

plot.title

Main title for histogram and QQ --- passed to hist(main=) and qqgl(main=). If you set this to "default", it will include the fitting method and gld type, for example ``Starship fit of FMKL type GLD''.

...

arguments passed to plot AND hist

Details

summary Gives the details of the starship.adaptivegrid and optim steps.

References

Freimer, M., Mudholkar, G. S., Kollia, G. & Lin, C. T. (1988), A study of the generalized tukey lambda family, Communications in Statistics - Theory and Methods 17, 3547--3567.

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

http://tolstoy.newcastle.edu.au/~rking/gld/

Examples

Run this code

data <- rgl(100,0,1,.2,.2)
starship.result <- starship(data,optim.method="Nelder-Mead",initgrid=list(lcvect=(0:4)/10,
ldvect=(0:4)/10),return.data=TRUE)
plot(starship.result)