Malaria (by microscopy) in Lagos State [message #27260] |
Sat, 08 July 2023 20:31 |
cleonard
Messages: 3 Registered: November 2017
|
Member |
|
|
Hello,
I am calculating the weighted number of children 6- 59 months tested for malaria by microscopy and the number of children positive for malaria by microscopy. The dataset I am using is NGPR81FL. I am using the exact guidelines as shown in the DHS Guide to Statistics:
Number of de facto children tested using microscopy (hv042 = 1 & hv103 = 1 & hc1 in 6:59 & hml32 in 0,1,6)
Number of de facto children tested using microscopy who are positive for malaria (hml32 = 1)
I am also taking into account the complex sample design which include strata, cluster and weights. I am coding in R and have attached my code sample to this post.
I am comparing my results with the 2021 Nigeria MIS Report, specifically pg. 147 (Table B.46), which does not match up. I am getting 205 for weighted number of children tested by microscopy, and 2 children who tested positive.
Has anyone else had the issue? I wonder if something is wrong with my code. Thank you very much.
Colleen
|
|
|
Re: Malaria (by microscopy) in Lagos State [message #27313 is a reply to message #27260] |
Mon, 17 July 2023 12:15 |
Janet-DHS
Messages: 895 Registered: April 2022
|
Senior Member |
|
|
Following is a response from DHS staff member, Tom Pullum:
I use Stata and am not referring to your R code. To get the tabulation it is only necessary to use the weights, not the full svyset command. I opened NGPR81FL.dta and entered the following Stata line:
tab hml32 if hv103==1 & hv024==36 & hc1>=6 & hc1<=59 [iweight=hv005/1000000]
You do not need to specify the value of hv042; it is 1 for all cases. The Stata line produces the following output:
final result of malaria from |
blood smear test | Freq. Percent Cum.
---------------------------------+-------------------------- ---------
negative | 338.41388 97.42 97.42
positive | 8.95048 2.58 100.00
---------------------------------+-------------------------- ---------
Total | 347.36436 100.00
This matches the weighted denominator of 347 and the proportion .026 in Table B.46 on page 147. That is, I match the report and do not get the numbers you got. Can you see where the difference is? Hope so. Good luck and thanks for using DHS data!
|
|
|