PracTools (version 1.2.5)

HMT: Generate an HMT population

Description

Generate a population that follows the model in Hansen, Madow, and Tepping (1983)

Usage

HMT(N=5000, H=10)

Arguments

N

population size

H

number of strata

Value

N x 3 matrix with columns:

strat

stratum ID

x

auxiliary variable x

y

analysis variable y

Details

HMT generates a population based on the model: \(E(y)= \alpha + \beta x\), \(Var(y)=\sigma^2x^{3/2}\). Both x and y have gamma distributions. Strata are formed to have approximately the same total of \(x\).

References

Hansen,M.H., Madow,W.G., and Tepping,B.J. (1983). An evaluation of model-dependent and probability sampling inferences in sample surveys. Journal of the American Statistical Association, 78, 776-793.

Examples

# NOT RUN {
    # generate HMT population with 1000 units and 5 strata and plot results
pop <- HMT(N=1000, H=5)
plot(pop[,"x"],pop[,"y"])
# }