BalancedSampling (version 1.6.3)

lcubelandingphase: Landing phase for the local cube method

Description

Landing is done by dropping balancing variables (from rightmost column). Euclidean distance is used in the Xspread space.

Usage

lcubelandingphase(prob,probflight,Xspread,Xbal)

Value

Returns a vector of length N with indicators.

Arguments

prob

vector of length N with inclusion probabilities

probflight

vector of length N with probabilities from flightphase

Xspread

matrix of (standardized) auxiliary variables of N rows and q columns

Xbal

matrix of balancing auxiliary variables of N rows and r columns

References

Grafström, A. and Tillé, Y. (2013). Doubly balanced spatial sampling with spreading and restitution of auxiliary totals. Environmetrics, 24(2), 120-131.

Examples

if (FALSE) {
# Example 1
# Select sample
set.seed(12345);
N = 100; # population size
n = 10; # sample size
p = rep(n/N,N); # inclusion probabilities
X = cbind(runif(N),runif(N)); # matrix of auxiliary variables
pflight = lcubeflightphase(p,X,cbind(p,X)); # flight
indicators = lcubelandingphase(p,pflight,X,cbind(p,X)); # landing
# final sample
s = (1:N)[indicators==1];
}