Figures
Abstract
Background
Childhood diarrhea continues to be a public health problem in developing countries, including Ethiopia. Detecting clusters and trends of childhood diarrhea is important to designing effective interventions. Therefore, this study aimed to investigate spatiotemporal clustering and seasonal variability of childhood diarrhea in northwest Ethiopia.
Methods
Retrospective record review of childhood diarrhea was conducted using quarterly reported data to the district health office for the seven years period beginning July 1, 2007. Thirty three districts were included and geo-coded in this study. Spatial, temporal and space-time scan spatial statistics were employed to identify clusters of childhood diarrhea. Smoothing using a moving average was applied to visualize the trends and seasonal pattern of childhood diarrhea. Statistical analyses were performed using Excel® and SaTScan programs. The maps were plotted using ArcGIS 10.0.
Results
Childhood diarrhea in northwest Ethiopia exhibits statistical evidence of spatial, temporal, and spatiotemporal clustering, with seasonal patterns and decreasing temporal trends observed in the study area. A most likely purely spatial cluster was found in the East Gojjam administrative zone of Gozamin district (LLR = 7123.89, p <0.001). The most likely spatiotemporal cluster was detected in all districts of East Gojjam zone and a few districts of the West Gojjam zone (LLR = 24929.90, p<0.001), appearing from July 1, 2009 to June 30, 2011. One high risk period from July 1, 2008 to June 30, 2010 (LLR = 9655.86, p = 0.001) was observed in all districts. Peak childhood diarrhea cases showed a seasonal trend, occurring more frequently from January to March and April to June.
Citation: Azage M, Kumie A, Worku A, Bagtzoglou AC (2015) Childhood Diarrhea Exhibits Spatiotemporal Variation in Northwest Ethiopia: A SaTScan Spatial Statistical Analysis. PLoS ONE 10(12): e0144690. https://doi.org/10.1371/journal.pone.0144690
Editor: Agricola Odoi, University of South Africa (South Africa) and University of Tennessee (USA), UNITED STATES
Received: April 18, 2015; Accepted: November 22, 2015; Published: December 21, 2015
This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication
Data Availability: Due to third party restrictions data are available from the Amhara Regional Health Woreda office. Interested researchers may contact the core process owner to get permission in order to access the data: Endalkachew Desalegn Amhara Regional Health Bureau Research and Technology Transfer Core Process Owner endalk111@gmail.com Mobile phone: +251935466857 Office phone: +251582220191 Bahir Dar, Ethiopia.
Funding: The Addis Ababa University, Bahir Dar University, the University of Connecticut and USAID provided funding to collect data for the study by covering expenses during field visits and for software to analyze the data. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Abbreviations: CSA, Central Statistics Agency; LLA, Log Likelihood Ratio; RR, Relative risk; SSA, Sub-Saharan African; WHO, World Health Organization; WASH, Water, Sanitation and Hygiene
Introduction
The burden of communicable disease continues to be a public health problem in developing countries. The distribution of disease in the population is known to differ among groups, with children being the most vulnerable. Diarrhea, a gastrointestinal infection, is caused by a wide range of pathogens, including bacteria, viruses and protozoa [1, 2], and is often transmitted to humans via food or water that has been contaminated with fecal matter [2]. The disease is a major public health burden throughout the world, particularly in developing countries [3, 4], and children are severely and disproportionately affected. A 2012 review article on low and middle income countries indicated that childhood diarrhea had been reduced only modestly from nearly 1.9 billion episodes in 1990 to nearly 1.7 billion episodes in 2010 [5]. Diarrhea is also a major public health burden in Sub-Saharan Africa (SSA) countries [6], and was responsible for more than 50% of childhood morbidity and 50–80% of childhood mortality in the region [3, 4].
As a Sub-Saharan African country [7], Ethiopia has experienced childhood diarrhea and the associated high morbidity and mortality [4, 8]. Although studies conducted on childhood diarrhea morbidity have shown a slight decline in the past decade, the national prevalence of diarrhea has remained high in comparison to other developed countries. In the two weeks preceding the Ethiopian Demographic Health Survey conducted in 2000, childhood diarrhea prevalence was 24%; in 2005, it was 18%; and in 2011, it was 13.5% [9–11]. Moreover, other community-based studies conducted in different parts of the country at various times have shown the prevalence of diarrheal disease in the range of 11–35% [12–16]. Facility-based reports indicate it has been responsible for 11–18% of hospital admissions from all childhood diseases in the period between 2010 to 2012 [17, 18], and a community-based study in southern Ethiopia revealed 30% of child mortality was caused by diarrhea [19].
Childhood diarrhea incidence varies among countries as well as among the different regions of a given country. These differences can be attributed to climatic, environmental, behavioral, socio-demographic, and economic factors [20–22]. In particular, the disease has been widespread in areas where water resources are scarce and unsafe drinking water supplies, poor hygiene, and lack of sanitation are prevalent [23].
These factors have also favored the occurrence of clusters of childhood diarrhea cases in space and time. The understanding of the epidemiology and the disease pattern in time and space [24] is needed to manage diarrheal diseases. Thus, the identification of geographical areas with ongoing disease transmission, using geographic information systems and spatiotemporal statistical analyses, has become indispensable. However, research findings on diarrhea related diseases in Ethiopia have focused solely on prevalence, characteristics of the individuals affected and factors associated with diarrhea [12–16, 25–30]. The findings of these studies are insufficient and fail to capture transmission dynamics in time and space. For effective disease management, knowing when and where disease prevalence reaches a maximum is essential.
The current study investigated spatial, temporal and spatiotemporal patterns of childhood diarrhea in the Amhara region of northwest Ethiopia using spatial SaTScan statistical analysis. Several techniques for spatial analysis had been used previously in different countries to describe geographical patterns and identify clusters of enteric diseases, including diarrhea, in time and space [31–36]. Spatial SaTScan method was used since it is widely suggested that it performs very well in identifying clusters compared to other spatial analysis techniques [37–39], for two reasons. The first is that SaTScan identifies a cluster at any location of any size up to a set maximum, which limits the problem of multiple statistical tests. Second, spatial scan statistics have higher power than other available methods for detecting local clusters [38, 40].
Materials and Methods
Retrospective data on childhood diarrhea were collected at the district level by reviewing reports made quarterly to the routine surveillance system of the district health office over for a seven year period from July 1, 2007 to June 30, 2014 at the district level.
Study area
The study area was located in Awi, East Gojjam and West Gojjam administrative zones of the Amhara regional state of Ethiopia, which is in the northwestern part of the country between 9°20' and 14°20' North latitude and 36° 20' and 40° 20' East longitude (Fig 1). Spatiotemporal analysis was conducted in 33 districts in this area. Based on the 2007 National Population and Housing Census of Ethiopia, the projected population for 2013/14 in the study area is about 6,102,870, of which 2,998,149 (49.1%) were male and 826,246 (13.5%) lived in urban areas. Of the total population, 708,430 (11.6%) were children under the age of five years [41]. The study area has two rainy seasons: the main one is from June to September, followed by a shorter one from March to May. The dry season ranges from October to February. The mean maximum and minimum temperatures in the study area are 26°C and 11°C, respectively [42].
Data collection and data quality control
The number of childhood diarrhea cases was obtained retrospectively from the surveillance system of district health offices over the period of 1 July 2007 to 30 June 2014 (7 years). From July 1, 2007 to June 30, 2013 the data were reported by health institutions to the district health office quarterly (July to September, October to December, January to March and April to June) and from July 1, 2013 to June 30, 2014 they were reported monthly. The health institutions recorded childhood diarrhea cases based on the recommended World Health Organization (WHO) guidelines that define childhood diarrhea as the passage of loose or watery stools at least three times in a 24-hour period [3].
The projected population data of under-five children in a district for each year were obtained from the Amhara Regional State, Bureau of Finance and Economic Development. The population estimate was made based on the results of the 2007 National Population and Housing Census of Ethiopia from the Central Statistics Agency (CSA) [41]. The population density varied across the studied districts over the course of the study period. The projected population data of under-five children were used as the known underlying population at risk to calculate the annual incidence rate of childhood diarrhea. The district level polygon map was obtained from the CSA and spatial data layers were created in ArcGIS 10.0 [43].
A checklist was prepared based on the health management information system registry format to collect the data from the district health management records. To ensure its consistency, the checklist was first prepared in English, translated into Amharic (local language), and then back into English. Health management information system officers were recruited as data collectors. The data collectors were informed about the project’s objectives and trained in data collection procedures. Every day, after collection, the data were checked for completeness and consistency by the principal investigators of the study.
Data analysis
Trend and seasonal distribution analysis.
The incidence rate for each quarter for each district was calculated and childhood diarrhea incidences by year and by quarter were plotted to observe seasonal, annual fluctuations and trends in the study area. A smoothing technique was used to reduce data variability that might occur in time series as a result of irregularities and seasonality to observe trends and seasonal distribution of childhood diarrhea. Smoothing was performed based on the moving average window method (using a four-quarter window) [44]. The incidence rate of childhood diarrhea was transformed to log childhood diarrhea rate to satisfy the normality assumption for the linear regression model. Each value (log childhood diarrhea rate) in the chronological series was replaced by a weighted average of this value and neighboring ones to serve as a base for visualizing the trend and the seasonal variation. To estimate the contribution of each season, the weighted average was calculated across all seasons.
Spatial distribution of childhood diarrhea.
The distribution of excess hazard was described using the standard morbidity ratio, which was calculated as the ratio of the observed cases over the expected cases of a district. The expected cases were derived from the application of the discrete Poisson model on the purely spatial clusters.
Spatial, temporal, and space–time cluster analysis.
Purely spatial, temporal and spatiotemporal clusters were analyzed using the SaTScan™ spatial statistic developed by Kulldorff [45]. The discrete Poisson model was used as the number of cases of childhood diarrhea in each district was Poisson distributed. Children with diarrhea were taken as cases, and the projected population was the combined number of person-years lived used for fitting the Poisson model [45]. Then, Poisson data were analyzed with the purely spatial, the purely temporal, and the space-time scan statistics.
Purely spatial clusters.
For detection and analysis of spatial clusters of childhood diarrhea, a purely spatial analysis was conducted using the spatial scan statistic without taking time into account. This spatial statistical analysis method employs the creation of a circular window that scans the entire study area. The radius of the circle varies continuously from zero to a specified maximum size. The maximum-size specified the percentage of the maximum total population at risk within the scanning window. Since researchers recommend the maximum-size be no greater than 50%, that is a reported cluster can contain at most 50% of the total population at risk [40, 46], the maximum cluster size was set to 50% of the population at risk. Thus, numerous overlapping windows of different sizes are generated, which together cover the entire study area. Each circular window is considered a possible candidate cluster. The null hypothesis is that disease risk remains the same inside and outside the scanning window in space, while the alternative hypothesis is that the risk within the window is different from that outside the window. For each circle, the observed cases inside and outside the window are counted and compared to the number of expected cases, as calculated using the Poisson distribution. On this basis, the likelihood ratio within each circle is then calculated. Under the Poisson assumption, the likelihood function for a specific window is proportional to
Where:
- n is the observed number of cases inside the window (district),
- μ is the indirectly age-adjusted expected number of cases within the window (in that district) under the null hypothesis,
- N is the total number of cases in the study area, and
- I is an indicator function that is equal to 1 when the window has more new cases than expected under the null hypothesis and 0 otherwise. It should be noted that n/μ and (N-n)/(N-μ) are proportional to the age-standardized incidence ratios within and outside the window, respectively.
The circle with the maximum likelihood ratio and containing more cases than expected is identified as the most likely (primary) cluster that is least likely to have occurred by chance [45]. The likelihood ratio for this window comprises the maximum likelihood ratio test statistic. The p-value was estimated using Monte Carlo simulations (999) [47], by comparing the rank of the maximum likelihood from the real dataset with the maximum likelihoods from the random datasets [40]. A significance level of alpha < 0.05 was used to test whether the cluster was significant.
Spatiotemporal and temporal clusters.
The space-time scan statistic was employed to detect clusters in both space and time. This helped us to detect clusters that were not detected by the purely spatial statistic. To detect spatiotemporal clusters a cylindrical window was used. The base of the cylinder represents space, as in the purely spatial scan statistic, while the height reflects time. This statistic behaves as before, except the cylindrical window is then moved in both space and time, so that for each possible geographical location and size, it also visits each possible time interval. For this study, the time window was restricted to annual and the spatial window to 50 kilometers.
As in the previous case, a great many overlapping cylinders of different sizes were thus generated, jointly covering the entire study area and the period studied, with each cylinder reflecting a possible spatiotemporal cluster. Districts with a significant number of cases within the corresponding time were identified using a p-value that was determined using Monte Carlo simulations. Secondary clusters in addition to the most likely cluster (primary) were identified using an iterative manner as detailed in Kulldorff [47] for each purely spatial and space-time scan statistics. In the first iteration, only the most likely cluster was reported in the datasets, then this cluster was removed from the datasets, followed by a second iteration and a completely new analysis using the remaining data. This procedure was then repeated until no clusters remained with p-values of less than 0.05. The maximum cluster size was set to 50% of the population at risk. Purely temporal scan statistics used a window that moved in one dimension only using the height of the cylindrical window as the time dimension. A p-value was generated using Monte Carlo simulations in a manner similar to the spatiotemporal clusters. A significance level of alpha < 0.05 was used to identify a significant district.
Statistical analyses were performed and reported using Excel, SaTScan [48] and ArcGIS 10.0. Excel spreadsheet was used to describe data, draw line graphs and apply smoothing techniques. Spatial and spatiotemporal clusters were analyzed using SaTScan programs. The maps were plotted using ArcGis 10.0.
Ethical Considerations
Ethical clearance was approved by the research ethics committee of College of Medicine and Health Sciences, Bahir Dar University. A support letter from this institution was given to the Amhara Regional Health Bureau, which in turn issued a support letter to district health offices and requested cooperation for retrieving retrospective data on childhood diarrhea from records. All the information was kept confidential, and no individual identifiers were collected.
Results
Of the total of thirty six districts in the study area, data were obtained from thirty three districts. Three districts were excluded because of the data incompleteness. Fifteen of the districts were in the East Gojjam administrative zone, eleven in West Gojjam, and the remaining seven in Awi. About 740,764 childhood diarrhea cases were reported from July 1, 2007 to June 30, 2014.
Distribution and seasonal variation of childhood diarrhea
Childhood diarrhea during the study period showed a declining trend. A seasonal variation of childhood diarrhea, which peaked from January to March and April to June, was also observed. The result of smoothing using moving average revealed that, in the period between January to March and April to June, 6% and 3% excess childhood diarrhea cases, respectively, were reported from the base of moving average (Fig 2 and S1 File).
The overall average cumulative annual incidence rate during the study period was 155.3 per 1,000 population at risk. The highest incidence (546.8 per 1,000 population at risk) was reported in Dejen from July 1, 2008 to June 30, 2009, while the lowest occurred in Wenberma district (40.5 per 1,000 population at risk) from July 1, 2012 to June 30, 2013. The distribution pattern of childhood diarrhea varied greatly across the study districts. Most had higher incidence rates, 100 or more cases per 1,000 population at risk (Table 1).
Spatial distribution of childhood diarrhea
The excess hazard map shows the distribution of excess risk, which was defined as the ratio of the number of observed over the number of expected cases. Areas in blue (2 districts) and yellow (14 districts) had a lower risk than expected, as indicated by the excess risk of values of <1. In contrast, districts colored dark red (3) and red (14) had higher than expected or excess risk (>1). Of the 17 districts with excess risk, 11 were located in the East Gojjam administrative zone, 2 in Awi zone, and 4 were in West Gojjam (Fig 3 and S2 File).
Dark red color indicates districts with standard morbidity ratio (SMR) between 1.51 and 2.00, light red indicates SMR between 1.01 and 1.50, yellow color indicate SMR between 0.51 and 1.00, and blue color indicate SMR < = 0.50.
Purely spatial clusters of childhood diarrhea
The spatial pattern of childhood diarrhea was found to be nonrandom. A total of 12 clusters (1 most likely cluster and 11 secondary clusters) were identified throughout the study period, 8 of which were located in East Gojjam, and 3 in West Gojjam. The remaining 1 cluster was in the Awi administrative zone. The most likely cluster was observed in the East Gojjam zone of Gozamin district. The cluster window was centered at 10.349433 N, 37.729800 E. The base diameter of the window was 0 km, with a Relative Risk (RR) of 1.85. The Log Likelihood Ratio (LLR) for the most likely cluster was 7123.89, p <0.001. Secondary clusters were located in Aneded, Debere Elias/Mechkel, Dejen, Enarj Enawuga/Enemay, Enebise SarMidir, Gonchasiso Enese, Huletej Enese, and, all of which are found in the East Gojjam zone; Bure, Jebitenan, and South Achefer in the West Gojjam zone; and Gunagua/Jawi in the Awi zone (Table 2 and Fig 4).
Numerical identification of the clusters are in order of their likelihood ratio; Dark red color indicates the cluster with the highest likelihood ratio and labeled cluster 1 (most likely cluster or primary cluster) while cluster 2 (light red color), cluster 3 (flame red color) and cluster 4 (green), cluster 5 (light blue), cluster 6 (dark blue), cluster 7 (violet), cluster 8 (pink), cluster 9 (light apple), cluster 10 (gold), cluster 11 (rose) and cluster 12 (yellow) are secondary clusters from the highest to lowest likelihood ratio. Beige color indicates none cluster districs. For more detailed cluster information, refer to Table 2.
Spatiotemporal clusters of childhood diarrhea
The most likely spatiotemporal cluster was located in all districts of the East Gojjam administrative zone (Aneded, Awable, Basoliben, Bibugn, Debre Elias, Dejen, Enarj Enawuga, Enebise SarMidir, Enemay, Gonchasiso Enese, Gozamin, Huletej Enese, Mechkel, Senan, and Shebel Berenta,) and three districts of the West Gojjam zone (Dega Damot, Dembecha, and Jebitenan) (LLR = 24929.90, p<0.001) from July 1, 2009 to June 30, 2011. Secondary clusters were identified in Jawi (LLR = 1015.63, p<0.001), and Dangila/South Achefer (LLR = 685.18, p<0.001), from July 1, 2009 to June 30, 2010 and Quarit (LLR = 815.70, p<0.001), from July 1, 2008 to June 30, 2009 (Table 3).
Discussion
The trend of childhood diarrhea showed a declining pattern in space and time in northwest Ethiopia during the study period. This might be due to Water, Sanitation and Hygiene (WASH) interventions carried out by the government and non-governmental organizations [49, 50]. Interventions with WASH seem well progressing in Ethiopia in all these areas in the past decade. The majority of the population, however, still did not have access to improved water supply (51%) [9] and did not use improved sanitation (79%) in 2011 [51]. In spite of the slight decline, the annual incidence rate of childhood diarrhea relatively remained to be high in most districts. Such a high rate may be attributable to many factors in addition to the use of unimproved water supply and unimproved sanitation, such as poor hand hygiene, poor food hygiene, mothers’ limited education, improper house hygiene, and malnutrition [13, 22, 52–55].
A trend of seasonality in peak childhood diarrhea cases, which occurred more frequently from January to March and April to June, was also observed. In the study area, the period between January and March is a dry season and April to June is the beginning of the rainy season, both of which exhibit high temperatures compared to other months [56]. This finding is consistent with previous studies [21, 57–62] that showed climate factors significantly affected seasonal childhood diarrhea, with more cases observed during the periods of high temperature and drier months of the year, for a number of reasons. Increased ambient temperature favors the growth of bacterial and parasitic diarrhea and extends the survival of enterogastritis-causing bacteria in the external environment [63–66]. Environmental conditions such as surface run off and flooding increase the rate of contamination of drinking water sources at the beginning of the rainy season, as the first rainfall events after the dry season wash contaminants from the surface into the water sources [67]. The practice of open defecation is a common practice in Ethiopia, including the study area [51]. A growing concern of accessing to unsafe drinking water was found that 60% of springs and 45% of well water had excessive E.coli in a national study [68]. The dry season indirectly creates temporal changes in human behavior, such as increased consumption of water because of higher ambient temperatures. Drinking water scarcity may lead to a decline in hygiene practices, which in turn promotes the transmission of diarrhea [57, 67]. During the dry season, most hand dug wells and springs are out of commission, which forces households to use unimproved water sources such as rivers and ponds [67].
Another important finding of this study was that the spatial, temporal, and space-time distribution patterns of childhood diarrhea cases in the study districts were nonrandom. The purely spatial model identified one most likely cluster in the East Gojjam Zone. The spatiotemporal models indicated that the most likely cluster was located in all districts of East Gojjam zone and in three districts of the West Gojjam zone that neighbor East Gojjam. Moreover, the East Gojjam zone also contained seven secondary clusters identified by the spatial model and two secondary clusters in neighboring West Gojjam. In short, the study identified clusters that are closely related to a specific geographical area in the same region and share similar geographical parameters, such as altitude and weather conditions, and economic characteristics. Since the importance of cluster analysis in epidemiology is to detect aggregation of disease cases and ultimately to find evidence of risk factors on which prevention and control activities can be focused, the area in which the identified clusters were found should be given priority for childhood prevention and control interventions [69].
The spatiotemporal analysis provided further evidence for a higher than expected number of childhood diarrhea cases arising within a defined place and time period. The most likely spatiotemporal cluster was found in all districts of the East Gojjam zone, and in the Dega Damot, Dembecha, and Jebitenan districts of the West Gojjam zone between July 1, 2009 and June 30, 2011. Secondary space-time clusters were also identified between July 1, 2009 and June 30, 2010 in Dangila, Jawi, and South Achefer districts and in Quarit district on July 1, 2008 to June 30, 2009. The occurrence of clusters before 2012 may be explained in part by the commencement in 2011 of water, sanitation and hygiene interventions with community mobilization by the government [70].
Finally, this study detected significantly higher rates of purely temporal childhood diarrhea cluster from July 1, 2008 to June 30, 2010, in all districts. The reduction of childhood diarrhea following this occurrence of high incidence could also be due to the government interventions addressing drinking water supply and sanitation and the impact of health extension interventions using different approaches to behavioral changes and communication. For instance, recent evidence showed the provision of improved water supply had been increased from 29% in 2000 to 54% in 2011 [9, 10]; the use of improved sanitation had increased from 8% in 2000 to 22% in 2011 [51]; and the prevalence of open defecation had decreased from 61% in 2005 to 39% in 2010 [71]. Recent studies showed significantly greater decreases in childhood diarrhea in model households than in non-model households following intervention with health extension packages [12, 14]. One strategy to implement such packages is to provide every household with a minimum of 72 hours of training to bring positive behavioral change on hygiene and sanitation. This will be an important step toward the maintenance of a healthy environment through the raising of house to house health awareness with active community participation.
A limitation of this study is the incompleteness and non-representativeness of the data on childhood diarrhea. This limitation is believed to be minimal since the drawback of secondary data across all districts is more likely the same. Due to unavailability of data, this study could not test the association between incidence of childhood diarrhea and temperature, rainfall, improved water supply or sanitation. The spatial scan statistics method is highly efficient in detecting local clusters with good accuracy and can help epidemiologists or public health practitioners to evaluate disease clusters and early detection of disease outbreaks. This technique has been used to detect, analyze, and characterize the spatial, temporal, and spatiotemporal patterns of diarrhea clusters in recent years [69, 72–75]. Clusters in this method are assumed to be circular and dissimilar in size to areas within given district boundaries. This could result in the exclusion or inclusion of districts that register excess or less risks. The cluster analysis can, however, provide valuable information about spatial disparity of childhood diarrhea that may be relevant for further study in the study area. The results of this study will allow researchers to generate a hypothesis as to the contextual environmental or socioeconomic factors that might be influencing the distribution of childhood diarrhea. The identification and understanding of clusters in space and time at the village and agro-ecology are also areas for future research.
Conclusions
This study showed statistical evidence of spatial, temporal and spatiotemporal clustering, as well as seasonal patterns and decreasing temporal trends, in the distribution of childhood diarrhea cases in the study area. The advanced spatial statistical methods used constitute complementary tools for epidemiological surveillance of infectious diseases to guide prevention and control efforts. Thus, analyzing surveillance data using such techniques can provide useful information for public health planners in decision making by identifying underserved populations and allocating scarce resources for appropriate locations to implement specific interventions such as provision of safe drinking water supply, promotion of proper environmental and house hygienic practices, proper food hygienic and water handling practices. Taking into account the spatiotemporal variation in the study areas is highly recommended when striving to design an effective and efficient childhood diarrheal disease intervention strategies.
Supporting Information
S1 File. Smoothing using moving averages to visualize the trend of childhood diarrhea.
https://doi.org/10.1371/journal.pone.0144690.s001
(XLS)
S2 File. Observed, expected cases and excess risk in each district.
https://doi.org/10.1371/journal.pone.0144690.s002
(XLS)
Acknowledgments
We would like to thank the Awi and East and West Gojjam zonal health offices and the district health office for the dedication and commitment in providing information. We also thank the Ethiopian Central Statistics Agency for providing shaped files for the study areas. We are very grateful to Addis Ababa University, Bahir Dar University, the University of Connecticut and USAID for their technical support.
Author Contributions
Conceived and designed the experiments: MA AK AW AB. Performed the experiments: MA AK AW AB. Analyzed the data: MA AK AW AB. Contributed reagents/materials/analysis tools: MA AK AW AB. Wrote the paper: MA AK AW AB. Was involved in proposal development, data collection, data management, analysis, and write up of the article: MA. Were involved in proposal development, data analysis, and writing and editing the manuscript: AK AW AB. Read and approved the final manuscript: MA AK AW AB.
References
- 1. Fischer Walker CL, Sack D, Black RE. Etiology of diarrhea in older children, adolescents and adults: a systematic review. PLoS neglected tropical diseases. 2010;4(8):e768. Epub 2010/08/07. pmid:20689809; PubMed Central PMCID: PMC2914743.
- 2. Hasan KZ, Pathela P, Alam K, Podder G, Faruque SM, Roy E, et al. Aetiology of diarrhoea in a birth cohort of children aged 0–2 year(s) in rural Mirzapur, Bangladesh. J Health Popul Nutr. 2006;24(1):25–35. pmid:16796147
- 3.
UNICEF/WHO. Diarrhoea:Why children are still dying and what can be done. World Health Organization; 2009.
- 4. Walker CL, Rudan I, Liu L, Nair H, Theodoratou E, Bhutta ZA, et al. Global burden of childhood pneumonia and diarrhoea. Lancet. 2013;381(9875):1405–16. pmid:23582727
- 5. Fischer Walker CL, Perin J, Aryee MJ, Boschi-Pinto C, Black RE. Diarrhea incidence in low- and middle-income countries in 1990 and 2010: a systematic review. BMC public health. 2012;12:220. Epub 2012/03/23. pmid:22436130; PubMed Central PMCID: PMC3323412.
- 6. Lozano R, Naghavi M, Foreman K, Lim S, Shibuya K, Aboyans V, et al. Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012;380(9859):2095–128. pmid:23245604
- 7.
Collins RO, Burns JM. A History of Sub-Saharan Africa. New York City: Cambridge University Press; 2007.
- 8. Liu L, Johnson HL, Cousens S, Perin J, Scott S, Lawn JE, et al. Global, regional, and national causes of child mortality: an updated systematic analysis for 2010 with time trends since 2000. Lancet. 2012;379(9832):2151–61. pmid:22579125
- 9.
International CSAEaI. Ethiopia Demographic and Health Survey 2011. Addis Ababa, Ethiopia and Calverton, Maryland, USA: Central Statistical Agency and ICF International.: Central Statistical Agency [Ethiopia] and ORC Macro; 2012.
- 10.
Macro CSAEaO. Ethiopian Demographic Health Survey. Addis Ababa, Ethiopia, and Calverton, Maryland,USA:Centeral Statistical Agency and ORC Marco, 2000. Centeral Statistics Authority (Ethiopia) and ORC Macro; 2001.
- 11.
Macro CSAEaO. Ethiopia Demographic and Health Survey 2005. Addis Ababa, Ethiopia, and Calverton, Maryland, USA: Central Statistical Agency and ORC Macro. Central Statistical Agency [Ethiopia] and ORC Macro.; 2006.
- 12. Berhe F, Berhane Y. Under five diarrhea among model household and non model households in Hawassa, South Ethiopia: a comparative cross-sectional community based survey. BMC public health. 2014;14:187. Epub 2014/02/22. pmid:24555706; PubMed Central PMCID: PMCPmc3941798.
- 13. Dessalegn M KA, Tefera W. Predictors of under five childhood Diarrrhea: Mecha district, West Gojjam, Ethiopia. Ethiopian Journal Of Health development. 2011;25(3):192–200.
- 14. Gebru T, Taha M, Kassahun W. Risk factors of diarrhoeal disease in under-five children among health extension model and non-model families in Sheko district rural community, Southwest Ethiopia: comparative cross-sectional study. BMC public health. 2014;14:395. Epub 2014/04/25. pmid:24758243; PubMed Central PMCID: PMCPmc4031974.
- 15. Mengistie B, Berhane Y, Worku A. Prevalence of diarrhea and associated risk factors among children under-five years of age in Eastern Ethiopia: A cross-sectional study. Open Journal of Preventive Medicine. 2013;3(7):446–53.
- 16. Mohammed S.; Tilahun M.; Tamiru D. Morbidity and associated factors of diarrheal diseases among under five children in Arba-Minch district, Southern Ethiopia, 2012. Science Journal of Public Health. 2013;1(2):102–6.
- 17.
FDRE/MoH. Health and Health Related Indicators. Addis Ababa, Ethiopia: Federal Democratic of Republic Ethiopia, Ministry of Health 2010/2011.
- 18.
FDRE/MoH. Health and Health Related Indicators. Addis Ababa, Ethiopia: Federal Democratic of Republic Ethiopia, Ministry of Health; 2011/2012.
- 19. Deribew A TF, Girma B. Determinanets of under five mortality in Gilgel Gibe Field Research Centre, Southwest Ethiopia. Ethiopian Journal Of Health development. 2007;21(2):117–24.
- 20. Alexander KA, Carzolio M, Goodin D, Vance E. Climate change is likely to worsen the public health threat of diarrheal disease in Botswana. Int J Environ Res Public Health. 2013;10(4):1202–30. pmid:23531489
- 21. Bhandari GP, Gurung S, Dhimal M, Bhusal CL. Climate change and occurrence of diarrheal diseases: evolving facts from Nepal. J Nepal Health Res Counc. 2012;10(22):181–6. pmid:23281447
- 22. Cairncross S, Hunt C, Boisson S, Bostoen K, Curtis V, Fung IC, et al. Water, sanitation and hygiene for the prevention of diarrhoea. Int J Epidemiol. 2010;39(1).
- 23.
Prüss-Üstün A CC. Preventing disease through health environments Towards an estimate of the environmental burden of disease. Geneva: World Health Organization; 2006. 2006.
- 24.
Nykiforuk JC, Laura M. Flaman ML. Exploring the utilization of geographic information systems in health promotion and public health. Technical report; #08–001. In: the Centre for Health Promotion Studies, School of Public Health, Alberta Uo, editors. 2008.
- 25. Eshete WB. A stepwise regression analysis on under-five diarrhoael morbidity prevalence in Nekemte town, western Ethiopia: maternal care giving and hygiene behavioral determinants. East Afr J Public Health. 2008;5(3):193–8. pmid:19374323
- 26. Godana W, Mengistie B. Determinants of acute diarrhoea among children under five years of age in Derashe District, Southern Ethiopia. Rural Remote Health. 2013;13(2329):10.
- 27. Mediratta RP, Feleke A, Moulton LH, Yifru S, Sack RB. Risk factors and case management of acute diarrhoea in North Gondar Zone, Ethiopia. J Health Popul Nutr. 2010;28(3):253–63. pmid:20635636
- 28. Mekasha A, Tesfahun A. Determinants of diarrhoeal diseases: a community based study in urban south western Ethiopia. East Afr Med J. 2003;80(2):77–82. pmid:16167720
- 29. Mitike G. Prevalence of acute and persistent diarrhoea in north Gondar zone, Ethiopia. East Afr Med J. 2001;78(8):433–8. pmid:11921568
- 30. Teklemariam S, Getaneh T, Bekele F. Environmental determinants of diarrheal morbidity in under-five children, Keffa-Sheka zone, south west Ethiopia. Ethiop Med J. 2000;38(1):27–34. pmid:11144877
- 31. Borroto RJ, Martinez-Piedra R. Geographical patterns of cholera in Mexico, 1991–1996. Int J Epidemiol. 2000;29(4):764–72. pmid:10922357
- 32. Chaikaew N, Tripathi NK, Souris M. Exploring spatial patterns and hotspots of diarrhea in Chiang Mai, Thailand. Int J Health Geogr. 2009;8(36):8–36.
- 33. Cheng YJ, Tang FY, Bao CJ, Zhu YF, Liang Q, Hu JL, et al. Spatial analyses of typhoid fever in Jiangsu province, People's Republic of China. Geospat Health. 2013;7(2):279–88. pmid:23733290
- 34. Kelly-Hope LA, Alonso WJ, Thiem VD, Anh DD, Canh do G, Lee H, et al. Geographical distribution and risk factors associated with enteric diseases in Vietnam. Am J Trop Med Hyg. 2007;76(4):706–12. pmid:17426175
- 35. Micael Emch MA. Spatial and temporal patterns of diarrheal disease in Matlab,Bangladesh. Environment and Planning. 2001;33:339–50.
- 36. Pardhan-Ali A, Berke O, Wilson J, Edge VL, Furgal C, Reid-Smith R, et al. A spatial and temporal analysis of notifiable gastrointestinal illness in the Northwest Territories, Canada, 1991–2008. Int J Health Geogr. 2012;11(17):11–7.
- 37. Jackson MC, Huang L, Luo J, Hachey M, Feuer E. Comparison of tests for spatial heterogeneity on data with global clustering patterns and outliers. Int J Health Geogr. 2009;8(55):8–55.
- 38. Song C, Kulldorff M. Power evaluation of disease clustering tests. Int J Health Geogr. 2003;2(1):9. pmid:14687424
- 39. Kulldorff M, Tango T, Park PJ. Power comparisons for disease clustering tests. Computational Statistics & Data Analysis. 2003;42(4):665–84. doi: https://doi.org/http://dx.doi.org/10.1016/S0167-9473(02)00160-3.
- 40. Kulldorff M, Nagarwalla N. Spatial disease clusters: detection and inference. Stat Med 1995;14:799–810. pmid:7644860
- 41.
Development BoFaE. Population size of Amhara Region. Bureau of Finance and Economic Development 2014.
- 42. Mellander PE, Gebrehiwot SG, Gardenas AI, Bewket W, Bishop K. Summer rains and dry seasons in the upper Blue Nile Basin: the predictability of half a century of past and future spatiotemporal patterns. PloS one. 2013;8(7):e68461. Epub 2013/07/23. pmid:23869219; PubMed Central PMCID: PMCPmc3711817.
- 43.
ESRI. ArcGIS Desktop: Release 10. Redlands, CA, USA: Environmental Systems Research Institute. 2011.
- 44.
BOX GEP, Jenkins GM, Reninsel GC. Time Series Analysis, Forecasting and Control. Prentice Hall, Englewood Clifs, NJ; 1994. 3rd ed. ed1994.
- 45. Kulldorff M. A spatial scan statistic. Communications in Statistics—Theory and Methods. 1997;26(6):1481–96.
- 46. Elliott P, Wartenberg D. Spatial epidemiology: current approaches and future challenges. Environ Health Perspect. 2004;112(9):998–1006. pmid:15198920
- 47. Kulldorff M. An isotonic spatial scan statistic for geographical disease surveillance. Journal of the National Institute of Public Health. 1999;48:94–101.
- 48.
Kulldorff M. SaTScan TM:Ser Guide for version 9.2. 2013.
- 49.
FMOH. Health Sector Strategic Plan (HSDP-III), Planning and Programming Department. Addis Ababa, Ethiopia: Federal Ministry of Health(FMoH); 2005.
- 50. Alebachew A, Hatt L, Kukla M. Monitoring and evaluating progress towards Universal Health Coverage in Ethiopia. PLoS medicine. 2014;11(9):e1001696. Epub 2014/09/23. pmid:25244146; PubMed Central PMCID: PMCPmc4171462.
- 51.
WHO/UNICEF. Progress on sanitation and drinking-water—2014 update. Geneva2014.
- 52. Checkley W, Epstein LD, Gilman RH, Cabrera L, Black RE. Effects of acute diarrhea on linear growth in Peruvian children. Am J Epidemiol. 2003;157(2):166–75. pmid:12522024
- 53. Clasen TF, Bostoen K, Schmidt WP, Boisson S, Fung IC, Jenkins MW, et al. Interventions to improve disposal of human excreta for preventing diarrhoea. Cochrane Database Syst Rev. 2010;16(6).
- 54. Brooks JT, Shapiro RL, Kumar L, Wells JG, Phillips-Howard PA, Shi YP, et al. Epidemiology of sporadic bloody diarrhea in rural Western Kenya. Am J Trop Med Hyg. 2003;68(6):671–7. pmid:12892051
- 55. Hickman D, Jones MK, Zhu S, Kirkpatrick E, Ostrov DA, Wang X, et al. The effect of malnutrition on norovirus infection. mBio. 2014;5(2):e01032–13. Epub 2014/03/07. pmid:24595373; PubMed Central PMCID: PMCPmc3958801.
- 56. Ayalew D TK, Mamo G, Yitaferu B, Bayu W,. Outlook of future climate in northwestern Ethiopia. Agricultural Sciences 2012;3:608–24.
- 57. Chou WC, Wu JL, Wang YC, Huang H, Sung FC, Chuang CY. Modeling the impact of climate variability on diarrhea-associated diseases in Taiwan (1996–2007). Sci Total Environ. 2010;409(1):43–51. pmid:20947136
- 58. Jagai JS, Sarkar R, Castronovo D, Kattula D, McEntee J, Ward H, et al. Seasonality of rotavirus in South Asia: a meta-analysis approach assessing associations with temperature, precipitation, and vegetation index. PloS one. 2012;7(5):31.
- 59. Hashizume M, Armstrong B, Hajat S, Wagatsuma Y, Faruque AS, Hayashi T, et al. Association between climate variability and hospital visits for non-cholera diarrhoea in Bangladesh: effects and vulnerable groups. Int J Epidemiol. 2007;36(5):1030–7. pmid:17664224
- 60. D'Souza RM, Hall G, Becker NG. Climatic factors associated with hospitalizations for rotavirus diarrhoea in children under 5 years of age. Epidemiol Infect. 2008;136(1):56–64. pmid:17352836
- 61. Checkley W, Epstein LD, Gilman RH, Figueroa D, Cama RI, Patz JA, et al. Effect of El Nino and ambient temperature on hospital admissions for diarrhoeal diseases in Peruvian children. Lancet. 2000;355(9202):442–50. pmid:10841124
- 62. Thomas KM, Charron DF, Waltner-Toews D, Schuster C, Maarouf AR, Holt JD. A role of high impact weather events in waterborne disease outbreaks in Canada, 1975–2001. Int J Environ Health Res. 2006;16(3):167–80. pmid:16611562
- 63. Bentham G, Langford IH. Environmental temperatures and the incidence of food poisoning in England and Wales. International journal of biometeorology. 2001;45(1):22–6. Epub 2001/06/20. pmid:11411411.
- 64. El-Fadel M, Maroun R, Quba'a R, Mawla D, Sayess R, Massoud MA, et al. Determinants of diarrhea prevalence in urban slums: a comparative assessment towards enhanced environmental management. Environ Monit Assess. 2013;29:29.
- 65. Islam MS, Sharker MA, Rheman S, Hossain S, Mahmud ZH, Islam MS, et al. Effects of local climate variability on transmission dynamics of cholera in Matlab, Bangladesh. Transactions of the Royal Society of Tropical Medicine and Hygiene. 2009;103(11):1165–70. Epub 2009/05/30. pmid:19477477.
- 66. Huang D, Guan P, Guo J, Wang P, Zhou B. Investigating the effects of climate variations on bacillary dysentery incidence in northeast China using ridge regression and hierarchical cluster analysis. BMC infectious diseases. 2008;8:130. Epub 2008/09/26. pmid:18816415; PubMed Central PMCID: PMCPmc2569063.
- 67. Hoque BA, Hallman K, Levy J, Bouis H, Ali N, Khan F, et al. Rural drinking water at supply and household levels: quality and management. International journal of hygiene and environmental health. 2006;209(5):451–60. Epub 2006/06/13. pmid:16765086.
- 68.
FMOH, WHO, UNICEF. Rapid assessment of drinking water quality in the Federal Republic of Ethiopia- Coutry Report of the pilot project implementation in 2004–2005. 2007.
- 69. Odoi A, Martin SW, Michel P, Middleton D, Holt J, Wilson J. Investigation of clusters of giardiasis using GIS and a spatial scan statistic. Int J Health Geogr. 2004;3(1):11. Epub 2004/06/05. pmid:15176979; PubMed Central PMCID: PMCPmc436063.
- 70.
FMOH. Health Sector Development Plan IV Midterm Review Report. Addis Ababa, Ethiopia: Federal Ministry of Health (FMOH); 2013.
- 71. Galan DI, Kim SS, Graham JP. Exploring changes in open defecation prevalence in sub-Saharan Africa based on national level indices. BMC public health. 2013;13(527):1471–2458.
- 72. Tang F, Cheng Y, Bao C, Hu J, Liu W, Liang Q, et al. Spatio-temporal trends and risk factors for Shigella from 2001 to 2011 in Jiangsu Province, People's Republic of China. PloS one. 2014;9(1):e83487. Epub 2014/01/15. pmid:24416167; PubMed Central PMCID: PMCPmc3885411.
- 73. Osei FB, Duker AA. Spatial dependency of V. cholera prevalence on open space refuse dumps in Kumasi, Ghana: a spatial statistical modelling. Int J Health Geogr. 2008;7:62. Epub 2008/12/18. pmid:19087235; PubMed Central PMCID: PMCPmc2628349.
- 74. Luquero FJ, Banga CN, Remartinez D, Palma PP, Baron E, Grais RF. Cholera epidemic in Guinea-Bissau (2008): the importance of "place". PloS one. 2011;6(5):e19005. Epub 2011/05/17. pmid:21572530; PubMed Central PMCID: PMCPmc3087718.
- 75. Bompangue Nkoko D, Giraudoux P, Plisnier PD, Tinda AM, Piarroux M, Sudre B, et al. Dynamics of cholera outbreaks in Great Lakes region of Africa, 1978–2008. Emerging infectious diseases. 2011;17(11):2026–34. Epub 2011/11/22. pmid:22099090; PubMed Central PMCID: PMCPmc3310557.