Next Article in Journal
Respiratory Condition of Family Farmers Exposed to Pesticides in the State of Rio de Janeiro, Brazil
Previous Article in Journal
Resident Perceptions of Neighborhood Conditions, Food Access, Transportation Usage, and Obesity in a Rapidly Changing Central City
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 15
Issue 6
10.3390/ijerph15061202
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatiotemporal Clustering Analysis of Malaria Infection in Pakistan

by
Muhammad Farooq Umer
1,
Shumaila Zofeen
1,
Abdul Majeed
2,
Wenbiao Hu
3,
Xin Qi
1,4,* and
Guihua Zhuang
1,*
1
School of Public Health, Xi’an Jiaotong University Health Science Center, Xi’an 710061, China
2
Directorate of Malaria Control Program, Islamabad 44000, Pakistan
3
School of Public Health and Social Work, Queensland University of Technology, Kelvin Grove, QLD 4059, Australia
4
Global Health Institute, Xi’an Jiaotong University Health Science Center, Xi’an 710061, China
*
Authors to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2018, 15(6), 1202; https://doi.org/10.3390/ijerph15061202
Submission received: 17 April 2018 / Revised: 24 May 2018 / Accepted: 4 June 2018 / Published: 7 June 2018
Browse Figures
Versions Notes

Abstract

:
Despite tremendous progress, malaria remains a serious public health problem in Pakistan. Very few studies have been done on spatiotemporal evaluation of malaria infection in Pakistan. The study aimed to detect the spatiotemporal pattern of malaria infection at the district level in Pakistan, and to identify the clusters of high-risk disease areas in the country. Annual data on malaria for two dominant species (Plasmodium falciparum, Plasmodium vivax) and mixed infections from 2011 to 2016 were obtained from the Directorate of Malaria Control Program, Pakistan. Population data were collected from the Pakistan Bureau of Statistics. A geographical information system was used to display the spatial distribution of malaria at the district level throughout Pakistan. Purely spatiotemporal clustering analysis was performed to identify the high-risk areas of malaria infection in Pakistan. A total of 1,593,409 positive cases were included in this study over a period of 6 years (2011–2016). The maximum number of P. vivax cases (474,478) were reported in Khyber Pakhtunkhwa (KPK). The highest burden of P. falciparum (145,445) was in Balochistan, while the highest counts of mixed Plasmodium cases were reported in Sindh (22,421) and Balochistan (22,229), respectively. In Balochistan, incidence of all three types of malaria was very high. Cluster analysis showed that primary clusters of P. vivax malaria were in the same districts in 2014, 2015 and 2016 (total 24 districts, 12 in Federally Administered Tribal Areas (FATA), 9 in KPK, 2 in Punjab and 1 in Balochistan); those of P. falciparum malaria were unchanged in 2012 and 2013 (total 18 districts, all in Balochistan), and mixed infections remained the same in 2014 and 2015 (total 7 districts, 6 in Balochistan and 1 in FATA). This study indicated that the transmission cycles of malaria infection vary in different spatiotemporal settings in Pakistan. Efforts in controlling P. vivax malaria in particular need to be enhanced in high-risk areas. Based on these findings, further research is needed to investigate the impact of risk factors on transmission of malaria in Pakistan.

1. Introduction

Malaria is a deadly vector-borne infectious disease, rendering around 42% of the world’s population at risk [1]. In the year 2016, malaria resulted in 216 million cases and 445,000 deaths worldwide [2]. Malaria in humans is predominantly due to Plasmodium falciparum and Plasmodium vivax, which are transmitted mostly by the same species of the genus Anopheles [2,3]. Around 99% of malarial deaths were attributed to P. falciparum, with the majority of them occurring in Africa. P. vivax, on the other hand, is the dominant cause of malaria in many countries of the Eastern Mediterranean Regional Office (EMRO, Cairo, Egypt) and in the American region [2,4], with 2.5 billion people globally at endemic risk of its morbidity and mortality [5]. Experts believe that in countries having disease load of both P. falciparum and P. vivax malaria, the real challenge would be the elimination of P. vivax, as it has the ability to remain dormant, evading blood-stage treatment and thus causes relapsing infections because of its prolonged neglect based on mistaken assumptions of its harmlessness [6,7]. Mixed Plasmodium infection is another significant type of malaria which is usually under-reported, because of over-reporting of the more expressed P. falciparum or P. vivax in conventional diagnostic methods [8,9].
Pakistan is classified epidemiologically as a moderate malaria-endemic country with an average national annual parasitic incidence of 1.69%, having diversity within and between its provinces and districts [10]. Pakistan remains one of the highest malaria burden sharing countries in the EMRO region of the World Health Organization (WHO, Geneva, Switzerland), which also includes Sudan, Yemen, Somalia and Afghanistan [2,11]. It is pertinent to note that both P. falciparum and P. vivax malaria are prevalent in Pakistan, and incidence of P. vivax is over 4 times that of P. falciparum [12]. In 2015, 81% of the total deaths due to P. vivax outside sub-Saharan Africa took place in only four countries: Ethiopia, India, Indonesia and Pakistan [13]. The coexistence of two species of malaria (P. falciparum and P. vivax) in the same community increases the likelihood of mixed Plasmodium infections in a single individual [14] and Pakistan is no exception to this, as mixed plasmodium infections have being regularly reported in many parts of the country, making situation more serious [15,16,17].
Malaria has significant impact on developing countries like Pakistan, despite all the disease control strategies designed and executed by the concerned authorities [18,19]. One of the key reasons for this plight is that, epidemiologically, there are not sufficient data and/or robust studies available in Pakistan to evaluate the incidence of malaria [20]. This implies that, to make our efforts more fruitful in controlling malaria, effective surveillance of the disease, as well as information on incidence and/or prevalence of infections or incidence of clinical cases and deaths, are crucial for identifying the most affected population groups to direct resources to these areas for malaria control and prevention [21,22].
Malaria incidence, like many infectious diseases, exhibits spatial and temporal variations with respect to climatic and environmental conditions [1,23]. Thus, the spatial and temporal distribution of malaria—the most basic aspect of its epidemiology—needs to be studied more in the areas of its occurrence [21,24]. Malaria transmission nationwide in Pakistan is seasonal for both P. vivax (peaks usually in April) and P. falciparum (peaks usually in October) [25]. This characteristic of malaria transmission and distribution may be explored by geographical information system (GIS) mapping and spatial analysis, which are helpful in planning, implementing and measuring impact of disease interventions and then monitoring at subnational, national and regional levels for eventual control of malaria [26,27].
Spatial and spatiotemporal studies on malaria, including its incidence, prevalence, mortality, risk measurement and many other dimensions, have been helpful in many ways in successful design and implementation of programs on malaria control and prevention [24,27]. However, very few studies on the spatial pattern of malaria have been done in Pakistan before. The current study was designed to find out the spatiotemporal pattern of the annual incidence of malaria in Pakistan at the district level, and to detect clusters of high-risk areas of the disease in the country. The results in this study may provide robust evidence for local malaria control and prevention strategy in Pakistan.

2. Methods

2.1. Study Area

This study aimed to explore malaria cases and incidence at the district level in Pakistan. Geographically, Pakistan is located within the coordinates of 30°00′ N, 70°00′ E, with total land area approximately 796,096 km2. According to Pakistan Bureau of Statistics (Census Report 2017), the population of Pakistan is 207.8 million (excluding Azad Kashmir and Gilgit Baltistan) with an average annual growth rate of 2.4% since its last census completed in 1998 [28].
Pakistan is comprised of the federal capital territory (Islamabad, Pakistan), five provinces (Punjab, Sindh, Khyber Pakhtunkhwa (KPK, Pakistan), Balochistan and Gilgit Baltistan), a group of Federally Administered Tribal Areas (FATA, Peshawar, Pakistan) and an autonomous state (Azad Kashmir, Pakistan). Each of these administrative units is further divided into districts and there is a total of 146 districts in Pakistan. Figure 1 shows the administrative units in Pakistan. Pakistan has diverse climatic, geographic and topographic features, ranging from high snow-covered mountains to hot and dry deserts and moderate plains and coastline [29]. The Human Development Index (HDI) of Pakistan is 0.515, ranking 146 out of 187 countries [10].

2.2. Data Collection

The Directorate of Malaria Control Program (DoMC), Islamabad, Pakistan, is the sole federal organization for the national-level planning and policy making to control the disease in the country, and it is implemented by the respective provinces. The aggregated, district-wise annual malaria data by species from 2011 to 2016 were obtained from the DoMC, through a formal request. More details of incidence time (e.g., day and month) and home address (e.g., street number) of malaria cases were not available. All the malaria cases diagnosed at each public sector healthcare facility from individual districts in Pakistan are supposed to report digitally to the Directorate of Malaria Control Program, where data are compiled into a summary sheet. This summary was done at the district level and was used in the current study, which included P. vivax, P. falciparum, mixed Plasmodium and total positive malaria cases, respectively. The national policy with respect to malaria diagnosis and treatment during the study period has been consistent, with improved accessibility of laboratory testing in general. The diagnosis on malaria is made on the basis of WHO-recommended methods of microscopy or rapid diagnostic tests (RDT) [2]. Ten out of 146 districts in Pakistan (mainly Islamabad and districts from Gilgit Baltistan) are not entitled to report cases to DoMC. Thus, the current study represents data from 136 of the total 146 districts in the country. As Gilgit Baltistan has recently been designated as a separate province and the malaria control program was not operational in this part of the country, no data collection or reporting was done there and Gilgit Baltistan was not included in our study. Islamabad, although the capital of the country, is not included in the malaria control program and that is why no data were available on its malaria infection in our study. There are three independent sources which report to DoMC via District Health Information System (DHIS), Malaria Information System (MIS) and Disease Early Warning System (DEWS), but they do not work with the same indicators of malaria. Malaria cases are reported mainly by DHIS, whose reports confirmed cases based on microscopy or RDT. MIS, which was developed recently for the electronic recording of disease statistics, can be used for planning and policy making and contains a few more indicators related to program management and administration than DHIS. DEWS reports the data on suspected malaria cases (not included in our study) and not the confirmed cases diagnosed by microscopy or RDT. The confirmed cases of malaria are mainly reported to DoMC by DHIS, which has its limitations. Firstly, a few of the public sector healthcare facilities do not report to DHIS (security-compromised areas), and secondly, it does not include private healthcare facilities (80% of the patients utilize private healthcare services, which are not regulated and do not have data records [30]), implying that the malaria burden in the country is expected to be 5 to 6-fold more than reported [31,32]. Despite all its limitations, data from DoMC is still the best available, as it is systematically collected by reliable sources and covers almost all parts of the country.
District-wise population data were downloaded from Pakistan Bureau of Statistics, which is based on the latest census completed in 2017. This district-wise population of 2017 was taken as a reference and it was retrospectively extrapolated, depreciating by 2.4% annually, which is the average annual growth rate of the country between 1998 and 2017.

2.3. Data Analyses

Statistical analyses were performed by steps. Firstly, malarial cases data and population data were compiled in Microsoft Excel according to the requirements of the GIS software. Then, the annual incidence rate (IR) and relative risk of each species of malaria cases (P. vivax, P. falciparum, mixed Plasmodium and total positive) at the district level were calculated using the mean malaria incidence rate of the whole country by species in each year as the reference. Using the relative risk at the district level can help quantify district IR compared to the national level. The descriptive statistics of IR and sequence plots of annual cases for each species of malaria were performed using statistical package SPSS 13.0 [33].
The aligned annual malaria cases data by species, along with annual population and latitude/longitude coordinates of centroids of districts, respectively, were input in SaTScan 9.4.4 for the spatial cluster analysis. SaTScan is used to perform spatial, temporal and space–time analysis using scan statistics. This software has widespread application in fields like epidemiology, ecology, economics, geography and many others [34]. In the field of health and epidemiology, SaTScan has been used to detect spatial and space–time disease clusters and to evaluate the statistical significance of these clusters for geographical surveillance of the disease [35,36]. In our study, purely spatial scanning analysis was performed for clusters with high rates using the discrete Poisson model. The clusters included primary clusters (the cluster with highest risk among all clusters in the map) and secondary clusters (other statistically significant high-risk clusters) in each year from 2011 to 2016, identified as a result of SaTScan analysis. The selection of “No geographic overlap” was used as a default setting so that secondary clusters may not overlap the most significant (primary) cluster. The maximum spatial cluster size on the spatial window was set to be at 10% of the population at risk, and the maximum radii of the reported clusters was set to be 200 km from the SaTScan setting, based on the total land area and population distribution in Pakistan. A similar method of setting maximum population size and cluster radii has also been used in another spatial cluster study [36]. Relative risk of high-risk clusters of each species of malaria for each year was estimated using the Poisson model and by SaTScan, considering case number and population size of districts within each cluster. The relative risk was then plotted in an Excel sheet containing the aligned data.
The spatial patterns of malaria IR, the relative risk and the information on primary and secondary clusters by species at the district level (2011–2016) were displayed in maps using ArcGIS 10.2 [37]. This study obtained ethical approval from all the relevant authorities. The approval letter for malaria by species data was obtained from Directorate of Malaria Control Program Pakistan, then an Exemption approval letter was obtained from National Bioethical Committee Pakistan, letter number No.4-87/NBC-279-Exempt./17/1139. The research protocol was approved by Institutional Review Board of Xi’an Jiaotong University, China.

3. Results

A large number of malaria cases by species were reported across the country during the period of this study (2011–2016). We found 1,593,409 total positive cases in the 6-year study period in Pakistan, with highest burden in KPK. A total of 1,221,026 cases of P. vivax were reported, with the majority of cases being found in KPK (474,478). There were 306,172 P. falciparum cases, with the highest number of cases in Balochistan (145,445). The mixed Plasmodium cases totaled 66,211, with Sindh (22,421) and Balochistan (22,229) being the top contributors (both provinces reported almost equal numbers of mixed Plasmodium cases).
Table 1 explains the IR (per 100,000) of malaria by species at the district level during the study period. The incidence of malaria by species varies considerably among the districts and over the years. Standard deviation is higher than the respective mean of IR for each species of malaria in each year of the study period, showing disparities in disease incidence among the districts in Pakistan, and this difference is further evidence of the wide gaps amid all the minimum and maximum values of IR in this table. IR of P. vivax malaria was higher than the other two species (P. falciparum and mixed Plasmodium malaria) at the district level in each year of the study period, except in 2011, when the IR of P. falciparum remained the highest of all three malaria species at the district level.
Figure 2 shows the IR of aggregated total positive cases from 2011 to 2016 and represents a summary of the disease burden in the country at the district level during the whole study period (Figure 2A). We divided IR into five categories to demarcate the districts with no data and the ones with high to very high and low to very low IR of the disease. It also shows the sequence plot of annual cases of malaria by species, 2011–2016. This figure upholds the findings in Table 1, as it confirms the annual variations in burden of malaria in the country. Figure 2B exhibits that the case counts for total positive and P. vivax are in close proximity to one another and look alike, implying that P. vivax is the principal burden bearer of malaria in the country and is responsible for annual variation in the disease load. The figure also demonstrates that P. falciparum infection is declining with the passage of time, while mixed Plasmodium has a potential rising trend.
Figure 3 illustrates the IR (per 100,000) of malaria by species at the district level, 2011–2016. In both Figure 2 and Figure 3, the darker shades in the geographical maps indicate districts with higher IRs of malaria, while the lighter shades represent the ones with lower IRs. Figure 3 expounds on Figure 2 in detail, thus conveying that IR of P. vivax malaria in Pakistan is much higher than that of P. falciparum and mixed Plasmodium malaria, and the patterns of distribution of total positive and P. vivax malaria in the corresponding years were similar. This figure further illustrates that total positive and P. vivax malaria with high IRs are concentrated in the south, west and northwest parts of Pakistan and that P. falciparum malaria is shrinking over the years, while mixed Plasmodium is gradually expanding.
Figure 4 displays the relative risk of malaria by species at the district level from 2011 to 2016. The distribution patterns of relative risk (RR) of total positive and P. vivax malaria in the corresponding years look alike and are largely in medium categories (1–19.99), covering most districts of the country. RR values of P. falciparum malaria are mostly in medium to high categories (1–19.99 and 20–39.99) and RR values of mixed Plasmodium malaria attained peak values (RR 40 and above) in a few districts.
Figure 5 exhibits the high-risk clusters of malaria by species at the district level from 2011 to 2016. This figure displays most significant (primary) and secondary clusters of the disease and the “no cluster areas” in the geographical maps. Once again, there is a similarity in the patterns of distribution of total positive and P. vivax disease clusters.
Table 2 shows details about most significant high-risk districts (primary clusters) of the disease during the study period. It gives yearly information on number of districts, coordinates of centroid, radii, population, number of cases, IR, relative risk and statistical significance of the primary clusters of malaria by species. It was established that the number of cases and IRs of P. vivax malaria were much higher than those of the other two species of malaria (P. falciparum and mixed Plasmodium) in the respective years. Another important finding is that the districts (number and locations) in primary clusters of (i) P. vivax are same for 2014, 2015 and 2016 (total 24 districts, 12 in FATA, 9 in KPK, 2 in Punjab and 1 in Balochistan); (ii) P. falciparum are same for 2012 and 2013 (total 18 districts, all in Balochistan); (iii) mixed Plasmodium are same for 2014 and 2015 (total 7 districts, 6 in Balochistan and 1 in FATA). We also tested spatial cluster analysis using other parameters (100 km, 300 km and 400 km for maximum radius of cluster, and 5%, 15% for maximum population at risk) and did not find significant differences from the above results.

4. Discussion

In this study, the incidence of malaria by species in Pakistan fluctuated among districts, which is in accordance with previous findings [38]. Our study confirmed the occurrence of three species of malaria (P. vivax, P. falciparum and mixed Plasmodium) in the country, and incidence of P. vivax malaria was significantly higher than for the other two species [39,40]. The increase in the incidence of P. vivax in recent years may be due to the focused efforts of the malaria control program on case detection. Most of the malaria high endemic districts were found in Balochistan, FATA and KPK, with Sindh estimated as a moderately high endemic area. These findings are similar to those in past studies conducted in Pakistan [11,41]. According to our study, spatial distribution of the most significant high-risk clusters of malaria by species exhibited that P. vivax primary clusters comprised districts in FATA and KPK, P. falciparum primary clusters were mainly in Balochistan, while mixed Plasmodium primary clusters included districts in Sindh and Balochistan. In general, the pattern of malaria incidence remained consistent over the study period. The finding in our study that the same districts persisted in primary clusters in different years can be useful for prompting policy makers to involve authorities from different administrative units to put combined efforts in combating the disease.
The high incidence of malaria by species in FATA, Balochistan and, to some extent, in KPK, is explained by the fact that most of the districts in these areas are lagging behind in their social and living standards compared to the rest of country [42], which is due to compromised the law and order situation in these areas, resulting in very low-performing health facilities (e.g., limited health care resources and unavailability to health facilities for diagnosis and treatment in some communities), therefore increasing malaria in particular, and other diseases in general, in these areas [31,43]. In Balochistan, alarmingly, incidence of all three species of malaria was very high, demanding contemplation by the concerned authorities. In some districts of Balochistan, P. vivax malaria was dominant, as supported by the evidence from other studies [44,45], whereas, in its eastern and southeastern districts, P. falciparum malaria burden was higher [46,47]. In the northern parts adjoining KPK, both P. vivax and P. falciparum malaria coexisted in large numbers [48]. Moreover, in this study, Balochistan, along with Sindh, reported the highest burden of mixed Plasmodium malaria, which reiterates the presence in abundance of both P. vivax and P. falciparum malaria in these areas. In most districts of KPK, P. vivax turned out to be higher than P. falciparum and mixed Plasmodium malaria, as was shown in the other studies [15,49,50]. Sindh was estimated as a moderately high malaria-endemic area in this study, with some high-risk districts adjoining Balochistan and a few near coastal areas [51].
The hypothesis that P. falciparum and P. vivax are endemic in separate regions can be attributed to the variability of climatic conditions, proximity to or remoteness from the coastal line and particularly due to antimalarial drug resistance [40,52]. Incidence of P. falciparum malaria is higher in Balochistan and Sindh, where chloroquine drug resistance may have contributed to this species’ dominance and, also, these two provinces have a coastal line in proximity to many of its districts [40,52]. The use of the latest antimalarial drugs, like artemisinin-based combination therapy, could help address the problem [40,52]. On the other hand, P. vivax is endemic throughout the country, including Sindh and Balochistan, and chloroquine is still considered the first-line drug therapy against P. vivax [40,52].
P. vivax and mixed Plasmodium malaria incidences were reasonably high in many districts of Sindh, which is a serious concern, whereas P. falciparum malaria showed a decline over the years in this province in the course of this study. Punjab was assessed as a low malaria-endemic area in this study, despite the fact that it covers 53% of the total population in Pakistan. The lower incidence of malaria by species in Punjab could be attributed to the fact that, firstly, during the last three decades, ecological changes (waterlogging and salinization) resulted in selection of the fittest yet less efficient vector species in this province [53] and, secondly, because of better healthcare delivery infrastructure than the rest of country due to better social and living standards here [52]. It is important to mention that, although overall incidence rate of malaria by species was significantly low in Punjab, southern districts had higher IRs than the central and northern ones within the province, thus requiring attention by the authorities [19].
In Pakistan, an estimated 20% of the population consult public sector healthcare facilities for their health needs. These health-seeking statistics are not uniform throughout the country but have been found varying in different provinces. According to Pakistan Social and Living Standards Measurement Survey 2014–2015, 17% of the people visited public sector health facilities for consultation and/or treatment in Punjab, 20% in Sindh, and 28% each in Khyber Pakhtunkhwa and Balochistan. These percentages suggest that the better the living standards in a province, the more their paying capacity and, thus, the more its people seek healthcare from private sector health facilities. The living standards and therefore paying capacity of Balochistan and KPK are relatively lesser as compared to other provinces, hence malaria burden in these areas is higher. Moreover, this also reflects the performance of public sector health facilities in KPK and Balochistan that despite these facilities are being used by the community the malaria burden is increasing. As no adjustments were made in this study for the malaria burden reported by public sector health facilities versus the private sector ones, so considering the higher percentages of population utilizing private sector healthcare facilities, it is important to mention here that estimates shown here are to be shown patterns of incidence rather than absolute values.
This study has several strengths to its credit. This is the first study in Pakistan exploring spatiotemporal pattern and high-risk clusters of malaria incidence at the district level. Latest data (malaria by species and population data) were used in this study to describe current patterns of disease. The study encompassed almost whole country at the district level and discussed all three species of malaria infections incidence in Pakistan. This study provides novel information on the malaria burden in Pakistan with respect to time and space discovering malaria hotspots at the district level. The results in this study may provide scientific evidence for malaria control and prevention strategies in Pakistan.
Our study had certain limitations as well that can be taken into account in designing the subsequent studies. This study caters only those patients who consulted public sector healthcare facilities (20% of the total population) and gives no information about the larger section of population going to private sector healthcare facilities. Thus, the actual malarial incidence may be underestimated. The data did not describe individual patients’ profiles (age, sex, socioeconomic and education status, etc.), as aggregated data was used in this study. The data we collected from DoMC was aggregated annual malaria data, and month and day of malaria incidence at the district level were not available. Thus, seasonality of malaria at the district level in Pakistan could not be established. Geographical boundaries of districts were consistent in the study period, but the aggregated malaria data at the district level have no more specific geographical details, e.g., home address information. Thus, the results at the district level (including using centroid of districts to detect spatial clusters by SaTScan) may have bias from a geographical perspective, e.g., masking the potential impact of socioeconomic disparity across different communities on malaria incidence within each district. Moreover, the exact neighborhood of each district, which is significant in influencing distribution of clustering districts, cannot be reflected in SaTScan analysis.
Some recommendations can be proposed based on the findings of this study. P. vivax malaria may no longer be considered as benign, keeping in view the morbidity and mortality associated with it [5]. Also, the number of hospitalizations with P. vivax infection have increased over time, and mortality rates attributable to both species of malaria (P. vivax and P. falciparum) in our part of the world are comparable [54,55]. As the diagnosis of mixed Plasmodium malaria is becoming more challenging now, owing to the conventional and sometimes poor microscopy practices, enhancing microscopy skills through trainings and, if feasible, bringing new diagnostic techniques, like polymerase chain reaction (PCR), into practice may improve the diagnosis [17,39]; although costly, PCR helps in precise detection of mixed infections. Irrational and high use of antimalarial drugs, as reported in different studies [40,52], has contributed to drug-resistant malaria in Pakistan, which should be noticed by authorities before the threat gets too serious to tackle. The knowledge and practices about malaria and its prevention are seriously deficient among the people of Pakistan, especially those living in rural settings, and this aspect requires work on a large scale for prospects of success in controlling malaria in Pakistan [56,57]. The high burden of malaria needs the attention of authorities to locate the high-risk areas and then direct resources and interventions (long-lasting insecticidal nets, sprays, drugs against malaria, community awareness campaigns) towards these areas by involving local community and local authorities for long-lasting control and, eventually, elimination of the disease. The control strategies are already in place in Pakistan, including in the areas of its high incidence, but due to limited supplies and sometimes inequitable distribution, disease control and prevention still needs to be improved. Public sector healthcare facilities need to deliver better services to the community, as they are utilized mostly by the underprivileged ones. The private sector healthcare facilities should be regularized by legislative bodies to streamline their performances by maintaining and sharing data with the malaria control program, thus controlling the disease in totality, including the disease-reporting system. More public health investigations and studies need to be conducted in hotspot areas explored in our study in Pakistan. Besides, seasonality of malaria and its environmental drivers (e.g., climate change and extreme weather) are important aspects in the pattern of disease development, yet very few detailed studies have been reported in Pakistan [58]. Thus, this should be explored for better understanding of the disease progression in the country.

5. Conclusions

This study described spatial and temporal variations in all three species of malaria widespread in Pakistan at the district level. As a result, high-risk districts of malaria by species were identified by conducting spatial cluster analysis, which has its significance in exploring area-specific evidence on disease progression. The study results suggested ramping up malaria control activities in Balochistan, FATA, KPK and Sindh, which are high burden disease areas. Diagnosis of mixed Plasmodium malaria is a challenge to be taken seriously, keeping in mind its high-risk areas (i.e., Balochistan and Sindh). Efforts in controlling P. vivax malaria need to be enhanced, since it is the major burden bearer in most parts of the country, while, simultaneously, not being complacent about P. falciparum malaria, which is still causing damages in the country, especially in Balochistan. More studies are required with detailed spatial and temporal analysis in the individual high-risk districts for further exploration of patterns in malaria progression and possible risk factors, for its control and prevention in Pakistan and its regions.

Author Contributions

M.F.U. and X.Q. designed the study. M.F.U. implemented all statistical analysis and drafted the manuscript. G.Z. and X.Q. conceptualized the idea, revised the study protocol and edited the manuscript. S.Z. assisted in data alignment and drafting. A.M. helped in data collection from DoMC. W.H. provided statistical assistance. All the authors read and approved the final manuscript.

Funding

This research was funded by Chinese Scholarship Council (CSC) for M.F.U.’s PhD studies at Xi’an Jiaotong University, China; and X.Q. was funded by China Postdoctoral Science Foundation (Approval No.: 2015M570842) and National Natural Science Foundation of China (Approval No.: 81502892).

Conflicts of Interest

The authors declare that they have no competing interests.

Abbreviations

PPlasmodium
EMROEastern Mediterranean Regional Office
WHOWorld Health Organization
GISGeographical Information System
KPKKhyber Pakhtunkhwa
FATAFederally Administered Tribal Areas
DoMCDirectorate of Malaria Control Program
DHISDistrict Health Information System
MISMalaria Information System
DEWSDisease Early Warning System
IRIncidence Rate
RRRelative Risk
PCRPolymerase Chain Reaction

References

  1. Faburay, B. The case for a “one health” approach to combating vector-borne diseases. Infect. Ecol. Epidemiol. 2015, 5, 28132. [Google Scholar] [CrossRef] [PubMed]
  2. World Health Organization (WHO). World Malaria Report 2017; WHO Press: Geneva, Switzerland, 2017; pp. 33–41. [Google Scholar]
  3. Hundessa, S.H.; Williams, G.; Li, S.; Guo, J.; Chen, L.; Zhang, W.; Guo, Y. Spatial and space–time distribution of Plasmodium vivax and Plasmodium falciparum malaria in China, 2005–2014. Malar. J. 2016, 15, 595. [Google Scholar] [CrossRef] [PubMed]
  4. Gunda, R.; Chimbari, M.J.; Shamu, S.; Sartorius, B.; Mukaratirwa, S. Malaria incidence trends and their association with climatic variables in rural Gwanda, Zimbabwe, 2005–2015. Malar. J. 2017, 16, 1–13. [Google Scholar] [CrossRef] [PubMed]
  5. Gething, P.W.; Elyazar, I.R.F.; Moyes, C.L.; Smith, D.L.; Battle, K.E.; Guerra, C.A.; Patil, A.P.; Tatem, A.J.; Howes, R.E.; Myers, M.F.; et al. A long neglected world malaria map: Plasmodium vivax endemicity in 2010. PLoS Negl. Trop. Dis. 2012, 6, e1814. [Google Scholar] [CrossRef] [PubMed]
  6. Cotter, C.; Sturrock, H.J.W.; Hsiang, M.S.; Liu, J.; Phillips, A.A.; Hwang, J.; Gueye, C.S.; Fullman, N.; Gosling, R.D.; Feachem, R.G.A. The changing epidemiology of malaria elimination: New strategies for new challenges. Lancet 2013, 382, 900–911. [Google Scholar] [CrossRef]
  7. Feachem, R.G.A.; Phillips, A.A.; Hwang, J.; Cotter, C.; Wielgosz, B.; Greenwood, B.M.; Sabot, O.; Rodriguez, M.H.; Abeyasinghe, R.R.; Ghebreyesus, T.A.; et al. Shrinking the malaria map: Progress and prospects. Lancet 2010, 376, 1566–1578. [Google Scholar] [CrossRef]
  8. Imwong, M.; Nakeesathit, S.; Day, N.P.J.; White, N.J. A review of mixed malaria species infections in anopheline mosquitoes. Malar. J. 2011, 10, 253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Ginouves, M.; Veron, V.; Musset, L.; Legrand, E.; Stefani, A.; Prevot, G.; Demar, M.; Djossou, F.; Brousse, P.; Nacher, M.; et al. Frequency and distribution of mixed Plasmodium falciparum-vivax infections in French Guiana between 2000 and 2008. Malar. J. 2015, 14, 10–15. [Google Scholar] [CrossRef] [PubMed]
  10. Directorate of Malaria Control Program Pakistan. National Malaria-Strategic Plan-Pakistan (2015–2020), 2014. Available online: http://www.dmc.gov.pk/documents/pdfs/1National Malaria-Strategic Plan-Pakistan.pdf (accessed on 10 November 2017).
  11. Program Management Unit DoMC. Malaria Annual Report 2016, 2017. Available online: http://dmc.gov.pk/documents/pdfs/Malaria%20Annual%20Report%202016.pdf (accessed on 15 November 2017).
  12. Rahman, S.; Jalil, F.; Khan, H.; Jadoon, M.A.; Ullah, I.; Rehman, M.; Khan, A.M.; Khan, A.; Hayat, A.; Iqbal, Z. Prevalence of malaria in district Shangla, Khyber Pakhtunkhwa, Pakistan. J. Entomol. Zool. Stud. 2017, 5, 678–682. [Google Scholar]
  13. Owens, S. Malaria and the millennium development goals. Arch. Dis. Child. 2015, 100, S53–S56. [Google Scholar] [CrossRef] [PubMed]
  14. Gupta, B.; Gupta, P.; Sharma, A.; Singh, V.; Dash, A.P.; Das, A. High proportion of mixed-species Plasmodium infections in India revealed by PCR diagnostic assay. Trop. Med. Int. Health 2010, 15, 819–824. [Google Scholar] [CrossRef] [PubMed]
  15. Khatoon, L.; Khan, I.U.; Shah, S.A.; Jan, M.I.; Ullah, F.; Malik, S.A. Genetic diversity of Plasmodium vivax and Plasmodium falciparum in Kohat District, Pakistan. Braz. J. Infect. Dis. 2012, 16, 184–187. [Google Scholar] [CrossRef]
  16. Fareed, N.; Ghaffar, A.; Malik, T.S. Spatio-temporal extension and spatial analyses of dengue from Rawalpindi, Islamabad and Swat during 2010–2014. Climate 2016, 4, 23. [Google Scholar] [CrossRef]
  17. Shahzadi, S.; Akhtar, T.; Hanif, A.; Sahar, S.; Niaz, S.; Bilal, H. Molecular detection of malaria in South Punjab with higher proportion of mixed infections. Iran. J. Parasitol. 2014, 9, 37–43. [Google Scholar] [PubMed]
  18. Oneeb, M.; Maqbool, A.; Lateef, M.; Babar, M.E. Detection of Plasmodium falciparum infection in Anopheles stephensi in Punjab, Pakistan. Pak. J. Zool. 2015, 47, 1192–1195. [Google Scholar]
  19. Dash, A.P.; Valecha, N.; Anvikar, A.R.; Kumar, A. Malaria in India: Challenges and opportunities. J. Biosci. 2008, 33, 583–592. [Google Scholar] [CrossRef] [PubMed]
  20. Rana, M.S.; Tanveer, A.; Latif, A.A.; Tahir, A.H. Significant reduction of malaria in the Punjab, Pakistan after introduction of roll back malaria strategy in 2003. Med. Forum Mon. 2011, 22, 54–59. [Google Scholar]
  21. Weiss, D.J.; Mappin, B.; Dalrymple, U.; Bhatt, S.; Cameron, E.; Hay, S.I.; Gething, P.W. Re-examining environmental correlates of Plasmodium falciparum malaria endemicity: A data-intensive variable selection approach. Malar. J. 2015, 14, 68. [Google Scholar] [CrossRef] [PubMed]
  22. World Health Organization (WHO). Global Malaria Programme. Global Technical Strategy for Malaria 2016–2030; WHO Press: Geneva, Switzerland, 2015; pp. 1–15. [Google Scholar]
  23. Ikeda, T.; Behera, S.K.; Morioka, Y.; Minakawa, N.; Hashizume, M.; Tsuzuki, A.; Maharaj, R.; Kruger, P. Seasonally lagged effects of climatic factors on malaria incidence in South Africa. Sci. Rep. 2017, 7, 2458. [Google Scholar] [CrossRef] [PubMed]
  24. Hay, S.I.; Snow, R.W. The malaria atlas project: Developing global maps of malaria risk. PLoS Med. 2006, 3, 2204–2208. [Google Scholar] [CrossRef] [PubMed]
  25. Directorate of Malaria Control Program, Pakistan Medical Research Council (PMRC). Malaria Indicator Survey in 38 High Risk Districts of Pakistan—2013–14; PMRC: Islamabad, Pakistan, 2014. Available online: http://dmc.gov.pk/documents/pdfs/3Final%20MIS%20Report-From%20Printe-30%20April%202015%20(2).pdf (accessed on 5 October 2017).
  26. Hay, S.I.; Guerra, C.A.; Gething, P.W.; Patil, A.P.; Tatem, A.J.; Noor, A.M.; Kabaria, C.W.; Manh, B.H.; Elyazar, I.R.F.; Brooker, S.; et al. A world malaria map: Plasmodium falciparum endemicity in 2007. PLoS Med. 2009, 6, 286–302. [Google Scholar] [CrossRef]
  27. Cohen, J.M.; Le Menach, A.; Pothin, E.; Eisele, T.P.; Gething, P.W.; Eckhoff, P.A.; Moonen, B.; Schapira, A.; Smith, D.L. Mapping multiple components of malaria risk for improved targeting of elimination interventions. Malar. J. 2017, 16, 1–12. [Google Scholar] [CrossRef] [PubMed]
  28. Pakistan Bureau of Statistics (PBS). Press Release on Provisional Summary Results of the 6th Population and Housing Census–2017. PBS: Islamabad, Pakistan, 2017; Government of Pakistan, Ministry of Statistics website. Available online: http://www.statistics.gov.pk/assets/publications/Population_Results.pdf (accessed on 3 October 2017).
  29. Ali, H.; Qamer, F.M.; Ahmed, M.S.; Khan, U.; Habib, A.H.; Chaudhry, A.A.; Ashraf, S.; Khan, B.A. Ecological ranking of districts of Pakistan: A geospatial approach. Pak. J. Bot. 2012, 44, 263–268. [Google Scholar]
  30. Meghani, S.T.; Sehar, S.; Punjani, N.S. Comparison and analysis of healthcare delivery system: Pakistan versus China. Int. J. Endorsing Health Sci. Res. 2014, 2, 46–49. [Google Scholar]
  31. Directorate of Malaria Control Program. Pakistan Malaria Program Review (MPR), 2014. Available online: dmc.gov.pk/documents/pdfs/4MPR%20Pakistan.pdf (accessed on 20 November 2017).
  32. Rowland, M.; Rab, M.A.; Freeman, T.; Durrani, N.; Rehman, N. Afghan refugees and the temporal and spatial distribution of malaria in Pakistan. Soc. Sci. Med. 2002, 55, 2061–2072. [Google Scholar] [CrossRef]
  33. SPSS Inc. SPSS for Windows, version 13.0; SPSS Inc.: Chicago, IL, USA, 2005. [Google Scholar]
  34. Kulldorff, B.M. SaTScanTM User Guide V9.4, 2015. Available online: https://www.satscan.org/cgi-bin/satscan/register.pl/SaTScan_Users_Guide.pdf?todo=process_userguide_download (accessed on 9 September 2017).
  35. Kulldorff, M. A spatial scan statistic. Commun. Stat. Theory Methods 1997, 26, 1481–1496. [Google Scholar] [CrossRef]
  36. Qi, X.; Hu, W.; Page, A.; Tong, S. Spatial clusters of suicide in Australia. BMC Psychiatry 2012, 12, 86. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Environmental Systems Research Institute (ESRI). ArcGIS Desktop, version 10.2; ESRI: Redlands, CA, USA, 2014. [Google Scholar]
  38. Ullah, Z.; Khattak, A.A.; Bano, R.; Hussain, J.; Awan, U.A.; Rahman, S.U.; Mahsud, M.A.J. High incidence of malaria along the Pak-Afghan bordering area. J. Pak. Med. Assoc. 2018, 68, 42–45. [Google Scholar] [PubMed]
  39. Kakar, Q.; Khan, M.A.; Bile, K.M. Malaria control in Pakistan: New tools at hand but challenging epidemiological realities. East Mediterr Health J. 2010, 16, S54–S60. [Google Scholar] [CrossRef] [PubMed]
  40. Mengal, M.H.; Mengal, M.A.; Mengal, M.A.; Rautio, A.; Rehman, F.; Mengal, M.A. Prevalence of drug resistance malaria in Pakistan (Plasmodium vivax and P. falciparum). J. Appl. Emerg. Sci. 2016, 5, 12–19. [Google Scholar]
  41. Khattak, A.A.; Venkatesan, M.; Nadeem, M.F.; Satti, H.S.; Yaqoob, A.; Strauss, K.; Khatoon, L.; Malik, S.A.; Plowe, C.V. Prevalence and distribution of human Plasmodium infection in Pakistan. Malar. J. 2013, 12, 297. [Google Scholar] [CrossRef] [PubMed]
  42. Gov Pakistan Stat Div, Pakistan Bureau of Statistics (PBS). Pakistan Social and Living Standards Measurement Survey (2014–15); PBS: Islamabad, Pakistan, 2016. Available online: http://www.pbs.gov.pk/content/pakistan-social-and-living-standards-measurement (accessed on 2 December 2017).
  43. Karim, A.M.; Hussain, I.; Malik, S.K.; Lee, J.H.; Cho, I.H.; Kim, Y.B.; Lee, S.H. Epidemiology and clinical burden of malaria in the war-torn area, Orakzai Agency in Pakistan. PLoS Negl. Trop. Dis. 2016, 10, e0004399. [Google Scholar] [CrossRef] [PubMed]
  44. Malik Tareen, A.; Tareen, A.M.; Rafique, M.; Wadood, A.; Qasim, M.; Rahman, H.; Shah, S.H.; Khan, K.; Pirkani, G.S. Malaria burden in human population of Quetta, Pakistan. Eur. J. Microbiol. Immunol. 2012, 2, 201–204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Yasinzai, M.I.; Kakarsulemankhel, J.K. Prevalence of human malaria infection in district Jhal Magsi of Balochistan. Rawal Med. J. 2012, 37, 9–12. [Google Scholar]
  46. Yasinzai, M.I.; Kakarsulemankhel, J.K. Frequency of human malaria infection in south-eastern areas of Balochistan: District Lasbella. Pak. J. Med. Sci. 2012, 28, 167–170. [Google Scholar]
  47. Farooq, M.A.; Salamat, A.; Iqbal, M.A. Malaria—An experience at CMH Khuzdar (Balochistan). J. Coll. Phys. Surg. Pak. 2008, 18, 257–258. [Google Scholar]
  48. Yasinzai, M.I.; Kakarsulemankhel, J.K. Incidence of human malaria infection in northern hilly region of Balochistan, adjoining with NWFP, Pakistan: District Zhob. Pak. J. Biol. Sci. 2008, 11, 1620–1624. [Google Scholar] [PubMed]
  49. Awan, Z.U.R.; Shah, H.T.A.; Shah, A.H.; Khan, M.A.; Suleman, M. Malaria among the students of religious schools of Bannu district, Khyber Pakhtunkhwa, Pakistan. Pak. J. Zool. 2012, 44, 959–962. [Google Scholar]
  50. Majid, A.; Rehman, M.U.; Ahmad, T.; Ali, A.; Ali, S. Prevalence of malaria in human population of district Mardan, Pakistan. World J. Zool. 2016, 11, 63–66. [Google Scholar]
  51. Nizamani, M.A.; Kalar, N.A.; Khushk, I.A. Burden of malaria in Sindh, Pakistan: A two years surveillance report. J. Liaquat Univ. Med. Health Sci. 2006, 5, 76–83. [Google Scholar]
  52. Khan, S.Y.; Khan, A.; Arshad, M.; Tahir, H.M.; Mukhtar, M.K.; Ahmad, K.R.; Arshad, N. Irrational use of antimalarial drugs in rural areas of eastern Pakistan: A random field study. BMC Public Health 2012, 12, 941. [Google Scholar] [CrossRef] [PubMed]
  53. Herrel, N.; Amerasinghe, F.P.; Ensink, J.; Mukhtar, M.; van der Hoek, W.; Konradsen, F. Adult anopheline ecology and malaria transmission in irrigated areas of south Punjab, Pakistan. Med. Vet. Entomol. 2004, 18, 141–152. [Google Scholar] [CrossRef] [PubMed]
  54. Zubairi, A.B.; Nizami, S.; Raza, A.; Mehraj, V.; Rasheed, A.F.; Ghanchi, N.K.; Khaled, Z.N.; Beg, M.A. Severe Plasmodium vivax malaria in Pakistan. Emerg. Infect. Dis. 2013, 19, 1851–1854. [Google Scholar] [CrossRef] [PubMed]
  55. Price, R.N.; Douglas, N.M.; Anstey, N.M. New developments in Plasmodium vivax malaria: Severe disease and the rise of chloroquine resistance. Curr. Opin. Infect. Dis. 2009, 22, 430–435. [Google Scholar] [CrossRef] [PubMed]
  56. Khan, M.P.; Ahmad, F.; Ashfaq, A.; Danish, S.H.; Warsi, S.M.A.; Khan, A.I. Malaria and dengue: Assessment of knowledge among households in rural Sindh, Pakistan. J. Dow Univ. Health Sci. 2015, 9, 106–111. [Google Scholar]
  57. Imtiaz, F.; Nisar, N.; Shafi, K.; Nawab, F. Knowledge and practices regarding malaria control and its treatment among patients visiting outpatient clinics of civil hospital Khairpur. Ann. Abbasi Shaheed Hosp. Karachi Med. Dent. Coll. 2015, 21, 221–227. [Google Scholar]
  58. Jamil, S.; Khan, M. Seasonal variations of vivax and falciparum malaria: An observation at a tertiary care hospital. J. Ayub Med. Coll. Abbottabad 2012, 24, 93–95. [Google Scholar] [PubMed]
Ijerph 15 01202 g001 550
Figure 1. The administrative units in Pakistan.
Figure 1. The administrative units in Pakistan.
Ijerph 15 01202 g001
Ijerph 15 01202 g002 550
Figure 2. (A) Raw incidence (per 100,000) of aggregated total positive cases at the district level in Pakistan, and (B) sequence plot of annual cases of malaria by type in Pakistan, 2011–2016.
Figure 2. (A) Raw incidence (per 100,000) of aggregated total positive cases at the district level in Pakistan, and (B) sequence plot of annual cases of malaria by type in Pakistan, 2011–2016.
Ijerph 15 01202 g002
Ijerph 15 01202 g003 550
Figure 3. Raw incidence (per 100,000) of malaria by type at the district level in Pakistan, 2011–2016.
Figure 3. Raw incidence (per 100,000) of malaria by type at the district level in Pakistan, 2011–2016.
Ijerph 15 01202 g003
Ijerph 15 01202 g004 550
Figure 4. Relative risk of malaria by type at the district level in Pakistan, 2011–2016.
Figure 4. Relative risk of malaria by type at the district level in Pakistan, 2011–2016.
Ijerph 15 01202 g004
Ijerph 15 01202 g005 550
Figure 5. High-risk clusters of malaria by type at the district level in Pakistan, 2011–2016.
Figure 5. High-risk clusters of malaria by type at the district level in Pakistan, 2011–2016.
Ijerph 15 01202 g005
Table
Table 1. Summary of descriptive statistics of malaria incidence rate (per 100,000) by type at the district level in Pakistan, 2011–2016.
Table 1. Summary of descriptive statistics of malaria incidence rate (per 100,000) by type at the district level in Pakistan, 2011–2016.
YearPlasmodium vivaxPlasmodium falciparumMixed PlasmodiumTotal Positive
Mean ± SDMin–MaxMean ± SDMin–MaxMean ± SDMin–MaxMean ± SDMin–Max
2011264.43 ± 523.150–3745.63181.17 ± 560.990–5322.227.26 ± 47.130–523.01452.85 ± 866.790–5903.79
2012259.12 ± 629.890–4913.11109.42 ± 265.420–1811.714.25 ± 19.660–189.07372.79 ± 786.090–5805.01
2013439.90 ± 1528.290–15,615.6399.04 ± 251.980–2289.9128.75 ± 105.690–1112.85568.09 ± 1826.100–19,018.40
2014430.25 ± 1202.550–11,509.6970.54 ± 169.640–1198.9226.18 ± 80.300–525.08526.98 ± 1338.770–12,442.91
2015278.85 ± 584.710–3813.3263.42 ± 155.330–1085.4121.47 ± 74.440–505.52363.76 ± 733.580–4335.27
2016617.27 ± 2096.720–19,935.0598.61 ± 245.610–1914.2870.93 ± 215.630–1687.26786.18 ± 2434.840–22,071.55
Mixed Plasmodium: Plasmodium vivax, Plasmodium falciparum. Total positive: P. vivax, P. falciparum and mixed Plasmodium.
Table
Table 2. Summary of primary high-risk clusters’ information of malaria by type at the district level in Pakistan, 2011–2016.
Table 2. Summary of primary high-risk clusters’ information of malaria by type at the district level in Pakistan, 2011–2016.
Type of MalariaYearNumber of Districts in Primary ClusterCoordinates of Primary Cluster CentroidPrimary Cluster Radius (km)Primary Cluster PopulationNumber of Cases in Primary ClusterIncidence Rate/100,000RRp Value
Plasmodium vivax2011533°00′ N, 70°44′ E59.532,175,18143,2111987.922.61<0.0001
20121733°69′ N, 70°31′ E142.4211,009,12863,356574.311.29<0.0001
2013832°65’N, 70°39′ E49.462,725,80559,4382182.023.45<0.0001
20142533°00′ N, 70°44′ E181.7618,172,987113,900627.28.76<0.0001
20152533°00′ N, 70°44′ E181.7618,705,11576,566409.67.98<0.0001
20162533°00′ N, 70°44′ E181.7619,253,991136,436707.110.96<0.0001
Plasmodium falciparum2011830°30′ N, 68°72′ E118.931,293,46321,4471659.233.96<0.0001
20121930°04′ N, 67°82′ E195.575,220,34618,683357.119.17<0.0001
20131930°04′ N, 67°82′ E195.575,412,26111,211207.310.33<0.0001
2014225°40′ N, 63°43′ E75.41646,9994947765.149.20<0.0001
2015134°04′ N, 71°20′ E0854,7774978582.844.10<0.0001
20161930°04′ N, 67°82′ E195.576,039,01813,472222.615.31<0.0001
Mixed Plasmodium2011829°98′ N, 69°54′ E153.447,719,21713,752178.322.21<0.0001
20121325°98′ N, 69°31′ E174.8714,988,376247416.533.96<0.0001
20131225°34′ N, 69°75′ E190.8412,975,991428933.17.63<0.0001
2014731°22′ N, 68°78′ E137.221,173,9363008256.472.87<0.0001
2015731°22′ N, 68°78′ E137.221,205,5871658137.640.95<0.0001
20161432°29′ N, 69°81′ E154.616,303,9839941157.417.96<0.0001
Total positive20111332°29′ N, 69°81′ E153.665,289,93476,4181445.69.96<0.0001
2012933°00′ N, 70°44′ E77.623,986,26649,8461247.914.19<0.0001
2013732°65′ N, 70°39′ E49.462,725,80566,0892426.219.85<0.0001
20141532°29‘ N, 69°81′ E158.426,447,60381,9911272.511.67<0.0001
20152533°00′ N, 70°44′ E181.7618,705,11588,869475.47.17<0.0001
20162533°00′ N, 70°44′ E181.7619,253,991157,613816.99.00<0.0001
RR: relative risk. Mixed Plasmodium: P. vivax, P. falciparum. Total positive: P. vivax, P. falciparum and mixed Plasmodium.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/ijerph15061202

AMA Style

Umer MF, Zofeen S, Majeed A, Hu W, Qi X, Zhuang G. Spatiotemporal Clustering Analysis of Malaria Infection in Pakistan. International Journal of Environmental Research and Public Health. 2018; 15(6):1202. https://doi.org/10.3390/ijerph15061202

Chicago/Turabian Style

Umer, Muhammad Farooq, Shumaila Zofeen, Abdul Majeed, Wenbiao Hu, Xin Qi, and Guihua Zhuang. 2018. "Spatiotemporal Clustering Analysis of Malaria Infection in Pakistan" International Journal of Environmental Research and Public Health 15, no. 6: 1202. https://doi.org/10.3390/ijerph15061202

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