![]() |
MCLPDFName:
![]()
with K The standard McLeish distribution can be generalized with location and scale parameters in the usual way.
<SUBSET/EXCEPT/FOR qualification> where <x> is a variable, a number, or a parameter; <alpha> is a positive number of parameter that specifies the value of the shape parameter; <loc> is an optional number or parameter that specifies the value of the location parameter; <scale> is an optional positive number or parameter that specifies the value of the scale parameter; <y> is a variable or a parameter (depending on what <x> is) where the computed McLeish pdf value is stored; and where the <SUBSET/EXCEPT/FOR qualification> is optional.
LET Y = MCLPDF(X1,ALPHA) PLOT MCLPDF(X,ALPHA) FOR X = -10 0.01 10
LET Y = MCLEISH RANDOM NUMBERS FOR I = 1 1 N To generate a McLeish probability plot or a McLeish Kolmogorov-Smirnov or chi-square goodness of fit test, enter the following commands
MCLEISH PROBABILITY PLOT Y MCLEISH KOLMOGOROV SMIRNOV GOODNESS OF FIT Y MCLEISH CHI-SQUARE GOODNESS OF FIT Y To generate a PPCC or Kolmogorov-Smirnov plot, enter the following commands
LET ALPHA2 = <value> MCLEISH PPCC PLOT Y MCLEISH KS PLOT Y The default values for ALPHA1 and ALPHA2 are 1.0 and 15.5. For the McLeish distribution, the shape parameter acts a bit like a scale parameter. For this reason, the KS PLOT with the location and scale parameters fixed will probably work better than the PPCC PLOT. One recommendation is to set the scale parameter to 1 and the location parameter to the mode of the data (the McLeish distribution is symmetric). You can use the relative histogram or the kernel density plot to determine a useful value for the location. For example,
LET KSSCALE = 1 MCLEISH KS PLOT Y
Y1LABEL Probability X1LABEL X LABEL CASE ASIS TITLE CASE ASIS CASE ASIS Y1LABEL DISPLACEMENT 12 MULTIPLOT 2 2 MULTIPLOT CORNER COORDINATES 0 0 100 95 MULTIPLOT SCALE FACTOR 2 TITLE Alpha = 1 PLOT MCLPDF(X,1) FOR X = -10 0.01 10 TITLE Alpha = 2 PLOT MCLPDF(X,2) FOR X = -10 0.01 10 TITLE Alpha = 5 PLOT MCLPDF(X,5) FOR X = -10 0.01 10 TITLE Alpha = 10 PLOT MCLPDF(X,10) FOR X = -10 0.01 10 END OF MULTIPLOT MOVE 50 97 JUSTIFICATION CENTER TEXT McLeish Distribution ![]()
Date created: 4/20/2005 |