



ORIGINAL ARTICLE 

Year : 2021  Volume
: 14
 Issue : 10  Page : 463470 

Prediction of malaria cases in the southeastern Iran using climatic variables: An 18year SARIMA time series analysis
Hamid Reza Tohidinik^{1}, Hossein Keshavarz^{2}, Mehdi Mohebali^{2}, Mandana Sanjar^{3}, Gholamreza Hassanpour^{2}
^{1} HIV/STI Surveillance Research Center, and WHO Collaborating Center for HIV Surveillance, Institute for Futures Studies in Health, Kerman University of Medical Sciences, Kerman, Iran ^{2} Center for Research of Endemic Parasites of Iran (CREPI); Department of Medical Parasitology & Mycology, School of Public Health, Tehran University of Medical Sciences, Tehran, Iran ^{3} Department of Medical Parasitology & Mycology, School of Public Health, Tehran University of Medical Sciences, Tehran, Iran
Date of Submission  01May2021 
Date of Decision  21Oct2021 
Date of Acceptance  26Oct2021 
Date of Web Publication  29Oct2021 
Correspondence Address: Gholamreza Hassanpour Center for Research of Endemic Parasites of Iran (CREPI); Department of Medical Parasitology & Mycology, School of Public Health, Tehran University of Medical Sciences, Tehran Iran Hossein Keshavarz Center for Research of Endemic Parasites of Iran (CREPI); Department of Medical Parasitology & Mycology, School of Public Health, Tehran University of Medical Sciences, Tehran Iran
Source of Support: None, Conflict of Interest: None
DOI: 10.4103/19957645.329008
Objective: To predict future trends in the incidence of malaria cases in the southeast of Iran as the most important area of malaria using Seasonal Autoregressive Integrated Moving Average (SARIMA) model, and to check the effect of meteorological variables on the disease incidence. Methods: SARIMA method was applied to fit a model on malaria incidence from April 2001 to March 2018 in Sistan and Baluchistan province in southeastern Iran. Climatic variables such as temperature, rainfall, rainy days, humidity, sunny hours and wind speed were also included in the multivariable model as covariates. Then, the best fitted model was adopted to predict the number of malaria cases for the next 12 months. Results: The bestfitted univariate model for the prediction of malaria in the southeast of Iran was SARIMA (1,0,0)(1,1,1)_{12} [Akaike Information Criterion (AIC)=307.4, validation root mean square error (RMSE)=0.43]. The occurrence of malaria in a given month was mostly related to the number of cases occurring in the previous 1 (p=1) and 12 (P=1) months. The inverse number of rainy days with 8month lag (β=0.329 2) and temperature with 3month lag (β=0.002 6) were the best predictors that could improve the predictive performance of the univariate model. Finally, SARIMA (1,0,0)(1,1,1)_{12} including mean temperature with a 3month lag (validation RMSE=0.414) was selected as the final multivariable model. Conclusions: The number of malaria cases in a given month can be predicted by the number of cases in the prior 1 and 12 months. The number of rainy days with an 8month lag and temperature with a 3month lag can improve the predictive power of the model. Keywords: Malaria; Time series; SARIMA; Forecasting; Climate; Iran
How to cite this article: Tohidinik HR, Keshavarz H, Mohebali M, Sanjar M, Hassanpour G. Prediction of malaria cases in the southeastern Iran using climatic variables: An 18year SARIMA time series analysis. Asian Pac J Trop Med 2021;14:46370 
How to cite this URL: Tohidinik HR, Keshavarz H, Mohebali M, Sanjar M, Hassanpour G. Prediction of malaria cases in the southeastern Iran using climatic variables: An 18year SARIMA time series analysis. Asian Pac J Trop Med [serial online] 2021 [cited 2022 May 26];14:46370. Available from: https://www.apjtm.org/text.asp?2021/14/10/463/329008 
1. Introduction   
Malaria is a serious vectorborne disease in the world, particularly in the tropical and subtropical regions^{[1]}. It is caused by protozoan parasites of the genus Plasmodium, which are transmitted across humans by infected female Anopheles mosquitoes^{[2]}. The disease is the fifth cause of disabilityadjusted life years in children younger than 10 years in 2019^{[3]} and remained an important cause of mortality and morbidity in many parts of the world where can have negative effects on the health and socioeconomic status of the population^{[4],[5]}. According to the World malaria report 2018, there were 219 million cases of malaria globally in 2017 and 435 000 malaria deaths. The majority of these cases (92%) were in the African Region, followed by the SouthEast Region with 5% of the cases and the Eastern Mediterranean Region with 2%^{[1]}.
Iran is one of the malariaendemic countries in the world. In spite of successful implementation of malaria control measures that led to a significant reduction in the incidence of the disease over the last decade^{[6],[7]}, it has remained a major public health concern in Iran. It can be explained by drug resistance^{[8],[9]} and the extremely genetic polymorphic nature of Plasmodium spp.^{[10],[11],[12]}. More than 90% of malaria cases occur in the southern and southeastern regions that share borders with Pakistan and Afghanistan^{[13]}. Sistan and Baluchistan province is an endemic area of malaria in Iran with an appropriate climate for breeding various malaria vectors^{[7]}. The province is susceptible to malaria epidemics due to several factors including climate diversity^{[14]}.
Temperature, precipitation, and humidity are important meteorological factors associated with the population dynamics of malaria vectors which could affect the disease spread^{[15],[16]}.
Early detection, prevention, and containment of malaria epidemics is one of the four principal components of the WHO’s global malaria control strategy^{[17]}. Therefore, developing predictive models is an essential part of malaria surveillance that enables policymakers and public health staff to predict future incidence of the disease and act proactively^{[18]}. Seasonal Integrated Moving Average (SARIMA) model^{[19]} is widely used to predict different infectious diseases including malaria^{[20],[21],[22],[23]}. Some statistical models have been adopted for malaria in some regions of Iran^{[14],[15],[16],[18],[19],[20],[21],[22],[23],[24],[25],[26]}, but to the best of our knowledge, no study had applied SARIMA time series to predict malaria incidence in Sistan and Baluchistan province. Therefore, the objective of this study was to provide a SARIMA time series model for the prediction of malaria incidence in the southeast of Iran, and to check if the inclusion of climatic variables enhances the predictive power of the model.
2. Subjects and methods   
2.1. Ethical approval
The study protocol was approved by ethical committee of Tehran University of Medical Sciences with registration number IR.TUMS. SPH.REC13970270.
2.2. Study area
Sistan and Baluchistan province, in the southeast of Iran (28.5° N, 60.5° E, [Figure 1]), is the second largest province of Iran with an area of 180 726 km^{2} and a population of 2.8 million. The climate is relatively hot and dry in most areas but has high humidity in the southern coastal regions. Mean temperature ranges from 17 °C in winter to 32 °C in summer, and mean annual precipitation of 58 mm range from 0.8 mm per month in summer to 12.3 mm in winter months.  Figure 1: Geographical location of the study area in the southeast of Iran.
Click here to view 
2.3. Data collection
We used the monthly number of confirmed cases of malaria from April 2001 to March 2019 that was provided by the Centers for Surveillance of Communicable Diseases in Zahedan, Iranshahr and Zabol Universities of Medical Sciences. We retrieved the monthly meteorological data for the same period from Meteorological Bureau of Sistan and Baluchistan province. We gathered monthly average temperature (°C), average minimum temperature (°C), average maximum temperature (°C), average wind speed, average number of rainy days, accumulated rainfall (millimeters), average relative humidity (%), and average sunny hours. We considered data from April 2001 to March 2018 as training data to fit the model, and kept the number of new malaria cases during April 2018March 2019 as the validation data.
2.4. Statistical analysis
We described number of malaria cases in different months/seasons as median (Q1, Q3). Spearman correlation coefficient was used to assess the correlation between number of malaria cases and meteorological variables. To find the best model fitted to our data and to predict the number of malaria cases, we applied SARIMA (p, d, q) (P, D, Q)s model using BoxJenkins approach in which p is nonseasonal autoregressive (AR) order, d is the number of nonseasonal differencing passes, q is nonseasonal moving average (MA) order, and P, D, Q are corresponding seasonal orders. In this model, s is the seasonal period of the data (12 months in this study).
To fit the model, we first used disease trend plot and BoxCox test to check the stationarity in the variance. Because of nonstationarity in the variance (θ=0), we used natural log transformed (Ln) of the number of malaria cases in the model. Second, DickyFuller test was used to evaluate stationarity in the means of series. As the test showed stationarity in means (P=0.07), we did not apply any nonseasonal differencing pass on data (d=0) but to adjust seasonality of malaria cases and meteorological variables, first order seasonally differencing was applied to the model (D=1).
Third, we plotted autocorrelation function (ACF) and partial autocorrelation function (PACF) plots on stationary data to identify possible values of MA (q, Q) and AR (p, P) components, respectively. Fourth, we applied likelihood ratio test to estimate parameters of SARIMA model. We checked the goodnessoffit of each model, by plotting ACF and PACF plots for residuals to check their normality, and using Ljungbox (Q) test for residuals to check whether they are white noise with mean of 0 and constant variance.
In the next step, we incorporated meteorological variables into the final SARIMA model to check their ability to improve the predictive power. To do so, we first removed autocorrelation within each individual series via prewhitening procedure in which SARIMA model was applied on each series of meteorological variables to remove their seasonal trend. Then, the crosscorrelation between residuals of models was computed and climatic variables significantly correlated with the number of malaria cases were checked as potential covariates in multiple SARIMA model. Variance inflation factor >5 was considered as collinearity between climatic variables. To find the best fitted model, we compared R^{2}, the Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) in different SARIMA models, where the highest R^{2} and lowest AIC and BIC showed better fitness in the models.
We finally predicted 12month number of malaria cases during April 2018March 2019 for selected SARIMA models and compared observed and outofsample predicted values. The root mean square error (RMSE) was computed for both training and validation data as an index for validity of predictions in the models. Smaller values of this index show better predictive power of the model. The root mean square error equals to:
with Y_{t} as the observed values, Yt as the predicted values at time t, and N as the number of observations. We used STATA 12 (Stata Corp, College Station, TX, USA) for all statistical analyses and considered a twosided P<0.05 as statistically significant.
3. Results   
3.1. Malaria cases
From April 2001 to March 2019, a total of 85 378 malaria cases including 60 589 indigenous and 24 789 imported cases were reported from Sistan and Baluchistan province. During this period, Plasmodium vivax and Plasmodium falciparum were the most common types of Plasmodium in the area [Appendix Figure 1]. In addition to seasonal pattern, the disease showed a rising trend since 2001 until 2009 with a peak in 2003, but it showed a downward trend from 2009 onward [Appendix Figure 2].
The monthly distribution showed that the highest number of malaria cases occurred in October with 684 (139, 887) [median (Q1, Q3)] cases followed by September with median 648 (117, 912) [median (Q1, Q3)] cases, and the lowest number happened in February with median 25 (8, 35) [median (Q1, Q3)] cases per month. Transmission of malaria is possible in all seasons of the year, but most cases have occurred in summer and then in autumn, respectively [Figure 2].  Figure 2: Boxplot of the number of malaria cases in different months from 2001 to 2018 in Sistan and Baluchistan province, Iran.
Click here to view 
3.2. The impact of climate on malaria
Spearman correlation coefficients show that number of malaria cases was directly associated with average temperature, maximum and minimum temperature, wind speed and sunny hours (P<0.01). It was also inversely associated with accumulated rainfall, rainy days and relative humidity (P<0.01) [Appendix Table 1].
[Figure 3]A and [Figure 3]B present the ACF and PACF plots using data from 2001 to 2018. ACF suggested q≤S2 and Q≤S1 while PACF revealed that p≤S1 and P≤S2.  Figure 3: Autocorrelation (A) and partial autocorrelation (B) functions based on the seasonally differenced, intransformed number of malaria cases from 2001 to 2018 in Sistan and Baluchistan province, Iran; Autocorrelation suggests q≤2, Q≤1, and partial autocorrelation suggests p≤1 and P≤2.
Click here to view 
Various SARIMA models were assessed and the bestfitted univariate model was SARIMA (1,0,0)(1, 1, 1)_{12} (AIC=307.4, BIC=323.7, training RMSE=0.52, validation RMSE=0.48). Some selected SARIMA models are depicted in [Table 1].  Table 1: Comparison of five candidate univariable SARIMA models, with different values of p, P, q, and Q, for the number of malaria cases in the southeast of Iran (April 2001March 2018).
Click here to view 
Coefficients and statistics of the parameters of SARIMA (1,0,0) (1, 1, 1)_{12} have been presented in [Table 2] where firstorder nonseasonal and seasonal AR coefficients, and the firstorder seasonal MA coefficient are significant (P<0.001).  Table 2: Coefficients and parameters of the SARIMA (1,0,0)(1,1,1)_{12} model to predict the number of malaria cases in the southeast of Iran (April 2001March 2018).
Click here to view 
The components of the best fitted SARIMA models for meteorological variables in prewhitening phase are shown in [Appendix Table 2].
After crosscorrelation analysis, we found average temperature with a 3month lag and inverse number of rainy days with an 8month lag as the best predictors that can relatively improve the predictive performance of the univariate model. [Table 3] revealed a nonsignificant reverse relation between mean temperature with a 3month lag and the number of malaria cases, and a significant direct association between inverse number of rainy days with an 8month lag and the occurrence of the disease. For identifying the best model for prediction, we compared the univariate model with multivariable models and SARIMA (1,0,0)(1, 1, 1)_{12} including mean temperature with a 3month lag, which had lower RMSE (0.414), was selected as the final multiple model. The distribution of residuals of the model was normal [Appendix Figure 3] and LjungBox test confirmed that residuals are white noise with zero mean and constant variance (P=0.54).  Table 3: Characteristics of different multivariable SARIMA models including meteorological variables to predict the incidence of malaria in the southeast of Iran (April 2001 to March 2018).
Click here to view 
The comparison between the predicted and observed monthly number of malaria cases of both univariate and multiple models are presented in [Table 4]. [Figure 4] shows the observed numbers and predicted values for malaria cases from 2010 to 2018 and outofsample predictions in 2019 applying SARIMA (1,0,0) (1, 1, 1)_{12}+temperature with the 3monthlag model. Finally, as illustrated in [Appendix Figure 4], we predicted the number of malaria cases for the next Iranian year (April 2019March 2020) in Sistan and Baluchistan province. We predicted that totally, about 174 cases of malaria would occur in this province during April 2019 to March 2020. Among them, about eight cases would occur in first three months of 2020.  Figure 4: (A) observed and predicted number of malaria cases from 2010 to 2019, (B) outofsample predictions for April 2018March 2019 in the southeast of Iran, applying SARIMA(1,0,0)(1,1,1)_{12}+temperature with 3month lag.
Click here to view 
 Table 4: Predicted numbers of malaria cases during April 2018March 2019 in the southeast of Iran, calculated from two different SARIMA (p,d,q)(P,D,Q)_{12} models.
Click here to view 
4. Discussion   
The results of this study showed that malaria in Sistan and Baluchistan province had a significant decreasing trend during 2001 to 2019, which could be due to improving preventive interventions such as the distribution of insecticidetreated mosquito nets and indoor and outdoor spraying^{[27]}, early detection tools, especially the use of rapid diagnostic tests, effective vector control, and proper implementation of modern malaria control programs^{[7]}.
Our study showed that the number of autochthonous malaria cases in Sistan and Baluchistan province starts to increase from May and reaches its peak in September and October, and then declines and reaches its lowest level in February. The main malaria vectors in Sistan and Baluchistan province are primarily Anopheles (An.) culicifacies and An. stephensi^{[28]}. An. culicifacies has two main activity peaks in Sistan and Baluchistan province; the main peak occurs in June and the second in October^{[29],[30]}. It takes two weeks for Anopheles mosquitoes to complete their life cycle^{[31]}. Incubation period of Plasmodium parasite may last 822 days in the mosquito’s body^{[32]} and 1214 days in the human host^{[33]}. It seems that during the first peak, Anopheles mosquitoes have ample time to generate a large population of infectious vectors, leading to an increase in the number of malaria cases in the following months.
This study also indicated that the occurrence of the disease was possible in all seasons of the year, but the least transmission occurred in winter. It confirms that the main vectors of malaria in this area are active almost all year long, and are able to transmit the disease^{[29],[31],[32],[33]}.
In this study, SARIMA (1,0,0)(1, 1, 1)_{12} model as univariate and SARIMA (1,0,0)(1, 1, 1)_{12} including mean temperature with a 3month lag as multivariable model were selected for the prediction of the number of malaria cases in Sistan and Baluchistan province while the latter had predictions closer to the observed number of malaria cases. It seems that the model could provide an acceptable prediction in Sistan and Baluchistan province.
Ostovar et al. in Hormozgan province, southern Iran^{[23]}, Breit et al. in Sri Lanka^{[34]} and Abeku et al. in Ethiopia^{[35]} used ARIMA model for predicting malaria cases. Ebhuoma et al. in KwazulaNatal, South Africa, suggested the SARIMA (0, 1, 1)(0, 1, 1)_{12} as the best model for predicting malaria^{[36]}. Wangdi et al. in Bhutan^{[31]} and Kumar et al. in New Delhi^{[16]} selected ARIMA (2, 1, 1)(0, 1, 1)_{12} and ARIMA (0, 1, 1) (0, 1, 0)_{12} to forecast malaria, respectively.
In our study, mean temperature with a 3month lag showed a reverse relationship with the number of malaria cases. The optimum temperature for Anopheles mosquitoes is 20 °C to 30 °C. At these temperatures, the Anopheles can live long enough to acquire and transmit the parasite^{[37]}. Increasing the temperature to above 30 °C has negative effects on the mosquitoes’ growth and reduces the propagation rate of Plasmodium in the vector’s body^{[38]}, leading to reduced ability of mosquitoes to transmit the disease, which decreases the number of malaria cases consequently.
In the study of Ostovar et al. in Hormozgan province, the temperature with a 2month lag had a negative relationship with malaria cases^{[23]}, while Mohammadkhani et al. in Sistan and Baluchistan province^{[26]} and Haghdoost et al. in Kahnooj in the southeast of Iran^{[39]} reported that the temperature with a 1month lag was directly related to malaria cases. This delay is reported as one month in China^{[38],[40]} and Bhutan^{[31]}, and two months in Tibet^{[41]}.
Applying a relatively long period of time and including a variety of climatic variables in the prediction model of the disease were among the strengths of the present study. On the other hand, we were unable to include some factors such as socioeconomic factors, population movement, immunity status, sanitation condition, implementation of malaria control programs and ecological factors in the model which could be considered in future studies about this disease.
Malaria showed a decreasing trend in the southeast of Iran as one of the traditional endemic areas of the disease in the region. SARIMA time series model including climatic variables was able to forecast the number of malaria cases with relatively good accuracy in this region.
Conflict of interest statement
The authors declare that they have no conflict of interest.
Acknowledgement
This study was financially supported by the Tehran University of Medical Sciences (project No: 970316040156).
Authors’ contributions
H.R.T., H.K. and Gh.H. developed the theoretical formalism, H.R.T., and M.S. performed the analytic calculations. All authors contributed to the final version of the manuscript. Gh.H. supervised the project.
References   
1.  
2.  Carter R, Mendis KN. Evolutionary and historical aspects of the burden of malaria. Clin Microbiol Rev 2002; 15(4): 564594. 
3.  The Global Burden of Disease. 2019 diseases and injuries collaborators. Global burden of 369 diseases and injuries in 204 countries and territories, 19902019: A systematic analysis for the global burden of disease study 2019. Lancet 2020; 396(10258): 12041222. 
4.  Zoghi S, Mehrizi AA, Raeisi A, Haghdoost AA, Turki H, Safari R, et al. Survey for asymptomatic malaria cases in low transmission settings of Iran under elimination programme. Malar J 2012; 11: 126. 
5.  Hassanpour G, Mohebali M, Mirhendi H, Zeraati H, Raeisi A, Keshavarz H. Asymptomatic malaria infections among immigrants in malariaelimination programmed areas of south eastern Iran may threaten malaria eradication. Travel Med Infect Dis 2019; 31: 101426. 
6.  Edrissian GH. Malaria in Iran: Past and present situation. Iranian J Parasitol 2006; 1(1):114. 
7.  Vatandoost H, Raeisi A, Saghafipour A, Nikpour F, Nejati J. Malaria situation in Iran: 20022017. Malar J 2019; 18(1): 200. 
8.  Eskandarian AA, Keshavarz H, Mahboudi F. Pyrimethamine and proguanil resistanceconferring key mutations in dihydrofolate reductase gene of Plasmodium falciparum isolates from Iran. Iran J Public Health 1970; 32(4): 812. 
9.  Heidari A, Dittrich S, Jelinek T, Kheirandish A, Banihashemi K, Keshavarz H. Genotypes and in vivo resistance of Plasmodium falciparum isolates in an endemic region of Iran. Parasitol Res 2007; 100(3): 589592. 
10.  Heidari A, Keshavarz H, Dittrich S, Ebrahimzadeh A, Jelinek T. Genotyping of Plasmodium falciparum field isolates in major endemic region of Iran and potential uses in identification of field strains. J Med Sci 2007; 7(2): 228232. 
11.  Ehtesham R, Fazaeli A, Raeisi A, Keshavarz H, Heidari A. Detection of mixedspecies infections of Plasmodium falciparum and Plasmodium vivax by nested PCR and rapid diagnostic tests in southeastern Iran. Am J Trop Med Hyg 2015; 93(1): 181185. 
12.  Heidari A, Keshavarz H, Rokni MB, Jelinek T. Genetic diversity in merozoite surface protein (MSP)1 and MSP2 genes of Plasmodium falciparum in a major endemic region of Iran. Korean J Parasitol 2007; 45(1): 5963. 
13.  Abbasi M, HanafiBojd AA, YaghoobiErshadi M, Vatandoost H, Oshaghi M, Hazratian T, et al. Resistance status of main malaria vector, Anopheles stephensi Liston (Diptera: Culicidae) to insecticides in a malaria Endemic Area, Southern Iran. Asian Pac J Trop Med 2019; 12(1): 4348. 
14.  Salehi M, Mohammad K, Farahani MM, Zeraati H, Nourijelyani K, Zayeri F. Spatial modeling of malaria incidence rates in Sistan and Baluchistan province, Islamic Republic of Iran. Saudi Med J 2008; 29(12): 17911796. 
15.  Ye Y, Louis VR, Simboro S, Sauerborn R. Effect of meteorological factors on clinical malaria risk among children: An assessment using villagebased meteorological stations and communitybased parasitological survey. BMC Public Health 2007; 7: 101. 
16.  Kumar V, Mangal A, Panesar S, Yadav G, Talwar R, Raut D, et al. Forecasting malaria cases using climatic factors in delhi, India: A time series analysis. Malar Res Treat 2014; 2014: 482851. 
17.  
18.  Cunha GB, LuitgardsMoura JF, Naves EL, Andrade AO, Pereira AA, Milagre ST. Use of an artificial neural network to predict the incidence of malaria in the city of Canta, state of Roraima. Rev Soc Bras Med Trop 2010; 43(5): 567570. 
19.  Ebhuoma O, Gebreslasie M, Magubane L. A seasonal autoregressive integrated moving average (SARIMA) forecasting model to predict monthly malaria cases in KwaZuluNatal, South Africa. S Afr Med J 2018; 108(7): 573578. 
20.  Tohidinik HR, Mohebali M, Mansournia MA, Niakan Kalhori SR, AliAkbarpour M, Yazdani K. Forecasting zoonotic cutaneous leishmaniasis using meteorological factors in eastern Fars province, Iran: A SARIMA analysis. Trop Med Int Health 2018; 23(8): 860869. 
21.  Jayaraj VJ, Avoi R, Gopalakrishnan N, Raja DB, Umasa Y. Developing a dengue prediction model based on climate in Tawau, Malaysia. Acta Trop 2019; 197: 105055. 
22.  Cong J, Ren M, Xie S, Wang P. Predicting seasonal influenza based on SARIMA model, in mainland China from 2005 to 2018. Int J Environ Res Public Health 2019; 16(23): 4760. 
23.  Ostovar A, Haghdoost AA, Rahimiforoushani A, Raeisi A, Majdzadeh R. Time series analysis of meteorological factors influencing malaria in South Eastern Iran. J Arthropod Borne Dis 2016; 10(2): 222236. 
24.  Raeisi A, Gouya MM, Nadim A, Ranjbar M, Hasanzehi A, Fallahnezhad M, et al. Determination of malaria epidemiological status in Iran’s malarious areas as baseline information for implementation of malaria elimination program in Iran. Iran J Public Health 2013; 42(3): 326333. 
25.  World Health Organization & Global Partnership to Roll Back Malaria. Using climate to predict infectious disease outbreaks: A review. 2004. [Online]. Available from: https://apps.who.int/iris/handle/10665/84175. [Accessed on 20 April 2021]. 
26.  Mohammadkhani M, Khanjani N, Bakhtiari B, Tabatabai SM, Sheikhzadeh K. The relation between climatic factors and malaria incidence in Sistan and Baluchestan, Iran. SAGE Open 2019; 9(3): 2158244019864205. 
27.  Sedaghat MM, Nadim A, Goudarzi A, Holakoei Naeini K, Nateghpour M, Ladoni H, et al. Entomological research and its specific course in Iran from 1935 to 2008. J School Public Health Institute Public Health Res 2013; 11(1): 99112. 
28.  HanafiBojd AA, Vatandoost H, YaghoobiErshadi MR. Climate change and the risk of malaria transmission in Iran. J Med Entomol 2020; 57(1): 5064. 
29.  Vatandoost H, Emami SN, Oshaghi MA, Abai MR, Raeisi A, Piazzak N, et al. Ecology of malaria vector Anopheles culicifacies in a malarious area of Sistan va Baluchestan province, southeast Islamic Republic of Iran. East Mediterr Health J 2011; 17(5): 439445. 
30.  SaneiDehkordi A, SoleimaniAhmadi M, Jaberhashemi SA, Zare M. Species composition, seasonal abundance and distribution of potential anopheline vectors in a malaria endemic area of Iran: Field assessment for malaria elimination. Malar J 2019; 18(1): 157. 
31.  Wangdi K, Singhasivanon P, Silawan T, Lawpoolsri S, White NJ, Kaewkungwal J. Development of temporal modelling for forecasting and prediction of malaria infections using timeseries and ARIMAX analyses: A case study in endemic districts of Bhutan. Malar J 2010; 9: 251. 
32.  Yang HM. Malaria transmission model for different levels of acquired immunity and temperaturedependent parameters (vector). Rev Saude Publica 2000; 34(3): 223231. 
33.  White NJ, Pukrittayakamee S, Hien TT, Faiz MA, Mokuolu OA, Dondorp AM. Malaria. Lancet 2014; 383(9918): 723735. 
34.  Briet OJ, Vounatsou P, Gunawardena DM, Galappaththy GN, Amerasinghe PH. Models for short term malaria prediction in Sri Lanka. Malar J 2008; 7: 76. 
35.  Abeku TA, de Vlas SJ, Borsboom G, Teklehaimanot A, Kebede A, Olana D, et al. Forecasting malaria incidence from historical morbidity patterns in epidemicprone areas of Ethiopia: A simple seasonal adjustment method performs best. Trop Med Int Health 2002; 7(10): 851857. 
36.  Ebhuoma O, Gebreslasie M, Magubane L. Modeling malaria control intervention effect in KwaZuluNatal, South Africa using intervention time series analysis. J Infect Public Health 2017; 10(3): 334338. 
37.  McMichael A. The health impacts of global climate change: Grappling with scenarios, predictive models, and multiple uncertainties. Ecosystem Health 1995; 1: 10. 
38.  Bi P, Tong S, Donald K, Parton KA, Ni J. Climatic variables and transmission of malaria: A 12year data analysis in Shuchen County, China. Public Health Rep 2003; 118(1): 6571. 
39.  Haghdoost AA, Alexander N, Cox J. Modelling of malaria temporal variations in Iran. Trop Med Int Health 2008; 13(12): 15011508. 
40.  Zhang Y, Bi P, Hiller JE. Meteorological variables and malaria in a Chinese temperate city: A twentyyear timeseries data analysis. Environ Int 2010; 36(5): 439445. 
41.  Huang F, Zhou S, Zhang S, Wang H, Tang L. Temporal correlation analysis between malaria and meteorological factors in Motuo County, Tibet. Malar J 2011; 10: 54. 
[Figure 1], [Figure 2], [Figure 3], [Figure 4]
[Table 1], [Table 2], [Table 3], [Table 4]
