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

RIGPPF

Name:
    RIGPPF (LET)
Type:
    Library Function
Purpose:
    Compute the reciprocal inverse Gaussian percent point function with shape parameters gamma and mu.
Description:
    The reciprocal inverse Gaussian distribution is the distribution of (1/X) when X has an inverse Gaussian distribution. It has the following cumulative distribution function:

      F(x;gamma,mu) = NORCDF[-(1/(mu*x) - 1)*SQRT(gamma*x)] - 
                 EXP(2/gamma**2)*
                 NORCDF[-(1/(mu*x) - 1)*SQRT(gamma*x)]
                 x >= 0; gamma, mu > 0

    with gamma and mu denoting the shape parameters and NORCDF denoting the cumulative distribution function of the standard normal distribution.

    The percent point function is the inverse of the cumulative distribution function. The percent point function for the reciprocal inverse Gaussian distribution does not exist in simple, closed form.

    The reciprocal inverse Gaussian percent point function can be computed in terms of the inverse Gaussian percent point function by

      RIGPPF(p;gamma,mu) = 1/IGPPF(1-p;gamma,mu)

    with IGPPF denoting the percent point function of the inverse Gaussian distribution.

    Dataplot uses this relationship to compute the reciprocal inverse gaussian percent point function.

    The reciprocal inverse Gaussian distribution can be generalized with location and scale parameters in the usual way.

Syntax:
    LET <y> = RIGCHAZ(<p>,<gamma>,<mu>,<loc>,<scale>)
                            <SUBSET/EXCEPT/FOR qualification>
    where <p> is a variable or a parameter;
                <gamma> is number or parameter that specifies the first shape parameter;
                <mu> is number or parameter that specifies the second shape parameter;
                <loc> is number or parameter that specifies the location parameter;
                <scale> is number or parameter that specifies the scale parameter;
                <y> is a variable or a parameter (depending on what <x> is) where the computed reciprocal inverse Gaussian ppf value is stored;
    and where the <SUBSET/EXCEPT/FOR qualification> is optional.

    Note that the location and scale parameters are optional.

Examples:
    LET A = RIGPPF(0.95,2,1)
    LET A = RIGPPF(P1,2,1)
    PLOT RIGPPF(P,2,1.5) FOR P = 0 0.01 0.99
Default:
    None
Synonyms:
    None
Related Commands:
    RIGCDF = Compute the reciprocal inverse Gaussian cumulative distribution function.
    RIGCHAZ = Compute the reciprocal inverse Gaussian cumulative hazard function.
    RIGHAZ = Compute the reciprocal inverse Gaussian hazard function.
    RIGPDF = Compute the reciprocal inverse Gaussian probability density function.
    RIGPPF = Compute the reciprocal inverse Gaussian percent point function.
    IGPDF = Compute the reciprocal inverse Gaussian probability density function.
    CHSPDF = Compute the chi-square probability density function.
    FPDF = Compute the F probability density function.
    NORPDF = Compute the normal probability density function.
    TPDF = Compute the t probability density function.
    WEIPDF = Compute the Weibull probability density function.
    WALPDF = Compute the Wald probability density function.
    FLPDF = Compute the fatigue life probability density function.
Reference:
    "Continuous Univariate Distributions--Volume 1", Second Edition, Johnson, Kotz, and Balakrishnan, Wiley, 1994, chapter 15.

    "Statistical Distributions", Third Edition, Evans, Hastings, and Peacock, Wiley, 2000, pp. 114-116.

Applications:
    Distributional Modeling
Implementation Date:
    1990/5: Original implementation
    2003/12: Modified to treat mu as a shape parameter instead of a location parameter
Program:
     
    X1LABEL Probability
    Y1LABEL X
    LABEL CASE ASIS
    X1LABEL DISPLACEMENT 12
    Y1LABEL DISPLACEMENT 12
    MULTIPLOT SCALE FACTOR 2
    MULTIPLOT 2 2
    MULTIPLOT CORNER COORDINATES 0 0 100 95
    TITLE GAMMA = 2, MU = 1
    PLOT RIGPPF(P,2,1) FOR P = 0.01  0.01  
    TITLE GAMMA = 5, MU = 1
    PLOT RIGPPF(P,5,1) FOR P = 0.01  0.01  
    TITLE GAMMA = 2, MU = 2
    PLOT RIGPPF(P,2,2) FOR P = 0.01  0.01  
    TITLE GAMMA = 5, MU = 2
    PLOT RIGPPF(P,5,2) FOR P = 0.01  0.01  
    END OF MULTIPLOT
    JUSTIFICATION CENTER
    MOVE 50 97
    CASE ASIS
    TEXT Reciprocal Inverse Gaussian Percent Point
        
    plot generated by sample program

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