Skip to main content

A time-series prediction model of acute myocardial infarction in northern of Iran: the risk of climate change and religious mourning

A Correction to this article was published on 26 May 2022

This article has been updated



Although various studies have been conducted on the effects of seasonal climate changes or emotional variables on the risk of AMI, many of them have limitations to determine the predictable model. The currents study is conducted to assess the effects of meteorological and emotional variables on the incidence and epidemiological occurrence of acute myocardial infarction (AMI) in Sari (capital of Mazandaran, Iran) during 2011–2018.


In this study, a time series analysis was used to determine the variation of variables over time. All series were seasonally adjusted and Poisson regression analysis was performed. In the analysis of meteorological data and emotional distress due to religious mourning events, the best results were obtained by autoregressive moving average (ARMA) (5,5) model.


It was determined that average temperature, sunshine, and rain variables had a significant effect on death. A total of 2375 AMI’s were enrolled. Average temperate (°C) and sunshine hours a day (h/day) had a statistically significant relationship with the number of AMI’s (β = 0.011, P = 0.014). For every extra degree of temperature increase, the risk of AMI rose [OR = 1.011 (95%CI 1.00, 1.02)]. For every extra hour of sunshine, a day a statistically significant increase [OR = 1.02 (95% CI 1.01, 1.04)] in AMI risk occurred (β = 0.025, P = 0.001). Religious mourning events increase the risk of AMI 1.05 times more. The other independent variables have no significant effects on AMI’s (P > 0.05).


Results demonstrate that sunshine hours and the average temperature had a significant effect on the risk of AMI. Moreover, emotional distress due to religious morning events increases AMI. More specific research on this topic is recommended.

Peer Review reports


Myocardial infarction (MI) as a result of acute coronary syndrome, refers to the permanent and irreversible cellular death of myocardium secondary to loss of blood flow and prolonged lack of oxygen supply (ischemia) [1]. ST-segment elevation myocardial infarction (STEMI), and non-ST-segment elevation are two types of MI. While the cause of mismatch between myocardial oxygen demand and myocardial oxygen consumption in STEMI is nearly always coronary plaque rupture resulting in thrombosis formation occluding a coronary artery, there are several potential causes of this mismatch in NSTEMI [2]. These life-threatening condition may appear suddenly and without any earlier symptoms, or after several angina attacks (chest pain) [3]. While only less than 10% of all deaths worldwide was attributed to cardiovascular disease at the beginning of the twentieth century, the current prevalence of this disease increased as the passage of the time and the impact of the industrial revolution on the social and economic dimensions of life [4]. AMI more than any other disease not only causes deaths but also results in disabilities and economic burdens [5].

While MI is most prevalent in older adults, it is also present in younger people [3]. The prevalence of this disease is increasing in developing countries; which means that 82% of the 16 million worldwide deaths related to AMI results in disability due to this disease [4]. Furthermore, AMI has increased twelve times in women and almost fourteen times in men in developing countries [6]. According to the Iranian Ministry of Health and Medical Education, 39.3% of annual death is due to heart disease and of these, 19.5% result from AMI [1]. The high incidence of AMI in developing countries is associated with less physical activity, weight gain, tobacco use, occupational and mental stress, and low health literacy [7].

Studies over the past decades have addressed various mechanisms and leading factors of AMI [8,9,10]. Risk factors affecting the occurrence of AMI are categorized into non-modifiable (age, gender, family history, weather) and modifiable (smoking, physical activity, psychological tension) factors [11,12,13,14]. Among non-modifiable factors, weather condition has not received much attention in the literature until recently. Few studies have been conducted for more than 50 years around the world [15] to investigate if weather conditions affect AMI. The non-modifiable factors described as “triggers” are defined as external stimuli triggering the internal mechanism that led to AMI. While acute triggers are classified into four categories of emotional, environmental, physical and chemical, in this study environmental and specially meteorological triggers have been further analyzed [16, 17].

Studies revealed that meteorological parameters and seasonal climate changes play an important role in the incidence of AMI [18, 19], so that approximately 4% of the incidence of AMI is associated with various types of these parameters [18,19,20]. Weather conditions such as average temperature, average humidity, wind speed, and wind pressure are identified parameters associated with AMI [7, 21, 22]. While the majority of studies revealed that the AMI incidence increase during winter and spring [15, 23,24,25], some studies showed that AMI occurs in the summer when the temperature and humidity are high and the atmospheric pressure is low [26]. It seems that physiological stressors such as sympathetic activation, high blood coagulation, and infection in cold climates (such as influenza and air pollution) are associated with AMI [24, 27, 28]. Hemodynamic changes are also intensified in the winter. Moreover, an increased immune response, high blood pressure, immobility and the risk of respiratory infections associated with low temperatures in winter can affect the incidence of AMI [29, 30]. The study findings by Versaci et al. [31] revealed that investigating the impact of climate change on coronary atherothrombotic events should be undertaken by considering season-specific patterns and the trends in climate features in the days preceding the acute event. Accordingly, a higher risk of acute myocardial infarction is expected with lower temperatures, lower minimum atmospheric pressure (ATM), and lower rainfall in Winter, greater changes in ATM and greater humidity in Spring, and higher temperatures in Summer [31].

The extensive literature has supported the exposition that pollutant agents such as fine particles (particulate particles with a diameter ≤ μm [PM2.5]) would increase the cardiovascular morbidity and mortality rate [32,33,34]. The yielded result from ‘Global Burden of Disease’ has ranked environmental pollutants among the top 10 determinants of cardiovascular disease out of the investigated 67 risk factors. The results have also highlighted the impact of PM2.5 on the ischemic cardiovascular disease [35].

Although various studies have been conducted on the effects of climate variables and seasonal climate changes on the risk of AMI, many of them have limitations to determine the predictable model. Epidemiological models can predict the prevalence of AMI before its onset using a combination of the variables. There is a balance between the accuracy and timeliness of the predictions made with this method; the closer we get to the time of the occurrence of AMI, the more accurate the predictions are, but less time is available for the necessary action [36]. It should be considered that the relationship between the predictive variables of AMI is very complex. Moreover, the biological, statistical, and mathematical models have been faced by some limitations due to their complexity for clinical applications.

Emotional triggers as the influencing factor of AMI incidence include anger, anxiety or depression secondary to hearing unexpected news, work-related stress, and war. Since triggers reduce sleep and appetite as well as increase in cortisol level, they are considered as factors contributing to cardiovascular disease [37, 38]. Steptoe et al. [39] showed that chronic stress at work, and in private life increases the occurrence of heart disease by 40–50%. Similarly, admission rates for AMI increased after the death of a loved ones or during the bereavement periods [40, 41]. Unlike anger or anxiety, the effect of bereavement is thought to persist for weeks and months [42]. While acknowledging the linkage between psychological stress and AMI. Brotman, Golden, Wittstein [43] discussed that the risk of AMI due to catastrophic emotional events is more prevalent in patients with pre-existing heart disease [43].

Iran is a country of four seasons and the temperature variation from summer to winter is significant. Iran seems an interesting context for the topic of this study as it is one of the only countries in the world that has the complete four seasons. A large part of the country suffers great extremes of heat and cold between summer and winter. Rainfall is mainly confined to winter and spring. Summers are hot with virtually continuous sunshine. Highly humid weather in summer especially along the Caspian Sea in the north and Persian Gulf increases the danger of heat exhaustion. Moreover, the official religion in Iran is Islam and the ruling regime that has tried to inject Islamic laws into every aspect of Iranian life. Therefore, the religious functions, and rituals are held every year in the form of mass folk flow and movement.

Considering the increased prevalence of AMI among in not only the elderly but also the younger generation and its impact on the quality of life requires a better understanding of all risk factors. Furthermore, religion and culture are intertwined with the lifestyle of Iranian people, and religious beliefs, functions, and rituals have an important impact on the lives of Iranian Muslims. The mournful religious activities, as well as nation-wide sadness during some months, are well-known as many people publicly display their mourning. The main objective of this study is to provide a prediction model for the occurrence of AMI using meteorological variables and emotional triggers.



The data of AMI patients was collected from the Mazandaran Province Heart Center, Iran located at the following coordinates (36.369 N, 52.270 W) because it offers the most comprehensive data of AMI patients in Iran [44]. Census sampling method was used between 25/03/2011 and 20/03/2018.

The following were considered as the definitive diagnostic criteria were: (1) existence of cardiac enzymes (CK or CK-MB) above the normal range; (2) ST-segment elevation or depression of more than 1; (3) abnormal Q waves; and (4) manifestation of Troponin enzyme in the blood [1].

The following variables were extracted: gender, the day, month, year and time of hospital admission and recovery situation. Also, weather variables were included daily temperature (Celsius) changes (minimum, maximum, and average), wind speed (meters per second) and its direction, rainfall (day), daily evaporation rate (mm), number of sunny days, and relative humidity (percent) between March 2011 to March 2018 were provided by the Meteorological Organization of Iran. Iran’s four climate seasons are: spring (April to June), summer (July to September), autumn (September to December) and winter (January to March). In this study, TRAMO/SEATS technique was used for all series in order to seasonally adjustment. The process of removing the seasonal component from a time series is known as deseasonalization or seasonal adjustment, and the time series thus obtained is called the deseasonalized, or seasonally adjusted, time series. TRAMO/SEATS method, developed by Gomez and Maravall [45]. TRAMO initially models the series with AR (1) and ARMA (1, 1) to determine the periodic and seasonal difference levels. The appropriate seasonal or non-seasonal ARMA model is selected according to BIC criterion; TRAMO also automatically identifies outliers and calculates other regression variables. Then, TRAMO passes the linearized series to SEATS, where the actual decomposition is done. In SEATS, first the spectral density function of the estimated model is decomposed into the spectral density function of the unobserved components, which are assumed orthogonal. SEATS then estimates the parameters of the two components (trend-cycle and seasonally adjusted component) [46, 47].

The emotional triggers in this study are the religious events in Iran. The ‘event’ variable in the model refers to mournful religious events in Iran.

Statistical method

Box–Jenkins (B–J) estimation models, one of the time series analysis methods, are widely used for modeling and predictions in the field of health and in different applied sciences. B–J method is a linear structure model that traces all the past values of the variable and its stochastic components to predict the values for all the future periods. Prior to have the details of B–J method of forecasting we need to see how a time series data of a particular variable is generated. There are three processes behind the cohort of a time series data.

  1. 1.

    AR (autoregressive) Process: Past values of the variable and error term generate the data

  2. 2.

    MA (moving average) Process: Only the errors or the disturbance term generate the data

  3. 3.

    ARMA (autoregressive and moving average) Process: Data is generated by the combination of AR and MA processes

Sometimes it is taken as ARIMA model where ‘I’ stands for the order of Integration of the series or how many differencing is done for making the time series of the variable to Stationary [48].

B–J methods are known as Autoregressive (AR), Moving Average (MA), Autoregressive Moving Average (ARMA) and Autoregressive Integrated Moving Average (ARIMA) [49, 50]. AR(p), MA (q) and ARMA(p,q) when the process is stationary in the time series that is being examined; ARIMA(p,d,q) can perform applications when it is not stationary [48].

Autoregressive model-AR(p)

The mathematical representation of AR model is given in Eq. 1. In Eq. 1, p is the number of autoregressive components. If there is no relationship between neighboring observation values, p value is zero. An AR(p) process is one where the current or present period’s value of a variable ‘y’ depends on only the past values plus an error term. If there are ‘p’ order in the process i.e. current value of y depends on the ‘p’ order of past (e.g. t − 1, t − 2, etc.) values and an error term of the current period then the AR(p) can be written as:

$$y_{t} = \mu + \varphi_{1} y_{t - 1} + \varphi_{2} y_{t - 2} + \varphi_{3} y_{t - 3} + \cdots + \varphi_{p} y_{t - p} + u_{t} = \mu + \sum \varphi_{i} y_{t - i} + t_{t}$$

In Eq. 1, μ is constant term, ut is White nose error term with zero mean, constant variance and zero auto-covariance [48].

Moving average-MA (q)

An MA(q) process, on the other hand, is the linear combination of all the q terms of the past values of the white noise terms depending on time. When there is no moving average component, q is zero. It is a white noise process in which the current value of yt depends on the current value of the white noise error term (ut) and all past values of the error terms. Because all the errors are white noise, so, an MA process is necessarily a stationary process. It is true further because it is the linear combination of all plus and minus values of the errors which hover around the value zero. The mathematical representation of MA model is given in Eq. 2.

$$y_{t} = u_{t} + \theta_{1} u_{t - 1} + \theta_{2} u_{t - 2} + \theta_{3} u_{t - 3} + \ldots + \theta_{q} u_{t - q} = u_{t} + \sum \theta_{i} u_{t - i}$$

An AR process is stationary if the characteristic root lies outside the unit circle or having values > 1. If it is so then then ϕ becomes less than 1. This means the condition ϕ < 1 lead to the values lying inside the unit circle representing stationarity of the AR process, the model is thus having stability property. The AR coefficients should then be less than unity or they should lie within the unit circle [48].

Autoregressive moving average-ARMA(p,q)

An ARIMA(p,q) process is the combination of AR and MA process, I being the order of integration which can be represented by ‘d’, number of differencing to convert the series from nonstationary to stationary. The model for ARMA(p,d,q) can then be written as,

$$y_{t} = \mu + \varphi_{1} y_{t - 1} + \varphi_{2} y_{t - 2} + \varphi_{3} y_{t - 3} + \cdots + \varphi_{p} y_{t - p} + u_{t} + \theta_{1} u_{t - 1} + \theta_{2} u_{t - 2} + \theta_{3} u_{t - 3} + \cdots + \theta_{q} u_{t - q}$$

Using lag operatör, we have

$$\left( {1 - \varphi_{1} L - \varphi_{2} L^{2} - \varphi_{3} L^{3} - \varphi_{4} L^{4} - \ldots \varphi_{p} L^{p} } \right)y_{t} = \mu + \left( {1 + \theta_{1} L + \theta_{2} L^{2} + \theta_{3} L^{3} + \ldots \theta_{q} L^{q} } \right)u_{t} or, \varphi \left( L \right)y_{t} = \mu + \theta \left( L \right)u_{t}$$

This relation stands (Eq. 4) for invertibility between the AR and MA process which means AR and MA processes can be made invertible from one to another [48, 51].

In this study, an AMI event was analyzed by using Autoregressive Moving Average (ARMA) model because of the time series’ structure is stationary. In order to determine whether the time series is stationary or not, a unit root test was applied to the data. Augmented Dickey Fuller (ADF) unit root test and Philips and Perron (PP) unit root test were used [52]. In the analysis, the optimal structure of ARMA model was discussed in various combinations of p and q values and the appropriate time series model was determined, the number of future AMI’s was estimated. The Ordinary Least Square (OLS) methods were used to estimate the parameters. Akaike Information Criterion (AIC) was used to evaluate the performance of the models [51, 53, 54].

Poisson regression analysis

The count model was applied to analyze the occurrence of AMI. Poisson Regression analysis is a count model method analysis that describes the relationship between independent (reason) variables and the dependent (response) variable obtained by counting. The underlying structure in the Poisson regression analysis is that the Y response variable is the discrete independent Poisson random variable. It is one of the methods shown as an alternative to classical linear regression analysis because of the failure of normality assumption [55].

The Poisson distribution has the following moments,

$${\text{E}}\left( {\text{X}} \right) = \upmu \;{\text{and}}\;{\text{Var}}({\text{X}}) = \upmu$$

The Poisson distribution assumes variance equal to the mean and, hence, it has limitations when dealing with overdispersal data, i.e., when the sample variance exceeds the sample mean. The opposite is called the under-dispersion problem (when the sample mean exceeds the sample variance). In this case negative binomial Poisson regression analysis is used [56].

The Poisson regression model determined the probability that the regressed variable (Y) occurs at a specific time interval.

The probability is modeled as:

$$P\left( {Y = y} \right) = \frac{{e^{ - \mu } \mu^{y} }}{y!}$$

where Y: dependent variable; \(\mu\) = mean parameter [56];

In our study, the dependent variable was number of patients with AMI per day. Explanatory variables were event (mournful religious activities), average temperature, average humidity, sunshine, raining, daily evaporation rate, wind direction, wind speed. The mean of dependent variable (number of patients with AMI) was 3.62 per day, which was similar to the variance result (3.64). If the mean and variance of the dependent variables are equal, the Poisson regression model would best fit the dependent variable of interest. There is not over dispersion and under dispersion problem in our data set. So, we used to Poisson regression analyze. In order to determine the most suitable model, the goodness of fit coefficients was taken into consideration. The analysis was performed with MATLAB R2016a and EViews 9.0.


The sample were 2375 patients with AMI. Table 1 shows the descriptive statistics for the study variables.

Table 1 Descriptive statistics for the study variable from 01 January to 16 March, 2013–2015

In order to perform the analysis of ARMA model, the correlograms of time series were examined with ACF and PACF graphs. Figure 1a and b show the ACF and PACF correlograms for the examined time series, respectively. Boundary lines around functions represent 95% confidence limits. There are no values exceeding the confidence limit in correlograms. In addition, the statistical significance of the autocorrelation coefficients was tested with the Ljung–Box Q statistic. Correlogram graphs and Q statistic values show that the series is stationary (Table 2).

Fig. 1
figure 1

ACF (a) ve PACF (b) figures

Table 2 AC, PAC and Q-stat results for heart attack series

Augmented Dickey Fuller (ADF) unit root test and Philips and Perron (PP) unit root test were used to determine whether the time series is stationary or not. The results of the ADF and PP tests are given in Table 3. Analysis results show that the time series does not contain unit roots and has a stationary structure.

Table 3 ADF and PP unit root tests

Various AR, MA and ARMA models were analyzed with different p and q degrees in order to model the number of heart attacks with time series. Since the correlogram graphs did not provide definite information about AR, MA and ARMA, the models were analyzed one by one until p and q values were in sixth degree. Table 4 shows the AIC results calculated with different p and q values. Accordingly, the model with the lowest AIC value determined the most suitable ARMA structure for the time series examined. Results show that ARMA (5,5) has the lowest AIC value in parameter combinations. The lowest calculated AIC value was calculated as 4.6058. Table 5 shows AR and MA parameter estimation values, standard errors, t-stats and probability values estimated by ARMA (5,5) model. Related numerical findings are given in Table 5.

Table 4 AIC results for ARMA model
Table 5 ARMA model parameter estimations

Autocorrelation and heteroscedasticity of residuals were examined in the process of examining whether the analyzed model was appropriate for the data set. The results of ARMA (5,5) model were investigated with Box–Ljung Q test by using model residuals. Here, H0 hypothesis: Model is appropriate; H1 hypothesis: Model is inappropriate. Ljung–Box Q test results show that ARMA (5,5) model is appropriate. For model residuals, ACF and PACF correlograms are shown in Fig. 2. As can be seen in the graphs, only a single value is out of the confidence limit with a relatively small margin.

Fig. 2
figure 2

Residual’s ACF (a) and PACF (b) graphs

The heteroscedasticity of residuals was examined by the ARCH test. According to ARCH results, there is no heteroscedasticity in the examined data set. As can be seen in the Table 6, the F-stat value was found to be 0.400478, prob. F (1.797) and Prob. Chi-Square 0.5264.

Table 6 ARCH results for ARMA

Omnibus test value is significant (Likelihood Ratio Chi-Square = 24.082, P < 0.001). In this case, the Poisson regression model is statistically significant. Based on Goodness of fit value (Deviance, scaled deviance, AIC, AICC, BIC) estimates, the Poisson regression model fit the data (Table7).

Table 7 Goodness of fit value of Poisson regression model

Table 8 shows the results from Poisson regression analysis. The model was run to estimate the relative risk (RR), which measures the rate at which the number of patients with AMI per day. Table 8 illustrates the RRs estimated across the explanatory variables. Average Temperature variable has significant effect on the AMI (β = 0.011, P = 0.014). For every extra degree of temperature increase 1.011 (95%CI 1.002–1.019) times risk of AMI. For every extra an hour sunshine 1.025 (95% CI 1.010–1.040) times more AMI risk was occurred, a statistically significant result, (β = 0.025, P = 0.001). The risk of AMI will be 0.948 times greater for Event (0) meaning the mournful religious activities can increase the risk of AMI. This means that the zero value of the event (Event = 0) reduces the risk of AMI. The fact that the event takes 1 (Event = 1) is increasing the risk of AMI (1/0.948 = 1.05 times). However, this risk is not significant (P > 0.05). Event meaning the mournful religious activities has not a statistically significant relationship with the number of patients with AMI per day (P > 0.05).

Table 8 Results of negative binomial Poisson regression analysis showing associations between AMI and suspected predictors


In this study, we have analysed the association between meteorological factors, emotional triggers, and epidemiological occurrence of AMI in Sari located in north of Iran during 2011–2018. It appeared that average temperature (°C) and sunshine (hour a day) are independent factors associated with an increased rate of hospital admission due to AMI during the specified period. Several studies have determined that ambient temperature is an important predictor of AMI indicating a U-shape correlation [57,58,59,60]. The findings of a systematic review and meta-analysis has shown that the association of heat exposure and heat wave is associated with increased risk of MI were immediate, whereas, the cold exposure has a delayed impact on MI incidence [61]. Contradictory findings suggest that warmer temperatures are not associated with a higher MI incidence, but the colder temperatures are associated with increased risk of MI [62]. Several underlying mechanisms are hypothesized to explain how ambient temperature is associated with increasing risk of MI that include: sympathetic nervous activity, higher blood pressure, heart rate, left ventricular end-diastolic pressure, and myocardial oxygen consumption, as well as reduced ischemia threshold, the change of hemodynamic, and coagulation [57, 63]. Furthermore, it is indicated that plaque rupture is more frequent in the colder ambient temperature that explains can cold temperatures influence on the pathogenesis of AMI [64].

Although the study has taken into account other meteorological factors such as minimum temperature, maximum temperature, wind speed, wind direction, average humidity, rain and evaporation, none had a significant impact on AMI hospital admission rates. The result yielded from the current study contradicting the findings of Madrigano et al. [7] and Radišauskas et al. [65] who showed a significant relationship between humidity, wind speed and wind pressure and AMI hospital admission. The positive association between average temperature and AMI admission is consistent with Amiya et al. [26]. Findings showing the increased in AMI incidence in summer. However, research findings in other countries have shown AMI admission in winter and spring is higher than other seasons [25, 66, 67]. According to Sharifnia et al. [68] due to geographical location of Sari (Capital of Mazandaran, Iran) a coastal city, there is a substantial change in meteorological parameters during 3 months of spring when the city undergoes a sudden rise in temperature. Higher sunshine duration was associated with higher risk of AMI. The current results provide useful information to alert health professionals to the danger of AMI incidence during this season and to take precautions when the hot season approaches.

The current study suggests an association between religious mourning and increased risk of AMI, but it was not statistically significant. Current evidences regarding the influence of religious practice on MI remains inconclusive. While some longitudinal studies that have been conducted in a large cohort of healthy people didn’t acknowledged the protective association between religious affiliation and practice [69, 70], the other study findings indicated that restful religiousness was longitudinally associated with higher odds of heart attack [71].

The context of the present study is Sari, Capital of Mazandaran state and the collected data covers the climate factors and admission numbers in this city. However, our findings can be generalized to predict the occurrence of AMI in other coastal cities with similar weather characteristics. To the best of the authors’ knowledge, the current results will be informative for people prone to AMI. This information can prepare then to take preventive measures before summer season and control the excessive emotions while attending the religious ceremonies.


Like many other research studies, our study is not free from limitations. Lack of access to more detailed health information of the patients such as past medical history, blood pressure and many other factors remains one of the primary limitations. Secondary, apart from meteorological and emotional factors, there are vast number of factors affecting AMI that were beyond the scope of our current research. Thirdly, the study is only limited to a specific province. Therefore, the results should be interpreted with caution.

Nevertheless, the large sample size remains as a unique feature of this study. The data for meteorological factors was extracted from the Central Station Monitoring for climate condition instead of measurements of exposure to environmental variables. This provided us the independent variables (meteorological data) that was unbiased with regards to outcome of this study.


Further studies are necessary in order to quantify the emotional variables with more accurate measures. Studies with samples from different regions and also longitudinal designs are suggested to verify the findings of this study. Future studies are recommended to incorporate more detailed patients’ information about all the recognized risk factors such as BMI, blood pressure, blood cholesterol, cigarette smoking, past medical history, BUN, Creatinine.


The results demonstrate that sunshine hours and average temperature had significant effects on risk of AMI.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Change history


  1. Sharif Nia H, Haghdoost AA, Nazari R, Bahrami N, Soleimani MA, Pormand K. Relationship of risk factors and ST segment changes with symptoms of acute coronary syndrome. Koomesh. 2013;15(1):46–53.

    Google Scholar 

  2. Vallabhajosyula S, Vallabhajosyula S, Burstein B, et al. Epidemiology of in-hospital cardiac arrest complicating non–ST-segment elevation myocardial infarction receiving early coronary angiography. Am Heart J. 2020;223:59–64.

    PubMed  Google Scholar 

  3. Mitka M. Heart disease a global health threat. JAMA. 2004;291(21):2533–2533.

    CAS  PubMed  Google Scholar 

  4. World Health Report. Reducing risks, promoting healthy life. Geneva: World Health Organization; 2002.

    Google Scholar 

  5. Raymond DBM. The early heart attack care program catheterization. J Cardiovasc Diagn. 2000;53(3):157.

    Google Scholar 

  6. Ahmadi A, Sajjadi H, Etemad K, Khaledifar A, Mobasherii M. Epidemiological characteristics and determinants of mortality in acute coronary syndrome in Iran. J Mazandaran Univ Med Sci. 2015;25(124):1–9.

    Google Scholar 

  7. Madrigano J, Mittleman MA, Baccarelli A, et al. Temperature, myocardial infarction, and mortality: effect modification by individual and area-level characteristics. Epidemiology. 2013;24(3):439–46.

    PubMed  PubMed Central  Google Scholar 

  8. Honda T, Fujimoto K, Miyao Y. Influence of weather conditions on the frequent onset of acute myocardial infarction. J Cardiol. 2016;67(1):42–50.

    PubMed  Google Scholar 

  9. Moschos N, Christoforaki M, Antonatos P. Seasonal distribution of acute myocardial infarction and its relation to acute infections in a mild climate. Int J Cardiol. 2004;93(1):39–44.

    PubMed  Google Scholar 

  10. Sen T, Astarcioglu MA, Asarcikli LD, et al. The effects of air pollution and weather conditions on the incidence of acute myocardial infarction. Am J Emerg Med. 2016;34(3):449–54.

    PubMed  Google Scholar 

  11. Čulić V. Acute risk factors for myocardial infarction. Int J Cardiol. 2007;117(2):260–9.

    PubMed  Google Scholar 

  12. Longo DL, Fauci AS, Kasper DL, Hauser SL, Jameson JL, Loscalzo J. Harrison’s manual of medicine. New York: McGraw-Hill Medical; 2013.

    Google Scholar 

  13. Urden LD, Stacy KM, Lough ME. Critical care nursing-e-book: diagnosis and management. Amsterdam: Elsevier Health Sciences; 2017.

    Google Scholar 

  14. Marx J, Walls R, Hockberger R. Rosen’s emergency medicine-concepts and clinical practice e-book. Amsterdam: Elsevier Health Sciences; 2013.

    Google Scholar 

  15. Kriszbacher I, Bódis J, Csoboth I, Boncz I. The occurrence of acute myocardial infarction in relation to weather conditions. Int J Cardiol. 2009;135(1):136–8.

    PubMed  Google Scholar 

  16. Edmondson D, Newman JD, Whang W, Davidson KW. Emotional triggers in myocardial infarction: Do they matter? Eur Heart J. 2012;34(4):300–6.

    PubMed  PubMed Central  Google Scholar 

  17. Glozier N, Tofler GH, Colquhoun DM, et al. Psychosocial risk factors for coronary heart disease. Med J Aust. 2013;199(3):179–80.

    PubMed  Google Scholar 

  18. Goggins WB, Chan EY, Yang C. Weather, pollution, and acute myocardial infarction in Hong Kong and Taiwan. Int J Cardiol. 2013;168(1):243–9.

    PubMed  Google Scholar 

  19. Patel NJ, Pant S, Deshmukh AJ, et al. Seasonal variation of acute myocardial infarction related hospitalizations in the United States: perspective over the last decade. Int J Cardiol. 2014;172(3):441–2.

    Google Scholar 

  20. Čulić V. Seasonal distribution of acute myocardial infarction: a need for a broader perspective. Int J Cardiol. 2006;109(2):265–6.

    PubMed  Google Scholar 

  21. Messner T, Lundberg V, Wikström B. A temperature rise is associated with an increase in the number of acute myocardial infarctions in the subarctic area. Int J Circumpolar Health. 2002;61(3):201–7.

    PubMed  Google Scholar 

  22. Radišauskas R, Bernotienė G, Bacevičienė M, Ustinavičienė R, Kirvaitienė J, Krančiukaitė-Butylkinienė D. Trends of myocardial infarction morbidity and its associations with weather conditions. Medicina. 2014;50:182–9.

    PubMed  Google Scholar 

  23. Morabito M, Modesti PA, Cecchi L, et al. Relationships between weather and myocardial infarction: a biometeorological approach. Int J Cardiol. 2005;105(3):288–93.

    PubMed  Google Scholar 

  24. Rumana N, Kita Y, Turin TC, et al. Seasonal pattern of incidence and case fatality of acute myocardial infarction in a Japanese population (from the Takashima AMI Registry, 1988 to 2003). Am J Cardiol. 2008;102(10):1307–11.

    PubMed  Google Scholar 

  25. Spencer FA, Goldberg RJ, Becker RC, Gore JM. Seasonal distribution of acute myocardial infarction in the second National Registry of Myocardial Infarction. J Am Coll Cardiol. 1998;31(6):1226–33.

    CAS  PubMed  Google Scholar 

  26. Amiya S, Nuruki N, Tanaka YT, et al. Relationship between weather and onset of acute myocardial infarction: Can days of frequent onset be predicted? Int J Cardiol. 2009;54(2):231–7.

    Google Scholar 

  27. Fares A. Winter cardiovascular diseases phenomenon. N Am J Med Sci. 2013;5(4):266–79.

    PubMed  PubMed Central  Google Scholar 

  28. Sheth T, Nair C, Muller J, Yusuf S. Increased winter mortality from acute myocardial infarction and stroke: the effect of age. J Am Coll Cardiol. 1999;33(7):1916–9.

    CAS  PubMed  Google Scholar 

  29. Jorge JEL, Cagy M, Mesquita ET, Costa TLM, Moscavitch SD, Rosa MLG. Seasonal variation in hospitalizations due to heart failure in Niterói city, Southeastern Brazil. Rev Saude Publ. 2009;43(3):555–7.

    Google Scholar 

  30. Lin GM, Li YH, Lin CL, Wang JH, Han C. Seasonal variation in cardiac death of patients with angiographic coronary artery disease from the ET-CHD registry, 1997–2006. Int J Cardiol. 2013;168(3):3090–2.

    PubMed  Google Scholar 

  31. Versaci F, Biondi-Zoccai G, Giudici AD, et al. Climate changes and ST-elevation myocardial infarction treated with primary percutaneous coronary angioplasty. Int J Cardiol. 2019;294:1–5.

    PubMed  Google Scholar 

  32. Brook RD, Rajagopalan S, Pope CA III, et al. Particulate matter air pollution and cardiovascular disease: an update to the scientific statement from the American Heart Association. Circulation. 2010;121(21):2331–78.

    CAS  PubMed  Google Scholar 

  33. Wang M, Hou Z-H, Xu H, et al. Association of estimated long-term exposure to air pollution and traffic proximity with a marker for coronary atherosclerosis in a nationwide study in China. JAMA Netw Open. 2019;2(6):e196553–e196553.

    PubMed  PubMed Central  Google Scholar 

  34. Stanley Young S, Kindzierski WB. Evaluation of a meta-analysis of air quality and heart attacks, a case study. Crit Rev Toxicol. 2019;49(1):85–94.

    CAS  PubMed  Google Scholar 

  35. Landrigan PJ, Fuller R, Acosta NJ, et al. The lancet commission on pollution and health. Lancet. 2018;391(10119):462–512.

    PubMed  Google Scholar 

  36. Matola Y, White G, Magayuka S. The changed pattern of malaria endemicity and transmission at Amani in the eastern Usambara mountains, north-eastern Tanzania. J Trop Med Hyg. 1987;90(3):127–34.

    CAS  PubMed  Google Scholar 

  37. Khan MS, Aouad R. The effects of insomnia and sleep loss on cardiovascular disease. Sleep Med Clin. 2017;12(2):167–77.

    PubMed  Google Scholar 

  38. von Känel R, Princip M, Schmid J-P, et al. Association of sleep problems with neuroendocrine hormones and coagulation factors in patients with acute myocardial infarction. BMC Cardiovasc Disord. 2018;18(1):1–12.

    Google Scholar 

  39. Steptoe A, Kivimäki M. Stress and cardiovascular disease. Nat Rev Cardiol. 2012;9(6):360.

    CAS  PubMed  Google Scholar 

  40. Schwartz B, French W, Mayeda G, et al. Emotional stressors trigger cardiovascular events. Int J Clin Pract. 2012;66(7):631–9.

    CAS  PubMed  Google Scholar 

  41. Mittleman M, Maclure M, Sherwood J, Kondo H, Tofler G, Muller J. Death of a significant person increases the risk of acute MI onset. Paper presented at: Circulation. 1996.

  42. Moon JR, Kondo N, Glymour MM, Subramanian S. Widowhood and mortality: a meta-analysis. PLoS ONE. 2011;6(8):e23465.

    CAS  PubMed  PubMed Central  Google Scholar 

  43. Brotman DJ, Golden SH, Wittstein IS. The cardiovascular toll of stress. The Lancet. 2007;370(9592):1089–100.

    Google Scholar 

  44. Goudarzian AH, Sharif Nia H, Jafari H, et al. Inpatient satisfaction with health system transformation project in Mazandaran Educational Hospitals, Iran. J Mazandaran Univ Med Sci. 2016;26(136):190–5.

    Google Scholar 

  45. Gómez V, Maravall Herrero A. Programs TRAMO and SEATS: instructions for the user (beta version: September 1996). Banco de España: Servicio de Estudios; 1995.

    Google Scholar 

  46. Barska M. Seasonality testing for macroeconomic time series–comparison of X-12-ARIMA and TRAMO/SEATS procedures. Śląski Przegląd Statystyczny. 2014;18(12):121–40.

    Google Scholar 

  47. Atuk O, Ural BP. Seasonal adjustment methods: an application to the Turkish monetary aggregates. Centr Bank Rev. 2012;2(1):21–37.

    Google Scholar 

  48. Das RC. Forecasting incidences of COVID-19 using Box–Jenkins method for the period July 12-Septembert 11, 2020: a study on highly affected countries. Chaos Solitons Fract. 2020;140:110248.

    Google Scholar 

  49. Makridakis S, Hibon M. ARMA models and the Box–Jenkins methodology. J Forecast. 1997;16(3):147–63.

    Google Scholar 

  50. Bustos OH, Yohai VJ. Robust estimates for ARMA models. J Am Stat Assoc. 1986;81(393):155–68.

    Google Scholar 

  51. Meko DM. Applications of Box–Jenkins methods of time series analysis to the reconstruction of drought from tree rings. 1981.

  52. Gujarati D, Porter D. Basic econometrics. London: McGraw-Hill Higher Education; 2003.

    Google Scholar 

  53. Gujarati DN. Basic econometrics. New Delhi: Tata McGraw-Hill Education; 2009.

    Google Scholar 

  54. Aygören H, Sarıtaş H, Moralı T. Forecasting ISE 100 indice using artifıcial neural networks and newton numerical search models. Int J Alanya Fac Bus. 2012;4(1):73–88.

    Google Scholar 

  55. Frome EL, Kutner MH, Beauchamp JJ. Regression analysis of Poisson-distributed data. J Am Stat Assoc. 1973;68(344):935–40.

    Google Scholar 

  56. Bermúdez L, Karlis D, Santolino M. A discrete mixture regression for modeling the duration of non-hospitalization medical leave of motor accident victims. Accid Anal Prev. 2018;121:157–65.

    PubMed  Google Scholar 

  57. Danet S, Richard F, Montaye M, et al. Unhealthy effects of atmospheric temperature and pressure on the occurrence of myocardial infarction and coronary deaths. A 10-year survey: the Lille-World Health Organization MONICA project (Monitoring trends and determinants in cardiovascular disease). Circulation. 1999;100(1):E1-7.

    CAS  PubMed  Google Scholar 

  58. Braga AL, Zanobetti A, Schwartz J. The effect of weather on respiratory and cardiovascular deaths in 12 U.S. cities. Environ Health Perspect. 2002;110(9):859–63.

    PubMed  PubMed Central  Google Scholar 

  59. Panagiotakos DB, Chrysohoou C, Pitsavos C, et al. Climatological variations in daily hospital admissions for acute coronary syndromes. Int J Cardiol. 2004;94(2–3):229–33.

    PubMed  Google Scholar 

  60. Dilaveris P, Synetos A, Giannopoulos G, Gialafos E, Pantazis A, Stefanadis C. CLimate impacts on myocardial infarction deaths in the Athens TErritory: the CLIMATE study. Heart. 2006;92(12):1747–51.

    CAS  PubMed  PubMed Central  Google Scholar 

  61. Sun Z, Chen C, Xu D, Li T. Effects of ambient temperature on myocardial infarction: a systematic review and meta-analysis. Environ Pollut. 2018;241:1106–14.

    CAS  PubMed  Google Scholar 

  62. García-Lledó A, Rodríguez-Martín S, Tobías A, Alonso-Martín J, Ansede-Cascudo JC, de Abajo FJ. Heat waves, ambient temperature, and risk of myocardial infarction: an ecological study in the Community of Madrid. Rev Española Cardiol (Engl Ed). 2020;73(4):300–6.

    Google Scholar 

  63. Marchant B, Ranjadayalan K, Stevenson R, Wilkinson P, Timmis AD. Circadian and seasonal factors in the pathogenesis of acute myocardial infarction: the influence of environmental temperature. Heart. 1993;69(5):385–7.

    CAS  Google Scholar 

  64. Katayama Y, Tanaka A, Taruya A, et al. Increased plaque rupture forms peak incidence of acute myocardial infarction in winter. Int J Cardiol. 2020;320:18–22.

    PubMed  Google Scholar 

  65. Radišauskas R, Bernotienė G, Bacevičienė M, Ustinavičienė R, Kirvaitienė J, Krančiukaitė-Butylkinienė D. Trends of myocardial infarction morbidity and its associations with weather conditions. Medicina. 2014;50(3):182–9.

    PubMed  Google Scholar 

  66. Hernández EG, O’Callaghan AC, Doménech JC, et al. Seasonal variations in admissions for acute myocardial infarction. The PRIMVAC study. Rev Española Cardiol. 2004;57(1):9–12.

    Google Scholar 

  67. Alijani B. Natural seasons in Iran. J Geograph Res Tehran Univ. 1998;35:21–33.

    Google Scholar 

  68. Nia HS, Chan YH, Froelicher ES, et al. Weather fluctuations: predictive factors in the prevalence of acute coronary syndrome. Health Promot Perspect. 2019;9(2):123.

    Google Scholar 

  69. Schnall E, Wassertheil-Smoller S, Swencionis C, et al. The relationship between religion and cardiovascular outcomes and all-cause mortality in the women’s health initiative observational study. Psychol Health. 2010;25(2):249–63.

    PubMed  Google Scholar 

  70. Salmoirago-Blotcher E, Fitchett G, Hovey KM, et al. Frequency of private spiritual activity and cardiovascular risk in postmenopausal women: the women’s health initiative. Ann Epidemiol. 2013;23(5):239–45.

    PubMed  PubMed Central  Google Scholar 

  71. Ahrenfeldt LJ, Hvidt NC, Kjøller ST, Möller S, Lindahl-Jacobsen R. Religiousness and diseases in Europe: findings from SHARE. J Relig Health. 2019;58(6):1925–37.

    CAS  PubMed  Google Scholar 

Download references


We thank all of the participants that contribute to completing this study. In addition, we acknowledge Miss Ameneh Yaghoobzadeh for her help. Furthermore, we express our gratitude to Mazandaran University of Medical Sciences for supporting.


Not applicable.

Author information

Authors and Affiliations



The authors confirm contribution to the paper as follows: Study conception and design: HSN, NN, AAH. Data collection: AHG, RP. Analysis and interpretation of results: OG, SPS, HSN. Draft manuscript preparation: ESF, OG, NN, FKF. All authors reviewed the results and approved the final version of the manuscript.

Corresponding author

Correspondence to Fatemeh Khoshnavay Fomani.

Ethics declarations

Ethics approval and consent to participate

The study was conducted in accordance with the Declaration of Helsinki and with Ethics Committee of Mazandaran University of Medical Sciences, Sari, Iran approvals and was registered as IR.MAZUMS.REC.96–10232. Written informed consent was obtained from participants prior to enrolment. All patients who referred to the Heart Center of Mazandaran Province signed an informed- consent form to use their data in research project. The informed-written consents are available.

Consent for publication

Not applicable. The manuscript does not contain any individual person’s data.

Competing interests

The authors didn’t have any financial conflicts to declare with any persons or organizations regarding to this study.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

The original online version of this article was revised: The name of the corresponding author has been corrected

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Sharif Nia, H., Gorgulu, O., Naghavi, N. et al. A time-series prediction model of acute myocardial infarction in northern of Iran: the risk of climate change and religious mourning. BMC Cardiovasc Disord 21, 563 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Weather
  • Acute myocardial infarction
  • Time series
  • Iran