The deviate,
associated with the lower tail probability,
, for the Normal distribution is defined as the solution to
where
The method used is an extension of that of
Wichura (1988).
is first replaced by
.
(a) |
If , is computed by a rational Chebyshev approximation
where and , are polynomials of degree . |
(b) |
If , is computed by a rational Chebyshev approximation
where and , are polynomials of degree . |
(c) |
If , is computed as
where and , are polynomials of degree . |
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.
Wichura (1988) Algorithm AS 241: the percentage points of the Normal distribution Appl. Statist. 37 477–484
- 1:
– IntegerInput
-
On entry: the length of the array
tail.
Constraint:
.
- 2:
– const Nag_TailProbabilityInput
-
On entry: indicates which tail the supplied probabilities represent. Letting
denote a variate from a standard Normal distribution, and
, then for
, for
:
- The lower tail probability, i.e., .
- The upper tail probability, i.e., .
- The two tail (confidence interval) probability, i.e., .
- The two tail (significance level) probability, i.e., .
Constraint:
, , or , for .
- 3:
– IntegerInput
-
On entry: the length of the array
p.
Constraint:
.
- 4:
– const doubleInput
-
On entry:
, the probabilities for the Normal distribution as defined by
tail with
,
.
Constraint:
, for .
- 5:
– IntegerInput
-
On entry: the length of the array
xmu.
Constraint:
.
- 6:
– const doubleInput
-
On entry: , the means with , .
- 7:
– IntegerInput
-
On entry: the length of the array
xstd.
Constraint:
.
- 8:
– const doubleInput
-
On entry: , the standard deviations with , .
Constraint:
, for .
- 9:
– doubleOutput
-
Note: the dimension,
dim, of the array
x
must be at least
.
On exit: , the deviates for the Normal 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 | . |
-
- 11:
– NagError *Input/Output
-
The NAG error argument (see
Section 3.7 in How to Use the NAG Library and its Documentation).
None.