Hajekestimator: The Hajek estimator

Description Usage Arguments See Also Examples

View source: R/Hajekestimator.r

Description

Computes the H<e1>jek estimator of the population total or population mean.

Usage

1
Hajekestimator(y,pik,N=NULL,type=c("total","mean"))

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.

N

population size; N is only used for the total estimator; for the mean estimator its value is NULL.

type

the estimator type: total or mean.

See Also

HTestimator

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
# 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 Hajek estimator of the population mean
Hajekestimator(y[s==1],pik[s==1],type="mean")
# Computes the Hajek estimator of the population total
Hajekestimator(y[s==1],pik[s==1],N=N,type="total")

sampling documentation built on May 30, 2017, 1:14 a.m.

Search within the sampling package
Search all R packages, documentation and source code