Skip to main content
  • Research article
  • Open access
  • Published:

Prevalence, patterns, and predictors of diarrhea: a spatial-temporal comprehensive evaluation in India



Spatial analysis has been vital in mapping the spread of diseases and assisting in policy making. Targeting diarrhea transmission hotspots is one of the potential strategies for reducing diarrhea cases. This study aimed to examine the spatial-temporal variations and to identify the modifiable determinants of diarrhea while controlling for the spatial dependence in the data.


An ecological study on diarrhea data from DLHS-3 and NFHS- 4 in India. Moran’s I and LISA were used to detect the spatial clustering of diarrhea cases and to test for clustering in the data. Spatial regression was used to identify the modifiable factors associated with the prevalence of diarrhea. The study comprised of the prevalence of diarrhea among the children below the age of five years (U-5 s) across different states in India. The determinants of diarrhea were obtained using spatial lag models. The software used were GeoDa 1.6.6 and QGIS 2.0.


The presence of spatial autocorrelation in DLHS-3 and NFHS-4 (Moron’s I = 0.577 and 0.369 respectively) enforces the usage of geographical properties while modeling the diarrhea data. The geographic clustering of high-prevalence districts was observed in the state of UP consistently. The spatial pattern of the percentage of children with diarrhea was persistently associated with the household with a sanitation facility (%) (p = 0.023 and p = 0.011). Compared to the diarrhea cases in the period 2007–2008, no much reduction was observed in the period 2015–2016. The prevalence of diarrhea and percentage of household with sanitation were ranging between 0.1–33.8% and 1.3–96.1% in the period 2007–2008 and 0.6–29.1% and 10.4–92.0% in the period 2015–2016 respectively. The least and highest prevalence of diarrhea being consistently from Assam and UP respectively.


Despite improvements in controlling spread of diarrheal disease, the burden remains high. Focus on widespread diarrheal disease control strategy by addressing the social determinants of health like basic sanitation is crucial to reduce the burden of diarrhea among U-5 s in India. The identification of hotspots will aid in the planning of control strategies for goal setting in the targeted regions.

Peer Review reports


The time between birth and age 5 of a child’s life represents a sensitive period for development [1, 2]. The developmental construct of the physical health and well-being of a U-5 includes physical readiness for the school day, physical liberation, and gross and fine motor skills. Diarrhea is a gastrointestinal infection caused by pathogenic microorganisms including E. coli, Rotavirus, Salmonella spp., Vibrio cholerae, Campylobacter, and Shigella [3]. Highly preventable morbidity, diarrhea, is a leading cause of malnutrition and the third leading cause of mortality among U-5’s [4]. Diarrheal illness remained a key threat to global U-5 health and was estimated to be a reason for 15% of all-cause mortality among U-5 s [5]. Globally 1.339 billion U-5 s and 79 thousand neonates aged 0–27 days were estimated to be dead due to diarrhea [6]. India, Nigeria, Afghanistan, Pakistan, and Ethiopia accounts for more than half of the deaths caused by diarrhea [7, 8].

To prevent diarrheal disease and to monitor and control U-5MR, the government of India introduced schemes including child survival and safe motherhood programme in 1992, target-free approach in 1996, reproductive and child health programme-1 in 1997, reproductive and child health programme-2 in 2005, national rural/urban health mission between 2005 and 2012, national health mission since 2013 to-date, the expanded program on immunization, and a program for the control of diarrheal diseases and acute respiratory infection. As a result of various efforts, India has recorded a decline of 52.73% in U-5MR with 55 and 29 deaths per 1000 live births in 2011 and 2015 respectively [9,10,11,12]. India has set its way to achieve the SDG target for under-five child mortality of 25 by 2030 [13, 14]. The average annual rate of decline in U-5MR at 3.1% in India was considered insufficient with the MDG- 4 target to minimize the U-5 mortality to 29/1000 live births by 2015 [15,16,17]. To prevent U-5MR and to achieve the various targets, special attention need to be put on most prevalent diseases among U-5 s including diarrhea. Prevalence of diarrhea is highly heterogeneous across the states of India, and it ranges from 0.1–33.8% in period 2007–2008 and 0.6–29.1% in period 2015–2016 [18, 19]. Though the reduction in the range of diarrhea prevalence over time is observed, generating more evidence through spatial analysis will help target the ones at the rear. Spatial analysis is shown to be important in mapping the extent of infectious diseases and assisting the control policies [20,21,22,23]. Studies examined cultural dimensions of neighborhoods, conveying critical factors that add to the ethnic divisions among spaces in the unequal distribution of developmental health [1]. The spatial techniques help us identify the hotspots and provide the information on significant clusters and the associated factors. The spatial evaluation will enable public health officers and policymakers for strategic planning to reduce the prevalence of diarrhea. More than 1 billion people worldwide are deprived of basic facilities including clean water sources and sanitation [20]. In light of recent developments in India, evidence on diarrheal diseases and its determinants need to be studied for improved planning and need-based organization of health services within the country. Targeting the hotspots may represent an efficacious strategy for further reduction in diarrhea cases. Developing countries like India, with highly heterogeneous diarrhea occurrence, may benefit the most from a targeted community-wide interventional approach that reaches all diarrhea-affected and diarrhea-exposed children. On a priority basis, interventions should be channeled to the most affected districts in India, as revealed by the hotspot analysis. Availability of effective interventions to the locations at the rear can prevent diarrheal occurrence [24].

To the best of our knowledge, the spatial-temporal analysis of diarrhea prevalence in India has not been previously reported. In the present study, the spatiotemporal mappings of regions affected by the pervasive illnesses are reported for instituting the relation between the disease and the dynamics associated with its evolution. The mapping will assist in the formulation and implementation of need-based preventive strategies for the diarrheal disease. This study aims at addressing the social determinants of health to reduce the burden of diarrhea among U-5 s in India.



This ecological study uses the aggregate level data from a large-scale multi-round survey conducted in a representative sample of households in India. The surveys were conducted by various regional agencies and coordinated by the IIPS, Mumbai. In the present study, we use the third round of the DLHS-3 [18] and the fourth round of the NFHS-4 [19] conducted in the period 2007–2008 and 2015–2016 respectively. Data was obtained (from under the National Data Sharing and Accessibility Policy (NDSAP) of Government of India and is freely available and accessible. The data on EAG states, Daman and Diu, Dadra and Nagar Haveli and Lakshadweep, were not available at both the time points and hence are not analyzed. The schematic presentation of data used in the analysis is shown in Fig. 1.

Fig. 1
figure 1

Schematic presentation of the data collected by the national agency and the data used for the analysis in the present study

In addition to the common variables data, DLHS-3 provides data on the percentage of the pucca house [pucca house(%)], women aware about diarrhea [awareness (%)] and households with BPL card [BPL (%)] [18]. The present study utilizes data on the above-listed variables for the analysis. The software used in this study was GeoDa [25] 1.6.6 (Arizona State University, Tempe, AZ, USA) and QGIS 2.0. [26]. The outcome variable considered in this study was the district level prevalence of diarrhea.

Exploratory spatial analysis

Spatial cluster detection was performed to identify the cluster of locations with a higher prevalence of diarrhea for the formulation of policies on MCH. Global spatial autocorrelation Moran’s I was used to measure the overall clustering of the data and to project the strength and pattern of spatial autocorrelation. LISA statistic was used to detect the hot spots of the district-level percentage of children suffering from diarrhea [27]. LISA was used to assess the influence of individual districts on the magnitude of the global statistic and to identify the clusters [28]. A significance map shows the areas with a statistically significant LISA statistic value [29]. Hotspots are the geographical units with high diarrhea prevalence and are surrounded by other geographic units with high diarrhea prevalence. Similarly, a cold spot is a geographic location with low diarrhea prevalence surrounded by other geographic units with low diarrhea prevalence. High-high and low-low districts suggest the clustering of geographies with similar values of diarrhea prevalence whereas the high-low and low-high districts indicate spatial outliers [30]. Queens first order contiguity matrix was used to provide the weights [31]. A p < 0.05 was considered significant throughout.

Confirmatory spatial analysis

To explore the relationship between the percentage of diarrhea and a set of explanatory variables regression techniques were used. The application of appropriate regression technique depends on the violation of assumptions required. The OLS estimation is built on assumptions that random error terms or independent variables are uncorrelated, have a constant variance and follow a normal distribution. The estimates obtained using OLS are biased and inefficient when the above assumptions are violated [32]. Multicollinearity condition number was used to investigate if the independent variables are correlated. The Breusch-Pagan test was used to test for heteroskedasticity. Normality of error distribution was assessed using the Jarque-Bera test of normality of regression disturbances. Spatial regression is a technique used to predict the value of an outcome variable based on values of a set of explanatory variables, taking into account the spatial dependence. Two forms of spatial dependence are a spatial error and spatial lag. The dependency is said to be on the spatial error if the error terms across different spatial units are correlated. The spatial dependence is accounted to spatial lag if the dependent variable y in place i is affected by the independent variables in both places i and j [33,34,35,36]. The spatial error model with the vector of spatially autocorrelated error term (ε), the matrix of weights indicating the spatial framework of neighborhood effects among the residuals (W), the spatial error coefficient (λ), the vector of spatially autocorrelated error term (u), and vector of uncorrelated error terms (ε), is mentioned in eq. (1).

$$ \mathrm{y}\kern0.5em =\kern0.5em \mathrm{X}\upbeta +\mathrm{u}\kern5em \mathrm{where}\kern0.5em \mathrm{u}\kern0.5em =\kern0.5em \uplambda \mathrm{Wu}+\kern0.5em \upvarepsilon $$

If the spatial error coefficient λ = 0, then it is said to exist no spatial correlation between the errors. The spatial lag model for the matrix of weights indicating the spatial framework of neighborhood effects among the dependent variable (W), the matrix of observations on the explanatory variable (X), the vector of error terms (ε), and the parameters β and ρ, is mentioned in eq. (2).

$$ \mathrm{y}\kern0.5em =\kern0.5em \mathrm{XB}\kern0.5em +\kern0.5em \uprho \mathrm{Wy}\kern0.5em +\kern0.5em \upvarepsilon $$

A spatial lag is a variable that averages the neighboring values of a location and controls spatial autocorrelation in the dependent variable. Robust LM(error) and Robust LM(lag) test were used to investigate the spatial dependence. The analysis mentioned above were performed using GeoDA software. The schematic for the choice of an appropriate model for the data is depicted in Fig. 2.

Fig. 2
figure 2

Schematic for the choice of the most suitable model among OLS, spatial lag & spatial error models

Bivariate LISA will be then used to examine the spatial relationship between the exposure and the response variable across the geographic locations of India.


Overall clustering of data

Moron’s scatter plot in Fig. 3 provides a visual representation of data attached to the shapefile and the spatial associations in the neighborhood of each district-level observation. The Global spatial autocorrelation Moran’s I value for the prevalence of diarrhea is observed to decrease significantly over time from 0.579 to 0.369 when investigated using the data from DLHS-3 and NFHS-4 respectively. The result indicates segregation of elements of the significant clusters over time and is suggestive of the need to include a component of geography in the model for identifying the determinants of diarrhea.

Fig. 3
figure 3

The shapefile of India and the Moran’s scatter plot for diarrhea prevalence in India in the period (a) 2007–2008 (b) 2015–2016


Table 1 highlights that the least and the most diarrhea prevalent districts were consistently from the states of Assam and UP, respectively.

Table 1 The prevalence of diarrhea and the most endemic and non-endemic states in India

Hotspot detection

This study revealed the presence of 32 hotspots and 63 coldspots of diarrhea in the study region using the DLHS 3 data. The hotspots were observed to exist in the districts of Haryana, North-west UP, North and East Maharashtra, and North-West Orissa. Coldspots existed in the districts of East Karnataka, North Tamil Nadu, Chhattisgarh, and Assam. NFHS4 hotspots were observed in the Uttarakhand, UP, Madhya Pradesh, Odisha, and Meghalaya. Coldspots were observed in the states of Assam, Andhra Pradesh, Sikkim, and Kerala. The significance map showed locations with significant local Moran statistics as depicted in Fig. 4. The observed blank places on the map of India are due to the unavailability of data for empowered action group states. The list of spatial outliers and the statistically significant clusters which contributed to the global autocorrelation measure are listed in Additional file 1: Tables S1 and S2.

Fig. 4
figure 4

Local indicators of spatial association cluster map and the significance map for diarrhea prevalence in the study region in the year (a) 2007–2008 (b) 2015–2016

Spatial regression

On the data derived from DLHS-3 and NFHS-4 respectively, the diagnostic tests for OLS assumptions suggest avoidable dependencies between the explanatory variable with the multicollinearity condition number 20.403 and 22.624. A violation of assumption on normally distributed errors surveys using the Jarque-Bera test was observed (p < 0.001 & p < 0.001). The Breusch-pagan test indicated the possible existence of heteroskedasticity (p < 0.001 & < 0.001), due to the spatial dependence in the data, justifying the use of spatial regression. The estimates obtained using various models are shown in Table 2.

Table 2 Estimated regression coefficients obtained from OLS, spatial lag and spatial error models

For the model from DLHS-3, it was observed that the spatial error model compared to the spatial lag model had Akaike Information Criterion (1489.62 vs. 1488.24) Schwarz criterion (1513.8 vs. 1515.88). The Robust LM(error) statistic was not significant (p = 0.465) whereas the Robust LM(lag) statistic was significant (p = 0.004). The spatial lag model was observed to have the best fit for the data derived from the DLHS-3 survey. For a unit increase in the awareness on diarrhea among women (%), improved water supply (%) and sanitation facility (%), the prevalence of diarrhea is expected to decrease by 0.089, 0.036 and 0.048% respectively. With a unit increase in illiterate (%) and pucca house, we expect 0.077 and 0.072% increase in the diarrhea prevalence. On the same lines, for the model obtained based on data derived from NFHS-4, it was observed that spatial error model compared to the spatial lag model, had Akaike Information Criterion (3507.93 vs. 3507.15) and Schwarz criterion (3525.52 vs 3529.13). The Robust LM (error) statistic was not significant (p = 0.764) whereas the Robust LM (lag) statistic was significant (p = 0.003). The spatial lag model was observed to have the best fit and is hence used for further discussion. It is to be noted that after a significant eight years gap, improved sanitation facility was observed to have a persistent significant association with the diarrhea prevalence in India. It was observed that with a unit increase in the improved sanitation facility (%), the diarrhea prevalence is expected to decrease by 0.031%.

Table 3 highlights that the districts with the lowest and the highest percentage of household with proper sanitation.

Table 3 The descriptive for sanitation facility in India over the specified period

Bivariate LISA is used to address if the geographic locations underprivileged in terms of sanitation facility were also lacking in terms of diarrhea prevention. The Bivariate LISA map shows the local patterns of spatial correlation between the diarrhea prevalence of a district and the average sanitation facility for its neighbors. The results are presented in Fig. 5.

Fig. 5
figure 5

Bivariate Local indicators of spatial association cluster map and the significance map for the prevalence of diarrhea with the sanitation (%) in the study region in the year (a) 2007–2008 (b) 2015–2016


The present study findings provide valuable insights on factors affecting diarrhea prevalence. The spatial evaluation of evidence on diarrhea conceals the fact that the situation has improved with time. Most of the states have shown good progress in reducing the prevalence of diarrhea among U-5 except UP and Odisha which remained persistent hotspot till 2015–2016. The state of Maharashtra, Punjab, and Haryana had several hotspots in period 2007–2008 which were no more significant in the period 2015–2016. The results signify that the strategic approach of the government is fruitful and the sweats of focusing on low performing states are paying off. These remarkable achievements are the result of various initiatives of the Government [14]. The implemented community intervention National Rural Health Mission between 2005 and 2012 might have assisted in reducing the spatial inequalities in diarrhea indicators. The observed reduction may also be a result of the state government schemes like Punjab rural water supply and sanitation project [37]. The solid waste disposal concessions awarded to 50 municipal corporations in India including Kolhapur district of Maharashtra would have contributed to decreasing diarrheal disease prevalence in the state of Maharashtra [38]. The impact on improved water and sanitation could be a fruit of several other national and international projects including the second phase of community-led infrastructure finance facility-2015 [39] and sulabh international [40]. Despite Nirmal Bharat Abhiyan and other schemes in the state of UP and Odisha, no significant improvement was noted. In the period 2007–2008, moderate spatial autocorrelation was noted, and the factors having significant impact on diarrhea prevalence were illiteracy and awareness among mothers, use of an improved source of drinking water, sanitation facility, living in the pucca house and having a BPL card. After eight significant years, in the period 2015–2016, a mild spatial autocorrelation was still present, and a household with improved sanitation facility (%) was observed to be solely contributing towards diarrhea prevalence. The results highlight the importance of having a better sanitation facility in every household. Studies [41, 42] with a similar finding report that the mode of water transportation, the existence of excess water in the street, domestic water storage conditions, feces disposal and the presence of vectors predispose the under-five children to diarrhea. About 2.4 billion people lacked access to proper sanitation facility, and four out of five of these underprivileged people lived in Asia alone [38, 43]. The millennium development goal to halve the proportion of people without access to basic sanitation by 2015 [38] seems unmet. To respond to these public health concerns, urban India needs to have universal access to toilets with safe collection and treatment of human excreta [44]. Indiscriminate stool disposal by the mothers, lack of hand-washing before feeding their children and hand-washing without soap need to be addressed [45, 46].

Poor sanitation not only has a direct impact on health but also has indirect contribution to the living situation, reduced education outcomes and poverty [38]. Studies [9] reveal that low socioeconomic status, poor maternal literacy, inadequate breastfeeding, malnutrition, poor sanitation and hygiene practices of the mother are associated with a higher incidence of diarrheal diseases in young children. Poor socioeconomic status is one of the factors influencing the diarrheal incidence [47]. Educational status of the mother showed a positive correlation with the incidence of diarrheal diseases [48].

The present study addressed the situation of diarrhea prevalence using the aggregate level data. The study is beneficial for region-specific allocation of resources. Despite efforts mentioned above, there are few limitations in the study. The unit of analysis is a district and not an individual. Analysing aggregate data lead to the ecological fallacy. The findings cannot be generalized at the individual level. A better identification of risk factors is expected if the meta-data with the names of finer blocks were available. Data on all the states were not available and hence were not included for the analysis. Analysis for the period 2007–2008 included data only on 237 districts compared to the period 2015–2016 data on 599 districts. This may have minor effect on the global measure of spatial autocorrelation Moran’s I, but the local measure LISA and BLISA remains unaffected. Due to the unavailability of data, factors like place of residence, wealth index, parent’s education, final say on health care, and exposure to mass media were not considered for this study. These factors might play an important role, and hence future work may utilize information on these aspect while explaining spatial dissimilarities in the prevalence of diarrhea.


Despite working towards improving the modifiable risk factors, the management of diarrhea remains crucial. Proper sanitation and hygiene will help in substantial reduction of diarrhea cases. The identification of hotspots will assist in the planning of control strategies in the locations at the rear rather than plan for the whole state or nation. In addition to available schemes, customized interventions, region-specific allocation of funds and awareness programs would decrease the diarrhea prevalence, ensuring better child health.



Below Poverty line


District Level Household and Facility Survey


Empowered action group


International Institute of Population Sciences


Local Indicators of Spatial Association


Lagrange Multiplier


Maternal and Child Health


Millennium Development Goal


National Family Health Survey


Ordinary Least Squares


Sustainable Development Goal


Under Five


Under Five mortality rate


Uttar Pradesh


  1. Minh A, Muhajarine N, Janus M, Brownell M, Guhn M. A review of neighborhood effects and early child development: how, where, and for whom, do neighborhoods matter? Health & place. 2017;46:155–74.

    Article  Google Scholar 

  2. Lahariya C, Paul VK. Burden, differentials, and causes of child deaths in India. The Indian J Pediatr. 2010;77(11):1312–21.

    Article  Google Scholar 

  3. Yongsi H. Pathogenic microorganisms associated with childhood diarrhea in low-and-middle income countries: case study of Yaoundé–Cameroon. Int J Environ Res Public Health. 2008;5(4):213–29.

    Article  Google Scholar 

  4. Diarrhoeal disease. WHO 2018.

  5. Cunliffe NA, Witte D, Ngwira BM, Todd S, Bostock NJ, Turner AM, Chimpeni P, Victor JC, Steele AD, Bouckenooghe A. Efficacy of human rotavirus vaccine against severe gastroenteritis in Malawian children in the first two years of life: a randomized, double-blind, placebo controlled trial. Vaccine. 2012;30:A36–43.

    Article  CAS  Google Scholar 

  6. Black RE, Cousens S, Johnson HL, Lawn JE, Rudan I, Bassani DG, Jha P, Campbell H, Walker CF, Cibulskis R. Global, regional, and national causes of child mortality in 2008: a systematic analysis. Lancet. 2010;375(9730):1969–87.

    Article  Google Scholar 

  7. UNICEF: Progress for children: a world fit for children statistical review. In.: Unicef; 2007.

  8. Kammerer K, Higginson A, Hickey P. Recommendations for diarrhea management during medical stability operations. Mil Med. 2012;177(7):870-6.

    Article  Google Scholar 

  9. Lakshminarayanan S, Jayalakshmy R. Diarrheal diseases among children in India: current scenario and future perspectives. J Natural Sci, Biol, Med. 2015;6(1):24.

    Article  Google Scholar 

  10. Bhan MK. Accelerated progress to reduce under-5 mortality in India. Lancet Glob Health. 2013;1(4):e172–3.

    Article  Google Scholar 

  11. Chaikaew N, Tripathi NK, Souris M. Exploring spatial patterns and hotspots of diarrhea in Chiang Mai, Thailand. Int J Health Geogr. 2009;8(1):36.

    Article  Google Scholar 

  12. George CM, Perin J, De Calani KJN, Norman WR, Perry H, Davis TP Jr, Lindquist ED. Risk factors for diarrhea in children under five years of age residing in peri-urban communities in Cochabamba, Bolivia. Am J Trop Med Hyg. 2014;91(6):1190–6.

    Article  Google Scholar 

  13. Goal SD: .

  14. Interventions E. Commodities and guidelines for reproductive, maternal, newborn and child health: a global review of the key interventions related to reproductive, maternal, newborn and child health (RMNCH). Geneva: WHO; 2011.

    Google Scholar 

  15. Greenwell F. MDG 4–5: reduce child mortality and improve maternal health; 2011.

    Google Scholar 

  16. Unies N: Beyond the Midpoint: Achieving the Millennium Development Goals. In.; 2010.

  17. Welfare MoHaF: In.; MoHFW 2018.

  18. Survey-3 DLHF: DLHS-3, 2007–08. In: IIPS/MoHFW. vol. 4; 2010.

  19. Sheet IF. National family health Survey-4. In: International Institute for Population Studies; 2017.

    Google Scholar 

  20. Bessong PO, Odiyo JO, Musekene JN, Tessema A. Spatial distribution of diarrhoea and microbial quality of domestic water during an outbreak of diarrhoea in the Tshikuwi community in Venda, South Africa. J Health Popul Nutr. 2009;27(5):652.

    Article  Google Scholar 

  21. Njemanze PC, Anozie J, Ihenacho JO, Russell MJ, Uwaeziozi AB. Application of risk analysis and geographic information system technologies to the prevention of diarrheal diseases in Nigeria. Am J Trop Med Hyg. 1999;61(3):356–60.

    Article  CAS  Google Scholar 

  22. Fang L, Yan L, Liang S, de Vlas SJ, Feng D, Han X, Zhao W, Xu B, Bian L, Yang H. Spatial analysis of hemorrhagic fever with renal syndrome in China. BMC Infect Dis. 2006;6(1):77.

    Article  Google Scholar 

  23. Sarkar R, Prabhakar AT, Manickam S, Selvapandian D, Raghava MV, Kang G, Balraj V. Epidemiological investigation of an outbreak of acute diarrhoeal disease using geographic information systems. Trans R Soc Trop Med Hyg. 2007;101(6):587–93.

    Article  Google Scholar 

  24. Barros AJ, Victora CG. Measuring coverage in MNCH: determining and interpreting inequalities in coverage of maternal, newborn, and child health interventions. PLoS Med. 2013;10(5):e1001390.

    Article  Google Scholar 

  25. Anselin L, Syabri I, Kho Y. GeoDa: an introduction to spatial data analysis. Geogr Anal. 2006;38(1):5–22.

    Article  Google Scholar 

  26. Graser A: Learning QGIS 2.0: Packt Publishing Ltd; 2013.

  27. Grubesic TH, Murray AT. Detecting hot spots using cluster analysis and GIS. In: Proceedings from the fifth annual international crime mapping research conference, Vol. 26. 2001.

  28. Anselin L. Local indicators of spatial association—LISA. Geogr Anal. 1995;27(2):93–115.

    Article  Google Scholar 

  29. Anselin L: The Moran scatterplot as an ESDA tool to assess local instability in spatial association: regional research institute, West Virginia University Morgantown, WV; 1993.

  30. Pfeiffer D, Robinson TP, Stevenson M, Stevens KB, Rogers DJ, Clements AC. Spatial analysis in epidemiology, vol. 142. New York: Oxford University Press; 2008.

    Book  Google Scholar 

  31. Anselin L, Rey S. Properties of tests for spatial dependence in linear regression models. Geogr Anal. 1991;23(2):112–31.

    Article  Google Scholar 

  32. Hession SL, Moore N. A spatial regression analysis of the influence of topography on monthly rainfall in East Africa. Int J Climatol. 2011;31(10):1440–56.

    Article  Google Scholar 

  33. Wismadi A, Brussel M, Zuidgeest M, Sutomo H, Nugroho LE, van Maarseveen M. Effect of neighbouring village conditions and infrastructure interdependency on economic opportunity: a case study of the Yogyakarta region, Indonesia. Comput Environ Urban Syst. 2012;36(5):371–85.

    Article  Google Scholar 

  34. Anselin L, Bera AK. Spatial dependence in linear regression models with an introduction to spatial econometrics. Statistics Textbooks and Monographs. 1998;155:237–90.

    Google Scholar 

  35. Anselin L. Under the hood: issues in the specification and interpretation of spatial regression models. Agric Econ. 2002;27(3):247–67.

    Article  Google Scholar 

  36. LeSage JP, Altman M, Gill J, McDonald MP. Spatial regression models. Numerical Issues Statist Comp Soc Scientist. 2004:219–37.

  37. Singh J. Punjab rural water supply and sanitation project: problems faced by rwss committees. The Social ION. 2017;6(1):76–87.

    Google Scholar 

  38. Varmani R, Singh AK: Urban Sanitation and Hygiene Promotion in Uttar Pradesh.

  39. Jack M, Morris R: The community led infrastructure finance facility (CLIFF). Article written for the International Association of Local and Regional Development Funds in Emerging markets, drawn from: http://www.citiesallianceorg/doc/resources/financing/cliff/articleIADF (FullVersionFinal)_Jun05.pdf (Last accessed 07/06/08) 2005.

  40. Pathak B. Sulabh sanitation and social reform movement. International NGO Journal. 2011;6(1):014–29.

    Google Scholar 

  41. Simiyu S. Water risk factors pre-disposing the under five children to diarrhoeal morbidity in Mandera district, Kenya. East Afr J Public Health. 2010;7(4).

  42. Kumar S, Subita L. Diarrhoeal diseases in developing countries: a situational analysis. Kathmandu Univ Med J. 2012;38(2):83–8.

    Google Scholar 

  43. WHO U. The global water supply and sanitation assessment. Geneva: World Health Organisation; 2000.

    Google Scholar 

  44. Wankhade K. Urban sanitation in India: key shifts in the national policy frame. Environ Urban. 2015;27(2):555–72.

    Article  Google Scholar 

  45. Jayalakshmy R, Roy G, Premarajan K. Incidence and risk factors of acute diarrheal disease among under-five children in urban slums in Pondicherry z–one year follow-up study. Indian J Matern Child Health. 2011;13:1–11.

    Google Scholar 

  46. Ghosh S, Sengupta P, Mondal S, Banu M, Gupta D, Sircar B. Risk behavioural practices of rural mothers as determinants of childhood diarrhoea. J Commun Dis. 1997;29(1):7–14.

    CAS  PubMed  Google Scholar 

  47. Gupta P, Murali M, Seth A. Epidemiology of diarrhea in urban slums. Indian Pediatr. 1998;35:147–50.

    CAS  PubMed  Google Scholar 

  48. Singh J, Gowriswari D, Chavan B, Patiat R, Debnath A, Jain D, Sharma R, Sharma R, Datta K. Diarrhoeal diseases amongst children under five. A study in rural Alwar. J Commun Dis. 1992;24(3):150–5.

    CAS  PubMed  Google Scholar 

Download references


The authors would like to acknowledge Ms. Amitha Puranik, MAHE, Manipal, for the fruitful discussions. The authors also recognize the support from Mr. Shreenidhi SM, MAHE, Manipal, for extracting the data and Ministry of Health and Family Welfare (MoHFW), Government of India, for making the data freely available for research purposes.



Availability of data and materials

The datasets analyzed for the conclusion made in this current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations



NN and KS conceived and designed the investigation. AK, SNR participated in the modification of the methodology. NN and KS gathered the freely available data and analyzed. NN, AK, SNR, SK, and UB participated in the writing of the manuscript. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Nilima.

Ethics declarations

Competing interest

Dr. Bhaskaran Unnikrishnan, one of the authors of the present study, is an editorial board member of BMC Public Health.

Ethics approval and consent to participate

The study used secondary data available and released under the National Data Sharing and Accessibility Policy (NDSAP) of Government of India. The data set had no identifiable information on the survey participants, therefore no ethical approval is required for this work.

Consent for publication

Not applicable.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Tables S1 and S2. A list of spatial outliers and the most significant clusters are presented in Tables S1 and S2. (DOCX 20 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Nilima, Kamath, A., Shetty, K. et al. Prevalence, patterns, and predictors of diarrhea: a spatial-temporal comprehensive evaluation in India. BMC Public Health 18, 1288 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: