lcubelandingphase: Landing phase for the local cube method

In jlisic/BalancedSampling: Balanced and Spatially Balanced Sampling

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)	

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

Value

Returns a vector of length N with indicators.

References

Grafstr<c3><b6>m, A. and Till<c3><a9>, Y. (2013). Doubly balanced spatial sampling with spreading and restitution of auxiliary totals. Environmetrics, 24(2), 120-131.

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(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];

## End(Not run)

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