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

DPUPPF

Name:
    DPUPPF (LET)
Type:
    Library Function
Purpose:
    Compute the doubly Pareto uniform percent point function
Description:
    The doubly Pareto uniform distribution has the following percent point function:

      G(p;m,n,alpha,beta) = lambda1*(beta-alpha) + alpha   0 < p < pi1;
  lambda2*(beta-alpha) + alpha   pi1 <= p <= pi2;
  lambda3*(beta-alpha) + alpha   pi2 < p < 1;
  alpha < beta; m, n > 0

    where

      pi1 = n/(m + m*n + n)
      pi2 = (m*n)/(m + m*n + n)
      pi3 = m/(m + m*n + n)

      lambda1 = 1 - (pi1/p)**(1/m)
      lambda2 = (p - pi1)/pi2
      lambda3 = (pi3/(1-p))**(1/n)

    and with m and n denoting the shape parameters, alpha denoting the location parameter, and (beta - alpha) denoting the scale parameter.

    This distribution is uniform between alpha and beta. It has Paretian tails for both the lower and upper tails. The m parameter controls the shape of the lower tail and the n parameter controls the shape of the upper tail.

    The case where alpha = 0 and beta = 1 is referred to as the standard doubly Pareto uniform distribution.

Syntax:
    LET <y> = DPUPPF(<p>,<m>,<n>,<alpha>,<beta>)
                            <SUBSET/EXCEPT/FOR qualification>
    where <p> is a number, parameter, or variable in the interval (0,1);
                <y> is a variable or a parameter (depending on what <p> is) where the computed double Pareto uniform ppf value is stored;
                <m> is a number, parameter, or variable that specifies the first shape parameter;
                <n> is a number, parameter, or variable that specifies the second shape paremeter;
                <alpha> is a number, parameter, or variable that specifies the location parameter;
                <beta> is a number, parameter, or variable (<beta> - <alpha> is the scale parameter);
    and where the <SUBSET/EXCEPT/FOR qualification> is optional.

    If <alpha> and <beta> are omitted, they default to 0 and 1, respectively.

Examples:
    LET A = DPUPPF(0.95,1.4,3.2,0,1)
    LET Y = DPUPPF(P,1.4,3.2,0,1)
    PLOT DPUPPF(P,1.4,3.2,-5,5) FOR X = 0.01 0.01 0.99
Default:
    None
Synonyms:
    None
Related Commands:
    DPUCDF = Compute the doubly Pareto uniform cumulative distribution function.
    DPUPDF = Compute the doubly Pareto uniform probability density function.
    TSPPDF = Compute the two-sided power probability density function.
    POWPDF = Compute the power probability density function.
    GTRPDF = Compute the generalized trapezoid probability density function.
    TSOPDF = Compute the two-sided ogive probability density function.
    OGIPDF = Compute the ogive probability density function.
    TSSPDF = Compute the two-sided slope probability density function.
    SLOPDF = Compute the slope probability density function.
    BETPDF = Compute the Beta probability density function.
    JSBPDF = Compute the Johnson SB probability density function.
Reference:
    Singh, Van Dorp, Mazzuchi "A Novel Asymmetric Distribution with Power Tails", Communications in Statistics, Theory and Methods, Vol. 36 (2), to appear.

    Van Dorp, Singh, and Mazzuchi "The Doubly-Pareto Uniform Distribution with Applications in Uncertainty Analysis and Econometrics", Mediterranean Journal of Mathematics, Vol. 3 (2), pp. 205-225.

Applications:
    Distributional Modeling
Implementation Date:
    2007/10
Program:
     
    LABEL CASE ASIS
    TITLE CASE ASIS
    TITLE OFFSET 2
    .
    MULTIPLOT 2 2
    MULTIPLOT CORNER COORDINATES 0 0 100 95
    MULTIPLOT SCALE FACTOR 2
    .
    LET ALPHA = 0
    LET BETA  = 1
    .
    LET M = 0.5
    LET N = 0.5
    TITLE M = ^m, N = ^n
    PLOT DPUPPF(P,M,N,ALPHA,BETA) FOR P = 0.01  0.01  0.99
    .
    LET M = 2
    LET N = 0.5
    TITLE M = ^m, N = ^n
    PLOT DPUPPF(P,M,N,ALPHA,BETA) FOR P = 0.01  0.01  0.99
    .
    LET M = 0.5
    LET N = 2
    TITLE M = ^m, N = ^n
    PLOT DPUPPF(P,M,N,ALPHA,BETA) FOR P = 0.01  0.01  0.99
    .
    LET M = 2
    LET N = 2
    TITLE M = ^m, N = ^n
    PLOT DPUPPF(P,M,N,ALPHA,BETA) FOR P = 0.01  0.01  0.99
    .
    END OF MULTIPLOT
    .
    CASE ASIS
    JUSTIFICATION CENTER
    MOVE 50 97
    TEXT Doubly Pareto Uniform Percent Point Functions
        
    plot generated by sample program

Date created: 1/8/2008
Last updated: 1/8/2008
Please email comments on this WWW page to alan.heckert@nist.gov.