The balanced repeated replication (BRR) method requires that the full sample be drawn by using a stratified sample design with two primary sampling units (PSUs) per stratum. Let H be the total number of strata. The total number of replicates R is the smallest multiple of 4 that is greater than H. However, if you prefer a larger number of replicates, you can specify the REPS=number option. If a number number Hadamard matrix cannot be constructed, the number of replicates is increased until a Hadamard matrix becomes available.
Each replicate is obtained by deleting one PSU per stratum according to the corresponding Hadamard matrix and adjusting the original weights for the remaining PSUs. The new weights are called replicate weights.
Replicates are constructed by using the first H columns of the Hadamard matrix. The rth () replicate is drawn from the full sample according to the rth row of the Hadamard matrix as follows:
If the element of the Hadamard matrix is 1, then the first PSU of stratum h is included in the rth replicate and the second PSU of stratum h is excluded.
If the element of the Hadamard matrix is –1, then the second PSU of stratum h is included in the rth replicate and the first PSU of stratum h is excluded.
Note that the “first” and “second” PSUs are determined by data order in the input data set. Thus, if you reorder the data set and perform the same analysis by using BRR method, you might get slightly different results, because the contents in each replicate sample might change.
The replicate weights of the remaining PSUs in each halfsample are then doubled to their original weights. For more details about the BRR method, see Wolter (2007) and Lohr (2010).
By default, an appropriate Hadamard matrix is generated automatically to create the replicates. You can request that the Hadamard matrix be displayed by specifying the VARMETHOD=BRR(PRINTH) methodoption. If you provide a Hadamard matrix by specifying the VARMETHOD=BRR(HADAMARD=) methodoption, then the replicates are generated according to the provided Hadamard matrix.
You can use the VARMETHOD=BRR(OUTWEIGHTS=) methodoption to save the replicate weights into a SAS data set.
Suppose that is a population parameter of interest. Let be the estimate from the full sample for . Let be the estimate from the rth replicate subsample by using replicate weights. PROC SURVEYMEANS estimates the variance of by

with H degrees of freedom, where H is the number of strata.
If a parameter cannot be computed from one or more replicates, then the variance estimate is computed by using those replicates from which the parameter can be estimated. For example, suppose the parameter is a ratio. If a replicate r contains observations such that the denominator of the ratio is zero, then the ratio cannot be computed from replicate r. In this case, the BRR variance estimate is computed as

where the summation is over the replicates where the parameter can be computed, and is the number of those replicates.