The Horvitz-Thompson estimator
Computes the Horvitz-Thompson estimator of the population total.
HTestimator(y,pik)
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. |
data(belgianmunicipalities) attach(belgianmunicipalities) # Computes the inclusion probabilities pik=inclusionprobabilities(Tot04,200) N=length(pik) n=sum(pik) # Defines the variable of interest y=TaxableIncome # Draws a Poisson sample of expected size 200 s=UPpoisson(pik) # Computes the Horvitz-Thompson estimator HTestimator(y[s==1],pik[s==1])
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.