I'm trying to reproduce some results from the American Community Survey
PUMS data using the "survey" package. I'm using the one-year 2012 estimates
for New Hampshire
(http://www2.census.gov/acs2012_1yr/pums/csv_pnh.zip) and comparing to the
estimates for user verification from
http://www.census.go
hi michael, you probably need
options( "survey.replicates.mse" = TRUE )
i also think you don't want type = "Fay" and you do want scale = 4/80 and
rscales = rep( 1 , 80 ) as well, but what you have might be equivalent
(not sure)
regardless, this blog post details how to precisely replicate the
Setting the options did the trick. The code using "Fay" came from another
post, but this works:
options( "survey.replicates.mse" = TRUE)
pums_p.rep <- svrepdesign(repweights = pums_p[207:286],
weights = ~PWGTP,
combined.weights = TRUE,
3 matches
Mail list logo