Generate an HMT population

Share:

Description

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

Usage

1
HMT(N=5000, H=10)

Arguments

N

population size

H

number of strata

Details

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

Value

N x 3 matrix with columns:

strat

stratum ID

x

auxiliary variable x

y

analysis variable y

Author(s)

Richard Valliant, Jill A. Dever, Frauke Kreuter

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

1
2
3
    # generate HMT population with 1000 units and 5 strata and plot results
pop <- HMT(N=1000, H=5)
plot(pop[,"x"],pop[,"y"])

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.