balstrat: Balanced Stratification

In RJauslin/StratifiedSampling: Different Methods for Stratified Sampling

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)

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.

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.

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

RJauslin/StratifiedSampling documentation built on Feb. 5, 2025, 4:18 a.m.