Good morning and I hope you are well. I am trying to analyse data from the COVID-19 survey and I have a question.
I am interested in home-schooling outcomes. When I run regressions with the SVY command in stata hundreds of strata are omitted from estimation (about 200-300 participants are therefore omitted). What I am doing for the moment is to treat the strata with single sampling units as certainty units. Do you think this is right? Are there any available corrections? Should I merge single units?
Updated by Alita Nandi over 1 year ago
- Status changed from New to Feedback
- % Done changed from 0 to 50
We advise users to use "scaled" option, as the "certainty" option means that these units do not contribut to the standard error. If you specify "scaled" then as per Stata Help, "the scaling factor comes from using the average of the variances from the strata with multiple sampling units for each stratum with one sampling unit."
But when you say hundreds of strata get dropped do you mean the statement in the output "# strata omitted because they contain no population members"? This will happen whichever singleunit option you use as this just refers to the fact that there are many strata where there are no observations used in your analyses so the strata get dropped.
Also, note that if you are subsetting the sample, then use "over()" rather than "if" with svy.
I am sorry, if I have misunderstood your problem. In that case you could you please attach this part of the syntax? If it is confidential, you don't need to specify the actual variables you have used, use depvar indepvar instead.