lcubestratified: Stratified doubly balanced sampling with pooling of landing phases
Description
This is a fast implementation of stratified doubly balanced sampling. To have a fixed sample size, include the inclusion probabilities as a balancing variable in Xbal and make sure the inclusion probabilities sum to a positive integer (within each stratum). Euclidean distance is used in the Xspread space.
Usage
lcubestratified(prob,Xspread,Xbal,integerStrata)
Arguments
prob
vector of length N with inclusion probabilities
Xspread
matrix of (standardized) auxiliary variables of N rows and q columns
Xbal
matrix of balancing auxiliary variables of N rows and r columns
integerStrata
vector of length N with stratum number
Value
Returns a vector of length N with sampling indicators.
References
Chauvet, G. (2009). Stratified balanced sampling. Survey Methodology, 35, 115-119.
Grafstr<U+00F6>m, A. and Till<U+00E9>, Y. (2013). Doubly balanced spatial sampling with spreading and restitution of auxiliary totals. Environmetrics, 24(2), 120-131.
# NOT RUN {# Example 1N = 10;
n = 5;
p = rep(n/N,N);
Xspread = cbind(1:N);
strata = c(1,1,1,1,1,1,2,2,2,2);
indicators = lcubestratified(p,Xspread,cbind(p),strata);
s = (1:N)[indicators==1];
# }