plot.starship | R Documentation |
Plots to compare a fitted Generalised Lambda Distribution to data. This works with fitted gld objects from starship
or fit.fkml
.
## S3 method for class 'starship' plot(x, data, ask = NULL, one.page = FALSE, breaks = "Sturges", plot.title = "default",...)
x |
An object of class |
data |
Data to which the gld was fitted. Leave this as NULL if the return.data argument was TRUE in the call that created ( |
ask |
Ask for user input before next plot. The default of NULL changes to TRUE if one.page is FALSE and |
one.page |
If TRUE, put the two plots on one page using |
breaks |
Control the number of histogram bins — passed to |
plot.title |
Main title for histogram and QQ — passed
to |
... |
arguments passed to |
This function produces a Q-Q plot of the data against the fitted generalised lambda distribution and a histogram of the data overlaid with the fitted probability density function.
Robert King, robert.king.newcastle@gmail.com, https://github.com/newystats/
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
https://github.com/newystats/gld/
starship
,
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.