StratifiedSampling (version 0.4.2)

balstrat: Balanced Stratification

Description

Select a stratified balanced sample. The function is similar to balancedstratification of the package sampling.

Usage

balstrat(X, strata, pik, rand = TRUE, landing = TRUE)

Value

A vector with elements equal to 0 or 1. The value 1 indicates that the unit is selected while the value 0 is for rejected units.

Arguments

X

A matrix of size (\(N\) x \(p\)) of auxiliary variables on which the sample must be balanced.

strata

A vector of integers that specifies the stratification.

pik

A vector of inclusion probabilities.

rand

if TRUE, the data are randomly arranged. Default TRUE

landing

if TRUE, landing by linear programming otherwise supression of variables. Default TRUE

Details

The function implements the method proposed by Chauvet (2009). Firstly, a flight phase is performed on each strata. Secondly, a flight phase is applied on the whole population by aggregating the strata. Finally, a landing phase is applied by suppression of variables.

References

Chauvet, G. (2009). Stratified balanced sampling. Survey Methodology, 35:115-119.

See Also

ffphase, landingRM

Examples

N <- 100
n <- 10
p <- 4
X <- matrix(rgamma(N*p,4,25),ncol = p)
strata <- as.matrix(rep(1:n,each = N/n))
pik <- rep(n/N,N)

s <- balstrat(X,strata,pik)

t(X/pik)%*%s
t(X/pik)%*%pik

Xcat <- disj(strata)

t(Xcat)%*%s
t(Xcat)%*%pik