The lower tail probability for the gamma distribution with parameters
and
,
, is defined by:
The mean of the distribution is
and its variance is
. The transformation
is applied to yield the following incomplete gamma function in normalized form,
This is then evaluated using
nag_incomplete_gamma (s14bac).
The input arrays to this function are designed to allow maximum flexibility in the supply of vector arguments by re-using elements of any arrays that are shorter than the total number of evaluations required. See
Section 2.6 in the g01 Chapter Introduction for further information.
- 1:
– IntegerInput
-
On entry: the length of the array
tail.
Constraint:
.
- 2:
– const Nag_TailProbabilityInput
-
On entry: indicates whether a lower or upper tail probability is required. For
, for
:
- The lower tail probability is returned, i.e., .
- The upper tail probability is returned, i.e., .
Constraint:
or , for .
- 3:
– IntegerInput
-
On entry: the length of the array
g.
Constraint:
.
- 4:
– const doubleInput
-
On entry: , the value of the gamma variate with , .
Constraint:
, for .
- 5:
– IntegerInput
-
On entry: the length of the array
a.
Constraint:
.
- 6:
– const doubleInput
-
On entry: the parameter of the gamma distribution with , .
Constraint:
, for .
- 7:
– IntegerInput
-
On entry: the length of the array
b.
Constraint:
.
- 8:
– const doubleInput
-
On entry: the parameter of the gamma distribution with , .
Constraint:
, for .
- 9:
– doubleOutput
-
Note: the dimension,
dim, of the array
p
must be at least
.
On exit: , the probabilities of the beta distribution.
- 10:
– IntegerOutput
-
Note: the dimension,
dim, of the array
ivalid
must be at least
.
On exit:
indicates any errors with the input arguments, with
- No error.
-
On entry, | invalid value supplied in tail when calculating . |
-
-
On entry, | , |
or | . |
- The solution did not converge in iterations, see nag_incomplete_gamma (s14bac) for details. The probability returned should be a reasonable approximation to the solution.
- 11:
– NagError *Input/Output
-
The NAG error argument (see
Section 3.7 in How to Use the NAG Library and its Documentation).
This example reads in values from a number of gamma distributions and computes the associated lower tail probabilities.