landingphase: Landing phase of the cube method

In jlisic/BalancedSampling: Balanced and Spatially Balanced Sampling

Description

Landing is done by dropping balancing variables (from rightmost column).

Usage

landingphase(prob,probflight,Xbal)	

Arguments

prob	vector of length N with inclusion probabilities
probflight	vector of length N obtained from the flightphase
Xbal	matrix of balancing auxiliary variables of N rows and q columns

Value

Returns a vector of length N with inclusion indicators.

References

Deville, J. C. and Till<c3><a9>, Y. (2004). Efficient balanced sampling: the cube method. Biometrika, 91(4), 893-912.

Chauvet, G. and Till<c3><a9>, Y. (2006). A fast algorithm for balanced sampling. Computational Statistics, 21(1), 53-62.

Examples

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

## End(Not run)

jlisic/BalancedSampling documentation built on May 19, 2019, 12:46 p.m.