Home » Topics » Child Health » Confidence Intervals ORS useNigeria 2013
Confidence Intervals ORS useNigeria 2013 [message #3138] 
Tue, 21 October 2014 23:46 
pr2423
Messages: 2 Registered: October 2014 Location: NYC

Member 


Hi,
I am trying to replicate % of children who received ORS among children who had diarrhea in the past 15 days from the Nigeria 2013 data. I am using SAS for my analysis. For country, region and urban/rural I get the correct point estimates and standard error only when I ignore the strata command in SAS and use only clusters. I know that the data was first stratified by state and then urban/rural, is it okay to ignore stratification while constructing confidence intervals? would this change if I wanted confidence intervals for each state? I would really appreciate any insight into this! am I missing something about the study/sample design?
Please find the code that gives me the correct answer below, adding the strata command in proc survey freq(strata = v022) to this I get the wrong answer
/*creating variable for having been given ORS*/
data ors;
set ors;
if h13 = 1 or h13 = 2 then ors =1;
else ors =0;
run;
data ors;
set ors;
wgt = v005/1000000;
run;
/*h11 = 2 when child had diarrhea in the past 15 days, v024 = region*/
proc surveyfreq
data = ors;
where h11 =2 and v024 =5;
cluster v021;
weight wgt;
tables ors;
run;



Re: Confidence Intervals ORS useNigeria 2013 [message #3150 is a reply to message #3138] 
Thu, 23 October 2014 13:02 
LizDHS
Messages: 1516 Registered: February 2013

Senior Member 


Dear User,
Below is a response from one of our sampling experts, Dr. Mahmoud Elkasabi:
Quote:The standard errors in the 2013 Nigeria DHS report were calculated using the region by urban/rural (V024 cross V025) as stratification variable, which is the case for most of the DHS surveys. However, the 2013 Nigeria DHS used a more detailed stratification which is sate cross urban/rural and the stratification is coded in V022, but attention must be paid when use V022 (for the 2013 Nigeria DHS) as it has one singlecluster stratum, the urban of AKWA IBOM state. You may need to collapse it with the rural stratum of the state because sampling error calculation will not allow singlecluster stratum. If you want to replicate the DHS results on standard errors, you need to use V024 cross V025 as stratification. But you can use V022 as stratification, this gives the exact stratification used in the sample selection and this may give you more accurate results for other analysis such as modeling.
By the way, if you just want to produce point estimation with confidence interval (CI), you would better use Proc Surveymeans instead of Proc Surveyfreq. Proc Surveymeans will give you the 95% CI in the output. However, note that your results may not exactly match the published results, since DHS uses mean minus/plus 2 times the estimated standard error as confidence limits. While the SAS output uses mean minus/plus t(d.f, 95%) as confidence limits. You can use "Class" statement to produce results for urban/rural and regions separately within the Proc Surveymeans statement, and use "Domain" to produce results based on other socialeconomic, education or age group variables. Hope this helps.




Re: Confidence Intervals ORS useNigeria 2013 [message #3168 is a reply to message #3160] 
Tue, 28 October 2014 16:26 
LizDHS
Messages: 1516 Registered: February 2013

Senior Member 


Dear User,
Here is a response from one of our sampling experts, Dr. Mahmoud Elkasabi,
Quote:To produce state level results, you can use V022 (sate cross urban/rural) as stratification variable. In Proc Surveymeans, you can add "Domain SHSTATE". It should produce the indicator results on the state level with the SE and 95% CI. However, we do not recommend you produce state level estimates, since the survey sample was not designed for this purpose, especially for the main indicators. The enough sample size to produce a DHS indicator differs based on the indicator type. For the TFR and the mortality rates, at least 800 to 1000 households is required for producing reliable indicators (with acceptable precision).
Regarding the mortality rates, if you are able to replicate the rate itself, it would be easy to apply the Jackknife variance estimation by applying the formula given in the DHS report appendix B; you need to calculate a series of UMR each time by dropping one cluster, then apply the formula.
Thank you for your post.



Goto Forum:
Current Time: Thu Mar 4 13:14:41 Coordinated Universal Time 2021
