fastflightcube {sampling} | R Documentation |
Execute the fast flight phase of the cube method (algorithm of Chauvet and Tillé, 2005, 2006). The data are sorted following the argument order. Inclusion probabilities equal to 0 or 1 are tolerated.
fastflightcube(X,pik,order=1,comment=TRUE)
X |
matrix of auxiliary variables on which the sample must be balanced. |
pik |
vector of inclusion probabilities. |
order |
1, the data are randomly arranged, 2, no change in data order, 3, the data are sorted in decreasing order. |
comment |
a comment is written during the execution if comment is equal to TRUE . |
Tillé, Y. (2006), Sampling Algorithms, Springer.
Chauvet, G. and Tillé, Y. (2006). A fast algorithm of balanced sampling. to appear in Computational Statistics.
Chauvet, G. and Tillé, Y. (2005). New SAS macros for balanced sampling. In INSEE, editor, Journées de Méthodologie Statistique, Paris.
Deville, J.-C. and Tillé, Y. (2004). Efficient balanced sampling: the cube method. Biometrika, 91, 893-912.
Deville, J.-C. and Tillé, Y. (2005). Variance approximation under balanced sampling. Journal of Statistical Planning and Inference, 128/2:411–425.
# Matrix of balancing variables X=cbind(c(1,1,1,1,1,1,1,1,1),c(1,2,3,4,5,6,7,8,9)) # Vector of inclusion probabilities. # The sample has the size equal to 3. pik=c(1/3,1/3,1/3,1/3,1/3,1/3,1/3,1/3,1/3) # pikstar is almost a balanced sample and is ready for the landing phase pikstar=fastflightcube(X,pik,order=1,comment=TRUE) pikstar