Next Article in Journal
Determination of Propane-1,3-sultone in Workplace Air for Occupational Exposure Assessment
Previous Article in Journal
Physical and Physiological Profiles of Aerobic and Anaerobic Capacities in Young Basketball Players
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Space-Time Clustering Characteristics of Tuberculosis in Khyber Pakhtunkhwa Province, Pakistan, 2015–2019

1
Department of Fundamental & Applied Sciences, Universiti Teknologi PETRONAS, Seri Iskandar 32610, Perak, Malaysia
2
School of Mathematical and Computer Sciences, Heriot-Watt University Malaysia, Putrajaya 62200, Malaysia
3
Center for Applied Intelligent Systems Research (CAISR), Halmstad University, SE-301 18 Halmstad, Sweden
4
Department of Geological Engineering, Universitas Islam Riau, Pekanbaru 28284, Indonesia
5
Department of Statistics, University of Peshawar, Peshawar 25120, Pakistan
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(4), 1413; https://doi.org/10.3390/ijerph17041413
Submission received: 8 December 2019 / Revised: 15 January 2020 / Accepted: 21 January 2020 / Published: 21 February 2020

Abstract

:
The number of tuberculosis (TB) cases in Pakistan ranks fifth in the world. The National TB Control Program (NTP) has recently reported more than 462,920 TB patients in Khyber Pakhtunkhwa province, Pakistan from 2002 to 2017. This study aims to identify spatial and space-time clusters of TB cases in Khyber Pakhtunkhwa province Pakistan during 2015–2019 to design effective interventions. The spatial and space-time cluster analyses were conducted at the district-level based on the reported TB cases from January 2015 to April 2019 using space-time scan statistics (SaTScan). The most likely spatial and space-time clusters were detected in the northern rural part of the province. Additionally, two districts in the west were detected as the secondary space-time clusters. The most likely space-time cluster shows a tendency of spread toward the neighboring districts in the central part, and the most likely spatial cluster shows a tendency of spread toward the neighboring districts in the south. Most of the space-time clusters were detected at the start of the study period 2015–2016. The potential TB clusters in the remote rural part might be associated to the dry–cool climate and lack of access to the healthcare centers in the remote areas.

1. Introduction

Tuberculosis (TB) is an aerosol transmittable disease that is caused by a bacillus mycobacterium and is mainly spread from a TB-infected person to other nearby people through coughing, sneezing, spitting, or talking. TB is still one of the major public health challenges at the global level, particularly in developing countries, and is listed in the top-ten major causes of mortality at the world-level. It is the second fatal disease After HIV/AIDS caused by a single infectious agent [1,2]. The developing countries such as Nigeria, Philippines, Indonesia, Pakistan, India, South Africa, and China highly contribute (above 60%) to total TB burden in the world, where the diagnosis and treatment of TB are difficult to access [3,4].
Pakistan is one of those developing countries where TB is a major public health challenge. Approximately half a million new TB cases are reported each year, and approximately 70,000 people die each year due to TB disease [5]. In terms of TB case burden, Pakistan ranks 5th in the world due to the high occurrence of multidrug-resistant TB (MDR-TB) cases [3,6,7]. Khyber Pakhtunkhwa is the northwestern province of Pakistan, where approximately 55,000 new cases of TB (all types) are reported each year [5]. During 2002–2017, more than 462,920 TB cases were reported in Khyber Pakhtunkhwa province [8].
Several epidemiological studies on TB have been previously conducted in different districts of Khyber Pakhtunkhwa [9,10,11,12,13,14,15,16,17]. These studies focused on the epidemiological characteristics and risk factors analysis of TB within a specific district. However, the space-time cluster analysis of TB incidence at the district-level in the province seems to be lacking. Space-time cluster detection can identify important space-time patterns in TB incidence in Khyber Pakhtunkhwa province, which could be useful for the evidence-based public health interventions to control the TB outbreak [18]. This study aims to identify the space-time clusters of TB cases at the district-level in Khyber Pakhtunkhwa province, Pakistan from January 2015 to April 2019.
A number of methods have been developed for the automatic detection of disease clusters in public health data. Generally, cluster detection methods are classified as either global or local tests. Global clustering methods are used to assess whether a global tendency for the disease to group together is apparent throughout the study region, but it does not identify the location of clusters [19,20,21,22]. These types of methods are appropriate, for example, for finding evidence of whether a disease is infectious or not. Local methods infer the locations and extent of clusters [23,24,25]. This group of methods (which infer the geographical location and time-period of the clusters) can further be divided into two categories [26]: clustering-based methods [27,28,29] and scan statistics methods [23,24,25]. Scan statistics methods are the most widely used methods in public health and epidemiology [30,31,32,33,34,35,36,37].
We used the scan statistics method [23,24,38] to detect the potential space-time clusters of TB cases in Khyber Pakhtunkhwa province, Pakistan. Detecting such clusters helps in identifying the important patterns in the TB incidence, and can thereby provide important information for public health managers to identify their targets of interest for interventions. Moreover, detecting space-time clusters of TB cases assist the epidemiologists in finding the possible environmental or social determinants of TB outbreak in the study area.

2. Materials and Methods

2.1. Study Area

The study was carried out in Khyber Pakhtunkhwa, the northwestern province of Pakistan. The total area of Khyber Pakhtunkhwa province is 74,521 km2 with a census population of 30.52 million according to the census 2017. The land area of the province has 25 districts of different sizes as shown in Figure 1. However, very recently in the year 2019, seven Federally Administered Tribal Areas (FATA) have been merged politically and administratively into the province of Khyber Pakhtunkhwa. The climate in this province is of the tropical monsoon type, but most of the districts are situated beyond the tropical zone with relatively high temperatures and a cool, dry winter that runs from December to February with March to June representing the hot and dry season. Summer extends from July to September and is generally rainy. October and November represent the receding monsoon period. The province has two rainy seasons: March to April and the summer monsoon from July to September.

2.2. Data Collection

The quarterly data on TB cases and population-at-risk in each district were collected from the quarterly reports (2015–2019) of the District Health Information System (DHIS) Khyber Pakhtunkhwa, Pakistan [39]. The District Health Information System (DHIS) collects data on the reported disease cases in each district on a monthly basis. All hospitals in a district report the registered disease cases to the respective DHIS office on a monthly basis. The district offices then send the monthly data to the provincial office, DHIS Cell, Khyber Pakhtunkhwa [40]. The collected dataset comprises the number of quarterly reported TB cases in each district and the corresponding population-at-risk. The collected dataset is available in supplementary information File S1.

2.3. Clusters Detection

The scan statistics software SaTScan™ ver.9.6 [41] with a discrete Poisson probability model was applied under retrospective analysis to detect the spatial and space-time clusters of TB cases at district-level in Khyber Pakhtunkhwa province, Pakistan from January 2015 to April 2019. SaTScan with the discrete Poisson model assumes the number of disease cases in each region is Poisson-distributed with the parameter λ as in Equation (1) [42].
λ = μ ( z ) n G μ ( G )
where µ(z) denotes the population count for area z, μ (G) denotes the total population of the whole study area, and nG is the total observed disease counts in the whole study area. Under the null hypothesis of no cluster exists, the expected number of disease cases in each region is proportional to the population size in that region. SaTScan is an open source software program that was originally developed in [43]. The space-time scan statistics is defined by a cylindrical shaped window with a circular base. The base of the cylindrical window corresponds to the geographical size of the cluster and the height corresponds to the temporal length of the cluster. The circular base is centered around one of the several location-centroids in the study area, with the radius varying continuously from zero to the maximum spatial window size. The height of the cylinder is varying from zero to the maximum temporal window size. The cylindrical window is then moved in space as well as time to check the clusters with all possible spatial and temporal sizes. As a result, a number of overlapping cylinders of various sizes are obtained that jointly cover the entire study region. Each cylindrical window reflects a possible space-time cluster. The likelihood function is maximized over all cylindrical windows, and the window with the maximum log-likelihood ratio (LLR) is assumed to be the most likely cluster, that is, the cluster least likely to be caused by chance, and other windows with a statistically significant LLR were measured as secondary clusters. The statistical significance of each cluster in the study area is based on comparing the likelihood ratio (LLR) against a null distribution achieved from the Monte Carlo simulation. The detail of how SaTScan works is given in the SaTScan user guide [42].
In our study, the maximum spatial window size was set to the default value, i.e.,  ≤ 50% of the total population of the study area and the maximum temporal window size was set to default value, i.e. , ≤ 50% of the study period [42]. A circular window shape was chosen. A Monte-Carlo approach with 9999 repetitions was performed to test the null hypothesis that there was no difference in Relative Risk (RR) between the TB clusters. The clusters with p-value < 0.001 were considered as statistically significant clusters.

3. Results and Discussion

The space-time cluster analysis by SaTScan identified a total of six space-time clusters of TB cases in Khyber Pakhtunkhwa, Pakistan from 1st quarter 2015 to 1st quarter 2019 (Table 1, Figure 2). The most likely cluster was seen in the northern district Upper Dir for the 1st quarter of 2017. The 1st secondary cluster was seen in two northern districts Buner and Swat for the whole year 2015. The 2nd secondary cluster occurred in the western district Bannu for two years 2015–2016. The 3rd secondary cluster appeared in western district Hangu for 2nd quarter 2015 to 1st quarter 2016. The 4th cluster covered the two central districts Noshera and Mardan for the period of one year (from 2nd quarter 2015 to 2nd quarter 2016). All these clusters are statistically significant (p-value < 0.001). The results show that five of the potential clusters have occurred in the years 2015–2016, one in 2017, and none of the clusters was seen in the years 2018–2019. The interval 2015–2016 exhibits the largest number of clusters in the study period. The high TB incidence in the start of the study period (2015–2016) might be due to the lack of access to the healthcare because the number of healthcare centers in Khyber Pakhtunkhwa province was lower (786) in the year 2015 as compared to (1445) in 2018 [44].
The geographical locations of the detected space-time clusters of TB cases in the study area are shown in Figure 2, which indicates some important patterns of TB incidence in the study area. It is obvious from Figure 2 that in the 1st quarter 2015, two potential TB clusters occurred. The one covered the two districts (Swat and Buner) in the north persisted for one year (2015), which shows a tendency of spread toward the neighboring districts in the central part of the province from 2015 to 2016 and then moved again to the north (Upper Dir) in 2017. The other cluster in 1st quarter 2015 appeared in one district in the west persisted for two years (2015–2016), showing no tendency of spread toward the neighboring districts. Such a space-time pattern of the potential cluster provides clues to disease etiology and risk behaviors, suggesting local environmental or social characteristics that promote increased risk. Moreover, the retrospective clusters guide the public health officials on the significance of the control strategies that have been implemented previously in the cluster-regions to control the prevalence.
To know the general pattern of spatial clusters of TB cases over time, we performed spatial cluster analysis for each year (2015–2018) using SaTScan. Five significant clusters were detected in the year 2015, seven in 2016, five in 2017 and four in 2018 (Table 2). The cluster with p-value < 0.001 is considered to be statistically significant. The largest number of spatial clusters were seen in 2016 compared to the other years. The interval 2015–2016 was also found to have a large number of clusters in space-time analysis (Table 1).
The geographical locations of the significant spatial clusters in each year were displayed in Figure 3. In 2015, the most likely cluster was seen in district Swat. In 2016, the same district with one neighboring district (Buner) appeared as the most likely cluster. In 2017, the most spatial cluster was moved from (Swat Buner) to the neighboring district Upper Dir in the west. The most likely cluster was moved back to districts Swat covering additional neighboring districts (Shangla and Battagram) in the South. The two districts (Swat and Buner) were found repeatedly in the potential cluster each year. The district Swat was seen in the most likely cluster in the year 2015, 2016, and 2018 showing it to be the main region of TB outbreak. Moreover, the two western districts Bannu and Hangu were found to be the spatial clusters throughout the study period (2016–2018). These results suggest three district Swat, Shangla, and Battagram for policymakers to be the most targeted districts for possible interventions because this part of the province was seen in the most likely cluster in each year (2016–2018). The two districts Bannu and Hangu are suggested to be the second targeted districts for possible interventions. These two secondary clusters show a tendency of spread toward the central districts in 2018. The TB case burden in the two western districts (Bannu and Hangu) might be due to the presence of the large proportion of Afghan refugees in these districts as evidence from the previous study [45]. Some previous studies on TB prevalence in the individual district have also been identified the evidence of TB outbreak in our resulted most likely cluster regions Swat and Buner [11,12].
It is evident from Figure 2 and Figure 3 that the most likely TB clusters occurred each year in rural hilly areas [46], which might be due to lack of access to healthcare. The population is scattered in these areas, i.e., the population is distributed in small villages that are very far from each other and hence most of the villages are very far from healthcare centers. In addition, the lack of awareness and misconceptions about TB might have caused TB clusters in these rural areas [47]. The snowfall and dry–cold climate in winter may also contribute to high TB cases in these districts.

4. Conclusions

This study provides a good understanding of the space-time clustering characteristics of TB in Khyber Pakhtunkhwa, Pakistan from 2015 to 2019. The most likely and the 1st secondary space-time clusters were seen in the northern part of the province such as (Dir Upper, Swat, and Buner), showing a tendency of spread toward the central part from 2015 to 2016 and then moved toward the north in 2017 (Figure 2). Most of these space-time clusters were seen in the years 2015–2016, which might be associated to the lack of healthcare centers in this period. In addition, the most likely spatial clusters were also seen each year in the northern districts (Dir Upper, Swat, and Shangla), showing a tendency of spread over the years from the north toward the neighboring districts in the south (Figure 3). This study suggests these northern districts for policymakers to be the most targeted regions for possible interventions. This targeted intervention may help control the TB epidemic in the province more effectively. Future research on finding the social and environmental determinants of TB outbreak in Khyber Pakhtunkhwa province is strongly recommended.

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-4601/17/4/1413/s1, File S1: Data on TB cases and population-at-risk in Khyber Pakhtunkhwa, Pakistan (2015–2019).

Author Contributions

H.D. and S.U. conceived the general concept of this research; S.U. designed the models and wrote the paper; S.C.D. and H.F.-T. contributed in the statistical assessment; H.K. and A. reviewed the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the International Grants: Universitas Islam Riau (UIR) Indonesia, (cost center: 015ME0-012) and Murata Science Foundation, Japan (cost center: 015ME0-107).

Acknowledgments

The authors wish to thank the Universiti Teknologi PETRONAS for providing the required facilities for conducting the research.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. World Health Organization. Global Tuberculosis Report. 2019. Available online: https://www.who.int/tb/publications/global_report/en/ (accessed on 13 January 2020).
  2. World Health Organization Media Centre 2019. Tuberculosis: Key Facts. 2019. Available online: https://www.who.int/news-room/fact-sheets/detail/tuberculosis (accessed on 13 January 2020).
  3. World Health Organization Media Centre. 2018. Available online: https://www.who.int/en/news-room/fact-sheets/detail/tuberculosis (accessed on 5 October 2019).
  4. World Health Organization. Tuberculosis Country Profles. 2018. Available online: https//www.who.int/tb/country/data/profiles/en/ (accessed on 11 October 2019).
  5. NTP. National TB Control Program-Pakistan. 2019. Available online: http://ntp.gov.pk/national_data.php (accessed on 1 September 2019).
  6. Pakistan Bureau of Statistics. Pakistan’s 6th Census: Population of Major Cities Census. Available online: http://www.pbscensusgov.pk/ (accessed on 7 September 2019).
  7. World Health Organization. WHO Country Cooperation Strategy at a Glance: Pakistan. 2018. Available online: https://apps.who.int/iris/bitstream/10665/136607/1/ccsbrief_pak_en.pdf (accessed on 7 September 2019).
  8. National TB Control Program Pakistan (NTP). Available online: http://www.ntp.gov.pk/webdatabase.php (accessed on 7 September 2019).
  9. Falahuddin, S.; Khan, L.; Iqbal, M.; Salahuddin, N. Statistical analysis of various risk factors of tuberculosis in district Mardan, Pakistan. J. Biostat. Epidemiol. 2016, 2, 47–51. [Google Scholar]
  10. Khattak, I.; Mushtaq, M.H.; Ahmad, M.U.D.; Khan, M.S.; Haider, J. Zoonotic tuberculosis in occupationally exposed groups in Pakistan. Occup. Med. 2016, 66, 371–376. [Google Scholar] [CrossRef] [PubMed]
  11. Ahmad, H.; Ali, I.; Ahmad, T.; Tufail, M.; Ahmad, K.; Murtaza, B.N. Prevalence of Brucellosis in Human Population of District Swat, Pakistan. Pak. J. Zool. 2017, 49, 371–376. [Google Scholar] [CrossRef]
  12. Akhtar, N.; Saeed, K.; Khan, S.; Rafiq, N. Prevalence of Tuberculosis in District Buner Khyber Pakhtunkhwa, Pakistan. World J. Med. Sci. 2015, 12, 292–296. [Google Scholar]
  13. Ahmad, T.; Jadoon, M.A. Cross Sectional Study of Pulmonary Tuberculosis at Civil Hospital Thana, District Malakand Khyber Pakhtunkhwa Pakistan. World J. Zool. 2015, 10, 161–167. [Google Scholar]
  14. Ahmad, T.; Jadoon, M.A.; Khattak, M.N.K. Prevalence of sputum smear positive pulmonary tuberculosis at Dargai, District Malakand, Pakistan: A four year retrospective study. Egypt. J. Chest Dis. Tuberc. 2016, 65, 461–464. [Google Scholar] [CrossRef]
  15. Noorrahim, M.S.K.; Shahid, M.; Shah, A.; Shah, M.; Rafiullah, H.A. Prevalence of Tuberculosis in Livestock Population of District Charsadda by Tuberculin Skin Test (TST); Akinik Publications: Delhi, India, 2015. [Google Scholar]
  16. Basit, A.; Hussain, M.; Shahid, M.; Ayaz, S.; Rahim, K.; Ahmad, I.; Rehman, A.U.; Hassan, M.F.; Ali, T.; Au, R.; et al. Occurrence and risk factors associated with Mycobacterium tuberculosis and Mycobacterium bovis in milk samples from North East of Pakistan. Pak. Vet. J. 2018, 38, 199–203. [Google Scholar] [CrossRef]
  17. Ahmad, T.; Ahmad, S.; Haroon, M.Z.; Khan, A.; Salman, S.; Khan, N.; Gul, G. Epidemiological study of tuberculosis. Eur. Acad. Res. 2013, 1, 1855–1858. [Google Scholar]
  18. Ullah, S.; Daud, H.; Dass, S.C.; Fanaee-T, H.; Khalil, A. An Eigenspace approach for detecting multiple space-time disease clusters: Application to measles hotspots detection in Khyber-Pakhtunkhwa, Pakistan. PLoS ONE 2018, 13, e0199176. [Google Scholar] [CrossRef]
  19. Kulldorff, M.; Hjalmars, U. The Knox method and other tests for space-time interaction. Biometrics 1999, 55, 544–552. [Google Scholar] [CrossRef]
  20. Mantel, N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27, 209–220. [Google Scholar] [PubMed]
  21. Roza, D.L.; Caccia-Bava, M.C.G.; Martinez, E.Z. Spatio-temporal patterns of tuberculosis incidence in Ribeirao Preto, state of Sao Paulo, southeast Brazil, and their relationship with social vulnerability: A Bayesian analysis. Rev. Soc. Bras. Med. Trop. 2012, 45, 607–615. [Google Scholar] [CrossRef] [PubMed]
  22. Liu, Y.; Li, X.; Wang, W.; Li, Z.; Hou, M.; He, Y.; Wu, W.; Wang, H.; Liang, H.; Guo, X. Investigation of space-time clusters and geospatial hot spots for the occurrence of tuberculosis in Beijing. Int. J. Tuberc. Lung Dis. 2012, 16, 486–491. [Google Scholar] [CrossRef] [PubMed]
  23. Kulldorff, M.; Athas, W.F.; Feuer, E.J.; Miller, B.A.; Key, C.R. Evaluating cluster alarms: A space-time scan statistic and brain cancer in Los Alamos, New Mexico. Am. J. Public Health 1998, 88, 1377–1380. [Google Scholar] [CrossRef]
  24. Kulldorff, M. Prospective time periodic geographical disease surveillance using a scan statistic. J. R. Stat. Soc. Ser. A-Stat. Soc. 2001, 164, 61–72. [Google Scholar] [CrossRef]
  25. Kulldorff, M.; Heffernan, R.; Hartman, J.; Assunção, R.; Mostashari, F. A space-time permutation scan statistic for disease outbreak detection. PLoS Med. 2005, 2, 0216–0224. [Google Scholar] [CrossRef]
  26. Fanaee-T, H.; Gama, J. An eigenvector-based hotspot detection. arXiv 2014, arXiv:1406.3191. Available online: https://arxiv.org/abs/1406.3191 (accessed on 7 September 2019).
  27. Xie, Z.; Yan, J. Kernel density estimation of traffic accidents in a network space. Comput. Environ. Urban Syst. 2008, 32, 396–406. [Google Scholar] [CrossRef]
  28. Birant, D.; Kut, A. ST-DBSCAN: An algorithm for clustering spatial–temporal data. Data Knowl. Eng. 2007, 60, 208–221. [Google Scholar] [CrossRef]
  29. Levine, N. Crime mapping and the Crimestat program. Geogr. Anal. 2006, 38, 41–56. [Google Scholar] [CrossRef]
  30. Wu, S.; Wu, F.; Hong, R.; He, J. Incidence analyses and space-time cluster detection of hepatitis C in Fujian province of China from 2006 to 2010. PLoS ONE 2012, 7, e40872. [Google Scholar] [CrossRef] [PubMed]
  31. Zhu, B.; Liu, J.; Fu, Y.; Zhang, B.; Mao, Y. Spatio-temporal epidemiology of viral hepatitis in China (2003–2015): Implications for prevention and control policies. Int. J. Environ. Res. Public Health 2018, 15, 661. [Google Scholar] [CrossRef] [PubMed]
  32. Tadesse, S.; Enqueselassie, F.; Hagos, S. Spatial and space-time clustering of tuberculosis in Gurage Zone, Southern Ethiopia. PLoS ONE 2018, 13, e0198353. [Google Scholar] [CrossRef] [PubMed]
  33. Qi, L.; Tang, W.; Zhao, H.; Ling, H.; Su, K.; Zhao, H.; Li, Q.; Shen, T. Epidemiological characteristics and spatial-temporal distribution of hand, foot, and mouth disease in Chongqing, China, 2009–2016. Int. J. Environ. Res. Public Health 2018, 15, 270. [Google Scholar] [CrossRef] [PubMed]
  34. Desjardins, M.R.; Whiteman, A.; Casas, I.; Delmelle, E. Space-time clusters and co-occurrence of chikungunya and dengue fever in Colombia from 2015 to 2016. Acta Trop. 2018, 185, 77–85. [Google Scholar] [CrossRef] [PubMed]
  35. Low, G.K.; Papapreponis, P.; Isa, R.M.; Gan, S.C.; Chee, H.Y.; Te, K.K.; Hatta, N.M. Geographical distribution and spatio-temporal patterns of hospitalization due to dengue infection at a leading specialist hospital in Malaysia. Geospat. Health 2018, 13, 642. [Google Scholar] [CrossRef] [PubMed]
  36. Umer, M.F.; Zofeen, S.; Majeed, A.; Hu, W.; Qi, X.; Zhuang, G. Spatiotemporal Clustering Analysis of Malaria Infection in Pakistan. Int. J. Environ. Res. Public Health 2018, 15, 1202. [Google Scholar] [CrossRef]
  37. Zhao, F.; Cheng, S.; He, G.; Huang, F.; Zhang, H.; Xu, B.; Murimwa, T.C.; Cheng, J.; Hu, N.; Wang, L. Space-time clustering characteristics of tuberculosis in China, 2005–2011. PLoS ONE 2013, 8, e83605. [Google Scholar] [CrossRef]
  38. Kulldorff, M. A spatial scan statistic. Commun. Stat.-Theory Methods 1997, 26, 1481–1496. [Google Scholar] [CrossRef]
  39. District heath Information System. Quarter Reports. Available online: http://www.dhiskp.gov.pk/reports.php (accessed on 7 September 2019).
  40. District Health Information Khyber Pakhtunkhwa, Pakistan. Available online: http://www.dhiskp.gov.pk (accessed on 1 September 2019).
  41. SaTScan v9.6 March 2018. Available online: https://www.satscan.org (accessed on 1 September 2019).
  42. Kulldorff, M. SaTScan User Guid. 2018. Available online: https://www.satscan.org/techdoc.html (accessed on 1 September 2019).
  43. Kulldorff, M.; Nagarwalla, N. Spatial disease clusters: Detection and inference. Stat. Med. 1995, 14, 799–810. [Google Scholar] [CrossRef]
  44. BOS Khyber Pakhtunwa. Development Statistics of Khyber Pakhtunkhwa. 2019. Available online: https://kpbos.gov.pk/allpublication/1 (accessed on 1 September 2019).
  45. Malik, M.S.; Afzal, M.; Farid, A.; Khan, F.U.; Mirza, B.; Waheed, M.T. Disease Status of Afghan Refugees and Migrants in Pakistan. Front. Public Health 2019, 7, 185. [Google Scholar] [CrossRef] [PubMed]
  46. Smeda.org. Districts Profiles Khyber Pakhtunkhwa. Available online: https://smeda.org/index.php?option=com_phocadownload&view=category&id=1:district-profiles (accessed on 13 January 2020).
  47. Khan, J.A.; Irfan, M.; Zaki, A.; Beg, M.; Hussain, S.F.; Rizvi, N. Knowledge, attitude and misconceptions regarding tuberculosis in Pakistani patients. J. Pak. Med. Assoc. 2006, 56, 211. [Google Scholar] [PubMed]
Figure 1. Map of Khyber Pakhtunkhwa province.
Figure 1. Map of Khyber Pakhtunkhwa province.
Ijerph 17 01413 g001
Figure 2. The locations of space-time clusters in the study area.
Figure 2. The locations of space-time clusters in the study area.
Ijerph 17 01413 g002
Figure 3. The locations of spatial clusters in the study area.
Figure 3. The locations of spatial clusters in the study area.
Ijerph 17 01413 g003
Table 1. The detected space-time clusters of TB cases in Khyber Pakhtunkhwa, Pakistan.
Table 1. The detected space-time clusters of TB cases in Khyber Pakhtunkhwa, Pakistan.
Cluster TypeDistrictsNumber of DistrictsTime-FrameObserved CountsExpected CountsLLRP-value
Most LikelyDir Upper11st quarter 201710686219.3731541.8<0.001
1st secondaryBuner, Swat21st quarter 2015 to 4th quarter 2015130192599.4611044.8<0.001
2nd secondaryBannu11st quarter 2015 to 4th quarter 201644952001.461170.532<0.001
3rd secondaryHangu12nd quarter 2015 to 4th quarter 20162216822.69810.864<0.001
4th secondaryNoshera, Mardan22nd quarter2015 to 2nd quarter 201663564203.86496.057<0.001
5th secondaryCharsadda14th quarter 2015505365.4523.865<0.001
Table 2. The detected spatial clusters of TB cases in Khyber Pakhtunkhwa, Pakistan.
Table 2. The detected spatial clusters of TB cases in Khyber Pakhtunkhwa, Pakistan.
YearCluster TypeDistrictsNumber of DistrictsObserved CountsExpected CountsLLRP-value
2015Most likelySwat1124062500.8111728.029<0.001
1st secondaryBannu124351292.4220.132<0.001
2nd secondaryHangu11022616.25113.73<0.001
3rd secondaryShangla11102852.7834.255<0.001
4th secondaryHaripur113811181.3716.602<0.001
2016Most likelyShangla, Swat247942363.251097.165<0.001
1st secondaryKarak, Hangu, Bannu, Kohat448142770.32714.454<0.001
2nd secondaryNoshera118911221.14166.741<0.001
3rd secondaryMansehra118721350.3595.685<0.001
4th secondaryTank1568332.6969.656<0.001
5th secondaryMardan123712039.1628.081<0.001
6th secondaryHaripur11133966.7614.12<0.001
7th secondaryDir Lower110401002.260.7310.997
2017Most likelyUpper Dir111025957.1118802.163<0.001
1st secondaryBannu120301181.08262.797<0.001
2nd secondarySwat130532335.65109.45<0.001
3rd secondaryTank1573396.3135.082<0.001
4th secondaryShangla1979766.3727.848<0.001
5th secondaryHangu1599524.665.1250.065
6th secondaryHaripur110691014.361.4960.859
2018Most likelyShangla, Swat, Battagram379972672.794206.852<0.001
1st secondaryBannu, Lakki Marwat, Karak, Hangu, Tank, Kohat757963510.64759.992<0.001
2nd secondaryChitral1605337.3987.288<0.001
3rd secondaryHaripur1946756.4622.787<0.001
4th secondaryMansehra112661173.853.720.183
5th secondaryDir Lower111531082.932.3312630.570

Share and Cite

MDPI and ACS Style

Ullah, S.; Daud, H.; Dass, S.C.; Fanaee-T, H.; Kausarian, H.; Alamgir. Space-Time Clustering Characteristics of Tuberculosis in Khyber Pakhtunkhwa Province, Pakistan, 2015–2019. Int. J. Environ. Res. Public Health 2020, 17, 1413. https://doi.org/10.3390/ijerph17041413

AMA Style

Ullah S, Daud H, Dass SC, Fanaee-T H, Kausarian H, Alamgir. Space-Time Clustering Characteristics of Tuberculosis in Khyber Pakhtunkhwa Province, Pakistan, 2015–2019. International Journal of Environmental Research and Public Health. 2020; 17(4):1413. https://doi.org/10.3390/ijerph17041413

Chicago/Turabian Style

Ullah, Sami, Hanita Daud, Sarat C. Dass, Hadi Fanaee-T, Husnul Kausarian, and Alamgir. 2020. "Space-Time Clustering Characteristics of Tuberculosis in Khyber Pakhtunkhwa Province, Pakistan, 2015–2019" International Journal of Environmental Research and Public Health 17, no. 4: 1413. https://doi.org/10.3390/ijerph17041413

APA Style

Ullah, S., Daud, H., Dass, S. C., Fanaee-T, H., Kausarian, H., & Alamgir. (2020). Space-Time Clustering Characteristics of Tuberculosis in Khyber Pakhtunkhwa Province, Pakistan, 2015–2019. International Journal of Environmental Research and Public Health, 17(4), 1413. https://doi.org/10.3390/ijerph17041413

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