Multilevel modeling with svy and domestic violence weights [message #23647] |
Tue, 26 October 2021 08:50 |
hbsheldon
Messages: 4 Registered: July 2021
|
Member |
|
|
Hello,
I am working on running a meprobit that looks at the determinants of domestic violence experience at the country-level (for several African countries, but with a regression for each country). I recently used this DHS Methodological Report 27 (https://www.dhsprogram.com/pubs/pdf/MR27/MR27.pdf) to estimate level 1 and level 2 weights for IPV acceptance (which is answered by all women, not just women in the domestic violence module). I am looking for recommendations on how to calculate level 1 and level 2 weights for domestic violence. Repeating the same procedure outlined in MR27 for the domestic violence weight yields weights that are inconsistent between groups for level 2.
1) Some other forums suggested using the geometric mean or a weight = 1 for level 2 weights. Which weight is preferred for level 2?
2) For level 1, would it be ok to use the domestic violence weighted version of the alpha values for level 1, or the domestic violence weight on its own? Based on my research, I am concerned that using the domestic violence weight without any sort of rescaling/manipulation would be incorrect.
Thank you in advance!
[Updated on: Tue, 26 October 2021 13:16] Report message to a moderator
|
|
|
|
Re: Multilevel modeling with svy and domestic violence weights [message #23657 is a reply to message #23656] |
Wed, 27 October 2021 17:16 |
hbsheldon
Messages: 4 Registered: July 2021
|
Member |
|
|
Thank you for your response! When I try the same technique as MR27 but only with replacing v005 with d005, I get the following error:
"weights in variable wt2_domv_4 not constant within groups defined by: v001"
where wt2_domv_4 is wt2_4 but for domestic violence.
Here's a comparison of weights by cluster below:
. list v001 wt2_1 wt1_1 wt2_domv_1 wt1_domv_1 v005_weight domv_weight in 1/50 if v044 == 1
v001 wt2_1 wt1_1 wt2_do~1 wt1_do~1 v005_w~t domv_w~t
-
2. 1 1.5 8.878671 1.5 6.819981 .118748 .091214
10. 1 1.5 8.878671 1.5 6.819981 .118748 .091214
13. 1 1.5 8.878671 1.5 6.819981 .118748 .091214
14. 1 1.5 8.878671 1.5 15.58844 .118748 .208488
19. 1 1.5 8.878671 1.5 6.819981 .118748 .091214
[Updated on: Wed, 27 October 2021 17:16] Report message to a moderator
|
|
|
Re: Multilevel modeling with svy and domestic violence weights [message #23670 is a reply to message #23657] |
Fri, 29 October 2021 15:42 |
Bridgette-DHS
Messages: 3195 Registered: February 2013
|
Senior Member |
|
|
Following is another response from DHS Senior Sampling Specialist, Mahmoud Elkasabi:
Yes, this is correct because the DV weight is calculated on the individual level and not on the cluster level as other weights. I would suggest approximating the wt2 using the women weight or the HH weight and then calculating wt1 based on d005 and the approximated wt2. See the Stata code below as an example.
gen DHSwt = v005 / 1000000
************************************************************************
* Stage B *** Approximate Levels-weight ***
************************************************************************
* Steps to approximate Level-1 and Level-2 weights from Household or Individual Weights
*Step 1. De-normalize the final weight, using approximated normalization factor
gen d_HH = DHSwt * (M/m_c)
gen dv_HH = (d005/1000000) * (M/m_c)
*Step 2. Approximate the Level-2 weight
* f the variation factor
gen f = d_HH / ((A_h/a_c_h) * (M_h/S_h))
* Calculating the levels-weight based on different values of alpha
local alphas 0 0.1 .25 .50 .75 0.90 1
local i = 1
foreach dom of local alphas{
gen wt2_`i' = (A_h/a_c_h)*(f^`dom')
gen wt1_`i' = d_HH/wt2_`i'
gen d1_`i' = dv_HH/wt2_`i'
local ++i
}
svyset v001, weight(wt2_2) strata(v022) , singleunit(centered) || _n, weight(d1_5)
svy: melogit cumodern age secondary urban i.children || v001:
estimates store model2
[Updated on: Fri, 18 March 2022 10:38] Report message to a moderator
|
|
|
|
|