SED navigation bar go to SED home page go to Dataplot home page go to NIST home page SED Home Page SED Staff SED Projects SED Products and Publications Search SED Pages
Dataplot Vol 2 Vol 1

GMCCDF

Name:
    GMCCDF (LET)
Type:
    Library Function
Purpose:
    Compute the generalized McLeish cumulative distribution function.
Description:
    The standard form of the generalized McLeish distribution has the following probability density function:

      f(X;ALPHA,A) = [1/(SQRT(PI)*GAMMA(ALPHA))]*
(ABS(X)/2)**(ALPHA-1/2)*K(X,ALPHA-1/2)*(1-A**2)**ALPHA*EXP(A*X)
   alpha > 0, |A| < 1

    with Klambda(.) denoting the modified Bessel function of the of the second kind of order lambda and GAMMA denoting the gamma function.

    The standard generalized McLeish distribution can be generalized with location and scale parameters in the usual way.

    The cumulative distribution function is computed by numerically integrating the probability density function. Dataplot performs the integration using the DQAG routine from the Slatec library.

Syntax:
    LET <y> = GMCCDF(<x>,<alpha>,<a>,<loc>,<scale>)
                            <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 first shape parameter;
                <a> is a positive number of parameter that specifies the value of the second 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 generalized McLeish cdf value is stored;
    and where the <SUBSET/EXCEPT/FOR qualification> is optional.
Examples:
    LET Y = GMCCDF(3,1.5,0.8)
    LET Y = GMCCDF(X1,ALPHA,A)
    PLOT GMCCDF(X,ALPHA,A) FOR X = -10 0.01 10
Note:
    DATAPLOT uses the routine DBESK from the SLATEC Common Mathematical Library to compute the modified Bessel function of the second kind. SLATEC is a large set of high quality, portable, public domain Fortran routines for various mathematical capabilities maintained by seven federal laboratories.
Default:
    None
Synonyms:
    None
Related Commands:
    GMCPDF = Compute the generalized McLeish probability density function.
    GMCPPF = Compute the generalized McLeish percent point function.
    MCLPDF = Compute the McLeish probability density function.
    GALPDF = Compute the generalized asymmetric Laplace cumulative distribution function.
    GIGPDF = Compute the generalized inverse Gaussian probability density function.
    BEIPDF = Compute the Bessel I-function probability density function.
    BEKPDF = Compute the Bessel K-function probability density function.
Reference:
    Johnson, Kotz, and Balakrisnan, "Continuous Univariate Distributions--Volume I", Second Edition, Wiley, 1994, pp. 50-53.
Applications:
    Distributional Modeling
Implementation Date:
    8/2004
Program:
     
    Y1LABEL Probability
    X1LABEL X
    LABEL CASE ASIS
    TITLE CASE ASIS
    CASE ASIS
    Y1LABEL DISPLACEMENT 16
    MULTIPLOT 2 2
    MULTIPLOT CORNER COORDINATES 0 0 100 95
    MULTIPLOT SCALE FACTOR 2
    LET A = 0.8
    TITLE Alpha = 1.5, A = 0.8
    PLOT GMCCDF(X,1.5,A) FOR X = -30  0.1 30
    LET A = -0.8
    TITLE Alpha = 1.5, A = -0.8
    PLOT GMCCDF(X,1.5,A) FOR X = -30  0.1 30
    LET A = 0.2
    TITLE Alpha = 1.5, A = 0.2
    PLOT GMCCDF(X,1.5,A) FOR X = -30  0.1 30
    LET A = -0.2
    TITLE Alpha = 1.5, A = -0.2
    PLOT GMCCDF(X,1.5,A) FOR X = -30  0.1 30
    END OF MULTIPLOT
    MOVE 50 97
    JUSTIFICATION CENTER
    TEXT Generalized McLeish Cumulative Distribution
        
    plot generated by sample program

Date created: 4/19/2005
Last updated: 4/19/2005
Please email comments on this WWW page to alan.heckert@nist.gov.