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

KUMPPF

Name:
    KUMPPF (LET)
Type:
    Library Function
Purpose:
    Compute the Kumaraswamy percent point function with shape parameters alpha and beta.
Description:
    The standard Kumaraswamy distribution has the following percent point function:

      G(p;alpha,beta) = [1 - (1-p)**((1/beta)]**(1/alpha)     
0 <= p <= 1; alpha, beta > 0

    with alpha and beta denoting the shape parameters.

    This distribution can be extended with lower and upper bound parameters. If a and b denote the lower and upper bounds, respectively, then the location and scale parameters are:

      location = a
      scale = b - a

    The general form of the distribution can then be found by using the relation

      G(p;alpha,beta,a,b) = loc + scale*G(p;alpha,beta,0,1)

    This distribution has been proposed as a more tractable alternative to the beta distribution.

Syntax:
    LET <y> = KUMPPF(<p>,<alpha>,<beta>,<a>,<b>)
                            <SUBSET/EXCEPT/FOR qualification>
    where <p> is a number, parameter, or variable containing values in the interval (0,1);
                <y> is a variable or a parameter (depending on what <p> is) where the computed Kumaraswamy ppf value is stored;
                <alpha> is a positive number, parameter, or variable that specifies the first shape parameter;
                <beta> is a positive number, parameter, or variable that specifies the second shape parameter;
                <a> is a number, parameter, or variable that specifies the lower limit;
                <b> is a number, parameter, or variable that specifies the upper limit;
    and where the <SUBSET/EXCEPT/FOR qualification> is optional.

    If <a> and <b> are omitted, they default to 0 and 1, respectively.

Examples:
    LET A = KUMPPF(0.95,0.2,1.2)
    LET Y = KUMPPF(P,0.5,2)
    PLOT KUMPPF(P,2,3) FOR P = 0 0.01 1
Default:
    None
Synonyms:
    None
Related Commands:
    KUMCDF = Compute the Kumaraswamy cumulative distribution function.
    KUMPDF = Compute the Kumaraswamy probability density function.
    RGTPDF = Compute the reflected generalized Topp and Leone probability density function.
    GTLPDF = Compute the generalized Topp and Leone probability density function.
    TOPPDF = Compute the Topp and Leone probability density function.
    TSPPDF = Compute the two-sided power probability density function.
    BETPDF = Compute the beta probability density function.
    TRIPDF = Compute the triangular probability density function.
    TRAPDF = Compute the trapezoid probability density function.
    UNIPDF = Compute the uniform probability density function.
    POWPDF = Compute the power probability density function.
    JSBPDF = Compute the Johnson SB probability density function.
Reference:
    Kumaraswamy (1980), "A Generalized Probability Density Function for Double-Bounded Random Processes", Journal of Hydrology, 46: 79-88.
Applications:
    Distributional Modeling
Implementation Date:
    2007/11
Program:
     
    CASE ASIS
    LABEL CASE ASIS
    TITLE CASE ASIS
    TITLE OFFSET 2
    .
    MULTIPLOT 3 3
    MULTIPLOT CORNER COORDINATES 0 0 100 95
    MULTIPLOT SCALE FACTOR 3
    .
    LET ALPHA = 2
    LET BETA  = 3
    TITLE Alpha = ^alpha, Beta = ^beta
    PLOT KUMPPF(P,ALPHA,BETA) FOR P = 0  0.01  1
    .
    LET ALPHA = 1.5
    LET BETA  = 6
    TITLE Alpha = ^alpha, Beta = ^beta
    PLOT KUMPPF(P,ALPHA,BETA) FOR P = 0  0.01  1
    .
    LET ALPHA = 1.5
    LET BETA  = 2
    TITLE Alpha = ^alpha, Beta = ^beta
    PLOT KUMPPF(P,ALPHA,BETA) FOR P = 0  0.01  1
    .
    LET ALPHA = 1.5
    LET BETA  = 1
    TITLE Alpha = ^alpha, Beta = ^beta
    PLOT KUMPPF(P,ALPHA,BETA) FOR P = 0  0.01  1
    .
    LET ALPHA = 0.5
    LET BETA  = 2
    TITLE Alpha = ^alpha, Beta = ^beta
    PLOT KUMPPF(P,ALPHA,BETA) FOR P = 0  0.01  1
    .
    LET ALPHA = 0.5
    LET BETA  = 1
    TITLE Alpha = ^alpha, Beta = ^beta
    PLOT KUMPPF(P,ALPHA,BETA) FOR P = 0  0.01  1
    .
    LET ALPHA = 0.5
    LET BETA  = 0.75
    TITLE Alpha = ^alpha, Beta = ^beta
    PLOT KUMPPF(P,ALPHA,BETA) FOR P = 0  0.01  1
    .
    LET ALPHA = 0.5
    LET BETA  = 0.25
    TITLE Alpha = ^alpha, Beta = ^beta
    PLOT KUMPPF(P,ALPHA,BETA) FOR P = 0  0.01  1
    .
    LET ALPHA = 1
    LET BETA  = 1
    TITLE Alpha = ^alpha, Beta = ^beta
    PLOT KUMPPF(P,ALPHA,BETA) FOR P = 0  0.01  1
    .
    END OF MULTIPLOT
    .
    JUSTIFICATION CENTER
    MOVE 50 97
    TEXT Kumaraswamy Percent Point Functions
     
        

    plot generated by sample program

Date created: 11/27/2007
Last updated: 11/27/2007
Please email comments on this WWW page to alan.heckert@nist.gov.