The DHS Program User Forum
Discussions regarding The DHS Program data and results
Home » Topics » Malaria » Multiple datasets (Malaria)
Multiple datasets [message #17635] Wed, 01 May 2019 12:42 Go to next message
annette.cassy is currently offline  annette.cassy
Messages: 4
Registered: February 2018
Location: Mozambique
Member
Hi
I am trying to do some analysis using more than one dataset.
I have appendend but I don't know how to use the weights with more than one dataset.
Can you help me, please?

I need to do:
1. Propottion of HH with universal covarage (I need to understand which variables to use) / year (survey) / province
2. % of people who sought care in Community health workers / year (survey) / province
3. % of under 5 that slept under ITN / year (survey) / province
4. % of population that slept under ITN by gender / year (survey) / province
5. % of HH with IRS / year (survey) / province


Annette Cassy
Re: Multiple datasets [message #17638 is a reply to message #17635] Wed, 01 May 2019 14:38 Go to previous messageGo to next message
Bridgette-DHS is currently offline  Bridgette-DHS
Messages: 3199
Registered: February 2013
Senior Member

Following is a response from DHS Lead Malaria Research Analyst, Cameron Taylor:

Thanks for your question. For appending datasets you do need to take into account weighting. You need to construct unique ID codes for the clusters (PSUs) in the separate surveys. I suggest the following way of creating a new cluster and stratum id that is unique taking into account both survey years. Construct a unique PSU ID with "egen clusterid=group(hv007 hv021)". Similarly, if hv024 is the stratum variable, you can construct a unique stratum ID with "egen stratumid=group(hv007 hv024)". Then you put clusterid and stratumid in the appropriate places in svyset.

For your questions about how to calculate different indicators I have included the code below. In each of the example code I limited the tabulation to the Mozambique 2015 data (hv007==2015), since I want to match the report to make sure my tabulation is correct before proceeding. I wasn't sure what you were wanting with request #2 (% of people who sought care in Community health workers). However, if you are wanting to look at malaria care seeking data you need to use the KR file since your unit of analysis is children with fever.

1. Proportion of HH with universal coverage (I need to understand which variables to use) / year (survey) / province
5. % of HH with IRS / year (survey) / province


append using MZHR62FL.DTA MZHR71FL.DTA

egen clusterid=group(hv007 hv021)
egen stratumid=group(hv007 hv024)

gen wgt=hv005/1000000

svyset [pw=wgt], psu(clusterid) strata(stratumid)

//Household ITN ownership
          gen mal_ITNinHH=0
                   forvalues x=1/7 {
                             replace mal_ITNinHH=1 if hml10_`x'==1
                   }
          lab var mal_ITNinHH "Household owns at least one ITN"

//Number of ITNs per household
          gen mal_numitnhh=0 
                   forvalues x=1/7 {
                             gen itnhh_0`x'=(hml10_`x'==1)
                   }
          replace mal_numitnhh=itnhh_01 + itnhh_02 + itnhh_03 + itnhh_04 + itnhh_05 + itnhh_06 + itnhh_07
          lab var mal_numitnhh "Number of ITNs per household"

//Potential ITN users in Household
          gen mal_potuse = mal_numitnhh*2
          lab var mal_potuse "Potential ITN users in household"

          
//Households with > 1 ITN per 2 members
//Potential users divided by defacto household members is greater or equal to one
          gen mal_hhaccess = ((mal_potuse/hv013)>=1)

          lab var mal_hhaccess "Households with >1 ITN per 2 household members"

//IRS Household Sprayed
//Was the household sprayed in the past 12 months by someone other than a household member
          g irs=0
          replace irs=1 if hv253==1 & (hv253a==1|hv253b==1|hv253c==1)
          lab var irs "Household sprayed in the past 12 months by someone other than a household member"

//TABULATIONS
svy: tab hv025 mal_hhaccess if hv013>=1 & hv007==2015, row ci obs

svy: tab hv025 irs if hv007==2015, row ci obs


3. % of under 5 that slept under ITN / year (survey) / province
4. % of population that slept under ITN by gender / year (survey) / province

append using MZPR62FL.DTA MZPR71FL.DTA

egen clusterid=group(hv007 hv021)
egen stratumid=group(hv007 hv024)

gen wgt=hv005/1000000

svyset [pw=wgt], psu(clusterid) strata(stratumid)


//Categorizing nets
          gen mal_netcat=0 if hml12==0
          replace mal_netcat=1 if hml12==1|hml12==2
          replace mal_netcat=2 if hml12==3
          lab var mal_netcat "Mosquito net categorization"

//ITN net variable
          gen mal_itn=(mal_netcat==1)
          lab var mal_itn "ITN"

//TABULATIONS
*Request #3
svy: tab hv024 mal_itn if hv103==1 & hml16<5 & hv007==2015 [iw=wgt] , row ci obs

*Request #4
svy: tab hv024 mal_itn if hv103==1 & hv007==2015, row ci obs

svy: tab hv104 mal_itn if hv103==1 & hv007==2015, row ci obs
Re: Multiple datasets [message #17640 is a reply to message #17638] Wed, 01 May 2019 15:04 Go to previous message
annette.cassy is currently offline  annette.cassy
Messages: 4
Registered: February 2018
Location: Mozambique
Member
Thanks a lot.
I am going to try this right now.


Annette Cassy
Previous Topic: Artemisin-based Combination Therapy (ACT) in Mozambique 2015
Next Topic: Sampling size determination
Goto Forum:
  


Current Time: Thu Nov 21 07:55:56 Coordinated Universal Time 2024