![]() |
PAPPPFName:
![]()
with The cumulative distribution function is computed using the following recurrence relation (from page 379 of Johnson, Kemp, and Kotz)
![]() The percent point function is computed by computing the cumulative distribution function until the appropriate probability is obtained.
<SUBSET/EXCEPT/FOR qualification> where <x> is a variable, number, or parameter in the interval (0,1); <theta> is a positive number or parameter that specifies the first shape parameter; <p> is a positive number or parameter that specifies the second shape parameter; <y> is a variable or a parameter where the computed Polya-Aeppli ppf value is stored; and where the <SUBSET/EXCEPT/FOR qualification> is optional.
LET Y = PAPPPF(P,2,0.3) PLOT PAPPPF(P,2,0.3) FOR P = 0 0.01 0.99
Evans (1953), "Experimental Evidence Concerning Contagious Distributions in Ecology", Biometrika, 40, pp. 186-211. Johnson, Kotz, and Kemp (1992), "Univariate Discrete Distributions", Second Edition, Wiley, pp. 378-382.
title size 3 tic label size 3 label size 3 legend size 3 height 3 multiplot scale factor 1.5 x1label displacement 12 y1label displacement 17 . multiplot corner coordinates 0 0 100 95 multiplot scale factor 2 label case asis title case asis case asis tic offset units screen tic offset 3 3 title displacement 2 x1label Probability y1label X . xlimits 0 1 major xtic mark number 6 minor xtic mark number 3 . multiplot 2 2 . title Theta = 0.5, P = 0.5 plot papcdf(x,0.5,0.5) for x = 0 0.01 0.99 . title Theta = 1, P = 0.5 plot papcdf(x,1,0.5) for x = 0 0.01 0.99 . title Theta = 2.5, P = 0.5 plot papcdf(x,2.5,0.5) for x = 0 0.01 0.99 . title Theta = 5, P = 0.5 plot papcdf(x,5,0.5) for x = 0 0.01 0.99 . end of multiplot . justification center move 50 97 text Percent Point for Polya-Aeppli ![]()
Date created: 6/20/2006 |