binsmooth (version 0.2.2)

sb_sample: Random sample from splinebins distribution

Description

Draw a random sample of points from a smoothed distribution obtained using splinebins.

Usage

sb_sample(splinebinFit, n = 1)

Arguments

splinebinFit

A list as returned by splinebins.

n

A positive integer giving the sample size.

Value

A vector of random deviates. Returns NA if an inaccurate fit is detected, as indicated by fitWarn.

Details

The approximate inverse of the CDF calculated by splinebins is used to generate random values of the smoothed distribution.

References

Paul T. von Hippel, David J. Hunter, McKalie Drown. Better Estimates from Binned Income Data: Interpolated CDFs and Mean-Matching, Sociological Science, November 15, 2017. https://www.sociologicalscience.com/articles-v4-26-641/

Examples

# NOT RUN {
# 2005 ACS data from Cook County, Illinois
binedges <- c(10000,15000,20000,25000,30000,35000,40000,45000,
              50000,60000,75000,100000,125000,150000,200000,NA)
bincounts <- c(157532,97369,102673,100888,90835,94191,87688,90481,
               79816,153581,195430,240948,155139,94527,92166,103217)
splinefit <- splinebins(binedges, bincounts, 76091)
sb_sample(splinefit, 5)
hist(sb_sample(splinefit, 3000))
# }