lcubestratified {BalancedSampling} | R Documentation |
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.
lcubestratified(prob,Xspread,Xbal,integerStrata)
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 |
Returns a vector of length N with sampling indicators.
Chauvet, G. (2009). Stratified balanced sampling. Survey Methodology, 35, 115-119.
Grafström, A. and Tillé, Y. (2013). Doubly balanced spatial sampling with spreading and restitution of auxiliary totals. Environmetrics, 24(2), 120-131.
## Not run: # Example 1 N = 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]; ## End(Not run)