Compute components of relvariance for a sample design where primary sampling units (PSUs) and elements are selected via srs. The input is an entire sampling frame.
1  BW2stageSRS(X, psuID)

X 
data vector; length is the number of elements in the population. 
psuID 
vector of PSU identification numbers. This vector must be as long as 
BW2stageSRS
computes the between and within population relvariance components
appropriate for a twostage sample in which PSUs are selected via srs (either with or without replacement). Elements within PSUs are assumed to be selected by srswor. The same number of elements is assumed to be selected within each sample PSU. The function requires that an entire frame of PSUs and elements be input.
(Use BW2stagePPSe
if only a sample of PSUs and elements are available.)
List with values:
B2 
between PSU unit relvariance 
W2 
within PSU unit relvariance 
unit relvar 
unit relvariance for population 
B2+W2 
B^2 + W^2 
k 
ratio of B^2 + W^2 to unit relvariance 
delta full 
intraclass correlation estimated as B^2/(B^2 + W^2) 
Richard Valliant, Jill A. Dever, Frauke Kreuter
Cochran, W.G. (1977, chap. 11). Sampling Techniques. New York: John Wiley & Sons.
Valliant, R., Dever, J., Kreuter, F. (2013, sect. 9.2.1). Practical Tools for Designing and Weighting Survey Samples. New York: Springer.
BW2stagePPS
, BW2stagePPSe
, BW3stagePPS
, BW3stagePPSe
1 2 3 4 5 6  data(MDarea.pop)
MDsub < MDarea.pop[1:100000,]
# psu's are defined by PSU variable
BW2stageSRS(abs(MDsub$Hispanic2), psuID=MDsub$PSU)
# psu's are defined by SSU variable
BW2stageSRS(abs(MDsub$Hispanic2), psuID=MDsub$SSU)

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
Please suggest features or report bugs with the GitHub issue tracker.
All documentation is copyright its authors; we didn't write any of that.