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

PLNHAZ

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

      h(x,sigma,p)=p*(1/(x*sigma))*NORPDF(ln(x)/sigma)/NORCDF(-ln(x)/sigma)
 for x, p, sigma > 0

    where sigma is a shape parameter, p is a shape (power) parameter, and NORCDF and NORPDF are the cumulative distribution function and the probability density function for the standard normal distribution respectively.

    The hazard function for the power-lognormal distribution is a multiple (p) of the hazard function of the lognormal distribution.

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

Syntax:
    LET <y> = PLNHAZ(<x>,<p>,<s>)        <SUBSET/EXCEPT/FOR qualification>
    where <x> is a non-negative number, parameter, or variable;
                  <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 <y1> is) where the computed power-lognormal hazard 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 = PLNHAZ(3,2,1)
    LET X2 = PLNHAZ(X1,POW,SD)
    LET X2 = PLNHAZ(X1,0.5,0.5)
Note:
    The general power-lognormal hazard function can be expressed in terms of the standard power-lognormal hazard function as follows:

      h(x,sigma,p,loc,scale)=h((x-loc)/scale,sigma,p)/scale
 for x > loc,  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.
    PLNPPF = Compute the power-lognormal percent point 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:
    1998/4
Program:
    LABEL CASE ASIS 
    X1LABEL X 
    Y1LABEL Hazard 
    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 PLNHAZ(X,0.5,0.2) FOR X = 0.01 .01 3 AND 
    PLOT PLNHAZ(X,0.5,0.4) FOR X = 0.01 .01 3 AND 
    PLOT PLNHAZ(X,0.5,0.7) FOR X = 0.01 .01 3 AND 
    PLOT PLNHAZ(X,0.5,1.0) FOR X = 0.01 .01 3 
    TITLE P = 1.0, SIGMA = 0.2, 0.4, 0.7, 1.0 
    PLOT PLNHAZ(X,1.0,0.2) FOR X = 0.01 .01 3 AND 
    PLOT PLNHAZ(X,1.0,0.4) FOR X = 0.01 .01 3 AND 
    PLOT PLNHAZ(X,1.0,0.7) FOR X = 0.01 .01 3 AND 
    PLOT PLNHAZ(X,1.0,1.0) FOR X = 0.01 .01 3 
    TITLE P = 5.0, SIGMA = 0.2, 0.4, 0.7, 1.0 
    PLOT PLNHAZ(X,5.0,0.2) FOR X = 0.01 .01 2 AND 
    PLOT PLNHAZ(X,5.0,0.4) FOR X = 0.01 .01 2 AND 
    PLOT PLNHAZ(X,5.0,0.7) FOR X = 0.01 .01 2 AND 
    PLOT PLNHAZ(X,5.0,1.0) FOR X = 0.01 .01 2 
    TITLE P = 20, SIGMA = 0.2, 0.4, 0.7, 1.0 
    PLOT PLNHAZ(X,20,0.2) FOR X = 0.01 .01 1 AND 
    PLOT PLNHAZ(X,20,0.4) FOR X = 0.01 .01 1 AND 
    PLOT PLNHAZ(X,20,0.7) FOR X = 0.01 .01 1 AND 
    PLOT PLNHAZ(X,20,1.0) FOR X = 0.01 .01 1 
    . 
    END OF MULTIPLOT 
    MOVE 50 97 
    JUSTIFICATION CENTER 
    TEXT Power Lognormal Hazard Functions 
        
    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.