[R] Using "survey" package with ACS PUMS

2014-09-30 Thread Michael . Laviolette
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

Re: [R] Using "survey" package with ACS PUMS

2014-09-30 Thread Anthony Damico
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

Re: [R] Using "survey" package with ACS PUMS

2014-09-30 Thread Michael . Laviolette
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,