Question

Status:Closed    Asked:Apr 14, 2017 - 11:39 AM

Should standard errors computed with ATUS replicate weights be expected to match BLS published results?

Hello again! I'm attempting to generate standard errors for an ATUS-X export using the RWT06 series of replicate weights. Using R, I've carefully reproduced the formula documented in the RWT06 description page (and in the ATUS reference documentation), however, my computed standard errors do not match those published by the BLS here:


https://data.bls.gov/cgi-bin/surveymo...


Should they match? If so, are you aware of a reference implementation I could compare to? If not, are there any published, bullet-proofed results from the ATUS-X data that I could try to replicate?


Thanks,

Chris

 
Do you have the same question? Follow this Question
 

Staff Answer

avatar

Joe_Grover

Staff

I was able to replicate the BLS table "Avg hrs per day - Watching TV - TUU10101AA01014236" using the rwt06 replicate weights from IPUMS ATUS and the BLS time use variable BLS_LEIS_TV divided by 60.

Is there a specific table you are unable to match?

 

Apr 17, 2017 - 01:47 PM

0
0
Report it

Answers

As it happens, that was exactly the series I was trying to replicate. I must just have a subtle error in my code. I'm going to take another crack at it. If you happen to have any reference code I could take a look at, that would be much appreciated.

 

Apr 18, 2017 - 04:47 PM

0
0
Report it

I initially replicated the table in Stata using the svyset command.

I was able to replicate the 2007 figuer in R as well with the following:
(assuming the dataframe is named atus2007)


atus2007$hrs = atus2007$BLS_LEIS_TV / 60


y = svrepdesign(

data = atus2007,

weights = ~WT06,

repweights = "RWT06_[1-9]+",

type = "JK1",

combined.weights = TRUE,

scale = .025

)


svymean(~hrs, design = y)

mean SE

hrs 2.6232 0.0295


 

Apr 20, 2017 - 12:27 PM

0
0
Report it

Okay, you're right that works perfectly. I didn't think I could do it that way because of this note in the RWT06 docs:


"The successive difference replication approach (SDR) is different from other methods for creating replicate weights such as balanced repeated replication (BRR) and jackknife estimation."


The R survey library doesn't support SDR, so I thought I had to compute these from scratch. However, after seeing your code I did some more Googling and found this sentence in the IPUMS-CPS replicate weights docs:


"Correspondence with StataCorp statisticians and IPUMS testing revealed that successive difference replicate weights can be treated as jackknife replicate weights if the options are specified correctly."


This still doesn't explain why my (really pretty simple) code doesn't give the correct answers, but I'm not going to worry about that for the moment. I might suggest adding a mention of that to the ATUS replicate weights docs.


Thanks for the help, and for sharing the example!

 

Apr 20, 2017 - 04:28 PM

0
0
Report it

Not an answer but instead, much more basic questions ...


I am using STATA 13 and creating a simple cross-tabulation of vetstat (veteran status) and empstat (employment status). I want the weighted observation counts of veterans to approximate 20 million using CPS/ASEC(March) for 2015.


a) What weight variable(s) should I use (repwtp, wtsupp, or repwtp1-repwtp160)?

b) If I want to compute standard errors and confidence intervals, what weight variable(s) should I use?


b) Exactly how should I specify the tabulate command?


I am very much a STATA rookie. So, any detailed help that you can provide will be much appreciated.


George (where georgegsheldon@gmail.com can also be used)

 

Aug 14, 2017 - 08:14 AM

0
0
Report it

Hi George,


For generating standard deviations using Stata and replicate weights, I would recommend starting with this documentation from the IPUMS USA website. While that page discusses using the IPUMS USA data, the syntax is basically the same You can then use the svy command to correctly apply the weights. Not all commands are available with svy in Stata. See this instructional page for an applied example of svy usage.


I hope this helps!

 

Aug 14, 2017 - 01:35 PM

0
0
Report it

OTHER QUESTIONS NEEDING ANSWERS

The highest one percent incwage earners are combined into the incwage"Top Cod...
To compute unemployment rates from Ipums-I, should I use perwt as pweights or...
Why is the income wage variable incwage-capped? and how is the cap determined...
USA IPUMS 1980 ancestry variable -- mapping to country names
Did the version of database change since 2011?
Are DATANUM, SERIAL, and PERNUM the same in single-year and 5-year ACS sample...
Login   |   Register


Share |