I am approaching with regard to Nepal DHS (2016) data that I am using to measure the prevalence of diarrhea in six regions. I'd like to calculate the random effects accounting for the differences across regions. Can anyone kindly mention the process or share the commands to execute this in STATA.

The variables are as follows:

(DV)Diarrhea

(IVs) mother's_age region religion education wealth child's_age sex

Thanks in advance for your help

Best regards

Ghose]]>

To account for variation across 6 regions, a fixed effects model would be preferable. A rule of thumb for random effects is that you have at least 30 cases at that level, i.e. at least 30 regions. For fixed effects, you only need to include something like "i.v024" on the right hand side of the model.]]>

Thanks indeed for the kind help!

Could you please also write a bit about the process of measuring random effects, lets say across the PSUs, or across countries e.g. when the study includes 5 countries in South Asia. How could I quantify the variations occurring at PSUs or country level.

Am I right in thinking that the difference in R2 for the models with and without PSUs/countries provides some indication of the PSU/country level variations.

]]>

DHS always advises that you make three adjustments for the survey design, with svyset and svy if you are using Stata. In order of importance, these are the sample weights, to eliminate a bias toward the over-sampled strata; an adjustment for the clusters, to account for intra-class correlation within the PSUs or clusters; and an adjustment for stratification. The full svyset command is this: "svyset v001 [pweight=v005], strata(stratmid) singleunit(centered)". Where I have put "centered" you could instead put "scaled" or "uncertainty" with little difference. Where I have put v001 you could put v021 (usually, v001=v021). "stratumid" is v022 or v023 in most surveys. There have been many postings about svyset.

The position for v001 within this command, which identifies it as the primary sampling unit, actually specifies that there is to be a random effect for v001. You do not need a separate specification of a random effect for clusters.

Because you have only 5 countries, you should have fixed effects, rather than random effects, for countries. On the right hand side of the model you would have something like "i.country", where country takes the values 1, 2, 3, 4, 5, say.

Having said this, I strongly advise against pooling countries this way. It makes much more sense to do a separate analysis for each country and then present the results side by side. If you pool the countries, you will be describing an "average" for the five countries, with complexities about whether the countries are weighted by sample size or population size, etc. Wouldn't be more interesting to do separate analyses and see how they are similar and how they are different. Even within a country, particularly within India, there is huge variation. You can do statistical tests of whether coefficients, means, proportions, etc. in one country are significantly different from the corresponding terms in another country.

]]>

I would like to learn SPSS code for the following stata command: "svyset v001 [pweight=v005], strata(stratmid) singleunit(centered)"

Your sincerely,

Dr. Tasnuva]]>

I am working with DHS data in STATA for my research and was wondering how to apply mother fixed effect in order to determine how many children are in the same household.

I know you can cluster by household number V002, however, I am unable to determine the total number of children per mother in my regression. I know that applying mother fixed effects will reduce the sample size.

I am very grateful for your assistance on how to apply mother fixed effect.

Regards,]]>

"There are probably several ways to do this, but I suggest the following, using b16 in the BR file. Open a BR file and then insert these lines:

tab b16 keep if b16>0 & b16<. gen nch=1 collapse (sum) nch, by(v001 v002 v003) tab nch

This will give you the number of the women's children (of any age) who are living in the same household as her. You can attach this number to the woman's record with a merge, but in your analysis you probably want to use a file in which the children are the units, i.e. the BR or KR file. If you are using outcomes that are only in the KR file, you can apply these lines to that file, getting the number of children under 5 who are living in the same household as the mother."

https://github.com/DHSProgram/DHS-Indicators-Stata

https://github.com/DHSProgram/DHS-Indicators-SPSS

]]>

Many thanks!

Regards]]>