The DHS Program User Forum
Discussions regarding The DHS Program data and results
Home » Topics » Domestic Violence » Calculation of prevalence using STATA command
Calculation of prevalence using STATA command [message #172] Wed, 20 March 2013 17:38 Go to next message
liladhar is currently offline  liladhar
Messages: 2
Registered: March 2013
Location: Denmark
Member
I want to calculate prevalence of "any sexual violence by husband/partner" among currently married women of 15-49 years. I am using data set "NPIR60FL". Following are the STATA command i tried.
keep if V044==1
keep if V501==1
gen rweight=D005/1000000
svyset V021 [pweight=rw], strata(V022) vce(linearized) singleunit(missing)
svy linearized : tabulate V013 D108, cell

# i got following result:

(running tabulate on estimation sample)

Number of strata = 25 Number of obs = 3373
Number of PSUs = 289 Population size = 3084.4815
Design df = 264

----------------------------------------
| Experienced any sexual
Age in | violence (D105H-I) by
5-year | husband/partner
groups | No Yes (D10 Total
----------+-----------------------------
15-19 | .0742 .0098 .084
20-24 | .1769 .0246 .2015
25-29 | .1677 .0243 .192
30-34 | .1339 .0258 .1597
35-39 | .1216 .0213 .1429
40-44 | .1025 .0188 .1213
45-49 | .0831 .0154 .0986
|
Total | .8599 .1401 1
----------------------------------------
Key: cell proportions

Pearson:
Uncorrected chi2(6) = 7.9120
Design-based F(5.10, 1346.15)= 0.6263 P = 0.6829

But in NDHS 2011 report table no 14.3 percentage of married women who have experienced sexual violence is 15.2 and in my result it is 14.01. Number of women is 3,084 in report and 3084.48 in my result.

So would you please let me know how can i get exact result as in report of NDHS 2011.
Re: Calculation of prevalence using STATA command [message #318 is a reply to message #172] Mon, 15 April 2013 14:29 Go to previous messageGo to next message
Bridgette-DHS is currently offline  Bridgette-DHS
Messages: 3017
Registered: February 2013
Senior Member
Here is a response from one of our DHS experts Kerry MacQuarrie, that should answer your question.

You are using the domestic violence weight correctly, with the difference between the n you estimated and that found in the final report being due to rounding only. The discrepancy between the 14.01% you arrived at and the 15.2% in the report is because Table 14.3 in the report is NOT restricted to sexual violence perpetuated by the husband/partner only (which also explains how the Table reports a figure for prevalence of sexual violence among never married women). Instead of relying on variable d108 (ever experienced any sexual violence (d105h-i) by husband/partner) alone, this table uses a combination of d108 and variables d124 and d125. D124 is ever forced to have sex by anyone other than husband/spouse in last 12 months and d125 is ever forced to perform unwanted sexual acts. I found that if I generated the following composite variable:

g anysv=.
replace anysv=0 if d108==0 | d124==0 | d125==0
replace anysv=1 if d108==1 | d124==1 | d125==1

that I was then able to replicate the figures in Table 14.3:

tab v501 anysv [iweight=dvwt], row

index.php?t=getfile&id=11&private=0

If you are interested in sexual violence perpetuated by husbands/partners and not by others, you are correctly using d108. Please note that there is another variable available in the dataset: d123, which captures whether women's first sexual initiation was forced and this variable, in conjunction with those that describe sexual initiation as occurring within (and only one) marriage is sometimes added to d108 to capture any sexual violence by husband/partner. This variable is not included in the composite variable producing the results for Table 14.3.

I hope this is helpful.
  • Attachment: capture.bmp
    (Size: 666.34KB, Downloaded 3003 times)

[Updated on: Mon, 15 April 2013 14:30]

Report message to a moderator

Re: Calculation of prevalence using STATA command [message #2055 is a reply to message #318] Fri, 18 April 2014 16:51 Go to previous messageGo to next message
malayaka is currently offline  malayaka
Messages: 18
Registered: March 2014
Member
Hello,

I am using STATA 11.2 and having difficulty generating the DHS Malawi 2010 domestic violence figures (see attached jpgs) for:

1) Women ages 15-19 who ever experienced sexual violence (%)
2) Women ages 15-19 who ever experienced physical violence since age 15 (%)

Below are the STATA codes I used for "sexual violence (please see attached pdf)." I then plan to crosstab by wealth, residence (urban/rural), and region. I am also at a loss at what to do for "ever experienced physical violence since age 15" -- which variables do I use (or need to combine - similar to "sexual violence") to get "ever experienced physical violence?"

Can someone please help with the correct code to use?

I used the STATA codes below, but am not able to generate the same results (see attached pdf). Guidance to generate both DV figures will be much appreciated.

STATA Code:

use v005 v012 v013 v021 v022 v023 v024 v025 v044 v101 v102 v106 v149 v155
v190 v201 v212 v313 v384a v384b v384c v463z v501 v502 v511 v536 v613 v714
v717 v731 v739 v763a d005 d108 d124 d125 using "/Users/Desktop/MWIR61FL.DTA


******************************
* DOMESTIC VIOLENCE (SEXUAL)
******************************
generate dweight = d005/1000000
svyset [iweight=dweight], psu(v021) strata(v023)

* Generate new variable (anysv = any sexual violence)
generate anysv=.
replace anysv=0 if d108==0|d124==0|d125==0
replace anysv=1 if d108==1|d124==1|d125==1

* Age group (15-19 yrs) and experience of sexual violence
svy: tabulate v013 anysv, row

* Age group (15-19 yrs) and experience of sexual violence
svy: tabulate anysv v190 if v013==1

* By age group (15-19 yrs), experience of sexual violence and place of
residence (urban/rural)
svy: tabulate anysv v102 if v013==1

* By age group (15-19 yrs), experience of sexual violence and region
(North, Central, South)
svy: tabulate anysv v101 if v013==1

Thank you in advance.
Re: Calculation of prevalence using STATA command [message #2056 is a reply to message #2055] Fri, 18 April 2014 18:26 Go to previous message
Liz-DHS
Messages: 1516
Registered: February 2013
Senior Member
Dear User,
I am not a STATA user or a programmer. However, if you can provide the table number and page number from the Malawi DHS 2010 report I might be able to provide some CSPro code that might be useful to you in figuring things out.
Previous Topic: Selection for domestic violence module V044
Next Topic: DHS Malawi - domestic violence among women 15-19
Goto Forum:
  


Current Time: Fri Mar 29 10:17:09 Coordinated Universal Time 2024