NAG Library Function Document
nag_rand_gen_discrete (g05tdc) generates a vector of pseudorandom integers from a discrete distribution with a given PDF (probability density function) or CDF (cumulative distribution function) .
||nag_rand_gen_discrete (Nag_ModeRNG mode,
const double p,
nag_rand_gen_discrete (g05tdc) generates a sequence of
, from a discrete distribution defined by information supplied in p
. This may either be the PDF or CDF of the distribution. A reference vector is first set up to contain the CDF of the distribution in its higher elements, followed by an index.
Setting up the reference vector and subsequent generation of variates can each be performed by separate calls to nag_rand_gen_discrete (g05tdc) or may be combined in a single call.
One of the initialization functions nag_rand_init_repeatable (g05kfc)
(for a repeatable sequence if computed sequentially) or nag_rand_init_nonrepeatable (g05kgc)
(for a non-repeatable sequence) must be called prior to the first call to nag_rand_gen_discrete (g05tdc).
Kendall M G and Stuart A (1969) The Advanced Theory of Statistics (Volume 1) (3rd Edition) Griffin
Knuth D E (1981) The Art of Computer Programming (Volume 2) (2nd Edition) Addison–Wesley
mode – Nag_ModeRNGInput
: a code for selecting the operation to be performed by the function.
- Set up reference vector only.
- Generate variates using reference vector set up in a prior call to nag_rand_gen_discrete (g05tdc).
- Set up reference vector and generate variates.
- Generate variates without using the reference vector.
, , or .
n – IntegerInput
On entry: , the number of pseudorandom numbers to be generated.
p[np] – const doubleInput
On entry: the PDF or CDF of the distribution.
- , for ;
- if , ;
- if , .
np – IntegerInput
the number of values supplied in p
defining the PDF or CDF of the discrete distribution.
ip1 – IntegerInput
On entry: the value of the variate, a whole number, to which the probability in corresponds.
itype – Nag_DiscreteDistribInput
: indicates the type of information contained in p
- p contains a probability distribution function (PDF).
- p contains a cumulative distribution function (CDF).
r[lr] – doubleCommunication Array
On entry: if , the reference vector from the previous call to nag_rand_gen_discrete (g05tdc).
On exit: the reference vector.
lr – IntegerInput
: the dimension of the array r
- if , approximately (for optimum efficiency in generating variates);
- otherwise .
- if or , ;
- if , lr should remain unchanged from the previous call to nag_rand_gen_discrete (g05tdc).
state – IntegerCommunication Array
the actual argument supplied must be the array state
supplied to the initialization functions nag_rand_init_repeatable (g05kfc)
or nag_rand_init_nonrepeatable (g05kgc)
On entry: contains information on the selected base generator and its current state.
On exit: contains updated information on the state of the generator.
x[n] – IntegerOutput
On exit: contains pseudorandom numbers from the specified discrete distribution.
fail – NagError *Input/Output
The NAG error argument (see Section 3.6
in the Essential Introduction).
6 Error Indicators and Warnings
Dynamic memory allocation failed.
On entry, argument had an illegal value.
On entry, lr
is too small when
, minimum length required
On entry, .
On entry, .
An internal error has occurred in this function. Check the function call and any array sizes. If the call is correct then please contact NAG
On entry, state
vector has been corrupted or not initialized.
The value of np
is not the same as when r
was set up
in a previous call.
Previous value of
Previous value of
On entry, at least one element of the vector p
is less than
or greater than
and the values of p
are not all in
stricly ascending order.
and the sum of the elements of p
do not equal one.
On entry, .
Constraint: if , .
On entry, some of the elements of the array r
have been corrupted or have not been initialized.
This example prints
pseudorandom variates from a discrete distribution whose PDF,
, is defined as follows:
The reference vector is set up and and the variates are generated by a single call to nag_rand_gen_discrete (g05tdc), after initialization by nag_rand_init_repeatable (g05kfc)
9.1 Program Text
Program Text (g05tdce.c)
9.2 Program Data
9.3 Program Results
Program Results (g05tdce.r)