landingRM: Landing by suppression of variables
In RJauslin/StratifiedSampling: Different Methods for Stratified Sampling

landingRM

R Documentation

Landing by suppression of variables

Description

This function performs the landing phase of the cube method using suppression of variables proposed by Chauvet and Tillé (2006).

Usage

landingRM(A, pikstar, EPS = 1e-07)

Arguments

`A`	matrix of auxiliary variables on which the sample must be balanced. (The matrix should be divided by the original inclusion probabilities.)
`pikstar`	vector of updated inclusion probabilities by the flight phase. See `ffphase`
`EPS`	epsilon value

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. and Tillé, Y. (2006). A fast algorithm of balanced sampling. Computational Statistics, 21/1:53-62

Examples

N <- 1000
n <- 10
p <- 4
pik <- rep(n/N,N)
X <- cbind(pik,matrix(rgamma(N*p,4,25),ncol= p))
pikstar <- ffphase(X,pik) 
s <- landingRM(X/pik*pikstar,pikstar)
sum(s)
t(X/pik)%*%pik
t(X/pik)%*%pikstar
t(X/pik)%*%s

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