The DHS Program User Forum
Discussions regarding The DHS Program data and results
Home » Data » Weighting data » SPA survey Tanzania: Setting survey design and using weights for health facilities and health care p
SPA survey Tanzania: Setting survey design and using weights for health facilities and health care p [message #10380] Fri, 22 July 2016 04:06 Go to next message
masjas
Messages: 4
Registered: July 2016
Location: Amsterdam
Member
Using data from the SPA survey Tanzania we are interested in exploring the association between health care providers characteristics and their response to question w167 (dichotomized in yes and no).

The data set contains facility data and health provider data and includes weights for facility and health provider.
We want to conduct multilevel mixed effect logistic regression analyses (MELOGIT in STATA). As health providers within facility are more likely to be more similar compared to health providers between facilities we want to include random effects for health facilities.
MELOGIT requires that weights for each level need to be included ( http://www.stata.com/new-in-stata/multilevel-models-survey-d ata/).
We have therefore set the survey design to indicate that health care providers (weight u005_2 which is transformed weight factor by dividing U005 by 100,000) are nested within health facilities (hfid with weight v005_2) as follows:
svyset hfid, weight(v005_2)|| _n, weight(u005_2)

Next: we run melogit as follows in STATA (example agecategorie, hfid is health facility unique ID):
svy: melogit w167 i.agecat || hfid;
We would like to hear your thoughts if this is the correct analyses?
STATA requires that weights from all levels are included in setting survey design; therefore we conducted above analyses including. However, the program R does not seem to have this requirement as within R the analyses run while only health care providers weights are included in survey design (while STATA does not run the melogit analyses then as it gives the following error: weights in variable u005_2 not constant within group defined by:hfid)
Re: SPA survey Tanzania: Setting survey design and using weights for health facilities and health care p [message #10387 is a reply to message #10380] Fri, 22 July 2016 10:43 Go to previous messageGo to next message
Bridgette-DHS is currently offline  Bridgette-DHS
Messages: 3214
Registered: February 2013
Senior Member
Following is a response from Senior DHS Stata Specialist, Tom Pullum:

I believe your svyset is correct. Providers nested within facilities is equivalent to households nested within clusters. However, strata should be included in svyset. I just tried this with the Tanzania 2006 SPA data (I merged the facility and provider data and used the outcome variable you indicated). You need to check the final report but usually the SPA strata are region by facility type.

egen strata = group(u001 u007)
gen v005_2= v005/1000000
gen u005_2= u005/1000000
*outcome variable
recode w167 (8/9=0), gen(w167b)

svyset inv_id, weight(v005_2) strata(strata) singleunit(centered) || _n, weight(u005_2)
svy: melogit w167b w102 || inv_id:
Re: SPA survey Tanzania: Setting survey design and using weights for health facilities and health care p [message #10396 is a reply to message #10387] Mon, 25 July 2016 09:29 Go to previous message
masjas
Messages: 4
Registered: July 2016
Location: Amsterdam
Member
Thank you for you clear answer, this will help me to continue.
Previous Topic: Domestic violence and HIV
Next Topic: Weighting variable in Pakistan 2012/13
Goto Forum:
  


Current Time: Thu Dec 26 05:55:59 Coordinated Universal Time 2024