Re: South Africa 2016 Domestic Violence [message #21992 is a reply to message #21971] |
Thu, 14 January 2021 10:33 |
Shireen-DHS
Messages: 140 Registered: August 2020 Location: USA
|
Senior Member |
|
|
Hello,
For the domestic violence indicators you will need to merge a DV and IR file. This is unique to the South Africa 2016 survey. Usually in other DHS surveys the domestic violence variables are only in the IR file. I provide Stata code below that merges the two files then creates the first domestic violence indicator in the final report (ever experienced physical violence). The denominator I obtain is 9 cases over what is reported in the final report but I have not been able to find a way to exclude these cases. However, the percentage matches since it's only 9 cases off.
Hope you find this useful.
Best,
Shireen Assaf
The DHS Program
Stata code:
*first open the DV file
use "ZADV71FL.DTA" , clear
*removing missing values from the first domestic violence indicator in the file
keep if d101a<.
gen in_DV=1
sort v001 v002 v003
save ZAtemp.dta, replace
*now open IR file
use "ZAIR71FL.DTA" , clear
*remove the same missing values
keep if d101a<.
gen in_IR=1
sort v001 v002 v003
*merge the files using the temporary DV file you saved.
merge 1:1 v001 v002 v003 using ZAtemp.dta
tab _merge
* this is to see which women are coming from each file
tab in_DV in_IR,m
*generate an age variable
gen age10=int((v012-5)/10)
replace age10=6 if age10>6
label define age10 1 "18-24" 2 "25-34" 3 "35-44" 4 "45-54" 5 "55-64" 6 "65+"
label values age10 age10
*tabulate age variable to see denominator
tab age10 [iweight=d005/1000000]
*generate variable for every experienced physical violence
gen PVever = (inrange(d105a,1,4) | inrange(d105d,1,4) | inrange(d105e,1,4) | inrange(d105f,1,4) | inrange(d130a,1,4) )
*the percentage matches the final report
tab PVever [iw=d005/1000000]
|
|
|