Calculation of prevalence using STATA command [message #172] |
Wed, 20 March 2013 17:38 |
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.
|
|
|