Sequential Random Sampling

If you specify the METHOD=SEQ option and do not include a SIZE statement, PROC SURVEYSELECT uses the equal probability version of Chromy’s method for sequential random sampling. This method selects units sequentially with equal probability and without replacement. For more information, see Chromy (1979) and Williams and Chromy (1980). For information about Chromy’s PPS selection method, see the section PPS Sequential Sampling.

Sequential random sampling controls the distribution of the sample by spreading it throughout the sampling frame or stratum, thus providing implicit stratification according to the order of units in the frame or stratum. You can use the CONTROL statement to sort the input data set by the CONTROL variables before sample selection. If you also use a STRATA statement, PROC SURVEYSELECT sorts by the CONTROL variables within strata. By default (or if you specify the SORT=SERP option), the procedure uses hierarchic serpentine ordering for sorting. If you specify the SORT=NEST option, the procedure uses nested sorting. See the section Sorting by CONTROL Variables for descriptions of serpentine and nested sorting. If you do not specify a CONTROL statement, PROC SURVEYSELECT applies sequential selection to the observations in the order in which they appear in the input data set.

Following Chromy’s method of sequential selection, PROC SURVEYSELECT randomly chooses a starting unit from the entire stratum (or frame, if the design is not stratified). With this unit as the first one, the procedure treats the stratum units as a closed loop. This is done so that all pairwise (joint) selection probabilities are positive and an unbiased variance estimator can be obtained. The procedure numbers units sequentially from the random start to the end of the stratum and then continues from the beginning of the stratum until all units are numbered.

Beginning with the randomly chosen starting unit, PROC SURVEYSELECT accumulates the expected number of selections (hits), where the expected number of selections $\mr{E}(S_{hi})$ is $n_ h / N_ h$ for all units i in stratum h. The procedure computes

\[  I_{hi} = \mr{Int} \Bigl ( \sum _{j=1}^ i \mr{E}(S_{hj}) \Bigr ) = \mr{Int} ( i n_ h / N_ h )  \]
\[  F_{hi} = \mr{Frac} \Bigl ( \sum _{j=1}^ i \mr{E}(S_{hj}) \Bigr ) = \mr{Frac} ( i n_ h / N_ h)  \]

where $\mr{Int}(\cdot )$ denotes the integer part of the number, and $\mr{Frac}(\cdot )$ denotes the fractional part.

Considering each unit sequentially, Chromy’s method determines whether unit i is selected by comparing the total number of selections for the first (i – 1) units,

\[  T_{h(i-1)} = \sum _{j=1}^{i-1} S_{hj}  \]

with the value of $I_{h(i-1)}$.

If $T_{h(i-1)} = I_{h(i-1)}$, Chromy’s method determines whether or not unit i is selected as follows. If $F_{hi} = 0$ or $F_{h(i-1)} > F_{hi}$, then unit i is selected with certainty. Otherwise, unit i is selected with probability

\[  (F_{hi} - F_{h(i-1)}) / (1 - F_{h(i-1)})  \]

If $T_{h(i-1)} = (I_{h(i-1)} + 1)$, Chromy’s method determines whether or not unit i is selected as follows. If $F_{hi} = 0$ or $F_{hi} > F_{h(i-1)}$, then the unit is not selected. Otherwise, unit i is selected with probability

\[  F_{hi} / F_{h(i-1)}  \]