Compute an approximate variance estimate using the approximation of
joint-inclusion probabilities proposed by Hájek (1964).
Estimator of class 2, it requires only first-order inclusion probabilities and
only for sample units.
Usage
var_Hajek(y, pik)
Value
a scalar, the estimated variance
Arguments
y
numeric vector of sample observations
pik
numeric vector of first-order inclusion probabilities for sample units