landingRM | R Documentation |
This function performs the landing phase of the cube method using suppression of variables proposed by Chauvet and Tillé (2006).
landingRM(A, pikstar, EPS = 1e-07)
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 |
EPS |
epsilon 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.
Raphaël Jauslin raphael.jauslin@unine.ch
Chauvet, G. and Tillé, Y. (2006). A fast algorithm of balanced sampling. Computational Statistics, 21/1:53-62
fbs
, balstrat
.
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
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.