Home » Topics » Domestic Violence » Multilevel Weights
Multilevel Weights [message #7019] |
Tue, 11 August 2015 13:16 ![Go to next message Go to next message](/theme/default/images/down.png) |
swinter
Messages: 2 Registered: August 2015
|
Member |
|
|
Greetings DHS experts!
I recently posted a question on the forum about the domestic violence weights that is somewhat related to this post, but is focused on different questions:
( http://userforum.dhsprogram.com/index.php?t=tree&th=4410 &goto=6918&S=85ccc984d9dd623f59e3088f0d1e73c9#msg_69 18)
This post is meant to focus on the multilevel nature of the DHS weights (or, seemingly, the lack thereof). I am conducting a number of two and three-level multilevel modelling using both individual-level country data (two-level models) and multi-country data (three-level models). I have been doing the modeling largely in Stata 14 (using melogit) and MLwin 2.34 (using PQL and MCMC). Unfortunately, both of these softwares provide very specific instructions in their manuals to use multilevel sampling weights instead of a single-level weight when conducting multilevel analyses (except MCMC in MLwin or crossed mixed effects models, which don't allow weights). For an example of these instructions, see following passage from Stata 14 manual:
"it is not sufficient to use the single sampling weight wij , because weights enter the log likelihood at both the group level and the individual level. Instead, what is required for a two-level model under this sampling design is wj , the inverse of the probability that group j is selected in the first stage, and wi|j, the inverse of the probability that individual i from group j is selected at the second stage conditional on group j already being selected. You cannot use wij without making any assumptions about wj. Given the rules of conditional probability, wij = wjwi|j. If your dataset has only wij , then you will need to either assume equal probability sampling at the first stage (wj = 1 for all j) or find some way to recover wj from other variables in your data; see Rabe-Hesketh and Skrondal (2006) and the references therein for some suggestions on how to do this, but realize that there is little yet known about how well these approximations perform in practice. What you really need to fit your two-level model are data that contain wj in addition to either wij or wi|j. If you have wij--that is, the unconditional inclusion weight for observation i, j--then you need to divide wij by wj to obtain wi|j" (Stata 14 Manual - "meglm -- Multilevel mixed-effects generalized linear model, p.21 available at: http://www.stata.com/manuals14/memeglm.pdf#memeglmMethodsand formulas)
From my reading of the DHS sampling literature, the multilevel nature of the DHS sampling is particularly important in the domestic violence sampling weights because, unlike the other weights (v005 and hv005), individual women sampled for the dv module do not have the same weight as the households. So, it seems to me, that for a two-level model, there should be, at a minimum, a PSU-level (level 2) weight and an individual-level (level 1) dv weight that incorporates the dv sample design and non-response. Does anyone have suggestions about how to tackle the multilevel weighting issue? Should individuals interested in multilevel modeling just assume the sampling probability at the first-stage (e.g. the PSU weight) is equal for all PSUs (e.g. wj=1 for all j)? Or, should we try to "recover" the multilevel weights using the technique cited in Rabe-Hesketh and Skrondal (2006) or another technique? Will DHS provide methodology for extracting the multilevel weights (I did try doing this based on the equations in the DHS sampling manual, but I have more unknowns than equations, particularly in regards to back-calculating the domestic violence weights). Any other thoughts?
Thanks so much!
|
|
|
Re: Multilevel Weights [message #7062 is a reply to message #7019] |
Tue, 18 August 2015 12:28 ![Go to previous message Go to previous message](/theme/default/images/up.png) ![Go to next message Go to next message](/theme/default/images/down.png) |
Liz-DHS
Messages: 1516 Registered: February 2013
|
Senior Member |
|
|
Dear User,
Here is a response from on of our technical experts, Dr. Tom Pullum:
Quote:We agree that multi-level analysis should include weights at each level. Until Stata 14, multi-level models in Stata could not use weights at all. Beginning with Stata 14, weights (multi-level weights) are allowed and should be used (as in MLwin). The problem is that DHS does not have cluster-level weights. The clusters are sampled with probability proportional to size, within strata. The sizes of the clusters, usually census enumeration areas, are part of the sampling frame, usually the most recent census. The sampling frame is not public information. DHS only has access to it within the country and is not allowed to make a copy of it. Statistical offices typically do not want to share it.
It is conceivable that there is some way to approximate the cluster-level weights or at least to get away from the invalid assumption that the first stage sampling probabilities are the same for all clusters. We will look into this. Thanks for raising the issue.
Thank you!
[Updated on: Wed, 19 August 2015 15:12] Report message to a moderator
|
|
|
|
|
|
Re: Multilevel Weights [message #10124 is a reply to message #10019] |
Wed, 29 June 2016 11:15 ![Go to previous message Go to previous message](/theme/default/images/up.png) ![Go to next message Go to next message](/theme/default/images/down.png) |
Liz-DHS
Messages: 1516 Registered: February 2013
|
Senior Member |
|
|
Dear User,
A response from sampling expert, Dr. Mahmoud Elkasabi:
Quote:
This means that V005 is the normalized version of the sampling weight. The main purpose of the normalization process is to avoid the large values for the number of weighted cases in the tables in DHS survey final reports. This applies for all the DHS standard weights, including weights for households, such as HV005, and individuals. The V005 was calculated by multiplying the individual sampling weight by a normalization factor at the national level. The normalization factor is the total number of completed cases divided by the total number of weighted cases. In case of the V005, it is the total number of completed women divided by the weighted total number of completed women. In case of the HV005, it is the total number of completed households divided by the weighted total number of completed households. Therefore the standard weights in the DHS data files are relative weights. Relative weights can be used to estimate means, proportions, rates and ratios because the normalization factor is cancelled out when used in both numerator and denominator, so it has no effect on the calculated indicator values. However, the standard weights are not valid for estimating totals. Also the normalized weight is not valid for pooled data, even for data pooled for women and men in the same survey, because the normalization factor is country and sex specific.
|
|
|
|
|
|
Re: Multilevel Weights [message #12659 is a reply to message #12416] |
Wed, 28 June 2017 16:55 ![Go to previous message Go to previous message](/theme/default/images/up.png) ![Go to next message Go to next message](/theme/default/images/down.png) |
Liz-DHS
Messages: 1516 Registered: February 2013
|
Senior Member |
|
|
A response from Senior Sampling Expert, Dr. Ruilin Ren:
Quote:
Unfortunately, there is nothing new. With the confidentiality requests from the DHS protocol, we cannot provide the selection probabilities (we cannot keep them) of different levels which are the components of the multi-level weights. So the question is not a technical one, but rather confidential obligations.
Thanks
Ruilin
|
|
|
|
Re: Multilevel Weights [message #13482 is a reply to message #13289] |
Wed, 08 November 2017 14:52 ![Go to previous message Go to previous message](/theme/default/images/up.png) ![Go to next message Go to next message](/theme/default/images/down.png) |
Liz-DHS
Messages: 1516 Registered: February 2013
|
Senior Member |
|
|
Dear User, A response from Dr. Tom Pullum:
Quote:
Individuals in households are not sampled. In the survey design, clusters are sampled and then households are sampled within clusters. After the household has been selected, all eligible respondents (based on age, sex, and de facto residence) are selected.
Ideally we would provide separate sampling fractions (or their inverse, the weights) for the clusters and then the households. At this time it is not possible to do a full multilevel model because we can only provide the product (hv005, etc.) after an adjustment for nonresponse. As has been stated in other responses, for privacy reasons we do not save the more detailed information.
|
|
|
|
|
Re: Multilevel Weights [message #14700 is a reply to message #14413] |
Mon, 30 April 2018 14:53 ![Go to previous message Go to previous message](/theme/default/images/up.png) ![Go to next message Go to next message](/theme/default/images/down.png) |
Hassen
Messages: 121 Registered: April 2018 Location: Ethiopia,Africa
|
Senior Member |
|
|
Dear DHS Experts, Can Multilevel Modeling be conducted Using DHS Data set by using SPSS? Can I make Weighting at each levels (community,household and individual levels)? I already planned to use Multilevel modeling to investigate factors associated with childhood nutritional status using 2016 Ethiopia DHS Data set. What are your recommendation to conduct it using SPSS?
I need your reply.
Thank you in advance!!
Hassen Ali(Chief Public Health Professional Specialist)
[Updated on: Tue, 01 May 2018 02:13] Report message to a moderator
|
|
|
Re: Multilevel Weights [message #14717 is a reply to message #14700] |
Tue, 01 May 2018 23:28 ![Go to previous message Go to previous message](/theme/default/images/up.png) ![Go to next message Go to next message](/theme/default/images/down.png) |
Liz-DHS
Messages: 1516 Registered: February 2013
|
Senior Member |
|
|
Dear User, A response from Dr. Tom Pullum:
Quote:
The DHS analysis group uses Stata, sometimes R. We do not use SPSS and we cannot help you with the syntax of multilevel commands in SPSS. Regarding weights, we do not have separate weights for the different levels, only the net weight, which is proportional to hv005 or v005. We intend to develop recommendations for how to partition this into cluster-level and household-level weights but are not yet prepared to suggest anything specific.
|
|
|
Re: Multilevel Weights [message #14718 is a reply to message #14717] |
Wed, 02 May 2018 03:14 ![Go to previous message Go to previous message](/theme/default/images/up.png) ![Go to next message Go to next message](/theme/default/images/down.png) |
Hassen
Messages: 121 Registered: April 2018 Location: Ethiopia,Africa
|
Senior Member |
|
|
Thank you very much for your Attractive Response!
I am waiting to get your reply on How to partition the samples into cluster-level and household-level!! I will come up with My Challenges after I have seen all issues regarding Multilevel analysis etc.
Again,Thank you in advance!!
Hassen Ali(Chief Public Health Professional Specialist)
|
|
|
|
Re: Multilevel Weights [message #15354 is a reply to message #15332] |
Tue, 03 July 2018 23:04 ![Go to previous message Go to previous message](/theme/default/images/up.png) ![Go to next message Go to next message](/theme/default/images/down.png) |
Liz-DHS
Messages: 1516 Registered: February 2013
|
Senior Member |
|
|
A response from Dr. Shireen Assaf:
Quote:
Dear user,
To be able to use multilevel modeling with DHS data and the svy command, a weight must be applied for each level. Since we only have one weight, we can make the assumption that all individuals in the household have the same weight. The Stata code to run a multilevel model for a binary outcome is below using the melogit command. However, you can use other mixed model commands if your outcome had a different distribution. The svyset should be the same.
gen wt=v005/1000000
gen wt2=1
svyset v001, weight(wt) strata(v023) , singleunit(centered) || _n, weight(wt2)
*this is a random intercept model
svy: melogit outcome var1 var2 || v001:
*random intercept and random slop for var3
svy: melogit outcome var1 var2 || v001: var3
*for this model you can also add the covariance(unstructured) option. Please read the stata documentation for the melogit command on this
* you can also check the Stata documentation for svyset, they discuss how to construct this for multistage sample design.
Hope this helps.
Best,
Shireen Assaf
Technical Specialist
The DHS Program
|
|
|
Re: Multilevel Weights [message #15885 is a reply to message #15354] |
Thu, 04 October 2018 06:35 ![Go to previous message Go to previous message](/theme/default/images/up.png) ![Go to next message Go to next message](/theme/default/images/down.png) |
teketo
Messages: 2 Registered: November 2016 Location: Debre Markos
|
Member |
|
|
Dear DHS,
I am doing analysis on maternal health service use (modern contraception among married women, antenatal care, health facility delivery and caesarean section). I am analysing the 2016 Ethiopia DHS data using SAS software.
I am just wondering how I can use Multilevel weights for this analysis. It will be great if I can get solutions how it will be done using SAS.
I am using a two level random effects model (Level 1: individuals, and Level 2: regions) for modern contraception among married women, antenatal care and health facility delivery. Moreover, for caesarean section, a three level random effects model (Level 1: individuals, Level 2: clusters and Level 3: regions) using MCMC method.
The other question I have is what should be the maximum number of a grouping (membership) variable? To make it clear, just have a look the following SAS Proc glimmix program:
Proc glimmix data = care;
Class region V001;
Model y (event = last) = a b c d e f g h i / solution;
Random intercept e f g / solution subject = region;
Random intercept a b c d / solution subject = V001 (region);
Run;
The grouping variables showing how one is nested within the other on the above SAS code are region level 3 (there are 11 regions) and V001 level 2 (there are 622 clusters). We are asking the program to produce random effects for each of the 11 regions and the 622 clusters. It will have a convergence issue it will take long time to process and even might stop processing.
Regards
Teketo
teketo
|
|
|
Re: Multilevel Weights [message #19021 is a reply to message #15354] |
Sat, 04 April 2020 12:00 ![Go to previous message Go to previous message](/theme/default/images/up.png) ![Go to next message Go to next message](/theme/default/images/down.png) |
aheto
Messages: 1 Registered: August 2014 Location: United Kingdom
|
Member |
![jascall12](/theme/default/images/skype.png)
|
|
Hi Dr. Shireen Assaf,
Thanks for this response.
How can I implement your suggestions below in R for the standard multilevel you presented, and also for Bayesian multilevel model? Any R codes to work with?
gen wt=v005/1000000
gen wt2=1
svyset v001, weight(wt) strata(v023) , singleunit(centered) || _n, weight(wt2)
*this is a random intercept model
svy: melogit outcome var1 var2 || v001:
Thanks
[Updated on: Sat, 04 April 2020 12:04] Report message to a moderator
|
|
|
Re: Multilevel Weights [message #19083 is a reply to message #19021] |
Fri, 17 April 2020 09:17 ![Go to previous message Go to previous message](/theme/default/images/up.png) |
Liz-DHS
Messages: 1516 Registered: February 2013
|
Senior Member |
|
|
A response from Dr. Shireen Assaf:
Quote:
Hello,
Unfortunately I do not have any R code for this. However, please study the survey package in R and you may find an answer.
Thank you.
Best,
Shireen
|
|
|
Goto Forum:
Current Time: Sat Aug 31 02:08:25 Coordinated Universal Time 2024
|