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

PLNPPF

Name:
    PLNPPF (LET)
Type:
    Library Function
Purpose:
    Compute the power-lognormal percent point function.
Description:
    The standard lognormal distribution has the following percent point function:

      G(f,sigma,p) = EXP(NORPPF(1-(1-f)**(1/p))/sigma)
 for 0 < f < 1,  sigma, p > 0

    where sigma is a shape parameter, p is a shape (power) parameter, and NORPPF is the percent point function of the standard normal distribution. The input value is a real number between 0 and 1 (since it corresponds to a probability).

    If p is 1, this distribution reduces to the lognormal distribution.

Syntax:
    LET <y> = PLNPPF(<x>,<p>,<s>)        <SUBSET/EXCEPT/FOR qualification>
    where <x> is a number, parameter, or variable in the range 0 to 1;
                  <p> is a positive number, parameter, or variable that specifies the power parameter;
                  <s> is an optional positive number, parameter, or variable that specifies the shape parameter;
                  <y> is a variable or a parameter (depending on what <x> is) where the computed power-lognormal ppf value is stored;
    and where the <SUBSET/EXCEPT/FOR qualification> is optional.

    If the <s> parameter is omitted, it defaults to 1.

Examples:
    LET A = PLNPPF(0.9,2,1)
    LET X2 = PLNPPF(X1,P,SD)
    LET X2 = PLNPPF(X1,0.5,3)
Note:
    The general power-lognormal percent point function can be expressed in terms of the standard power-lognormal percent point function as follows:

      G(f,sigma,p,loc,scale) = loc + scale*G(f,sigma,p)
 for 0 < f < 1,  sigma, p, scale > 0

    where loc and scale are the location and scales parameters, respectively.

Default:
    None
Synonyms:
    None
Related Commands:
    PLNCDF = Compute the power-lognormal cumulative distribution function.
    PLNPDF = Compute the power-lognormal probability density function.
    PLNHAZ = Compute the power-lognormal hazard function.
    PLNCHAZ = Compute the power-lognormal cumulative hazard function.
    PNRPDF = Compute the power-normal probability density function.
    LGNPDF = Compute the lognormal probability density function.
    HFNPDF = Compute the half-normal probability density function.
    NORPDF = Compute the normal probability density function.
Reference:
    "A Computer Program POWNOR for Fitting the Power-Normal and -Lognormal Models to Life or Strength Data from Specimens of Various Sizes", Nelson and Doganaksoy, NIST-IR 4760, March 1992.
Applications:
    Reliability
Implementation Date:
    1995/5
Program:
    LABEL CASE ASIS 
    Y1LABEL X 
    X1LABEL Probability 
    TITLE CASE ASIS 
    .  
    MULTIPLOT 2 2 
    MULTIPLOT CORNER COORDINATES 0 0 100 95 
    . 
    TITLE P = 0.5, SIGMA = 0.2, 0.4, 0.7, 1.0 
    PLOT PLNPPF(F,0.5,0.2) FOR F = 0.01 .01 0.99 AND 
    PLOT PLNPPF(F,0.5,0.4) FOR F = 0.01 .01 0.99 AND 
    PLOT PLNPPF(F,0.5,0.7) FOR F = 0.01 .01 0.99 AND 
    PLOT PLNPPF(F,0.5,1.0) FOR F = 0.01 .01 0.99 
    TITLE P = 1.0, SIGMA = 0.2, 0.4, 0.7, 1.0 
    PLOT PLNPPF(F,1.0,0.2) FOR F = 0.01 .01 0.99 AND 
    PLOT PLNPPF(F,1.0,0.4) FOR F = 0.01 .01 0.99 AND 
    PLOT PLNPPF(F,1.0,0.7) FOR F = 0.01 .01 0.99 AND 
    PLOT PLNPPF(F,1.0,1.0) FOR F = 0.01 .01 0.99 
    TITLE P = 5.0, SIGMA = 0.2, 0.4, 0.7, 1.0 
    PLOT PLNPPF(F,5.0,0.2) FOR F = 0.01 .01 0.99 AND 
    PLOT PLNPPF(F,5.0,0.4) FOR F = 0.01 .01 0.99 AND 
    PLOT PLNPPF(F,5.0,0.7) FOR F = 0.01 .01 0.99 AND 
    PLOT PLNPPF(F,5.0,1.0) FOR F = 0.01 .01 0.99 
    TITLE P = 20, SIGMA = 0.2, 0.4, 0.7, 1.0 
    PLOT PLNPPF(F,20,0.2) FOR F = 0.01 .01 0.99 AND 
    PLOT PLNPPF(F,20,0.4) FOR F = 0.01 .01 0.99 AND 
    PLOT PLNPPF(F,20,0.7) FOR F = 0.01 .01 0.99 AND 
    PLOT PLNPPF(F,20,1.0) FOR F = 0.01 .01 0.99 
    . 
    END OF MULTIPLOT 
    MOVE 50 97 
    JUSTIFICATION CENTER 
    TEXT Power Lognormal PPF's 
        
    plot generated by sample program

Date created: 11/13/2002
Last updated: 4/4/2003
Please email comments on this WWW page to alan.heckert@nist.gov.