We present an alternative way of evaluating the following integral, which takes advantage of the connection between Gamma distribution and Poisson process.
The integrand in is the probability density function (pdf) of a Gamma distribution with parameters (scale parameter) and (shape parameter) where and is a positive integer. Thus the integral is the right tail of a Gamma distribution. If the pdf is to model the lifetime of a biological life, the integral is the survival function (i.e. the probability that the life will survice beyond time ).
The integral can be evaluated using the method of integration by parts and the following provides the alternative method.
The goal is not to memorize , but rather to understand the thought process behind the formula. The thought process involves identifying the Poisson process that is associated with the Gamma pdf in the integral. The integral is equivalent to a probability calculation in the Poisson process. The right-hand-side of is the probability that there are less than changes in the Poisson process from time zero to time . The left-hand-side of is the probability that the change in the Poisson process occurrs after time . The description of the algorithm is followed by examples and a concluding remark.
- From the Gamma pdf in the left-hand-side of , pick out the three numbers , and . The first two are the parameters of the Gamma pdf and the third one is the lower limit of the integral.
- Consider the Poisson process with parameter . One thing to keep in mind here is that in this Poisson process, the number of changes in a unit time interval is modeled by and the number of changes in a time interval of length is modeled by
- The right-hand-side of is the probability that there are less than changes in the time interval (or at most changes).
The three parameters in the algorithm are , and . Consider the Poisson process with parameter . We need to find the probability that there are less than 2 changes in the time interval . The Poisson probability is . Thus we have:
The three parameters in the algorithm are , and . We are interested in the Poisson process with parameter and in finding the probability that there are less than 4 changes in time interval . The Poisson probability is . The integral, after adjusting for a multiplicative constant, is the probability that the fourth change occurs after time 1.5.
The reason that the summation in the right-hand-side of works as an answer to the integral on the left-hand-side is that the summation is how the Gamma pdf is derived (where the shape parameter is a positive integer). Let be the time until the occurrence of the change in the Poisson process that has parameter . Then the right-hand-side of is . The distribution function is then . Taking derivative of produces the integrand in the left-hand-side of . Thus the alternative approach indicated here is simply going back to the root of the Gamma pdf.