Uses the Midzuno's method to select a sample of units (unequal probabilities, without replacement, fixed sample size).
Usage
UPmidzuno(pik,eps=1e-6)
Arguments
pik
vector of the inclusion probabilities.
eps
the control value, by default equal to 1e-6.
Value
Returns a vector (with elements 0 and 1) of size N, the population size.
Each element k of this vector indicates the status of unit k
(1, unit k is selected in the sample; 0, otherwise). The value 'eps' is used to
control pik (pik>eps & pik < 1-eps).
References
Midzuno, H. (1952), On the sampling system with probability proportional to sum of size.
Annals of the Institute of Statistical Mathematics, 3:99-107.
Deville, J.-C. and Till<e9>, Y. (1998),
Unequal probability sampling without replacement through a splitting method,
Biometrika, 85:89-101.
# NOT RUN {#define the prescribed inclusion probabilitiespik=c(0.2,0.7,0.8,0.5,0.4,0.4)
#select a samples=UPmidzuno(pik)
#the sample is(1:length(pik))[s==1]
# }