The DHS Program User Forum
Discussions regarding The DHS Program data and results
Home » Data » Weighting data » Weight in Indonesia DHS 2002 and 2007 (svy commands did not work for poisson and logistic)
Weight in Indonesia DHS 2002 and 2007 [message #18589] Mon, 06 January 2020 19:58 Go to next message
lhuri is currently offline  lhuri
Messages: 16
Registered: December 2014
Location: Indonesia
Member
Dear DHS Team,

I'm wondering why svy commands do not work for Poisson and logistic regressions in Indonesia DHS 2002 and 2007.
I use these to apply weights

gen wgt = v005/1000000
svyset[pw=wgt], psu(v021) strata (v022)

It works for commands like svy: tab var, but when I use svy: poisson var1 i.var2, irr, it cannot provide CI. There's a note below the table saying "Note: Missing standard errors because of stratum with single sampling unit".

I use the same commands for DHS 2012 and 2017, and there's no problem. Is it because there are too many strata in 2002 and 2007 dataset? For comparison, IDHS 2017 has only 67 strata while 2002 has 680. Could you explain why this happens and what should I do?
Thank you.

Re: Weight in Indonesia DHS 2002 and 2007 [message #18605 is a reply to message #18589] Fri, 10 January 2020 09:03 Go to previous messageGo to next message
Bridgette-DHS is currently offline  Bridgette-DHS
Messages: 3016
Registered: February 2013
Senior Member

Following is a response from DHS Research & Data Analysis Director, Tom Pullum:

You need to add a "singleset" option to svyset. For example, use

svyset [pw=wgt], psu(v021) strata (v022) singleunit(centered)

The word in parentheses can be "centered" or "scaled" or "certainty". The results are virtually indistinguishable, whichever of the three you use.
Re: Weight in Indonesia DHS 2002 and 2007 [message #18614 is a reply to message #18605] Sat, 11 January 2020 06:31 Go to previous messageGo to next message
lhuri is currently offline  lhuri
Messages: 16
Registered: December 2014
Location: Indonesia
Member
Thanks a lot Bridgette and Tom,

It works now.

On a related note, I just want to make sure about the stratification used to weight the data in 2002, 2007, 2012, and 2017 Indonesia DHS.

According to the DHS manual ( https://dhsprogram.com/Data/Guide-to-DHS-Statistics/Analyzin g_DHS_Data.htm), it's said that "Typically, DHS samples are stratified by geographic region and by urban/rural areas within each region." I could follow this concept for 2012 and 2017 because, in these datasets, v022 consists of the number of all urban-rural area for each province in the survey year. For example, v022 in 2017 dataset contains 67 strata which represent urban and rural area in each of 34 provinces minus rural Jakarta, because Jakarta is the capital city and entirely considered as urban area). This also matches the explanation in the 2017 report "The sample of the 2017 IDHS was stratified by province and by urban and rural areas, and implicitly stratified by welfare concentration" (appendix B, p351).

However, in the 2002 and 2007 dataset, v022 shows that the number of strata is 680 and 751, respectively. I'm aware that there was a change of province number in Indonesia over time (30 provinces in 2002 survey and 33 in 2007 survey) and the appendix B in 2002 report (p267) does mention additional info about the five districts in East Java covered in the Safe Motherhood Project (SMP). However, I still don't quite understand where the 680 strata came from.

So there are several things I need to ask/clarify:

1. Does v022 in each survey year represent the number of strata by provinces and urban-rural residence?

2. If the sample of each survey year is stratified by provinces and urban-rural residence, why do 2002 and 2007 datasets have hundreds of sample? What do these numbers represent?

3. Why do v022 and v023 have the same unique values in the 2017 and 2012 datasets, but not in 2007 and 2002 datasets? For example, in 2017, v022 and v023 both have 67 unique values. Meanwhile, in the 2002, v022 (which is labeled as number of stratum) has 680 values and v023 (labeled as sample domain) has 26 values.

4. I've been following the guidelines including the Youtube tutorial and using "gen wgt = v005/1000000" followed by "svyset[pw=wgt], psu(v021) strata (v022)" to apply weights in the 2012 and 2017 dataset. To run analyses, I simply use svy commands like "svy: tab var", "svy: poisson var1 i.var2 i.var3, irr", etc. I haven't encountered any problem so far, but I'm just checking if this is the correct approach.

5. I tried adding 'singleset' option in dataset 2017 but it didn't change the results (no difference with the results generated by commands without 'singleset' option. Is it because there's no singleton PSU is 2017 strata?

Apologies if I have missed any guidelines, and for the long list of questions.
Thank you and have a nice weekend.
Re: Weight in Indonesia DHS 2002 and 2007 [message #18667 is a reply to message #18614] Fri, 24 January 2020 08:53 Go to previous messageGo to next message
Bridgette-DHS is currently offline  Bridgette-DHS
Messages: 3016
Registered: February 2013
Senior Member

Following is a response by DHS Lead Sampling Statistician, Ruilin Ren:

Your question has been discussed on a few posts - please search the forum. Basically, in the older surveys, the contents of the variable V022 was not the true stratification used in sampling, but the paired clusters which was used for sampling error calculation using the old ISSA software. Region by urban-rural stratification is true for most of the DHS surveys. Even in case where the stratification is slightly modified, using region by urban-rural as stratification will be good for all analysis.
Re: Weight in Indonesia DHS 2002 and 2007 [message #19453 is a reply to message #18667] Thu, 25 June 2020 08:32 Go to previous message
lhuri is currently offline  lhuri
Messages: 16
Registered: December 2014
Location: Indonesia
Member
Thank you for your explanation.
Previous Topic: Weighted result's total gives different total value with unweighted one.
Next Topic: Weights when aggregating survey results
Goto Forum:
  


Current Time: Thu Mar 28 07:00:38 Coordinated Universal Time 2024