Generates a simple pseudorandom sample from a finite population.
#include <imsls.h>
float *imsls_f_random_sample (int nrow, int nvar, float population[], int nsamp,..., 0)
The type double function is imsls_d_random_sample.
int nrow
(Input)
Number of rows of data in population.
int nvar
(Input)
Number of variables in the population and
in the sample.
float
population[] (Input)
nrow by nvar matrix containing the population to be sampled. If either
of the optional arguments IMSLS_FIRST_CALL or IMSLS_ADDITIONAL_CALL
are specified, then population contains a
different part of the population on each invocation, otherwise population contains
the entire population.
int nsamp
(Input)
The sample size desired.
nsamp by nvar matrix containing the sample. To release this space, use free.
#include <imsls.h>
float
*imsls_f_random_sample (int
nrow,
int nvar, float
population[],
int nsamp,
IMSLS_FIRST_CALL, int
**index,
int
*npop
IMSLS_FIRST_CALL_USER, int
index[],
int
*npop
IMSLS_ADDITIONAL_CALL, int
*index,
int
*npop,
float
*samp,
IMSLS_POPULATION_COL_DIM, int
population_col_dim,
IMSLS_RETURN_USER, int
samp[],
0)
IMSLS_FIRST_CALL, int **index,
int
*npop (Output)
This is the first invocation with this
data; additional calls to imsls_f_random_sample
may be made to add to the population. Additional calls should
be made using the optional argument IMSLS_ADDITIONAL_CALL
. Argument index is the address
of a pointer to an internally allocated array of
length nsamp containing the indices of the sample in the
population. Argument npop returns the number of items in the population.
If the population is input a few items at a time, the first call to imsls_f_random_sample should use IMSLS_FIRST_CALL, and subsequent calls should use IMSLS_ADDITIONAL_CALL. See example 2.
IMSLS_FIRST_CALL_USER, int index[], int *npop
(Output)
Storage for index is provided by
the user. See IMSLS_FIRST_CALL.
IMSLS_ADDITIONAL_CALL, int *index,
int
*npop,
float
*samp (Input/Output)
This is an additional invocation of
imsls_f_random_sample,
and updating for the subpopulation in population is
performed. Argument index is a pointer to
an array of length nsamp containing the indices of the sample in the population, as
returned using optional argument IMSLS_FIRST_CALL. Argument npop, also obtained
using optional argument IMSLS_FIRST_CALL, returns the number of items in the
population. It is not necessary to know the number of items in the
population in advance. npop is used to cumulate the population size and should not be
changed between calls to imsls_f_random_sample. Argument samp is a pointer to the array of size nsamp by nvar containing the sample. samp is the result of calling imsls_f_random_sample with optional argument IMSLS_FIRST_CALL. See Example 2
IMSLS_POPULATION_COL_DIM, int
population_col_dim (Input)
Column dimension of the matrix population.
Default:
x_col_dim = nvar
IMSLS_RETURN_USER, int samp[]
(Output)
User-supplied array of size nrow by nvar containing the sample. This option should not be
used if IMSLS_ADDITIONAL_CALL is used.
Routine imsls_f_random_sample generates a pseudorandom sample from a given population, without replacement, using an algorithm due to McLeod and Bellhouse (1983).
The first nsamp items in the population are included in the sample. Then, for each successive item from the population, a random item in the sample is replaced by that item from the population with probability equal to the sample size divided by the number of population items that have been encountered at that time.
In this example, imsls_f_random_sample is used to generate a sample of size 5 from a population stored in the matrix population.
#include <stdio.h>
#include <imsls.h>
int main()
{
int nrow = 176, nvar = 2, nsamp = 5;
float *population;
float *sample;
population = imsls_f_data_sets(2, 0);
imsls_random_seed_set(123457);
sample = imsls_f_random_sample(nrow, nvar, population, nsamp, 0);
imsls_f_write_matrix("The sample", nsamp, nvar, sample,
IMSLS_NO_ROW_LABELS,
IMSLS_NO_COL_LABELS,
0);
}
The sample
1764 36
1828 62
1923 6
1773 35
1769 106
Routine imsls_f_random_sample is now used to generate a sample of size 5 from the same population as in the example above except the data are input to RNSRS one observation at a time. This is the way imsls_f_random_sample may be used to sample from a file on disk or tape. Notice that the number of records need not be known in advance.
#include <stdio.h>
#include <imsls.h>
int main()
{
int i, nrow = 176, nvar = 2, nsamp = 5;
int *index, npop;
float *population;
float *sample;
population = imsls_f_data_sets(2, 0);
imsls_random_seed_set(123457);
sample = imsls_f_random_sample(1, 2, population, nsamp,
IMSLS_FIRST_CALL, &index, &npop,
0);
for (i = 1; i < 176; i++) {
imsls_f_random_sample(1, 2, &population[2*i], nsamp,
IMSLS_ADDITIONAL_CALL, index, &npop, sample,
0);
}
printf("The population size is %d\n", npop);
imsls_i_write_matrix("Indices of random sample", 5, 1, index, 0);
imsls_f_write_matrix("The sample", nsamp, nvar, sample,
IMSLS_NO_ROW_LABELS,
IMSLS_NO_COL_LABELS,
0);
}
The population size is 176
Indices of random sample
1 16
2 80
3 175
4 25
5 21
The sample
1764 36
1828 62
1923 6
1773 35
1769 106
|
Visual Numerics, Inc. PHONE: 713.784.3131 FAX:713.781.9260 |