HTestimator {sampling}R Documentation

The Horvitz-Thompson estimator

Description

Computes the Horvitz-Thompson estimator of the population total.

Usage

HTestimator(y,pik)

Arguments

y vector of the variable of interest; its length is equal to n, the sample size.
pik vector of the first-order inclusion probabilities; its length is equal to n, the sample size.

See Also

UPtille

Examples

# Belgian municipalities data base
data(belgianmunicipalities)
# Computes the inclusion probabilities
pik=inclusionprobabilities(belgianmunicipalities$Tot04,200)
N=length(pik)
n=sum(pik)
# Defines the variable of interest
y=belgianmunicipalities$TaxableIncome
# Draws a Poisson sample of expected size 200
s=UPpoisson(pik)
# Computes the Horvitz-Thompson estimator
HTestimator(y[s==1],pik[s==1])

[Package sampling version 2.1 Index]