Next Article in Journal
A Scoping Review of Capacity-Building Efforts to Address Environmental Justice Concerns
Previous Article in Journal
Risk Factors for Depressive Symptoms in Long-Haul Truck Drivers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Integrated Approach for Spatio-Temporal Cholera Disease Hotspot Relation Mining for Public Health Management in Punjab, Pakistan

Department of Computer and Software Engineering, National University of Sciences and Technology, Islamabad 44000, Pakistan
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(11), 3763; https://doi.org/10.3390/ijerph17113763
Submission received: 24 April 2020 / Revised: 14 May 2020 / Accepted: 18 May 2020 / Published: 26 May 2020

Abstract

:
Public health management can generate actionable results when diseases are studied in context with other candidate factors contributing to disease dynamics. In order to fully understand the interdependent relationships of multiple geospatial features involved in disease dynamics, it is important to construct an effective representation model that is able to reveal the relationship patterns and trends. The purpose of this work is to combine disease incidence spatio-temporal data with other features of interest in a mutlivariate spatio-temporal model for investigating characteristic disease and feature patterns over identified hotspots. We present an integrated approach in the form of a disease management model for analyzing spatio-temporal dynamics of disease in connection with other determinants. Our approach aligns spatio-temporal profiles of disease with other driving factors in public health context to identify hotspots and patterns of disease and features of interest in the identified locations. We evaluate our model against cholera disease outbreaks from 2015–2019 in Punjab province of Pakistan. The experimental results showed that the presented model effectively address the complex dynamics of disease incidences in the presence of other features of interest over a geographic area representing populations and sub populations during a given time. The presented methodology provides an effective mechanism for identifying disease hotspots in multiple dimensions and relation between the hotspots for cost-effective and optimal resource allocation as well as a sound reference for further predictive and forecasting analysis.

1. Introduction

With big data in healthcare as well as other domains being generously generated, it presents endless possibilities for creating meaningful public health analytical models. Big data in health comes from multiple health-related interdependent or independent sources including Electronic Health Records (EHR) or Electronic Medical Records (EMR), laboratory tests and results, pharmaceuticals, prescriptions, medical imaging, genomics, etc. [1]. The data acquired are examined for subsequent analysis using predictive modeling and analysis [2,3]. The research investigations and their results have tremendously affected the health care sectors including clinical settings, medical services, pharmaceuticals, public health and other related segments [4,5,6].
In order to meaningfully utilize and manipulate big data in health, it is important to channel its flow from the source to systems where it can be managed, analyzed and used for analytical purposes. With an increase in data analysis, storage and dissemination capabilities, a framework used to acquire, integrate, analyze and distribute data at multiple scales can be implemented. Such a framework is able to support top-down analytical approach for example, from the individual to the gene level, as well as the bottom-up approach, for example, data aggregation from individual to population level. In addition, using the breadth of evidence base will allow representing populations from individual cases and a shift from passive to active surveillance systems in real-time as presented in our prior work. The data are obtained from multiple heterogeneous resources in a standardized way ensuring privacy and security while making it accessible to multiple facilities supporting research and evidence-based public health policy and intervention program designing [7]. Furthermore, there are various types and sources on population data that can be effectively analyzed for meaningful patterns for public health intervention programs while staying privacy-compliant [8,9].
Traditionally, public health data acquisition and analysis use aggregated data. However, many infectious diseases, including cholera, have complex driving factors involved. Given the high variability of disease dynamics with respect to environmental and other driving factors, it is important to study individual case incidences as the evidence base for the population [10]. For instance, in the case of disease management, instead of analyzing a generalized view of disease distribution, we need to identify the time and location where incidences create a dense hotspot or a meaningful pattern.
In this paper, we present a mathematical model to represent diseases and associated factors and operations applied to them to study disease patterns. The purpose of our model is to elicit aspects of a disease in a comprehensible, quantitative and generalized manner. This generalized representation will help to investigate the quantitative dynamics and relationships among diseases and their features. In particular, we study the cholera dynamics in the Punjab province of Pakistan using point process representation and apply data mining techniques to identify cholera hotspots and relationships among the identified hotspots based on selected features.
The rest of the paper is organized as follows. Section 2 gives the related work conducted in a similar manner. Section 3 describes the disease management model, data representation and operations performed to obtain the desired results. The section also outlines the data used to validate the model in the declared study region. Section 4 gives the results and discussion of the results when the proposed methodology is applied to the dataset of cholera incidences.

2. Related Work

Mathematical models to represent diseases date back to 1909 and are still widely used to represent disease transmission [11,12,13,14]. Mathematical modeling represents elements of disease transmission dynamics and other candidate factors of interest. The mathematical model construction has allowed applications of solutions from other domains such as network theories [15,16,17], artificial intelligence [18,19], machine learning [20,21,22] and game theories [23,24,25] to disease dynamics. Due to a strong relationship between cholera and other environmental factors, cholera transmission models frequently include other features of interest [26,27,28].
Stochastic or point processes have been effectively applied in studies that involve disease analysis in the population or sub-population in a region. These processes use two or three-dimensional space to represent disease incidences and thus, constitute a spatial study of the incidences [29,30,31,32]. The spatial and spatial-temporal models are different from each other in terms of time. Introducing the element of time into spatial analysis constitutes a spatio-temporal analysis. Spatial analysis is performed using two or three-dimensional space represented as ( x , y ) or ( x , y , z ) respectively. The time element can be introduced as third dimension z, that is, ( x , y , t ) or as an additional data feature, that is, [ ( x , y ) + f ) ] where f is the datetime feature [33]. These representations are frequently utilized in multiple analytical modeling methods including spatial filtering, Bayesian method, cluster identification and regression modeling techniques [34,35,36].
Significant research has been conducted in spatio-temporal analysis for multiple domains including sociology, criminology, economics, biology, health and environmental sciences etc. [37,38,39,40,41,42,43,44]. In the context of health data clustering, using the above-mentioned mechanisms can greatly affect the resulting clusters since the process of clustering is sensitive to space and time scales [45]. Different space and temporal data ranges can result in a different set of clusters. Creating spatial hotspots is a relatively straightforward task with many successful statistical techniques available including Getis–Ord Gi* [46], local Moran’s I [47], score statistic [48], and independent component analysis [49] etc. In addition to statistical techniques, data mining attempting to uncover previously unknown interesting patterns in the disease data are also applied. These techniques, although based on statistical foundations, integrate data mining including Machine Learning (ML) and Artificial Intelligence (AI) based methods [50,51,52]. However, when there are other factors affecting the disease incidence, it becomes a complex phenomenon [53]. The candidate driving factors can be incorporated into the model as a third feature or a set of features. For example, factors such as population, socio-economic status, education, water sanitation index, pollution index can be incorporated in the framework one at a time. Combining all candidate features of interest in disease incidence representation for spatio-temoral analysis is not as straightforward and requires a systematic approach towards modeling and analyzing the complex phenomenon.
We introduce a layer-based approach for similar hotspots cluster detection or hotspot relationship mining based on correlation identified factors. We use a multi-step process involving partition-based hotspot detection and distance-based similarity approach for mining relationships between the hotspots. Our approach allows identification of similar hotspots based on multiple determinants or factors that are related to each other in a spatio-temporal region. The model presented contributes towards disease modeling and understanding its spatial, temporal dynamics as well as dynamics of other candidate factors driving or contributing to disease spread and their patterns over time. Specifically, the case study presented describes the mechanism for hotspot detection and their evolution over time given a geographic area. In addition, the work contributes towards assessing the similarity of different hotspots based on selected features of interest. This can lead to actionable information including simulation of hotspots trajectories to evaluate the impact of public health interventions given different spatio-tempoiral scales.

3. Methods and Materials

3.1. Study Area

This study was conducted for the cholera outbreak during 2015 to 2019 in the Punjab province of Pakistan. Punjab is the most populated province of Pakistan with a population of over 110 million and a population density of 536 persons per square kilometers shown in Figure 1. The urban population is about 40 million and rural up to 70 million according to Pakistan Bureau of the Statistics Census Results 2017 [54]. The varying dynamics of climatic conditions and an increasing urbanization rate of the province have put a burden on its water resources and has resulted in increased water-related and water-borne diseases including cholera. In addition, Punjab has a complex network of rivers, lakes and water bodies and has a history of flood incidences spawned over years [55].

3.2. Dataset

Cholera remains a global threat for public health. It is caused due to contaminated water or food that results in an acute diarrhea and can lead to death within hours if not taken seriously. It has been estimated that up to four million of the world population suffer from cholera annually and up to 143 thousand people die [56]. The data set for cholera incidents in Punjab, Pakistan is obtained from the passive disease surveillance system operated by Punjab Information Technology Board (PITB) [57] where data is reported through all levels of health authorities in Punjab as well as from online portal, mobile application and WHO reports. The data consist of disease name, patient location, date and time for disease incidences from 2015–2019 over multiple locations over all the districts in Punjab province. There were 40,332 cases reported during the time frame. We applied pre-processing steps to the data including, removal of records that had any of the three attributes missing, since they were necessary for spatio-temporal analysis. 35,720 records are obtained based on the pre-processing steps.
The population data were obtained from Federal Bureau of Statistics in Pakistan [54], whereas features of interest including weather attributes are obtained from World Weather Online [58] based on one reading per day, that is, an interval of 24 h using WorldWeatherOnline historical weather data API wrapper v 0.0.4 into pandas v.0.25.1 frame and csv files in python version 5.7.8. The data is downloaded through the API using latitude and longitude of each geocoded disease incidence. Data from the year 2015 to 2019 were analyzed. The presented modeling framework described in Section 3.3 is applied to the case reports from 2015-2019 to underlying incidence cases from the study region of 38 districts.

3.3. Disease Management Model

Mathematical models for disease and epidemiology representation help to study disease infection and transmission dynamics. A sub-population with a single disease infection transmitting to a secondary group of individuals can be represented using Poisson distribution [59,60,61]. However, with increasingly changing infectious disease dynamics, it is apparent that many external factors including environmental attributes, climatic conditions, population demographics, industrialization, living standards and conditions, education and migration etc affect the disease dynamics in a population and make up the determinants of diseases. The relationship between disease incidences and other features of interest or determinants is not, however, very straight forward. Active research is being conducted to evaluate the causal and effectual relationship between the two entities. An effective model for disease dynamics must be able to represent and link, in a meaningful way, all possible factors of interest in addition to clinically available disease data. We present a modeling and analytical framework for public health disease management based on spatio-temporal disease incidence integrated with candidate features for evidence-based public health decision support. An overview of the proposed model is provided in Figure 2. We present here the mathematical definitions used in the model. Formally, we define F , as a set of all factor types important in disease transmission dynamics ranged over by f and D as a set of all diseases ranged over by d. Table 1 lists variables and their descriptions defined in the model. Each incidence of a disease d is defined as a tuple
d i k = ( x i , y i , t k )
where d i is an instance of disease incidence at a location, defined by a pair of coordinates ( x , y ) , and t is the time of the occurrence incident.

3.3.1. Hotspot Detection

Clustering works as a high performance algorithm for hotspot detection in spatio-temporal data [62,63]. We employ a spatial algorithm for finding dense disease incidence areas called hotspots. We use a partition-based algorithm to find the disease spatial hotspots where a distance function is used to find locations with a high occurrence of diseases called hotspots. The time is scaled based on time step values within the clusters to identifying patterns among clusters. A disease incidence data stream is clustered into n clusters where each cluster C j with j = 1 , , n is defined by a tuple as in Equation (2).
C j = ( i d j , H j , d j , m j ) ,
where i d is the cluster identification, H is the centroid or hotspot in the cluster, d is the subset of d that belongs to cluster C and m j is the number of total disease incidences in cluster C j . Each hotspot is represented by a latitude and longitude, that is,
H = ( x , y ) ,
where x, y are represented in the form of longitude and latitude representing location of the hotspot.

3.3.2. Correlation Based Factor Selection

To find the potential causal relationship between the disease incidence and external factors, we measure the relationship between the factors f F and disease incidence using correlation. The sequence d ( i k ) is a time series sequence at location i and is a mapping function from a set of dates with a unit interval, for example, one day. In this case, the time step L is defined as
L k ( t ) = d ( t k )   o r   L k ( t ) = d ( t + k ) ,
where k = 1 for a day interval and k = 30 for a month interval. This allows finding auto-correlation of d given by (5), as well as cross-correlation of d with F represented as (6) and calculated using (5).
r k = t = k + 1 T ( d t a d ) ( d t k a d ) / ( d t a d ) 2 ,
where T is the time boundary and a d is the mean of the sequence. The factor values are agrregated over time based on the selected k and over space based on the defined hotspot locations. Our goal is to select most significant factors that are affecting the disease burden. For this purpose we find the correlation of each factor independently with the disease burden given by (6)
ρ k = d f i ,
where i = 1 , . . | F | and ⊕ is used represent correlation function. We then select factors contributing towards disease burden based on threshold Δ ρ . Therefore, all factors with correlation greater than Δ ρ k are contributing towards disease incidences as given by (7).
f = σ c F ,
where f is a set of all those factors that satisfy the condition c = ρ k < Δ ρ k

3.3.3. Intgrating Spatio-Temporal and Determinant Factors

Each hotspot H, given by (8), is defined by its location, in terms of ( x , y ) and values of set of features f satisfying a condition given by (7).
H = [ ( x , y ) + ( f 1 , f 2 , . . , f n ) ] ,
where n is the total number of correlated external factors. For example, (8) can represent a cholera disease hotspot at a certain location, ( l o n g , l a t ) with values of temperature, population density, pollution index averaged over L where f 1 = t e m p e r a t u r e _ v a l u e , f 2 = p o p u l a t i o n d e n s i t y _ v a l u e , f 3 = p o l l u t i o n i n d e x _ v a l u e .

3.3.4. Feature Based Hotspot Relation Mining

Hotspot relation mining is based on neighborhood tagging that involves identifying a pattern in the disease data that is based on multiple determinants. Based on L d k ( t ) in (4), we define a sequence of consecutive time slices s 1 , s 2 , , s T . Now that each H is defined in terms of multi-dimensional factors, in addition to space and time, we explore each cluster for the defined time slice and tag neighbors. Note that the definition of the neighborhood includes, x , y , t , + f 1 , f 2 , , f n , that is, the similar intensity of disease cases occurring close to each other geographically and temporally and the features of interest of that area for the given time slice are also similar. In this way, we are able to identify if any given location within a cluster is different from the rest of the areas given all the factor values. This allows us to create a distance matrix and heatmaps for hotspots combining multiple factors of interest that can be analyzed for trend. This constitutes a traditional hotspot analysis but only it is based on more factors than the traditional space and time. For example, at any given location and time, we combine the number of cases, population density, humidity, rainfall and temperature for cholera incidents. We find the similarity between two hotspots based on the distance among the selected factors. In order to avoid skewed results because of different units involved, we calculate the z-score f z for all factor readings as (9)
f z = w f A / S ,
where w f is the value for feature f, A and S are the mean and standard deviation for the factor f. Given two spots H i and H j , we find the distance between the locations based the corresponding features in a time slice
d i s t ( H i , H j ) = p = 1 n ( f p i z f p j z ) 2
This distance allows us to cluster similar locations based on selected features. (10) creates a distance matrix M between all the hotspots with elements M i j representing d i s t ( H i , H j ) . We define a binary distance matrix M b by defining a threshold ρ t h over M to determine whether H i and H j are in the same cluster based on ρ t h .
M b = 0 , if M i j ρ t h or i = j 1 , otherwise
M b is determined for all time slices and added over the sequence to determine the frequency of each feature over the study period. We then apply threshold over the frequency based matrix to identify four categories of similarity among the hotspots. This allows us to define a graph G consisting of non empty set of hotspots represented as node as H ( G ) and a set of edges E represented as pair e = u v where an edge e connects u , v H ( G ) if they are similar and fall in the same similarity group. Therefore, we define the graph G as a set containing an ordered pair H i and H j and a relation R
G = ( ( H i , H j ) , R ) ,
where R = { v e r y s i m i l a r , s i m i l a r , d i f f e r e n t , v e r y d i f f e r e n t }
Therefore, we are able to identify locations that have similar feature patterns when disease intensity was higher. This allows us to find trends of features in similar locations for the predictive analysis. This also enables us to identify if the neighboring hotspots are similar or different from each other when more determinants are taken into account. The model also gives an important reference for optimized resource allocation and predictive analysis for any disease outbreak based on identified features.

4. Results and Discussion

We applied our disease management model to cholera incidence data in 38 districts of Punjab, Pakistan from 2015–2019. We geocode the cholera incident locations to obtain disease incidence dataset using reported incidence location. Each disease incidence is converted to georeferenced point on the map. For identification of high disease incidence hotspots in the province, we used scikit-learn machine learning in Python, a partition-based k-means clustering algorithm and identified 15 locations as cholera hotspots. The weekly distribution of these clusters over time is analyzed. Figure 3 illustrates the incidence rate within each hotspot across the weeks. Figure 4 shows the distribution of disease incidence data in each cluster based on the minimum, first quartile, median, third quartile, and maximum values. We find the cluster centroid location as the hotspot and identify the corresponding tehsil and district based on the coordinates. Table 2 gives each H location’s corresponding tehsil and district with a number of cases, that is, m in the corresponding cluster. Since the analysis is based on the location of reported incidence data, the administrative area can be identified to any granularity level.
In addition, in order to study the trend of each hotspot independently and in relation to other hotspots, each hotspot is analyzed with respect to other hotspots for every week. This allows us to identify the most significant hotspot, significant hotspot, significant cold spot and most significant coldspot every weak among the different hotspots locations. The cholera dataset ranged over 180 weeks. For each week, we classify each hotspot as belonging to the given categories. For example, Figure 5a shows a significant coldspot at week starting from 1 December 2016 with 24 number of cases, at an identified lat long falling under Faisalabad district. Similarly, Figure 5b,c show a significant hotspot in the first week of August in 2016 with 19 number of cases and a most significant cold spot in the first week of April in 2018 with 10 number of cases respectively at Toba Tek Singh district while Figure 5d shows a most significant hotspot at Faisalabad district with 92 number of cases starting in the first week of September in 2017. This creates a pattern of all hotspots over time with a time step equal to a week and defines temporal hotspots for a spatial region as shown in Figure 6. Since the interpretations are more subjective and may not be objectively true, therefore, in order to make maps more meaningful and comprehensible for the decision making personnel, our model must represent the data so that it is objectively true and unbiased to the underlying methodology. The results of the analysis of the weekly hotspots pattern over time in the presented model show that each weak has a different threshold for the number of cases for a hotspot to be assigned to a certain category. Similarly, each hotspot over a period of time, belongs to different categories based on different incidence cases. For example, the hotspot in the Faisalabad district had 92 cases in the first week of September, was a very significant hotspot while the same location also became the most significant hotspot with only 34 cases in the first week of October in 2016. This information is valuable in investigating a location’s trend over time for resource allocation and location dynamics for disease patterns. Furthermore, the analysis can be applied to any spatial and temporal resolutions.
As additional factors of interest in studying cholera disease patterns, we used nine weather attributes including maximum temperature, minimum temperature, temperature, sun hour, ultraviolet (UV) index, humidity, cloud cover, precipitation and heat index. However, any other related features such as distance of hotspot from water bodies, flood-affected areas and sanitation index can also be included in the model to identify similar hotspots. We calculated z-scores corresponding to the probability of distribution across weeks of data among all clusters. We found a symmetrical set of disese incidences. The points in the normal quantile plot lie reasonably close to a straight line. We used the Pearson correlation coefficient since the data followed the normal distribution. We find the weekly correlation of each weather attribute with the disease incidence at each hotspot. The features with correlation confidence of 95% and above are selected for inclusion in further analysis. Precipitation, humidity and cloud cover were selected as three climatic features having a p-value less than 0.05 and correlation confidence greater than 95% (see Figure 7). We also include population density as an external feature for further analysis. We then integrate the determinant factors with disease incidence rate at each hotspot location and find the distance among these locations based on Equations (10) and (11), followed by the summation of the correlation matrices over the weekly range. Figure 8 gives the heat map for distances among the hotspots based on all the candidate factors. Based on these distances, we categorize each pair of locations as different, similar, very similar, and weak similar. Figure 9 shows the final connected pair of hotspots for each category. For example, Table 3 gives a list of 23 location pairs that are very similar to each other in terms of selected features. The tehsil level of these locations was identified, and it comes under the district level in the province. Similarly, 34 location pairs were identified as different, 10 location pairs as similar and 29 location pairs as weak similar based on hotspot relationship mining algorithm.
The analysis of results represented by Figure 7 and Figure 9 show an interesting perspective of the hotspots. There are cities where weather may have played a key role in the spread of Cholera, however, there are others where we find weak or very weak correlation among the disease hot spots with similar weather conditions. This indicates other factors besides weather that must have contributed towards its spread. We have found that many of these hotspots are at locations that are flood-stricken over the selected time-period. The analysis of Figure 9. gives us a reason to look for other factors that may also be included along with traditionally studied factors like floods, construction of major projects, the spread of other diseases, healthcare facilities etc. While other potential factors were not included in this study due to the unavailability of data, the model can take multiple parameters and see their strong and weak-correlation to infer similarities to address the spread. In our case, we have currently taken weather information, but we can easily extend this by considering other factors like floods in the area, socioeconomic conditions and the spread of other diseases. The locations and their corresponding features can be further explored for optimized resource allocation or as a foundation for predictive analysis for forecasting disease patterns in other locations with similar feature patterns.
It is important to allow multiple representations of the analysis to aid the decision-making process. Hotspot analysis presented in this research minimizes the subjectivity present in density-based methods. In addition, the hotspots are statistically significant with respect to the higher number of disease incidences that are based on the z-value. The patterns are identified based on the quantification of confidence. We use z-score for meaningful pattern identification confidence. The quantifiable pattern mining allows important decision-making based on the analysis, such as resource allocation to high disease clusters having a positive impact. When more factors are thrown into the model, the spatial representation does not remain as straightforward as a traditional spatial analysis. The neighboring hotspots in this case may not be in the same spatial neighborhood. The presented model allows a seamless introduction of any number of related features of interest selected through correlation and allows their representation in the model for their inclusion in further analysis.
The presented model framework, however, has its own set of challenges and limitations. First, the identification of hotspots is dependent upon the availability of accurate spatial and temporal data. The data collected in this study were obtained based on hospitals and health care entities under the direct administration of the provincial government. There are many incidences in the region that go unreported, thus affecting the accuracy of the disease distribution model. Second, the incidence addresses and date–time entry method is not standard across all healthcare entities. The obtained data collected needs to be preprocessed for typing or human errors, duplicates and missing information. More than 98% of the dataset needed minor preprocessing since it was obtained through hospital EHR systems and was geocoded successfully to the tehsil and district level. The patient privacy was protected by using minimum tehsil level resolution when identifying hotspots. However, the hotspots identification process used each incidence’s unabstracted location for precise and accurate results.
Infectious public health diseases, including cholera have complex dynamics involving multiple climatic, ecological, social and environmental factors. Our presented work serves as a reference for further investigations involving other determinants of public health interests in the region. For example, with extensive rivers network in the Punjab province, there are multiple flood-stricken areas during the selected time period. Flood data, if supplemented in the model can significantly affect the disease characteristic and pattern mining process results. The availability of relevant data on the required scale can be used by the presented model to accurately provide information for predictive analysis on real-time data to characterize the dynamics of disease in the region.

5. Conclusions

The integration of multiple data sources combined with big data technologies enables the exploration of new data relationships. Healthcare data, including disease incidence data combined with data from other sources such as epidemiology data, geographic data, census data, climate data and social determinants data etc., allows conversion of big health data into smart data using data mining and other artificial intelligence techniques. We present a model combining disease incidence data with features of interest in the public-health context for enriching evidence for smart decision support systems. The data elements and operations performed on them reveal disease hotspots, disease trends and hotspot locations categorized based on similarity on the basis of selected features. In addition, one of the major goals of public health policies is to manage cost by optimally utilizing available resources to high disease density areas. The model allows exploration of identified disease hotspots in multiple dimensions that can be continuously enriched by introducing new factors in the model.

Author Contributions

Conceptualization, F.K., S.A.K.; Methodology, F.K., S.A.K; Formal Analysis, F.K.; Data Curation, F.K.; Writing—Original Draft Preparation, F.K.; Writing—Review and Editing, F.K., W.H.B., I.M.; Visualization, F.K.; Investigation, F.K.; Software, F.K.; Validation F.K.;Administration; S.A.K.; Supervision, S.A.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dash, S.; Shakyawar, S.; Sharma, M.; Kaushik, S. Big data in healthcare: Management, analysis and future prospects. J. Big Data 2019, 6, 54. [Google Scholar] [CrossRef] [Green Version]
  2. Bansal, S.; Chowell, G.; Simonsen, L.; Vespignani, A.; Viboud, C. Big Data for Infectious Disease Surveillance and Modeling. J. Infect. Dis. 2016, 214, S375–S379. [Google Scholar] [CrossRef] [PubMed]
  3. Hay, S.I.; George, D.B.; Moyes, C.L.; Brownstein, J.S. Big Data Opportunities for Global Infectious Disease Surveillance. PLoS Med. 2013, 10, e1001413. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Shin, D.; Biocca, F. Health Experience Model of Personal Informatics: The Case of a Quantified Self. Comput. Hum. Behav. 2017, 69, 62–74. [Google Scholar] [CrossRef]
  5. Raghupathi, W.; Raghupathi, V. Big data analytics in healthcare: Promise and potential. Health Inf. Sci. Syst. 2014, 2, 3. [Google Scholar] [CrossRef]
  6. Shin, D.; Lee, S.; Hwang, Y. How do credibility and utility play in the user experience of health informatics services? Comput. Hum. Behav. 2017, 67, 292–302. [Google Scholar] [CrossRef]
  7. Khalique, F.; Khan, S.; Nosheen, I. PHF-A Framework for Public Health Monitoring, Analytics and Research. IEEE Access 2019, 7, 101309–101326. [Google Scholar] [CrossRef]
  8. Khalique, F.; Khan, S.A.; ul Ain Mubarak, Q.; Safdar, H. Decision Tree-Based Anonymized Electronic Health Record Fusion for Public Health Informatics. In Intelligent Computing; Springer: Cham, Switzerland, 2018. [Google Scholar]
  9. Park, Y.J.; Shin, D. Contextualizing privacy on health-related use of information technology. Comput. Hum. Behav. 2019, 105, 106204. [Google Scholar] [CrossRef]
  10. Pickle, L.W. Spatial Analysis of Disease. In Biostatistical Applications in Cancer Research; Beam, C., Ed.; Springer: Boston, MA, USA, 2002; pp. 113–150. [Google Scholar] [CrossRef]
  11. Chen, D. Modeling the Spread of Infectious Diseases: A Review. In Analyzing and Modeling Spatial and Temporal Dynamics of Infectious Diseases; Chen, D., Moulin, B., Wu, J., Eds.; Wiley: Hoboken, NJ, USA, 2014; pp. 19–42. [Google Scholar] [CrossRef]
  12. Bailey, N.T. The Biomathematics of Malaria; Charles Griffin & Company Ltd.: London, UK, 1982; pp. 93–95. [Google Scholar]
  13. Singh, S.; Shukla, J.B.; Chandra, P. Mathematical modeling and analysis of the spread of carrier dependent infectious diseases: Effects of cumulative density of environmental factors. Int. J. Biomath. 2009, 2, 213–228. [Google Scholar] [CrossRef]
  14. Shin, D.; Zhong, B.; Biocca, F.A. Beyond user experience: What constitutes algorithmic experiences? Int. J. Inf. Manag. 2020, 52, 102061. [Google Scholar] [CrossRef]
  15. Keeling, M.; Eames, K. Networks and Epidemic Models. J. R. Soc. Interface/R. Soc. 2005, 2, 295–307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Schaposnik, L.P.; Zhang, A. Modelling epidemics on d-cliqued graphs. Lett. Biomath. 2018, 5, 49–69. [Google Scholar] [CrossRef] [Green Version]
  17. Seibold, C.; Callender, H.L. Modeling epidemics on a regular tree graph. Lett. Biomath. 2016, 3, 59–74. [Google Scholar] [CrossRef] [Green Version]
  18. Ganasegeran, K.; Abdulrahman, S. Artificial Intelligence Applications in Tracking Health Behaviors During Disease Epidemics. In Human Behaviour Analysis Using Intelligent Systems; Springer: Cham, Switzerland, 2020; pp. 141–155. [Google Scholar] [CrossRef]
  19. Shaban-Nejad, A.; Michalowski, M.; Buckeridge, D. Health intelligence: How Artificial Intelligence Transforms Population and Personalized Health. Nat. Med. 2018. [Google Scholar] [CrossRef]
  20. Uddin, S.; Khan, A.; Hossain, M.E.; Moni, M.A. Comparing different supervised machine learning algorithms for disease prediction. BMC Med. Inform. Decis. Mak. 2019, 19, 1–16. [Google Scholar] [CrossRef]
  21. Tai, A.; Albuquerque, A.F.; Carmona, N.E.; Subramanieapillai, M.; Cha, D.S.; Sheko, M.; Lee, Y.; Mansur, R.B.; McIntyre, R.S. Machine learning and big data: Implications for disease modeling and therapeutic discovery in psychiatry. Artif. Intell. Med. 2019, 99, 101704. [Google Scholar] [CrossRef]
  22. Peiffer-Smadja, N.; Rawson, T.; Ahmad, R.; Buchard, A.; Pantelis, G.; Lescure, F.X.; Birgand, G.; Holmes, A. Machine learning for clinical decision support in infectious diseases: A narrative review of current applications. Clin. Microbiol. Infect. 2020, 26, 584–595. [Google Scholar] [CrossRef]
  23. Reluga, T.C. Game Theory of Social Distancing in Response to an Epidemic. PLoS Comput. Biol. 2010, 6, e1000793. [Google Scholar] [CrossRef]
  24. Poletti, P.; Ajelli, M.; Merler, S. Risk perception and effectiveness of uncoordinated behavioral responses in an emerging epidemic. Math. Biosci. 2012, 238, 80–89. [Google Scholar] [CrossRef]
  25. Chang, S.L.; Piraveenan, M.; Pattison, P.; Prokopenko, M. Game theoretic modelling of infectious disease dynamics and intervention methods: A review. J. Biol. Dyn. 2020, 14, 57–89. [Google Scholar] [CrossRef] [Green Version]
  26. Gaudart, J.; Moore, S.; Rebaudet, S.; Piarroux, M.; Barrais, R.; Boncy, J.; Piarroux, R. Environmental Factors Influencing Epidemic Cholera. Am. J. Trop. Med. Hyg. 2013, 89, 1228–1230. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Xu, M.; Cao, C.; Wang, D.; Kan, B. Identifying Environmental Risk Factors of Cholera in a Coastal Area with Geospatial Technologies. Int. J. Environ. Res. Public Health 2014, 12, 354–370. [Google Scholar] [CrossRef] [PubMed]
  28. Sedas, V. Influence of environmental factors on the presence of Vibrio cholerae in the marine environment: A climate link. J. Infect. Dev. Ctries. 2007, 1, 224–241. [Google Scholar] [CrossRef] [PubMed]
  29. Chen, K.; Liu, X.; Ding, L.; Huang, G.; Li, Z. Spatial Characteristics and Driving Factors of Provincial Wastewater Discharge in China. Int. J. Environ. Res. Public Health 2016, 13, 1221. [Google Scholar] [CrossRef] [Green Version]
  30. Chen, B.; Zhu, Z.; Chen, F.; Zhao, Y.; Qiu, X. Strategically Patrolling in a Chemical Cluster Addressing Gas Pollutants’ Releases through a Game-Theoretic Model. Int. J. Environ. Res. Public Health 2019, 16, 612. [Google Scholar] [CrossRef] [Green Version]
  31. Benes, V.; Bodlák, K.; Jesper, M.; Rasmus, W. A case study on point process modelling in disease mapping. Image Anal. Stereol. 2005, 24, 159–168. [Google Scholar] [CrossRef]
  32. Spatial Point Pattern Analysis. In Applied Spatial Data Analysis with R; Springer: New York, NY, USA, 2008; pp. 155–190. [CrossRef]
  33. Shi, Z.; Pun-Cheng, L.S. Spatiotemporal Data Clustering: A Survey of Methods. ISPRS Int. J. Geo-Inf. 2019, 8, 112. [Google Scholar] [CrossRef] [Green Version]
  34. Ahn, J.; Johnson, T.D.; Bhavnani, D.; Eisenberg, J.N.; Mukherjee, B. A space-time point process model for analyzing and predicting case patterns of diarrheal disease in northwestern Ecuador. Spat. Spatiotemporal Epidemiol. 2014, 9, 23–35. [Google Scholar] [CrossRef] [Green Version]
  35. Barro, A.S.; Kracalik, I.T.; Malania, L.; Tsertsvadze, N.; Manvelyan, J.; Imnadze, P.; Blackburn, J.K. Identifying hotspots of human anthrax transmission using three local clustering techniques. Appl. Geogr. 2015, 60, 29–36. [Google Scholar] [CrossRef]
  36. Chandola, V.; Vatsavai, R.; Kumar, D.; Ganguly, A. Analyzing Big Spatial and Big Spatiotemporal Data: A Case Study of Methods and Applications. Handb. Stat. 2015, 33, 239–258. [Google Scholar] [CrossRef]
  37. Kim, Y.L. Data-driven approach to characterize urban vitality: How spatiotemporal context dynamically defines Seoul’s nighttime. Int. J. Geogr. Inf. Sci. 2020, 34, 1235–1256. [Google Scholar] [CrossRef]
  38. Chen, Y.; Ong, J.H.Y.; Rajarethinam, J.; Yap, G.; Ng, L.C.; Cook, A.R. Neighbourhood level real-time forecasting of dengue cases in tropical urban Singapore. BMC Med. 2018, 16, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Lee, I.; Estivill-Castro, V. Exploration of Massive Crime Data Sets through Data Mining Techniques. Appl. Artif. Intell. 2011, 25, 362–379. [Google Scholar] [CrossRef]
  40. Hu, Y.; Wang, F.; Guin, C.; Zhu, H. A spatio-temporal kernel density estimation framework for predictive crime hotspot mapping and evaluation. Appl. Geogr. 2018, 99, 89–97. [Google Scholar] [CrossRef]
  41. Nakaya, T.; Yano, K. Visualising Crime Clusters in a Space-time Cube: An Exploratory Data-analysis Approach Using Space-time Kernel Density Estimation and Scan Statistics. Trans. GIS 2010, 14, 223–239. [Google Scholar] [CrossRef]
  42. Marco, M.; López-Quílez, A.; Conesa, D.; Gracia, E.; Lila, M. Spatio-Temporal Analysis of Suicide-Related Emergency Calls. Int. J. Environ. Res. Public Health 2017, 14, 735. [Google Scholar] [CrossRef] [Green Version]
  43. Kisa, A.; Network, G.; Alipour, V.; Gad, M.; Rabiee, N.; El Tantawi, M.; Cevik, M.; Banach, M.; Ayanore, M.; Skryabin, V.; et al. Health sector spending and spending on HIV/AIDS, tuberculosis, and malaria, and development assistance for health: Progress towards Sustainable Development Goal 3. Lancet 2020. [Google Scholar] [CrossRef]
  44. Rawassizadeh, R.; Dobbins, C.; Akbari, M.; Pazzani, M. Indexing Multivariate Mobile Data through Spatio-Temporal Event Detection and Clustering. Sensors 2019, 19, 448. [Google Scholar] [CrossRef] [Green Version]
  45. Feng, Y.; Chen, L.; Chen, X. The impact of spatial scale on local Moran’s I clustering of annual fishing effort for Dosidicus gigas offshore Peru. Chin. J. Oceanol. Limnol. 2019, 37, 330–343. [Google Scholar] [CrossRef]
  46. Getis, A.; Ord, J.K. The analysis of spatial association by use of distance statistics. In Perspectives on Spatial Data Analysis; Springer: New York, NY, USA, 2010; pp. 127–145. [Google Scholar]
  47. Goovaerts, P.; Jacquez, G.M. Detection of temporal changes in the spatial distribution of cancer rates using local Moran’s I and geostatistically simulated spatial neutral models. J. Geogr. Syst. 2005, 7, 137–159. [Google Scholar] [CrossRef] [Green Version]
  48. Tang, Z.Z.; Lin, D.Y. MASS: Meta-analysis of score statistics for sequencing studies. Bioinformatics 2013, 29, 1803–1805. [Google Scholar] [CrossRef] [Green Version]
  49. Khedher, L.; Ramírez, J.; Górriz, J.M.; Brahim, A.; Illán, I.A. Independent Component Analysis-Based Classification of Alzheimer’s Disease from Segmented MRI Data. In Artificial Computation in Biology and Medicine; Ferrández Vicente, J.M., Álvarez-Sánchez, J.R., de la Paz López, F., Toledo-Moreo, F.J., Adeli, H., Eds.; Springer: Cham, Switzerland, 2015; pp. 78–87. [Google Scholar]
  50. Shekhar, S.; Jiang, Z.; Ali, R.; Eftelioglu, E.; Tang, X.; Gunturi, V.; Zhou, X. Spatiotemporal Data Mining: A Computational Perspective. ISPRS Int. J. Geo-Inf. 2015, 4, 2306–2338. [Google Scholar] [CrossRef]
  51. Atluri, G.; Karpatne, A.; Kumar, V. Spatio-Temporal Data Mining: A Survey of Problems and Methods. ACM Comput. Surv. 2018, 51, 1–41. [Google Scholar] [CrossRef]
  52. Chretien, J.P.; Swedlow, D.; Eckstrand, I.; George, D.; Johansson, M.; Huffman, R.; Hebbeler, A. Advancing epidemic prediction and forecasting: A new US government initiative. Online J. Public Health Inform. 2015, 7, e13. [Google Scholar] [CrossRef]
  53. He, Z.; Deng, M.; Cai, J.; Xie, Z.; Guan, Q.; Yang, C. Mining spatiotemporal association patterns from complex geographic phenomena. Int. J. Geogr. Inf. Sci. 2020, 34, 1162–1187. [Google Scholar] [CrossRef]
  54. Statistics Division. Government of Pakistan Federal Bureau of Statistics. Pakistan. Available online: http://www.pbs.gov.pk (accessed on 12 September 2019).
  55. The Provincial Disaster Management Authority (PDMA). Available online: http://pdma.gop.pk/floods (accessed on 17 February 2020).
  56. Global Health Observatory (GHO) Data. Available online: https://apps.who.int/gho/data/node.main.175?lang=en (accessed on 15 December 2019).
  57. Punjab Information Technology Board. Digital Punjab: Disease Surveillance System. Available online: https://www.pitb.gov.pk/dss (accessed on 15 September 2019).
  58. WorldWeatherOnline Historical Weather Data API Wrapper. 2019. Available online: https://www.worldweatheronline.com/developer/api/historical-weather-api.aspx (accessed on 15 October 2019).
  59. Ali, M.; Goovaerts, P.; Nazia, N.; Haq, M.Z.; Yunus, M.; Emch, M. Application of Poisson Kriging to the Mapping of Cholera and Dysentery Incidence in an Endemic Area of Bangladesh. Int. J. Health Geogr. 2006, 5, 45. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Yan, P. Distribution Theory, Stochastic Processes and Infectious Disease Modelling. Math. Epidemiol. 2008, 1945, 229–293. [Google Scholar]
  61. Kim, M.; Paini, D.; Jurdak, R. Modeling stochastic processes in disease spread across a heterogeneous social system. Proc. Natl. Acad. Sci. USA 2019, 116, 401–406. [Google Scholar] [CrossRef] [Green Version]
  62. Lawson, A. Hotspot detection and clustering: Ways and means. Environ. Ecol. Stat. 2010, 17, 231–245. [Google Scholar] [CrossRef]
  63. Zheng, C.; Fu, J.; Li, Z.; Lin, G.; Jiang, D.; Zhou, X.n. Spatiotemporal Variation and Hot Spot Detection of Visceral Leishmaniasis Disease in Kashi Prefecture, China. Int. J. Environ. Res. Public Health 2018, 15, 2784. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Punjab, a province of Pakistan, is identified as the study area with population density for different districts. Punjab is the highest populus province of Pakistan with 38 districts.
Figure 1. Punjab, a province of Pakistan, is identified as the study area with population density for different districts. Punjab is the highest populus province of Pakistan with 38 districts.
Ijerph 17 03763 g001
Figure 2. The conceptual model for disease management involving complex interdependent driving factors integrated over spatial and temporal regions. At Layer 1, the data were divided into a series of time slices and each disease incidence is represented in space by (x,y) in corresponding time slice. Layer 2 represents a single time slice, where additional features of interest are added in a given time frame for particular hotspots. Layer 3 presents the trend of hotspots in each hotspot location over all time slices or identification of similar and different hotspots based on all space and time data analysis.
Figure 2. The conceptual model for disease management involving complex interdependent driving factors integrated over spatial and temporal regions. At Layer 1, the data were divided into a series of time slices and each disease incidence is represented in space by (x,y) in corresponding time slice. Layer 2 represents a single time slice, where additional features of interest are added in a given time frame for particular hotspots. Layer 3 presents the trend of hotspots in each hotspot location over all time slices or identification of similar and different hotspots based on all space and time data analysis.
Ijerph 17 03763 g002
Figure 3. Weekly disease incidence for 15 hotspots.
Figure 3. Weekly disease incidence for 15 hotspots.
Ijerph 17 03763 g003
Figure 4. Cholera incidence data distribution in the 15 clusters.
Figure 4. Cholera incidence data distribution in the 15 clusters.
Ijerph 17 03763 g004
Figure 5. Weekly pattern for four categories of hot and coldspots in at all cluster locations. (a) Significant coldspot. (b) Significant hotspot. (c) Most significant coldspot. (d) Most significant hotspot.
Figure 5. Weekly pattern for four categories of hot and coldspots in at all cluster locations. (a) Significant coldspot. (b) Significant hotspot. (c) Most significant coldspot. (d) Most significant hotspot.
Ijerph 17 03763 g005
Figure 6. Hotspot trend in an interactive 3D plot over 180 weeks.
Figure 6. Hotspot trend in an interactive 3D plot over 180 weeks.
Ijerph 17 03763 g006
Figure 7. Weekly correlation confidence for nine weather attributes with cholera disease incidence.
Figure 7. Weekly correlation confidence for nine weather attributes with cholera disease incidence.
Ijerph 17 03763 g007
Figure 8. Heatmap for feature-based distance among all hotspot locations.
Figure 8. Heatmap for feature-based distance among all hotspot locations.
Ijerph 17 03763 g008
Figure 9. Hotspot relation mining based on features and incidence rate. Width of the line shows similar to different areas. Graph shows location pair connected for: (a) Very similar hotspots. (b) Similar hotspots. (c) Weak similar hotspots. (d) Different hotspots.
Figure 9. Hotspot relation mining based on features and incidence rate. Width of the line shows similar to different areas. Graph shows location pair connected for: (a) Very similar hotspots. (b) Similar hotspots. (c) Weak similar hotspots. (d) Different hotspots.
Ijerph 17 03763 g009aIjerph 17 03763 g009b
Table 1. Variables and their descriptions defined in the mathematical model for disease management.
Table 1. Variables and their descriptions defined in the mathematical model for disease management.
VariableDescriptionVariableDescriptionVariableDescription
ddisease incidencemcluster size f z z-score of feature
Ccluster d subset of d that
belongs to C
Mdistance matrix
HhotspotFset of all features M b binary distance matrix
fset of features selected
based on correlation
Ltime stepRrelationship between
hostspots
kinterval lengthAmean of sequenceGgraph connecting hotspots
ρ correlation coeffecientSStandard Varitaioncross correlation function
rautocorrelation σ conditional operator w f value for feature f
Table 2. The 15 hotspots identified with their corresponding tehsil and district as well as the number of cases within each hotspot cluster.
Table 2. The 15 hotspots identified with their corresponding tehsil and district as well as the number of cases within each hotspot cluster.
DistrictTehsilNo. of Cases
BahawalnagarChishtian220
BahawalpurKhairpur Tamewali3266
FaisalabadTandlianwala7257
FaisalabadSammundri2699
JhelumChoa Saidan Shah406
KasurKasur3462
KhanewalAhmadpur Sial321
KhushabKhushab3664
LahoreFerozewala1756
LayyahLayyah416
OkaraOkara481
PakpattanPakpattan123
RawalpindiKallar Sayyedan1936
RawalpindiFateh Jang1091
Toba Tek SinghToba Tek Singh6306
Table 3. Location pairs with very similar conditions identified scaled to the tehsil level.
Table 3. Location pairs with very similar conditions identified scaled to the tehsil level.
S. NoLocation ALocation B
1ChishtianLayyah
2Khairpur TamewaliAhmadpur Sial
3Ahmadpur SialChoa Saidan Shah
4ChishtianKallar Sayyedan
5Toba Tek SinghKhushab
6Toba Tek SinghKallar Sayyedan
7Toba Tek SinghLayyah
8Toba Tek SinghKasur
9KasurKallar Sayyedan
10ChishtianToba Tek Singh
11ChishtianOkara
12LayyahKhushab
13ChishtianKhushab
14KhushabKallar Sayyedan
15LayyahKallar Sayyedan
16PakpattanKhushab
17LayyahKasur
18FerozewalaKallar Sayyedan
19PakpattanLayyah
20ChishtianFateh Jang
21PakpattanToba Tek Singh
22Kallar SayyedanFateh Jang
23OkaraFateh Jang

Share and Cite

MDPI and ACS Style

Khalique, F.; Khan, S.A.; Butt, W.H.; Matloob, I. An Integrated Approach for Spatio-Temporal Cholera Disease Hotspot Relation Mining for Public Health Management in Punjab, Pakistan. Int. J. Environ. Res. Public Health 2020, 17, 3763. https://doi.org/10.3390/ijerph17113763

AMA Style

Khalique F, Khan SA, Butt WH, Matloob I. An Integrated Approach for Spatio-Temporal Cholera Disease Hotspot Relation Mining for Public Health Management in Punjab, Pakistan. International Journal of Environmental Research and Public Health. 2020; 17(11):3763. https://doi.org/10.3390/ijerph17113763

Chicago/Turabian Style

Khalique, Fatima, Shoab Ahmed Khan, Wasi Haider Butt, and Irum Matloob. 2020. "An Integrated Approach for Spatio-Temporal Cholera Disease Hotspot Relation Mining for Public Health Management in Punjab, Pakistan" International Journal of Environmental Research and Public Health 17, no. 11: 3763. https://doi.org/10.3390/ijerph17113763

APA Style

Khalique, F., Khan, S. A., Butt, W. H., & Matloob, I. (2020). An Integrated Approach for Spatio-Temporal Cholera Disease Hotspot Relation Mining for Public Health Management in Punjab, Pakistan. International Journal of Environmental Research and Public Health, 17(11), 3763. https://doi.org/10.3390/ijerph17113763

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop