ffphase: Fast flight phase of the cube method

In RJauslin/StratifiedSampling: Different Methods for Stratified Sampling

Fast flight phase of the cube method

Description

This function computes the flight phase of the cube method proposed by Chauvet and Tillé (2006).

Usage

ffphase(Xbal, prob, order = TRUE)

Arguments

Xbal	A matrix of size (N x p) of auxiliary variables on which the sample must be balanced.
prob	A vector of inclusion probabilities.
order	if the units are reordered, Default TRUE.

Details

This function implements the method proposed by (Chauvet and Tillé 2006). It recursively transforms the vector of inclusion probabilities pik into a sample that respects the balancing equations. The algorithm stops when the null space of the sub-matrix B is empty. For more information see (Chauvet and Tillé 2006).

Value

Updated vector of pik that contains 0 and 1 for unit that are rejected or selected.

See Also

fastflightphase, cube.

Examples

N <- 100
n <- 10
p <- 4
pik <- rep(n/N,N)
X <- cbind(pik,matrix(rgamma(N*p,4,25),ncol= p))

pikstar <- ffphase(X,pik) 
t(X/pik)%*%pikstar
t(X/pik)%*%pik
pikstar

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