 Research
 Open Access
 Published:
The Ebola contagion and forecasting virus: evidence from four African countries
Health Economics Reviewvolume 5, Article number: 16 (2015)
Abstract
This paper is focused on examining the number of deaths’ increases participation in the propagating the Ebola virus during the period ranging from March to October 2014. An application of the MGARCHDCC model regressions on four countries has led to discover that the finding that human contact play a significant role in transmitting the Ebola virus. Our findings also reveal that Guinea has already suffered from a spreadlike virus originating from Sierra Lione and Liberia. Noteworthy also, other countries are now liable to such a risk; for instance, Nigeria is a country vulnerable to the propagation of this virus. Consequently, we undertake to conduct our forecasts for EGARCH model estimates implements; which has estimated a decrease in the Ebola virus incurred number of deadly Ebola virus over the two months following the November and December.
Background
The Ebola virus belongs to a family of viruses called filovirus^{a}, commonly carrying the Marburg virus. Discovered in 1976 in the western equatorial province of Sudan and in a nearby region of northern Zaire (now Democratic Republic of Congo (DRC)) [1,2], the Ebola virus had affected several African countries between 1976 and 2008 (Centers for Disease Control and Prevention, 2014; [3,5]; Leroy et al., 2004; [612] and [13]. Previous Ebola virus modeling papers have examined its 1995 outbreak in Kikwit, DRC [1416] and [17,18]), and the in Uganda in 2000–2001 (see [16,19] and [17]).
Currently, however, the Ebola has also had certain implications in Guinea, Sierra Leone and Liberia. Indeed, a small number of cases are reported in Nigeria (21 deaths) and a single case in Senegal (travelers arriving from Guinea). This consists, actually, in is the contagious Ebola virus spread among countries. Until September 8, 2014, 4366 cases had been reported, including 2218 deaths. The Ebola virus is discovered to be the most dangerous epidemic, recently resulting in a remark able increasing in the number of patients and deaths it’s since the outbreak out of this virus. The first Ebola outbreak among human beings occurred in 1976 [20], with a spread contagion due mainly to travelling. An infected person travelled for instance from Guinea to Liberia may be carrying the infection, which would in turn infect some other people there (The institute of tropical medicine, 2013 [21]).
Based on the Ministry Health reported information (September 2014), Guinea recorded 1040 confirmed deaths highly exceeding the confirmed cases. Figures relevant to the Ebola cases and deaths have been provided by WHO, with the most prominent number being registered in Liberia, (more than 1760), while, 2450 probable deaths have been reported with respect to Sierra Leone. These data are based on official information reported by the Ministry of Health up the October 5th for both Guinea and Sierra Leone and October 4th for Liberia.
Concerning Guinea, one might well note that the number of confirmed cases is important in respect of the probable and suspected cases. Still, the Liberian probable cases remain the most significant in relation to the entirety of: confirmed and suspected cases. As for Sierra Lione confirmed cases, they constitute the highest proportion as compared to the other three remaining countries. Finally, Nigeria exhibits the weakest level of cases regarding the three indicators involving confirmed, probable and suspected cases.
Actually, five species of Ebola have been detected to prevail in Bundibugyo, Côte d’Ivoire, Reston^{b} [22], Sudan and Zaire. The outbreaks of Ebola virus associated of hemorrhagic fever has taken place mainly in Africa, with a death rate comprised 25% and 90%.
As a matter of fact, the Ebola virus is transmitted to human beings from wild animals, and then spreads from humantohuman. The 2014 Ebola epidemic is the most prominent in history, affecting aerial several countries. Since the beginning of 2014, some cases of Ebola virus infection have been reported in several subSaharan African countries. These countries, health authority’s line with the World Health Organization (WHO) and its different partners, have all been mobilized to prevent the virus transmission to other countries. Travelers to or from these countries can find relevant health care related information on the Ministry of Foreign Affairs and International Development website.
In effect, Ebola is a communicative disease, as it can be transmitted by either contact with blood, body fluids or tissues of infected subjects, or sick or dead animals. So for airborne transmission has been indicated to occur. The incubation period (the time between the infection and symptoms’ onset) ranges from 2 to 21 days. A person who bears no symptoms is not contagious. Most often, the risk degree of the Ebola virus contamination is considered too low during the disease the early stages but contagion of the Ebola virus makes the symptoms worsen or improve. People are infectious as long as their blood and secretions bear traits of the virus.
The outbreak of the Ebola virus tends to occur at irregular intervals in a medical care environment. The first trigger was reported in the Forest region of Guinea. The Ebola outbreaks in Liberia, Sierra Leone and Nigeria were contained [23]. Until recently, the Ebola virus has been raging mainly in central and eastern Africa. Contagion occurred between 1976 and 2012, in the Democratic Republic of Congo and Sudan along with the three other countries which had remarkably been affected by serious epidemics’, namely; Gabon, Uganda and Congo. Although Ebola has resulted in a high mortality rate, it’s the real varies between 25% and 90%, depending on the strain. According to the WHO (October 2014), Ebola currently stands within an average of 54% of total deaths [13].
Based on Figure 1, and for an effective combating of the epidemics some conditions and resources need be considered and ensured, namely case management, contact tracing and monitoring, laboratory quality, safe burials along with social mobilization. The community is involved in the outbreak control as earlystage awareness of the Ebola virus infection and the possible protective measures are likely to be successful. In fact, they have greatly helped in reducing contamination among humans. Messages about risk reduction should rather focus on social mobilization and communitybased awareness campaigns to help reduce the disease propagation risk [24].
Indeed, the WHO intends to reduce the Ebola virus propagation risk between animals and human beings, about all contact with monkeys. In fact animals should be handled with gloves as well as other appropriate protective clothes. More importantly, animal resource based foods must be thoroughly cooked before eating them.
Other various factors may well contribute greatly in increasing or decreasing contamination risk [25]. In this respect, the present study attempts to provide certain suggestions whereby contamination risk to human beings can be reduced, a risk occurring mainly through direct contact with similar subjects, and, predominately with their body fluids. Normally, people need to wear gloves and appropriate personal protective equipment when taking care of an affected patient at home. To prevent such a contagion, people also have to wash out hands after visiting patients in hospitals. The various measures necessary for containing the Ebola outbreak include corpses ‘rapid burial and identifying subjects likely to be in contact with an infected person.
Concerning fight against infection in health care facilities, health workers should always apply the standard precautions on taking care of patients, regardless of the presumed diagnosis. These precautions pertain to the basic rules regarding hygiene including, wearing personal protective equipment (to protect oneself against splashing or direct contact with infected materials), safety injections as well as funeral rites.
Health workers responsible for taking care of the Ebola virus infection suspected or confirmed cases must avoid any contact with the patient’s blood or body fluids and with such contaminated surfaces as clothing and bed linen materials. On standing too close too sick person (say, i.e. within a distance of less than a meter), one has to should wear face protection (face shield or surgical mask and goggles), a clean sterile gown and gloves (sterile medical procedures). In him laboratory staff might well face the same risk. The Ebola samples taken from human beings or animals should be handled by specialized trained staff and processed in suitable equipped laboratories.
The WHO aims to prevent Ebola outbreaks by providing disease control and helping countries atrisk to develop special safe guarding measures and preparing plans. In this context, this present paper is dedicated to provide some general guidance regarding the fight against outbreaks caused by these two viruses. Both of the Ebola disease and Marburg virus outbreak and curbing stages are: preparation, warning, control and evaluation. Whenever an outbreak is detected, the WHO is involved in providing assistance through monitoring, supporting patients, mobilizing laboratory services, research contacts, infection control, logistic support training and assistance with safe burial practices. In This regards, the WHO has developed detailed preventive information to curb prevent and fight against with Ebola virus contamination.
Since December 2013, an Ebola virus originating hemorrhagic fever epidemic has predominantly prevailed in West Africa (Guinea, Liberia, Sierra Leone, and Nigeria are the major affected countries in addition to a confirmed case in Senegal). For the first time, such in West African countries as Congo, DRC and Uganda have been affected by with the Ebola virus. On the interest, international SOS organization is monitoring the situation closely and updates the site with the latest available information. In this respect, the present work constitutes a case study concerning of the Ebola virus transmission among the most affected countries by the virus (Guinea, Sierra Leone, Liberia and Nigeria).
This paper is designed to check whether these countries can be actually constitute a source of recent contagion judging by the number of Ebola virus incurred deaths. More particularly, this work is focused on showing whether the Ebola virus can produce contagious effects within the studied countries’ in conditional volatilities over the period ranging between March and October 2014. Previous contagion related studies have somewhat failed to take account for the important distinction between the concepts of interdependence and contagion^{c}. Indeed, the conditional test results’ dynamic correlation reveal well that the contagion effects appear to be multidirectional since the Ebola virus emanating return shocks can actually sweeping across the entirely of the markets. Yet, the contagion involatility effects are mostly driven by the Ebola virus originating negative return shocks. Such an empirical finding reveals well that the number of death shocks can turn out to be contagious not only at the volatility level, but also at the mean level, indicating that Ebola outbreak can be a major source of contagion over the period MarchOctober 2014.
In fact, the objective of the paper lies in highlighting that deaths’ number does actually constitute a contagion channel. In a first step, this study is concerned to test the persistence of shocks and the stylized facts of these countries’ considered Deaths returns number through EGARCH process [4]. In effect, the presence of structural changes in the series of conditional variances as detected via the ICSS algorithm heave motivated us to study Conditional Correlation Dynamics [26] in a second stage as a procedure to determine the contagion effect across markets. Finally, this study serves to help predict the number of Deaths’ volatility across the EGARCH process.
In this study, the cited process involves four major steps. On a first stage, the Iterated Cumulative Sums of Squares algorithm (ICSS) of [27] is applied to detect the presence of structural breaks in deaths number relevant to the studied countries. On a second stage, and structural breaks and asymmetry to be taken into estimation, the study undertakes to analyze the univariate EGARCH model and bring the structural breaks’ dummy variables into variance equation. Compared to the standard specification, the EGARCH framework exhibits several advantages namely: no need is imposed to artificially implement a nonnegative constraint on the model parameters, while asymmetries are allowed under the EGARCH formulation. Dubbed dynamic conditional correlations such a methodology has actually been developed by [26]. In reality, the major advantage lying implementing the DCCMGARCH models consist in detecting of plausible changes noticeable in the relationships governing the variables remains prevalent in the underlying studied data.
The remainder of this paper is structured as follow. Section Methods is devoted to highflying the surveyed econometric methodology. The relevant data and major empirical results are discussed in Section Empirical results. As for the ultimate section it comprises provides the pertinent concluding remarks, along with the retained economic and political implications.
Methods
For the purpose of the determining whether the number of deaths’ volatility dynamics does actually differ from that of the countries, we, firstly, undertakes to expose the exact local Whittle and structural break tests, and, then, discuss the GARCH specifications so as to be able to capture the possible conditional dynamic dependencies likely to be noted in the number of death return volatilities.
The Exact Local Whittle method
The classes of semiparametric frequency domain estimators follow the local Whittle approach as suggested by [28] and analyzed by [29] (dubbing it Gaussian semiparametric estimator). The analysis applied process is the following:
As for the Local Whittle estimator, it defined as the maximization of the local Whittle likelihood purpose, such as:
where: m = m(T) denotes a bandwidth number tending to infinity T → ∞ except at a slower speed than T;\( I\left(\lambda \right)=\frac{1}{2\pi T}{\left{\displaystyle \sum_{t=1}^T{e}^{it\lambda }}\right}^2 \), represents the periodogram of X _{ t }, and g_{ x }(λ) the spectral density of X _{ t }, \( {\lambda}_j=\frac{2\pi j}{n} \), and j = 1,…,n.
A notable disadvantage as compared to logperiodogram estimation is that a statistical optimization is highly. Still, this estimator underlying assumptions are weaker than those pertaining to the logperiodogram regression (LPR) estimator. In this regard, [29] have show that while \( d\in \left(\frac{1}{2},\frac{1}{2}\right) \);
Hence, the asymptotic distribution turns about to be extremely simple, which facilitates easy asymptotic inference. More particularly, this estimator is discovered to be more efficient than the LPR one. The reliability and asymptotic normality ranges concerning the Local Whittle estimator have explicitly been demonstrated by [30] and [31] to equate those associated the LPR estimator.
This exact LW procedure as frequency labelled, implies replacing \( {\lambda}_j^{2d}I\left({\lambda}_j\right) \) in (2.1) by \( {I}_{\Delta^{d_y}}\left({\lambda}_j\right) \), and is only valid if μ = 0 in (2.2). Since the relevant means are different from zero, [32] suggests demeaning {y _{ t }} with an appropriate estimator \( \widehat{\mu} \), and computing the exact LW estimator starting from the demeaned data. So the objective function to be minimized turns out to be:
Where: \( {I}_{\Delta^d\left(y\widehat{\mu}\right)}\left({\lambda}_j\right) \) is the periodogram of \( {\Delta}^d\left(y\widehat{\mu}\right) \). For fractional differences, to be determined, it is assumed that {y _{ t }} is given by a process similar to equation (2.1). It turns out that the first sample observation y _{1} is a reliable mean estimator in the case of large values of d, while the usual arithmetic mean ӯ helps ensure a significant task for small coefficient values of d. In this way, [32] suggests putting forward the subsequent weighted estimator, such as:
For the purpose of attaining, a feasible procedure, he considers two necessary steps, the first of which serves to determine an estimator of \( \widehat{d} \) independent from μ in order to get an estimator of the constant: \( \widehat{\mu}=\widehat{\mu}\left(\widehat{d}\right) \). As for the second step, the slope and Hessian of R _{ E }(m,d) are used to compute the feasible estimator as follows:
Besides, [32] demonstrates shows that the twostep ELW estimator (2ELW) proves to be consistent registering the same limiting distribution as the LW and ELW estimators under −0.5 < d < 2. Similarly, and as indicated as shown by [33], if an unknown mean (initial value) appears to undergo certain change by its sample average, simulations suggest that the ELW estimator is inconsistent for d > 1. It is actually for this reason that we undertake to apply the 2ELW. In addition, [33] resort to modify the ELW objective function in a bid to estimate the mean by means of combining two estimators: the sample average and the first observation. He indicates the resulting estimator as being a two Stage Exact Local Whittle (2ELW). Applying the tapered estimator of [30] in its first stage, the 2ELW estimator, bears the same \( N\left(0,\frac{1}{4}\right) \) limit distribution for \( N\left(\frac{1}{2},2\right) \) and is consistent when \( d>\frac{1}{2} \). Furthermore, the 2ELW estimator finite sample performance appears to inherit the 2ELW estimator, desirable properties. Moreover, it can also be computed with prior data detrending (2ELWd) as in [33].
Structural breakpoints detection
As developed by [27] the ICSS algorithm has been applied to detect the structural breakpoints on 4 series over the study period. As a starting point, the stock return for market i on day t can be written as:
where: (P _{ i,t }) is the closing number of deaths:
Next, we proceed by defining
where {a _{ i,t }} is with zero mean and unconditional variance \( {\sigma}_t^2 \), μ_{ i } denotes country i average return. Let \( {C}_k={\displaystyle \sum_{t=1}^k{a}_t^2},k=1,\dots, T \) be the cumulative sums of squares of {a _{ t }} series, then D _{ k } statistic can be calculated as follows:
The ICSS algorithm is adopted to detect the multiple breaks in the unconditional variance of {a _{ i,t }} series. Thus, the D _{ k } statistics based ICSS algorithm would initiate with testing the structural breaks over the whole sample. In doing so, the ICSS helps depict any significant break, by applying the new statistic to examine the break for each of the two subsamples (defined by the break). The algorithm proceeds in such a way as the statistics become insignificant for the entirety of subsamples defined by any significant break. On the last stage, a dummy variables set is set up for the normalized return volatility to be captured.
This section is devoted to provide a through description of the wavelet transform as applied to the number of the death data decompositions, together with the multivariate GARCH model as used in our proper analysis.
As countries comovement is but an outcome of transmissions emitted from each country, the global country transmission is most often represented by the number of deaths occurring in transferred to other countries. For the sake of accounting for such interdependencies, we consider appealing to an MGARCH model including the EGARCH structural changes associated with the variances introduced by [4]. In what follows is a presentation of the dynamic model EGARCH (1, 1).
The EGARCH model
The above GARCH specification helps restrict any shock effect on the conditional variance to be symmetric and within the same size, be it positive or negative. Noteworthy however, such a shock effect should disappear geometrically over time. Still it wellbecome a known fact notably with the help of death data, the different country reactions in accordance with to the shock size and sign. To overcome such limitations, several asymmetrical GARCH models, such as the Exponential GARCH model [34], have been introduced in this respect. Formally, an EGARCH (1, 1) model corresponds to the following:
With an EGARCH specification, positivity constraints on parameters are no longer required. Furthermore, it is henceforth possible to apprehend asymmetry in volatility reaction toward external shocks. Indeed, if γ_{0} > 0 (respectively γ_{0} < 0), a positive shock on the lagged conditional variance implies an increase (respectively a decrease) in current volatility.
Parameter γ_{1} helps capture the asymmetry effect associated with the shock size \( \frac{\left{\varepsilon}_{t1}\right}{h_{t1}} \). If γ_{1} = 0, a positive innovation would them have the same effect on conditional variance as would a negative innovation, while for γ_{1} > 0, a high size shock should have more effect on conditional variance than should a low size shock. So for the Ebola virus contagion to be determined, the MGARCHDCC model will be applied so that the virus propagation effect can be determined concerning four African countries (Guinea, Sierra Leone, Liberia and Nigeria).
The multivariate GARCH DCC model
In this section, the dynamic conditional correlations’ twostage model will be treated as proposed by [26]. Let us for instance, consider a vector consisting of any two variables Y _{ t } = [y _{1t } y _{2t }]’. Each variable constitutes a constant function along with its own past values. Thus, the autoregressive process reduced form is written as:
where: A(L) is the polynomial delay, and ε_{ t } = [ε_{1t }ε_{2t }]’is a vector of residuals training from the estimation auto regression process for each variable whose variancecovariance matrix is described by H _{ t } = {h _{ i }}_{ t } with s.
Actually, the DCCMGARCH model can be easily apprehended by rewriting the matrix of variancecovariance H _{ t } such as:
where: \( {D}_t= diag\left\{\sqrt{h_{it}}\right\} \) is the standard deviations diagonal matrix, a variable temporally different from the two previous equations’ estimation with a univariate GARCH process; R _{ t } = {ρ_{ ij,t }} representing the conditional correlation coefficients matrix. The D _{ t } contained elements are generated into a GARCH (P,Q) process, which can be formulated as:
In addition, [26] considers adopting a GARCHtype structure while modelling the correlations’ dynamics. Thus, an (M, N) order DCC process can be described by:
where: \( {\xi}_t=\left\{{\varepsilon}_{it}/\sqrt{h_{it}}\right\} \) is the vector englobing the standardized residuals derived from the univariate GARCH model estimation, as a matrix of these standardized residuals’ conditional variancecovariance, whereas Q _{ t } = {q _{ ij,t }} represents the unconditional variancecovariance matrix, which are temporally invariant. The parameters (a _{ m };b _{ n }) are supposed to respectively intercept the shock effects and delay the dynamic correlations at the level of recent contemporary. As for Q ^{*} _{ t }, it stands for a diagonal matrix containing the square root of the main diagonal elements of Q _{ t }. Regarding to our example, this matrix is written as:
with: \( {\rho_{12,}}_t=\frac{q_{{{}_{12,}}_t}}{\sqrt{q_{{{}_{11,}}_t}{q}_{{{}_{22,}}_t}}} \) denoting the dynamic conditional correlations i.e. the matrix elements R _{ t } whose main diagonal consists of 1.
The model parameters are estimated via the maximum likelihood DCC method. In this regard, [35] have show that the loglikelihood function can be expressed as:
The estimation process involves two steps, the first of which consists in substituting an identity matrix by a matrix R _{ t } within the loglikelihood the function. The advantage of such a method is that it allows for getting the likelihood function sum of the GARCH univariate models. In other words, this initial step, series to achieve we obtain the parameters of equation (2.13). As for the second step, it consists in estimating the equation (2.14) parameters via adopting the original likelihood function as described by equation (2.16). This procedure ensures maintaining the dynamic correlations among the studied variables.
Results
Data and descriptive statistics
The empirical data comprise total number of daily^{d} deaths caused by the Ebola virus as calculated by “OMS” of four African countries. Actually, the number of deaths recorded exclusively in Guinea, Sierra Leone, Liberia and Nigeria has been selected as subject of study as these four countries are the most affected with the Ebola virus. The sample ranges over the period comprised between March 2014 and the end of October 2014^{e}, yielding 206 observations for each series.
Table 1 presents a wide range of descriptive statistics concerning the four series under investigation during the period (MarchOctober 2014). The null hypothesis of no ARCH effects is rejected at a significance level of 1%, suggesting that the GARCH parameterization might well fit for the conditional variance processes.
In Table 1, the sample size and the unconditional mean indicate that the number of average daily deaths is fixed at 3.9 for Guinea, 5.7 for Sierra Leone, 11.9 for Liberia and 0.03 for Nigeria. With respect to Liberia and Sierra Leone, the sample provides a higher average return with a higher unconditional variance, highlighting the prevalence of a strong dispersion around the Ebola virus incurred deaths in these two countries. It is actually, this high dispersion feature which provides us with significant breaks through the conditional variance. Furthermore, distribution of the number of the virus related deaths seems to be symmetrical and leptokurtic, implying rejection of the normality null hypothesis. The mean and variance will, conditionally, be modelled. Therefore, the most interesting statistics, such as the skewness and the excess kurtosis, can be used for the purpose of testing whether the empirical distribution does have kurtosis and skewness just like a normal distribution. This undertaking has already been applied done with the JarqueBera test which rejected the null hypotheses, indicating that the sample distribution stems actually from a normal distribution set to a 1% significance level. To note, LB (24) statistics are employed for the take of testing the persistence of a high autocorrelation in the results’ first and second moments. The LB (24) statistics figured on Table 1 display just the existence of a highorder autocorrelation for Guinea, Sierra Leone, Liberia and Nigeria. These LB (24) statistics are supposed to be significant, reflecting the existence of a noticeable interdependence among second moments of returns. This also indicates that the returns’ heteroscedasticity should change with time. Such as result highlights well the use of estimation and variance of the autoregressive conditional heteroscedasticity (ARCH) model of [35]. As for the ARCH test, the null hypothesis of homoscedasticity is not accepted suggesting a further conditional timevarying in the number of death return dynamics.
Long memory dependency and the Ebola virus
In a first place, the Robinson estimators of long memory parameters relevant to Table 2 reported periods, prove to be lower than 0.5 for the entirely series. Such a result highlights well that long memory dependency within the period turns out to be critically important. This might well be due to the recurrent shocks and breaks recently occurring in the African countries. A possible explanation for this lies in the fact that the invasive occurrences in the number of Ebola virus related deaths respective to Guinea, Sierra Leone, Liberia and Nigeria, during the studied period corresponding have lasted for extended periods and increased the long memory property as a mean process which responded asymmetrically and gradually to such shocks and breaks as already pointed out by [36].
The plot of Figure 2 depicts the Sample Partial Autocorrelation Function clearly highlighting the prevalence of a significant autocorrelation. The long memory process captures a very low frequency cycle in the number of Ebola virus incurred deaths, by permitting a slowly declining autocorrelation of a hyperbolic shape in the long horizon (see Figures 1 and 3 relevant to Guinea and Sierra Leone). Inversely, however, the short memory process is characterized by a fast exponential decline in autocorrelations as depicted by Figures 2 and 4 pertaining to Liberia and Nigeria. For the purpose of verifying the persistence of the Ebola virus, pertinent real shocks, we consider applying the LW, 2ELW and 2 ELW with trending.
Based on Table 2, the Robinson estimators (GS) do actually prove that \( 0\le \widehat{d}\le 0.5 \). Indeed, this consists in a longmemory process case through still stationary, with a slow or smooth decay in the catchingup process. Concerning the ARFIMA(1,d,1) process, it has been demonstrated that \( 0\le \widehat{d}\le 0.5 \) with respect to Guinea and Sierra Leone, underling a long memory and stationary process. In regard of both Liberian and Nigerian, cases, it has been prove that \( \widehat{d}\ge 1 \). In effect, this corresponds to an explosive to process case and a situation marked by where there is a strong magnification effect, with any initial difference being unexpected to be potentially reversed in the future.
In reality, this coincides with the “stochastic divergence” case liable to comparison with the initial deterministic divergence case. Regarding the 2ELW estimator, it has been demonstrated that \( 0.5\le \widehat{d}\le 1 \), corresponding to a long memory process case, which is nonstationary through still reverting. In such a case, the process is featured with high persistence, whereby any distant past output difference would still have a longlasting present inference. With respect to the 2ELWd estimator, it has also been demonstrated that \( 0.5\le \widehat{d}\le 1 \) regarding the entirety of studied cases, except for Sierra Leone. What noting also, the number of deaths highest values (0.272 and 0.237, respectively for Guinea and Liberia), highlight well the persistent of shocks in the number of Ebola virus increased deaths and its contagion.
In fact, the Ebola virus can persist in the patient’s body for about 2 months and, the sick and dead people are contagious. Following the incubation period ranging between 2 and 21 days (a period with no symptoms and no contagious), the first symptoms^{f} begin to appear, and the affected person becomes contaminated on with the onset of the symptoms. The Ebola virus raised death occurs relatively quickly due to bleeding and unfortunately, there is no treatment or vaccine against the disease is yet available. This might well sound being plausible when the virus affects at fast countries within a small restricted geographic area, but spreads out to cover a wider geographical area and persists underway for several months (between 3 and 4 months). In addition, if the Ebola dynamics were to be modeled in real time, it would be critical to consider possible delays while report some cases and outcomes.
The Ebola virus transmission is marked by two major phases. The first started in July 2014, when Guinea, Liberia and Nigeria become affected with the shocks. As for, the second phase, it started in August 2014, when all the four countries become affected by at least a single shock during that month. Concerning Guinea and Nigeria, there were structural changes for two months, September (05/09, 24/09 and 26/09) and August (09/08, 16/08 and 27/08). This fact indicates an approximation of the Ebola virus quiet periods boundaries and times. This denotes that the number of fatal Ebola viruses appears to be is more dispersed during these two months owing to the increasing number of Ebola virus incurred deaths. As a matter of fact, 494 deaths were recorded in August reaching 739 deaths in late September in Guinea.
In Liberia, one might well note that the four points marking structural changes through the ICSS algorithm are sited during August and early September. During the end of the period, the number of deaths increased remarkably from 260 to 953. In regard of Nigeria, only two structural change points have been highlighted occurring between late August and early September. This might well has its explanation in the increase noticed in the number of deaths growing up from 1 to eight due to the Ebola virus. Unless it is imperative, the studied countries are recommended to suspend all their flights towards the countries where the Ebola haemorrhagic fever cases have proven to be remarkably serious, such as Guinea, Sierra Leone, Liberia and Nigeria. Indeed, for five months, the Ebola outbreak was the worst since the emergence of the initially high contagious hemorrhagic fever in 1976. Actually, it was at the origin of 1,145 reported deaths according to the World Health Organization (WHO), reporting 413 cases in Liberia, 380 in Guinea, Sierra Leone and 348 in Nigeria.
Such interesting findings appear to be unable to capture the asymmetry of volatility reaction to external shocks. For this reason, an EGARCH (1,1) model is estimated whose empirical results are reported in Table 3. Accordingly, the parameter γ_{0} is shown to be negative and significant except for Liberia and Nigeria. This may well further evidence in favor of variance reduction after a positive shock and vice versa.
Noteworthy, however, the estimated value of γ_{0} appears to differ greatly between (Guinea and Sierra Leone) on the one hand, and (Liberia and Nigeria) on the other hand. This suggests that the volatile countries’ exposure and liability to risks varies not only among (Guinea and Sierra Leone) and (Liberia and Nigeria), but also among the countries under review. Obviously, this might have an association with the number of deaths which seemed to be lower with respect to (Liberia and Nigeria) than to (Guinea and Sierra Leone), mainly on the onset of the Ebola virus disease. As for parameter δ_{1} which serves to measure the conditional volatility persistence, its estimators prove to be significant. The highest value related to this parameter has been recorded for Guinea and Liberia exceeding 0.9, in respect of its value registered in Sierra Leone and Nigeria. It also reveals volatility dynamics which appear to be more persistent with regard to Guinea. Besides, our results appear also to confirm the persistence of asymmetry in the volatility reaction to shock size. Indeed, the estimated value of γ_{1} is positive for all the countries except for Liberia, scoring its highest points in (Guinea and Sierra Leone), suggesting that the volatility dependency to market events is more important for both Guinea and Sierra Leone.
Regarding the number of deaths registered in Guinea, Sierra Leone Liberia and Nigeria, 260 daily data were pointed during the entire period from: 20140322 until 20141230. The insample period: from 20140322 until 20141014 consisted of 206 daily data points while the outofsample period from: 20110103 until 20140415 consisted of 54 daily data points. In Figure 2, the number of daily deaths recorded during the entire period is plotted. The insample period is plotted in black, and the green graphics indicate the lower graphics while and the upper graphics are highlights with the blue line.
A decrease in the number of Ebola virus incurred deaths was noticed during the following two months (November and December 2014). This is mainly due to a drastic reduction in the number of animals natural reservoirs of the virus. Indeed, the upper and lower figures appear to have a similar shape to the graph in blue. This indicates well a significant decrease in the number of the Ebola virus caused fatalities in Guinea. The virus transmission remains intense in Guinea with a clear exponential in the daily deaths number time series. In Sierra Leone, for instance, one may notice that the number of daily Ebola deaths remained constant during the following two months (November and December). This number marked some peaks from one day to another. As for Liberia, the number of the Ebola viruses is expected to increase over the next 60 days (November and December). Liberia contributed strongly in the high number of the whole outbreak cases registered during the last two months (November and December). Regarding, Nigeria, the number of forecast deaths registered via EGARCH process highlights this figure constancy over the coming months (November and December). For the sake of determining the causal relationship between Ebola and the studied countries, the Granger causality test is applied as introduced by [37].
The Granger causality test
This part is interested to test the presence of Granger causality relationship in a bivariate model among different series (Guinea, Sierra Leone Liberia and Nigeria) on the studied countries. To note, the Granger causality test is a statistical hypothesis test helpful in determining whether a certain time series is useful for forecasting another (see [7]). Ordinarily, regressions help reflect correlations, but Granger argued that causality in economics could well be reflected by measuring the ability of predicting a certain time series’ future values using another series’ past value. Since the “true causality” issue is deeply philosophical, econometricians assert that the Granger test finds only “predictive causality” [38].
In this section, the Granger causality tests are going to be administered, using. We use the unrestricted model through which causality will be investigated between the number of deaths respective to Guinea, Sierra Leone, Liberia and Nigeria. The relevant results are summarized in Table 4, below. The number occurring in each cell of the probability columns indicates the presence of a significant causal relationship among the four African countries under consideration. A superscript (*, **, ***) reflects the Ftest statistical significance for Granger causality at 10%, 5% and 1% significance levels respectively. The results indicate three patterns:

Group 1: it is composed of four countries; Guinea^{g} is the origin source of Ebola virus, i.e. the African countries are recently affected with the Ebola virus through Guinea. This group includes Guinea, Sierra Leone, Liberia and Nigeria.

Group2: is composed of such countries as Guinea, Liberia and Nigeria which are affected by Sierra Leone^{h} where a high number of deaths predominate.

Group 3: consists of such countries high as Guinea, Sierra Leone and Nigeria affected by the Ebola virus as transmitted from Liberia^{i}.

Group 4: includes Guinea, Sierra Leone and Liberia in which the virus is transmuted from Nigeria^{j}.
The number of causal relations hips is depicted in the following table. A causal relationship is persistent in each of the first three groups, but no relationship has been noticed to exist in the fourth group.
According to Table 4 pertaining to the Granger causality test, the disease appears to be transmitted from Guinea to Liberia (Fisher 4049 is calculated at a 5% significance) threshold. This implies that there is a high probability for the disease to have been spread from Guinea to Liberia.
Moreover, it is clear that the disease in Sierra Leone has resulted in the same disease in Liberia (Fisher calculated is 2.984 with a probability of 0.08). Accordingly, one might well infer that the Ebola virus may have been spread from Guinea through Sierra Leone to Liberia. It has also been shown that the Ebola virus has been spread from Liberia to Sierra Leone. These two findings several that the Ebola virus native countries are discovered to be Guinea and Sierra Leone. For the purpose checking the right transmission of the Ebola virus among these countries, we have made appeal to the DCCMGARCH technique. In fact, this framework enables us to estimate two main parameters, the first of which helps determine the persistence of the Ebola virus, while the second serves to demonstrate the Ebola virus shock contamination effect through the considered countries.
The two major main parameters considered are a _{ m } and b _{ n }, respectively, underlining the persistence of shock effect and contamination. Based on Table 5, one can notice that coefficient a _{ m } proves to be high and significant (0.583) highlighting the continuity of the Ebola virus impact. Consequently, a high probability remains as to the virus transmission from Guinea to Nigeria. As matter of fact, coefficient a _{ m } has been demonstrated and deemed important because the study period has covered only 8 months (MarchOctober 2014). During this fairly short period, several structural changes have been noticed to take place (5 in Guinea; 4 in Sierra Leone and Liberia and 2 in Nigeria) as figured in Table 6. This implies well that these dates help indicate structural changes of the peak rather than the Ebola virus shock persistence of in the surveyed countries. Regarding the second parameter, b _{ n }, the Ebola virus spread transmission is observed to occur between GuineaLiberia, GuineaSierra Leone and Sierra LeoneLiberia. Yet, parameter b _{ n } appears to be not significant owing mainly to of the lack of data (only 206 observations) and the GARCH family. Indeed, the latter entails the availability of a large amount of data, parameters’ cruciality. For this reason, we have reckoned it more useful to apply the DCCMGARCH rather than the BEKKMGARCH which admits further parameters than does the DCC version. As results, we consider it’s rather effective to study the Ebola virus contagion via DCC figures.
Figure 5 depicts the Dynamic Correlation Coefficients respective to four country pairs. The correlation variation among each pair can be well observed via DCC coefficients. Thus once positive and close to 1, correlation is discovered to indicate a similar trend in the number of deaths. However correlation proves to be negative with an absolute value close to 1, it would then indicate an opposite trend in the number of deaths. Actually, DCC has been discovered to be comprised between 0.4 and 0.6 during the first four months (MarchJune) for both pairs GuineaNigeria as well as GuineaLiberia, indicating the Ebola virus spread between Guinea and Nigeria. In regard of the other two pairs, contagion seems less prevalent in both GuineaSierra Leone and Sierra LeoneLiberia and has been proven to present exclusively during the last month (with DCC comprised between 0.3 and 0.5). So, it can be concluded that Guinea appears to be the source origin of contagion, as one may notice the Ebola virus shock spread from Guinea to Sierra Leone to Liberia. Indeed, on March 31, the WHO did confirm the virus spread to Liberia. On April the 17th, the number of deaths ranged between 131 and 209 in both countries.
On August 26, the African Development Bank estimated that Ebola outbreak could cost between 1 and 1.5 percentage points of the GDP in Liberia, Sierra Leone and Guinea, starting it a “country which began to recover difficult years of crisis, civil wars of the 60s, 80s and 90s”^{k}.
Discussion
On September 17th, a report released from the World Bank warned against a “catastrophic” economic impact the epidemic may have on Liberia, Sierra Leone and Guinea; due to the virus associated “fear factor” paralyze activities. According to the World Bank orientations, the combined GDPs of Liberia, Guinea and Sierra Leone could be cut by $ 359 million in 2014 and $ 809 million by 2015 if the epidemic were to remain uncontained. As a result, economic growth would then fall, in the following year, from 11.7 points to 8.9 points in Liberia and Sierra Leone, with a risk of plunging the two poorest countries in recession^{l}.
Criticized in an internal report for the delay to take the necessary measure relevant to curb Ebola outbreak, the World Health Organization announced that the two antiEbola vaccines would be tested as soon as possible in Switzerland, subject to authorization of the medical authorities. French researchers beginning to the Commission for Atomic Energy and Alternative Energies (CEA) have also announced that a rapid test has been developed to diagnose the virus within fifteen minutes. Since its introduction in early March, the virus has been responsible for killing more than 4,500 people, mainly in Liberia, Sierra Leone and Guinea, according to the latest WHO released report (October 2014).
Conclusions
Throughout, the scone of this study paper attempts have been made to test whether the Ebola virus can stand as an origin source of contagion during 2014 transmission process using the number of deaths data available from the World Health Organization. More particularly, we have undertaken to examine whether the Ebola virus can lead to contagious effects implicated on both conditional means and volatilities of its relevant number of deaths occurring during the recent period (MarchOctober 2014). In fact previous studies have failed to take account of the important distinction both of the interdependence and contagion concepts. In this paper, contagion is defined as being the significant spillover of the assetspecific idiosyncratic virus during the crisis and as well following the first affection by the Ebola virus [23]. For an effective control of the economic fundamentals, we have considered relying on an international number of death models, which have provided a theoretical basic background in selecting the pertinent health fundamentals.
The empirical findings have revealed that contagioninmean effects turn out to be multidirectional given the fact that the Ebola virus chocks emanating from any of the four studied countries can sweep across throughout the entirety of investigated countries; however, contagioninvolatility effects appear to be driven mainly by the Ebola epidemic stemming negative return shocks. This empirical result indicates well that shocks to countries’ return can eventually become contagious not only at the volatility level, but also at the mean level, significantly implying that the number of deaths communication can actually constitute a major cause of contagion throughout the recent period.
Endnotes
^{a}The virus family of which Ebola is a member, “Filoviruses” are far more ancient than previously thought.
^{b}EbolaReston virus was recognized in monkey export ability in the Philippines. No human virus was identified.
^{c}Specifically, contagion is defined in this paper as being significant spillovers of assetspecific idiosyncratic shocks during the crisis, after the Ebola virus or systematic risks have been accounted for [39].
^{d}Weekly data are used here to get meaningful statistical generalizations and obtain a better picture of the movements of Ebola virus cases’ number.
^{e}As a first step, stationarity in the time series is checked by applying the Augmented Dickey Fuller (ADF) test. The results allow us to reject the null hypothesis stipulating that the Ebola virus returns have a unit root in favour of the alternative hypothesis (even at 5% critical value).
^{f}Fever, weakness, muscle and joint pain, diarrhea, internal and external bleeding.
^{g}In Guinea, there are 1298 cases, of which 1062 are confirmed, with 768 deaths. The epidemic epicenter is Guéckédou, Guinea, but other affected provinces include also’: Beyla Conakry, Coyah, Dabola Dalaba, Forecariah Macenta Nzérékoré Kindia Kérouané, Kissidougou Dubréka and Lola Yomou. The conakry isolation center was set up by MSF. It seems that appearence of new cases in the Guéckédfou region is decreasing (World Health Organization).
^{h}In Sierra Leone, there are 2948 cases with 880 deaths, among 2596 confirmed cases. All regions are discovered to be affected. The isolation and diagnosis center is sited in Kenema. The MSF has also implemented an isolation Kailahun and Lakka structures (Lakka Hospital) on the outskirts of Freetown. The Choitrams Hospital can no longer accept hospitalization.
^{i}In Liberia, there are 3924 cases in total, 941 of which were confirmed, in 1795 with 1188 and 2210 suspicious deaths. It is likely that these figures are far below reality. The Diagnostic Center is located in the Liberian Institute of Biomedical Research (LIBR) near Monrovia. The hospitals fully devoted to the Ebola case are: ELWA, ELWA3, JFK in Monrovia, are hospitals which cannot cope with the influx of new patients.
^{j}In Nigeria, the last new case dates back to September 5 in Lagos and Port Harcourt September 1 (15 in Lagos and 4 in Port Harcourt). The epidemic appears under control in Nigeria but there are still 25 subject contacts under observation. The US CDC reduced its recommendation travel Level 1 (can proceed with caution).
^{k}Ebola could cost “a lot” to the African economy, Le Monde, August 26, 2014.
^{l}Ebola: an economic impact “catastrophic” envisaged by the World Bank, Le Monde, 17 October 2014.
References
 1.
Krünsch HR. Statistical aspects of selfsimilar processes. In: Proceedings of the First World Congress of the Bernoulli Society. Utrecht: VNU Science Press; 1987. p. 67–74.
 2.
World Health Organization. Ebola haemorrhagic fever in Zaire, 1976[958 KB, 24 pages]. Report of an International Convention. Bull World Health Organ. 1978;56(2):271–93.
 3.
Legrand J, Grais RF, Boelle PY, Valleron AJ, Flahault A. Understanding the dynamics of Ebola epidemics. Epidemiol Infect. 2007;135:610–21.
 4.
Nelson DB. Conditional Heteroscedasticity in asset returns: a new approach”. Econometrica. 1991;2:374–0.
 5.
Ndanguza D, Tchuenche J, Haario H. Statistical data analysis of the 1995Ebola outbreak in the democratic republic of Congo. Afr Mat. 2013;24:55–68.
 6.
Formenty P, Libama F, Epelboin A, Allarangar Y, Leroy E, Moudzeo H, et al. L’épidémie de fievrehémorragique avirus Ebola en Republique du Congo, 2003: une nouvellestratégie. Med Trop. 2003;63:291–5.
 7.
Granger CWJ. Investigating Causal Relations by Econometric Models and Crossspectral Methods. Econometrica. 1969;37(3):424–38.
 8.
Incla’n C, Tiao GC. Use of cumulative sums of squares for retrospective detection of changes of variance. J Am Stat Assoc. 1994;89:913–23.
 9.
Georges AJ, Leroy EM, Renaut AA, Benissan CT, Nabias RJ, Ngoc MT, et al. Ebola hemorrhagic feveroutbreaks in Gabon, 1994–1997:epidemiologic and health control issues. J Infect Dis. 1999;179:S65–75.
 10.
Heymann D, Weisfeld J, Webb P, Johnson K, Cairns T, Berquist H. Ebola hemorrhagic fever: Tandala, Zaire, 1977–1978. J Infect Dis. 1980;142:372–6.
 11.
Breman J, Piot P, Johnson K, White M, Mbuyi M, Sureau P, et al. The epidemiology of Ebolahemorrhagic fever in Zaire 1976. In: Pattyn SR, editor. Ebola Virus HaemorrhagicFever. The Netherlands: Elsevier, Amsterdam; 1978. p. 85–97.
 12.
Nkoghe D, Kone ML, Yada A, Leroy EA. Limited outbreak of Ebola haemorrhagic fever in Etoumbi Republic of Congo 2005. Trans R Soc Trop Med Hyg. 2011;105:466–72.
 13.
Pattyn S. Ebola Virus Haemorrhagic Fever. Elsevier, Amsterdam: The Netherlands; 1978.
 14.
Lekone PE, Finkenstädt BF. Statistical inference in a stochastic epidemic SEIR model with control intervention: Ebola as a case study. Biometrics. 2006;62:1170–7.
 15.
Velasco C. Gaussian semiparametric estimation of nonstationary time series. J Time Ser Anal. 1999;20:87–127.
 16.
Lee YN, Gao Y, Wang HY. Differential mediation of the Wnt canonical pathway by mammalian Dishevelleds1,2,and3. CellSignalb20, 443452. 2008 CrossRefWeb of Science.
 17.
Chowell G, Hengartner NW, CastilloChavez C, Fenimore PW, Hyman JM. The basic reproductive number of Ebola and the effects of public healthmeasures: the cases of Congo and Uganda. J Theor Biol. 2004;229:119–26.
 18.
Miranda ME, Ksiazek TG, Retuya TJ. Epidemiology of Ebola (subtype Reston) virus in the Philippines, 1996. J Infect Dis. 1999;179 suppl 1:S115–9.
 19.
Ferrari MJ, Bjørnstad ON, Dobson AP. Estimation and inference of R0of aninfectious pathogen by a removal method. Math Biosci. 2005;198:14–26.
 20.
White LF, Pagano M. A likelihoodbased method for realtime estimationof the serial interval and reproductive number of an epidemic. Stat Med. 2008;27:2999–3016.
 21.
Institute of topical medicine in 2013. Remco van de Pasa and Sara van Belle, Ebola, the epidemic that should never have happened, 2015, 95100.
 22.
McKinley T, Cook AR, Deardon R. Inference in epidemic models withoutlikelihoods. Int J Biostat. 2009;5:24.
 23.
World Health Organisation. Disease Outbreak News. Geneva: WHO; 2014.
 24.
Carol R. Ebola awareness and prevention activities under way for refugees in Liberia, new stories in bahn Refugee, camp Liberia the UN Refugee Agency. Centers for Disease Control and Prevention, 2014. Outbreaks Chronology: Ebola Hemorrhagic Fever.
 25.
Bausch DG, Towner JS, Dowell SF. Assessment of the risk of Ebola virus transmission from bodily fluids and fomites. J Infect Dis. 2007;196(2):S142–7.
 26.
Engle R. Dynamic conditional correlation: A simple class of multivariate generalized autoregressive conditional heteroskedasticity models. J Bus Econ Stat. 2002;20(3):339–50.
 27.
International Convention. Burke et al. Bulletin of the World Health Organization. 1978;56(2):271–293.
 28.
Khan AS, Tshioko FK, Heymann DL, Le Guenno B, Nabeth P, Kerstiëns B, et al. The reemergence of Ebola hemorrhagic fever, Democratic Republic of the Congo 1995. J Infect Dis. 1999;179:S76–86.
 29.
Report of an International Commission. Ebola haemorrhagic fever in Zaire,1976. Bull World Health Organ. 1978;56:271–93.
 30.
Shimotsu K, Phillips PCB. Exact local Whittle estimation of fractional integration. Annals Statistics. 2006;33:1890–933.
 31.
Shimotsu K. Exact local Whittle estimation of fractional integration with unknown mean and trend. Econometric Theory. 2010;26:501–40.
 32.
Shimotsu K. Exact Local Whittle Estimation of Fractional Integration with Unknown Mean and Time Trend. In: Econometric Theory, 2008, forthcoming. 2008.
 33.
Robinson PM. Logperiodogram regression of time series with long range dependence. Ann Stat. 1995;23:1048–72.
 34.
National Center for Emerging and Zoonotic Infectious Diseases. Oyok T et al. 2001.
 35.
Engle RF. Autoregressive conditional heteroskedasticity with estimates of the variance of United Kingdom inflation”. J Monet Econ. 1982;30:987–1007.
 36.
Andersen TG, Bollerslev T, Diebold FX. Parametric and nonparametric volatility Measurement. In: Hansen ad LP, AitSahalia Y, editors. Handbook of Financial Econometrics. Amsterdam: NorthHolland, forthcoming; 2002.
 37.
Engle RF, Granger CW. Cointegration and error correction: representation, estimation and testing. Econometrica. 1987;55:251–76.
 38.
Diebold F. Elements of Forecasting. 2nd ed. South Western: Cincinnati; 2001. p. 254. ISBN 0324023936.
 39.
Forbes K, Rigobon R. No contagion, only interdependence: measuring stock market comovements”. J Financ. 2002;5:2223–61.
Author information
Additional information
Competing interests
Sebastian Braun is the national winner of the 2010/2011 HTAcademy scholarship, which was sponsored by Pfizer Deutschland GmbH.
Authors’ contributions
Nadhem Selmi PhD in Qantitatives Methods at the University of Sfax in Tunisia. His primary research interests are in contagion of financial crisis, long memory process, GARCH process and Ebola virus in West Africa. Selmi read and approved the final manuscript.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0), which permits use, duplication, 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 license, and indicate if changes were made.
About this article
Received
Accepted
Published
DOI
Keywords
 Ebola 2014
 Contagion
 Virus
 EGARCH
 MGARCHDCC