 Research note
 Open Access
 Published:
Predicting the incidence of brucellosis in Western Iran using Markov switching model
BMC Research Notes volume 14, Article number: 79 (2021)
Abstract
Objective
Brucellosis is a zoonosis almost chronic disease. Brucellosis bacteria can remain in the environment for a long time. Thus, climate irregularities could pave the way for the survival of the bacterium brucellosis. Brucellosis is more common in men 25 to 29 years of age, in the western provinces, and in the spring months. The aim of this study is to investigate the effect of climatic factors as well as predicting the incidence of brucellosis in Qazvin province using the Markov switching model (MSM). This study is a secondary study of data collected from 2010 to 2019 in Qazvin province. The data include brucellosis cases and climatic parameters. Two state MSM with time lags of 0, 1 and 2 was fitted to the data. The Bayesian information criterion (BIC) was used to evaluate the models.
Results
According to the BIC, the twostate MSM with a 1month lag is a suitable model. The month, the averagewindspeed, the minimumtemperature have a positive effect on the number of brucellosis, the age and rainfall have a negative effect. The results show that the probability of an outbreak for the third month of 2019 is 0.30%.
Introduction
Brucellosis is a zoonosis almost chronic disease which is transmitted by direct or indirect contact with infected animals or products [1]. B. melitensis is the most common and acute pathogen that can cause various symptoms such as sudden tremors, general body aches [2,3,4]. Despite a good health care system in Iran, brucellosis is still an important endemic disease. Iran ranks fourth in the world of the incidence of brucellosis [5,6,7]. According to the National Guidelines, epidemiological and clinical data of patients in Iran are recorded online in the health monitoring system. Patients with epidemiological and clinical signs of brucellosis, i.e. fever, myalgia, results of two laboratory tests are considered as cases. Laboratory tests include the Wright test (> 1.8 indicating infection) and the 2ME (Mercaptoethanol Brucella agglutination test).
In Iran, the incidence of brucellosis is reported to be between 50 and 500 per 100,000 people, often of the B. melitensis type and is known as a local disease. The incidence rate of brucellosis in Qazvin is 27.43 per 100,000 during 2010–2019 [8].
In epidemiology, a sudden increase in the incidence of a disease at a specific time or place is called an outbreak. In this study, the outbreak of the disease means a sudden increase in the disease in Qazvin. Climatic factors directly and indirectly affect the growth and reproduction dynamics of pets, sheep/goat human interactions, pathogen replication and population immunity affect brucellosis ecology [9,10,11,12]. Temperature changes may have a greater impact on the brucellosis epidemic than other climatic factors [13].
Despite the significant effect of climatic factors on brucellosis, these factors have been overlooked in most studies. The primary purpose of this study was to investigate the effect of climatic factors on the incidence of brucellosis. The second purpose of this study is to use MSM for prediction. A number of researchers have used the MSM in disease prediction [14] but so far the performance of this method has not been evaluated in brucellosis data.
Main text
Method
This study is a secondary study of data collected from April 2010 to March 2019 in Qazvin, which is extracted from the database of the Qazvin University of Medical Sciences and the meteorological system of the province [8]. No individual data were used. The available information was cumulative without mentioning personal information.
The response variable is the number of brucellosis cases. Month, ruralratio, age, menratio, the ratio of contact with livestock, Nonpasteurized dairy and factors in the meteorological system include averagemonthlytemperature, totalmonthlyrainfall, averagewindspeed, maximummonthlytemperature, minimummonthlytemperature and windspeed are considered as explanatory variables. The number of cases of brucellosis is very low in some months and if we consider daily or weekly cases, the number of cases will be zero on some times. Also, because brucellosis bacteria survive in the environment for a long time, considering monthly cases will help to better understand the factors affecting brucellosis.
In this study, the number of monthly cases of human brucellosis in Qazvin is a dependent variable (Y_{t}). Independent variables (X_{t}s) with time lags of 0, − 1, and − 2 were considered in twostates and threestates models. \({{\varvec{\sigma}}}_{{\varvec{s}}{\varvec{t}}}\) is the standard deviation of the error in the model.
MSM including time series {Y_{t}}_{t=1,…,T} and a sequence of related variables x_{1},…,x_{T} is introduced with the relation between x_{t} and Y_{t} as follows:
where \({{\varvec{\epsilon}}}_{{\varvec{t}}}\) has a standard normal distribution and S_{t} is the state at time t of a nonobservable N state Markov chain.
According to the fitted MSM, beta is the estimation of the effect of each variable on the response in each state. SE is the standard deviation of beta and pvalue shows the significance of each variable in the relevant state.
In twostate models, if we define state 1 as the disease outbreak period and 2 state as the nonoutbreak period, the probability of an outbreak in period t + 1 can be as follows [14]:
A twostate MSM with switching all effects is considered. The data were analyzed through SPSS software version 26 and MSwM package of R software 3.6.3 [15].
Results
From 2010 to 2019, 3194 people were infected with brucellosis. Of these, the highest incidence related to the year 2015 with 512 (16%), and the lowest incidence related to the year 2010 with 192 (6%). In matching the seasons with the Gregorian months, spring is related to March, April and May, summer is related to June, July and August, fall is related to September, October, and November and finally, winter is related to December, January and February. Among seasons, the summer with 961 (30.1%) and winter with 805 (25.2%) have the highest number of infected, respectively. The highest number of patients is related to temperatures 26 ℃ (6.6%), 25 ℃ (6.4%). The highest number of infected people is related to zero rainfall of 444 cases (13.9%) and the lowest number of patients is related to the averagemonthlyrainfall of 17.8 (0.2%) and 64.5 (0.2%). The lowest number of patients is related to the averagewindspeed of 0.7 with 6 cases (0.2%).
Fitting MSM
The MSM was fitted with two and three states and both models were fitted with 0, − 1, and − 2 lags for climate variables. Temporal lag is defined as the time interval between climatic characteristics and the incidence of brucellosis. Based on a comparison between models, based on BIC, twostate MSM with a time lag of − 1 is suitable. For this reason, only this model is offered to provide more results (Fig. 1).
Age, month, ruralratio, menratio, Nonpasteurized dairy, averagemonthlytemperature, totalmonthlyrainfall, averagewindspeed, maximummonthlytemperature, Minimummonthlytemperature, averagewindspeed, ratio of contact with livestock were recognized as significant variables (Table 1).
Autocorrelation and partial autocorrelation of residual and squared residual for model is confirmed lack of autocorrelation in the residual and the model seems to fit logically and there's no serial dependency on the residual.
Figure 2 shows the smoothed and filtered probabilities for state of one and two. Smoothed probabilities are used to determine peaks and depressions and 0.5 is determined as the cutoff value for 1 and 2 states. The filtered probabilities are calculated using the first observation up to t and the smoothed probabilities are calculated using the total observations.
The Q–Q plot shows where the normality hypothesis is questionable for series. Transition probability matrix in MSM as follows:
where the probability of nonoutbreak state in both t and t + 1 periods is 0.72, the probability of changing series from nonoutbreak state in period t to outbreak state in period t + 1 is 0.28 and the probability of the series changing from the outbreak state in period t to nonoutbreak state in period t + 1 is 0.65. When we are in a nonoutbreak state (0.72), the process tends to stay the same state and the process is transferred to the outbreak state with a probability of 0.28.
The probability of an outbreak in t + 1 is as follows:
Since the data is up to the second month of 2019, the probability of an outbreak for the third month of 2019 (1 month later) is very low and is equal to 0.30%.
The biggest difference between the coefficients of the variables in two states is related to the averagewindspeed. Therefore, the averagewindspeed is the most important factor in incidence brucellosis.
Month, averagewindspeed and minimumtemperature coefficients are positive which indicate a positive effect on the number of brucellosis. The age and totalmonthlyrainfall coefficients are negative, indicating a negative effect on the number of brucellosis.
The temporal changes of the observed cases of brucellosis and the values estimated by the MSM are illustrated in Additional file 1: Figure S1. The model has a relatively good performance in identifying peaks incidence of brucellosis.
Discussion
The highest incidence of brucellosis is related to 2015 with 512 cases (16%). The average wind speed was 1.89. The incidence is the highest at zero totalmonthlyrainfall. Totalmonthlyrainfall was 0 at 444 days (13.9%), which includes most days of study. The minimumtemperature was − 1 and the maximumtemperature was 30 ℃. The mean age of the patients was about 38 years.
There is no clear pattern in the number of cases of brucellosis in the 8 years studied and fluctuations in the incidence of this disease can be seen with three peaks 2015 December, July 2015, June 2014. These results are inconsistent with the results of Lee's study. In the Lee’s study, the incidence of human brucellosis in South Korea peaked in September 2006 and has dropped dramatically which indicates effective eradication [16]. In Rafiemanesh’s study, the incidence of brucellosis decreased from 2007 to 2016 which indicates an increase in the coverage of prevention programs, especially livestock vaccination [17]. These results are inconsistent with the results of the present study.
In the 2010 year, the lowest number of cases has been reported, followed by an upward trend until 2011. The reason for the rising trend of the disease from 2010 to 2011 may be related to the improvement of the data registration in the country's health system. This result is consistent with Hashtarkhani’s study [18]. From 2014 to 2017, there is an upward trend in the number of brucellosis cases which is inconsistent with Hashtarkhani's study. In the study of Hashtarkhani after 1990, we see a decreasing trend in the incidence of the disease [18]. The results of the study show that there are the highest number of infected people in summer and winter seasons. These results are consistent with the results of the Tapak’s study. The results of the Tapak’s study show that hot summers and cold winters make the disease less common while climate moderation in these seasons exacerbates the disease [19]. Therefore, the temperate climate of these seasons in Qazvin increases the number of patients with brucellosis. Model fit results indicate the negative effect of age and totalmonthlyrainfall on the number of brucellosis. This result is consistent with the results of Entezari's study. The results of the Entezari's study indicate a negative relationship between rainfall and brucellosis. In fact, as the rainfall decreases, the number of infected people increases [20].
In zenoses, changes in climatic factors naturally affect the contamination rate and dynamics of animal hosts as well as human exposure to infected animals [21].
High levels of evaporation and sun exposure cause drought and limit the germination of the plant, while dry environments may cause human skin dryness and cracked skin and increase the risk [13].
Comparing the coefficients of explanatory variables in two state, the averagewindspeed is the most important factor in incidence brucellosis. This result is consistent with the results of Ahmadkhani’s study. The results of Ahmadkhani's study indicate a positive correlation between windspeed, temperature, greenness and incidence of brucellosis [22]. The results of the present study are inconsistent with those of Tapak. According to the Tapak’s study, the wind at high speeds reduces the disease. This is because the bacterium has a shorter lifespan in the air.
The mean age of the patients was about 38 years. This result has been confirmed in other studies [23, 24]. The results of this study indicate that P_{11} and P_{22} are larger than P_{21} and P_{12}, respectively. In other words, states do not tend to change. That is, when we are in a nonoutbreak state, the process tends to remain the same state. Also, P_{21} is higher than P_{12}, which indicates that prolonged periods of nonoutbreak lead to a reduction in the probability of outbreak during the year.
As the country's health progresses in many areas, the incidence of brucellosis is expected to decline. However, the results of the present study indicate a sharp increase in the disease between 2014 and 2017, which requires a lot of health attention in Qazvin province. The necessary learning for highrisk age and occupational groups and not consuming unpasteurized dairy products and not having contact with suspicious animals and cooperating with livestock vaccination should be on the agenda of the region's health institutions.
Conclusion
The MSM can be used to detect factors related to the incidence of brucellosis as well as to predict the incidence of brucellosis. Most climatic parameters were effective in incidence the disease, and the most influential factor was the averagewindspeed. The probability of disease outbreak in the third month of 2019 was predicted to be 0.30%.
Limitation
One of the limitations of this study is the limited period of the time series data and lack of daily information. Another limitation is the lack of comparison between different time series models.
Availability of data and materials
The dataset used and/or analysed during the current study are available from the Zahra Hosseinkhani on reasonable request.
Abbreviations
 MSM:

Markov switching model
References
Corbel M. Brucellosis in humans and animals. Geneva: World Health Organization; 2006.
Franco M, Mulder M, Gilman R, Smits H. Human brucellosis. Lancent Infect Dis. 2007;7(12):775–86.
Silva T, Costa E, Paixão T, Tsolis R, Santos R. Laboratory animal models for brucellosis research. J Biomed Biotechnol. 2011;2011:518323.
Zeinali M, Shirzadi M, Sharifian J. National guideline for Brucellosis control. Tehran: Ministry of Health and Medical Education; 2009. p. 10–7.
Haran M, Agarwal A, Kupfer Y, Seneviratne C, Chawla K, Tessler S. Brucellosis presenting as septic shock. BMJ Case Rep. 2011;2011.
Pappas G, Papadimitriou P, Akritidis N, Christou L, Tsianos E. The new global map of human brucellosis. Lancet Infect Dis. 2006;6(2):91–9.
Ramin B, MacPherson P. Human brucellosis. BMJ. 2010;341:c4545.
Bagheri H, Tapak L, Karami M, Hosseinkhani Z, Najari H, Karimi S, et al. Forecasting the monthly incidence rate of brucellosis in west of Iran using time series and data mining from 2010 to 2019. PLoS ONE. 2020;15(5):e0232910.
RodriguezMorales JA. Climate change, climate variability and brucellosis. Recent Pat Antiinfect Drug Discov. 2013;8(1):4–12.
Li Y, Li X, Liang S, Fang L, Cao W. Epidemiological features and risk factors associated with the spatial and temporal distribution of human brucellosis in China. BMC Infect Dis. 2013;13:547.
Yang Z, Li X, Shao Z, Ma W, Yuan X, Wu K, et al. Characteristics on spatial and temporal distribution as well as the driving effect of meteorological factors on brucellosis in Datong city, Shanxi province, 2005–2015. Zhonghua Liu Xing Bing Xue Za Zhi. 2018;39(9):1165–71.
Zhu H, Wang L, Lin D, Hong R, Ou J, Chen W, et al. Analysis on epidemiology and spatialtemporal clustering of human brucellosis in Fujian province, 2011–2016. Chin J Epidemiol. 2017;38(9):1212–7.
Liu K, Yang Z, Liang W, Guo T, Long Y, Shao Z. Effect of climatic factors on the seasonal fluctuation of human brucellosis in Yulin, northern China. BMC Public Health. 2020;20(1):506.
Ansari H, Mansournia M, Izadi S, Zeinali M, Mahmoodi M, HolakouieNaieni K. Predicting CCHF incidence and its related factors using timeseries analysis in the southeast of Iran: comparison of SARIMA and Markov switching models. Epidemiol Infect. 2015;143(4):839–50.
SanchezEspigares J, LpezMoreno A. MSwM: Fitting Markov switching models. R package version. 2014;1.4.
Lee H, Her M, Levine M, Moore G. Time series analysis of human and bovine brucellosis in South Korea from 2005 to 2010. Prev Vet Med. 2013;110(2):190–7.
Rafiemanesh H, Alimohamadi Y, Aghdam S, Safarzadeh A, Shokri A, Zemestani A. Time series and trend analysis of brucellosis in Oskou county, East Azerbaijan: 2007–2016. Health Promot Perspect. 2019;9(4):285–90.
Hashtarkhani S, Akbari M, Jarahi L, Etminani K. Epidemiological characteristics and trend of incidence of Human brucellosis in Razavi Khorasan province. Med J Mashhad Univ Med Sci. 2015;58(9):531–8.
Tapak L, ShirmohammadiKhorram N, Hamidi O, Maryanaji Z. Predicting the frequency of human brucellosis using climatic indices by three data mining techniques of radial basis function, multilayer perceptron and nearest Neighbor: a comparative study. IRJE. 2018;14(2):153–65.
Entezari M, Moradpour S, Amiri M. Spatial distribution and the impact of geographical factors on brucellosis in Chaharmahal and Bakhtiari Province, Iran. Int J Epidemiol Res. 2016;3(2):98–105.
AlTalafhah A, Lafi S, AlTarazi Y. Epidemiology of ovine brucellosis in Awassi sheep in Northern Jordan. Prev Vet Med. 2003;60(4):297–306.
Ahmadkhani M, Alesheikh A. Spacetime analysis of human brucellosis considering environmental factors in Iran. Asian Pac J Trop Dis. 2017;7(5):257–65.
AlTawfiq J, AbuKhamsin A. A 24year study of the epidemiology of human brucellosis in a healthcare system in Eastern Saudi Arabia. J Infect Public Health. 2009;2(2):81–5.
Buzgan T, Karahocagil M, Irmak H, Baran A, Karsen H, Evirgen O, et al. Clinical manifestations and complications in 1028 cases of brucellosis: a retrospective evaluation and review of the literature. Int J Infect Dis. 2010;14(6):e469–78.
Acknowledgements
We would like to thank the staffs of the Department of Zoonotic Diseases at the Center for Communicable Disease Control, Ministry of Health and Medical Education of Iran (MOHME). Moreover, we thank the ViceChancellor of Research and Technology of Hamadan University of Medical Sciences.
Funding
This study (980217943, IR.UMSHA.REC.1398.085) was funded by the ViceChancellor of Research and Technology of Hamadan University of Medical Sciences. The funders had no role in study design, data collection, and analysis, decision to publish, or preparation of the manuscript.
Author information
Authors and Affiliations
Contributions
Conceptualization: MMKh, MS, ZCh, ZH. Data curation: ZCh, ZH. Formal analysis: MMKh. Funding acquisition: ZCh. Methodology: MMKh, MS. Writing—original draft: MMKh, MS. Writing—review and editing: MMKh, MS, ZCh, ZH. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
We declare that we have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1.
Figure S1: Prediction values obtained using Markov switching model along with the observed values.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
MohammadianKhoshnoud, M., Sadeghifar, M., Cheraghi, Z. et al. Predicting the incidence of brucellosis in Western Iran using Markov switching model. BMC Res Notes 14, 79 (2021). https://doi.org/10.1186/s13104020054155
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13104020054155