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

PLNCHAZ

Name:
    PLNCHAZ (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)=-LOG[(NORCDF(-LOG(x)/sigma))**p) for x, p, sigma > 0

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

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

Syntax:
    LET <y> = PLNCHAZ(<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 cumulative 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 = PLNCHAZ(3,2,1)
    LET X2 = PLNCHAZ(X1,POW,SD)
    LET X2 = PLNCHAZ(X1,0.5,0.5)
Note:
    The general power-lognormal cumulative hazard function can be expressed in terms of the standard power-lognormal cumulative hazard function as follows:

      H(x,sigma,p,loc,scale)=H((x-loc)/scale,sigma,p)
 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.
    PLNHAZ = Compute the power-lognormal 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 Cumulative 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 PLNCHAZ(X,0.5,0.2) FOR X = 0.01 .01 3 AND 
    PLOT PLNCHAZ(X,0.5,0.4) FOR X = 0.01 .01 3 AND 
    PLOT PLNCHAZ(X,0.5,0.7) FOR X = 0.01 .01 3 AND 
    PLOT PLNCHAZ(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 PLNCHAZ(X,1.0,0.2) FOR X = 0.01 .01 3 AND 
    PLOT PLNCHAZ(X,1.0,0.4) FOR X = 0.01 .01 3 AND 
    PLOT PLNCHAZ(X,1.0,0.7) FOR X = 0.01 .01 3 AND 
    PLOT PLNCHAZ(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 PLNCHAZ(X,5.0,0.2) FOR X = 0.01 .01 2 AND 
    PLOT PLNCHAZ(X,5.0,0.4) FOR X = 0.01 .01 2 AND 
    PLOT PLNCHAZ(X,5.0,0.7) FOR X = 0.01 .01 2 AND 
    PLOT PLNCHAZ(X,5.0,1.0) FOR X = 0.01 .01 2 
    TITLE P = 20, SIGMA = 0.2, 0.4, 0.7, 1.0 
    PLOT PLNCHAZ(X,20,0.2) FOR X = 0.01 .01 1 AND 
    PLOT PLNCHAZ(X,20,0.4) FOR X = 0.01 .01 1 AND 
    PLOT PLNCHAZ(X,20,0.7) FOR X = 0.01 .01 1 AND 
    PLOT PLNCHAZ(X,20,1.0) FOR X = 0.01 .01 1 
    . 
    END OF MULTIPLOT 
    MOVE 50 97 
    JUSTIFICATION CENTER 
    TEXT Power Lognormal Cumulative 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.