The DHS Program User Forum
Discussions regarding The DHS Program data and results
Home » Topics » Domestic Violence » Pakistan 2014 IPV analysis -- multi-level weights
Pakistan 2014 IPV analysis -- multi-level weights [message #8507] Fri, 06 November 2015 12:04 Go to next message
gspek is currently offline  gspek
Messages: 4
Registered: November 2015
I am using the Pakistan 2014 women's dataset for analysis of associations between women's empowerment (and individual and community level) and domestic violence experiences (at individual level). I am currently struggling with a couple of issues:

1. How do I correctly weight community-level independent variables? I tried to create community-level weighted means for average percentage of women working in a community, average decision making index value, and percentage in community accepting of any domestic violence justification scenarios (with participant's own values excluded). However, I'm not sure how to approach the weights in regressions where both the community and individual level independent variables are included.

foreach var of varlist works dm_sum dvaccept_any{
gen `var'_w = `var'*wgt if !missing(`var') & nocomm == 0
bys v001: gen `var'_agg_w = sum(`var'_w) if nocomm == 0
gen `var'_comm_w = (`var'_agg_w - `var'_w)/(clustersz - 1) if nocomm == 0

Furthermore, the number of respondents in each community is only 5-10. What is the best way to test if we are too underpowered to make any statements about the significance of the relationship?

2. If we choose to include men's community-level variables as well, how would that be weighted?

Thank you!
Re: Pakistan 2014 IPV analysis -- multi-level weights [message #8521 is a reply to message #8507] Sun, 08 November 2015 16:40 Go to previous message
Messages: 292
Registered: March 2013
Senior Member

I'm not sure there is a "right" answer here, although surely some wrong ones. I'll give a couple of thoughts:

1 - a weight is just a person's probability of selection (re-scaled in a DHS-manner, but still basically a probability weight). This does not vary by community, because there is random sampling within a cluster. So everyone in each PSU will have the same weight, and so you can just use the individual weights.

2 - Even adding in information from the Men, if the unit of analysis is still the individual woman you would want the women's weights.

3 - if you are trying to estimate "causal effects" of the community level variables, and you think those relationships are independent of the criteria used for sampling villages (so that, say, all women respond the same to a 1sd change in their "decision making" then you don't need to weight. Weighting buys you population-level mean estimates, but if there are "constant treatment effects" regardless of covariates, you may not need to weight*.

4 - cluster-level measurements are based on too few observations to be meaningful in and of themselves - as you say, there are wildly under-powered. A couple of things you could do: a) by averaging over many clusters, you can still get good estimates of community level variables, but each individual cluster-level point-estimate would be very, very noisy. But they may still mostly "agree" in some sense; b) so if in your hierarchical model you allow each cluster an unconstrained cluster-specific effect (like treating each cluster as a mini-experiment), you could look at those individual point-estimates on a scatter plot (say Beta across some variable you think would affect Beta); c) and then you could start restricting those Betas to have some particular distribution (a random slope model) and see how that changes your overall point estimate as you make your priors on the distribution of Beta more/less informative. I think this makes sense as a kind of model-checking or informal/additional inference procedure. A leave-one-out cross-validation approach might make sense too, depending on how you end up thinking about each of these within-cluster estimates.

I hope something there helped.
Previous Topic: Multilevel weights
Next Topic: Cluster identifiers and finite population correction
Goto Forum:

Current Time: Sat Aug 13 13:46:58 Coordinated Universal Time 2022