Next Article in Journal
Towards a More Robust Non-Rigid Robotic Joint
Next Article in Special Issue
An Empirical Algorithm for COVID-19 Nowcasting and Short-Term Forecast in Spain: A Kinematic Approach
Previous Article in Journal
Comparative Analysis of Classification Algorithms Using CNN Transferable Features: A Case Study Using Burn Datasets from Black Africans
Previous Article in Special Issue
Proposed Design of Walk-Through Gate (WTG): Mitigating the Effect of COVID-19
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of the Evolution of the COVID-19 Disease between Romania and Italy

1
Department of Exact Sciences, The University of Agricultural Sciences and Veterinary Medicine of Iasi, 700490 Iasi, Romania
2
Faculty of Energy Technology, University of Maribor, 8270 Krško, Slovenia
3
Center for Applied Mathematics and Theoretical Physics, University of Maribor, 2000 Maribor, Slovenia
4
Faculty of Civil Engineering, Transportation Engineering and Architecture, University of Maribor, 2000 Maribor, Slovenia
5
Institute of Mathematics, Physics and Mechanics, 1000 Ljubljana, Slovenia
*
Author to whom correspondence should be addressed.
Appl. Syst. Innov. 2020, 3(4), 44; https://0-doi-org.brum.beds.ac.uk/10.3390/asi3040044
Submission received: 14 September 2020 / Revised: 2 October 2020 / Accepted: 12 October 2020 / Published: 14 October 2020

Abstract

:
After the outbreak of COVID-19 in Italy, thousands of Romanian citizens who worked in Northern Italy, Spain or Germany returned to Romania. Based on the time-dependent susceptible–infected–recovered—SIR model, this paper compares the evolution of the COVID-19 disease between Romania and Italy, assuming that the parameter value of R0 in the time-dependent SIR model decreases to R1 < R0 after publicly announced restrictions by the government, and increases to a value of R2 < R1 when the restrictions are lifted. Among other things, we answer the questions about the date and extent of the second peak in Italy and Romania with respect to different values of R2 and the duration of the restrictions.

1. Introduction

In December 2019, several cases of pneumonia were detected in the Chinese province of Wuhan. Many common pneumonia and respiratory tract infections were excluded from the investigations, and a new coronavirus infection was discovered and confirmed. The new coronavirus was called SARS-CoV-2, and the disease caused by the virus was COVID-19. The virus is highly contagious and could be very dangerous, especially for the elderly and people with chronic diseases. To reduce the number of infected people countries have approached the problem in very different ways, mostly taking relatively strict measures to prevent the spread, some countries have started implementing measures relatively late, while others had a completely different approach of solving the problem, namely considering heard immunization, without setting restrictions.
At the beginning of 2020, information channels appeared more and more frequently in Romania about a new highly contagious virus that wreaked havoc shortly after its appearance in China. Given the fact that China is geographically located at a considerable distance from Romania, the majority of the population at that time did not consider that it could reach the territory of their country. Unfortunately, the global mobility of the population in the world yield a fast spreading of the virus and the first infections in Romania occurred within just one month. It was soon declared a pandemic by World Health Organization (WHO) officials, and the entire academic community mobilized to stop the spread of the disease. Several factors influence the spread of COVID-19. Predicting the course of the spread of the virus is very important for dealing with the disease and its consequences. To predict endemic and pandemic processes, researchers have developed various models.
Zhang and all [1] used the Poisson model to analyze the evolution of the disease in six Western European countries, and the statistical analysis allowed them to make predictions about the duration, maximum point and total percentage of the population that will be infected. Fanelli and Piazza [2] analyzed the dynamics of the new coronavirus in China, Italy and France. The authors simulated the effects of drastic isolation measures taken in Italy to reduce the infection rate.
Qianying et al. [3] proposed a conceptual model for Covid-19 outbreak in Wuhan, Kucharski et al. [4] estimated how the transmission of Covid19 had varied overtime during January 2020 and February 2020 in Wuhan. Rafiq et al. [5] applied a deterministic model to forecast the spread of COVID-19 in India and Okuonghae and Omame [6] analyzed the impact of non-pharmaceutical control measures on the population dynamics of the new Covid-19 in Lagos, Nigeria using a mathematical model.
Since there are many links between Romania and Italy, we try to analyze in this paper whether the evolution of the virus spread in Romania was similar to that in Italy [7]. The occurrence of the COVID-19 disease in Romania was not due to tourism or direct trade, cultural and economic exchange with China, but to the mostly economic and work related mobility of an important part of the Romanian population throughout Europe, especially to Italy, Spain, Germany or Great Britain. We chose a comparative analysis of the evolution of the spread in Italy and Romania, as Italy was the first European country to fight the major outbreak of Covid-19, it is notorious that many Romanian citizens work, live and travel to Italy, and a considerable amount of statistical data was available for the analysis.
According to the Italian National Institute statistics website [8], many Romanian citizens worked in the areas Emilia-Romagna, Lombardy Trentino Alto-Adige, Friuli Venezia Giulia, that were, in April 2020, heavily affected by the disease COVID-19.
The fact that the first documented case of Romanian citizen infected with COVID-19 was a man of Italian origin travelling to Romania and the next two cases were Romanian citizens traveling to areas in Italy affected by the COVID-19 disease comes in support of this theory of causality
The large number of isolated persons, with whom the statistical data start at 5600 [9,10], is due to the Romanian citizens who entered the country from the time of the first travel restrictions in Lombardy. The number of infected persons continued to increase, either because of the lack of responsibility of the infected persons who did not tell their doctors or friends about their living in the dangerous area of Italy, or because they did not know that they had been infected with the virus. The incubation period (i.e., the time between exposure to the virus and the onset of symptoms) for COVID-19 infected people is estimated to be on average five to six days, but it can be up to fourteen days [11]. The fact that many people, who came into contact with the virus showed no specific symptoms also facilitated the spread of the virus. Consequently, many people were hospitalized and a large number of medical staff was infected, too.
Government measures in Romania were taken as soon as the first cases of the disease were confirmed. The first case was registered on 27 February 2020 and two more the next day. On the 1st March 2020 the terms “quarantine” and “isolated person” were defined. The first data show that the number of people quarantined or in isolation was about 5600 [9,10], after medical check points were set in place at the border for citizens coming from the restricted areas of Italy.
The next restrictions were: a ban on visiting family members in hospitals, a ban on public or private events in open or enclosed spaces with more than 1000 people (8 March 2020), the suspension of the pre-school education process and the suspension of car, air and rail transport with Italy (10 March 2020).
When the number of 200 infected was exceeded (16 March 2020), a state of emergency was declared and other restrictive conditions were imposed: restaurants, hotels, bars would be closed, cultural, artistic, scientific, sports and entertainment events would be suspended, flights with Spain would be suspended and the restriction for travelling to/from Italy extended, universities would go online.
On 21 March 2020, when there were 367 confirmed infected patients, 4207 people in official quarantine and 55198 people in isolation, the restrictions were as follows: the borders were closed, shopping centers were closed, movement of people outside the home was prohibited between 22.00–06.00, dental surgeries, beauty centers were also closed. For private events where liturgical/religious acts are celebrated (weddings, baptisms, funerals), the number of participants had been limited to eight persons.
After statistical analysis of the data from [9,10] using variance and mean value analysis using the tests of Student and Fisher, it was found that the evolution of the COVID-19 infection in Romania differs significantly from that in Italy (p-value 0.001).
The national reports of Romania [12] and Italy [8] for the age distribution of the population reveals that in Romania (Italy) the older population, considered vulnerable to the COVID-19 infection, represents 36.4% (58.8%), respectively.
Italy imposed a restriction on national quarantine on 9 March 2020, when the number of infected persons was 9172, the number of daily infected persons was over 1000 and the number of deaths was 463 per day. In Romania, the maximum restrictions were imposed from 22 March 2020, when the number of infected persons was 433, the number of infected persons per day was 66 and the number of deaths on that day was three. With a population that represents one third of the Italian population, Romania reported 1/15 infected compared to Italy and the number of deaths is less than 1/35 compared to Italy.
Based on real data (see Figure 1) related to the first wave of the COVID -19 pandemic in Italy and Romania, we estimate the peak date, extent and duration of the second wave. For the calculations we use the Mathematica program to interpolate the functions from the time-dependent SIR model. We also study how the duration of the restrictions and the date when the restrictions were introduced in Italy and Romania affect the first wave and how they affect the second wave of the COVID -19 pandemic.
The paper is organized as follows: in the next section we consider the modified SIR model and two interesting theoretical examples. According to the presented model we adjust the parameters to the real-life evolution for Italy and Romania. In Section 3, we present the results of the corresponding SIR models for Italy and Romania. In Section 4, we present some remarks on immunity and finally, we answer in conclusions the questions asked in Section 2.

2. Time-Dependent SIR Model

In this section, the classical SIR model and the time-dependent SIR model are briefly described [13]. In the classical SIR model:
  d S d τ = β I S , d I d τ = β I S γ I , d R d τ = γ I ,  
where S ( τ ) , I ( τ ) and R ( τ ) denotes the number of susceptible persons, infected persons and persons recovered at time τ , respectively, the population of size N is closed, i.e., N = S ( τ ) + I ( τ ) + R ( τ ) . The constants β N and γ stand for: β N denotes the (average) number of susceptible persons an infected person comes into contact with (and is possible to spreads the disease), while 1 γ denotes the average time an infected person takes to recover (i.e., the time of being infectious). Therefore, the quotient R 0 = β N γ is constant and represents the number of new infections caused by an infected person before he/she recovers. In this paper R 0 is considered as time-dependent: R 0 = R 0 ( t ) ; changing after restrictions are announced.
Rewriting Equation (1) in a non-dimensional form, we define:
    u = S N , v = I N , w = R N , t = γ τ
to obtain
  d u d t = R 0 u v , d v d t = v ( R 0 u 1 ) , d w d t = v .
From the dynamical point of view [14,15], the ( u , v ) plane is invariant, v = 0 is the line of singular points, and ( u , v ) = ( 0 , 1 ) is unstable if R 0 > 1 , while all other singular points v [ 0 , 1 ) are stable. Note that in Model (3) all variables are less than (or equal to) one (since representing the percentage of the total population). Note that the exact formula characterizing R 0 is β S 0 γ where S 0 is the starting number of susceptible persons. However, if the population N is quite big and there are only few initially infected, the density S 0 / N is near 1, thus the approximation β S 0 γ β N γ is possible.
We assume (see Figure 2) the following time dependence for R 0 = R 0 ( t ) :
  • After introduction of the restrictions/interventions, at time t 1 , the value of R 0 decreases to a certain (lower constant) value R 1 (i.e., R 0 > R 1 ),
  • While the restrictions are in force this value remains constant (for time t [ t 1 , t 2 ] ),
  • After the end of the restrictions/interventions, after time t 2 , the value of R 0 increases to a certain (higher) value R 2 (i.e., R 2 > R 1 ),
  • The total duration of restrictions is denoted by Δ t = t 2 t 1 ,
  • The transitions R 0 R 1 and R 1 R 2 appear in a time comparable to 1 γ , which is much less than Δ t = t 2 t 1 .
This type of time dependent R 0 in a SIR-model causes the evolution of the infected-fraction v = v ( t ) , as shown in Figure 3. The maximum of the 1st wave is v 1 * and appears at time T 1 . The maximum of the 2nd wave is v 2 * and appears at time T 2 .
By R 3 a third value for R 0 is denoted (see Figure 4), as one considers again some restrictions to avoid the in-avoidable 2nd peak. Note also that the travels were one of the first things that were blocked by the restrictions, so closed population assumption is legitimate under the regime R 1 . It is known that even in the restricted period a large number of people was traveling from Romania to Germany, for agriculture. The traveling between countries with the same risk-factor was possible after the end of the 1st wave. Even the traveling from a country with higher risk factor to a country with lower risk factor was possible, but migrations were quite limited at this time and a quarantine was mandatory. We assumed the traveling persons were healthy and tested for COVID-19 and they return in the period close to some (low) multiple of 1 γ . This justifies the assumption of closed population also under the regime R 2 .
First we adjusted the value of R 0 at the beginning from the data [16,17] by simply taking logarithms of I = I ( t ) and adjusting the data ( t , log I ) for the linear regression line. For Italy we found R 0 [ 3.72 , 5.01 ] , while for Romania we found R 0 [ 3.56 , 4.8 ] . The lower value corresponds to value 1 γ = 11 days, while the higher value corresponds to 1 γ = 14 days [11]. Next, we read from the graph I = I ( t ) [16,17] the date of the beginning of the restrictions in Italy and Romania (i.e., we read t 1 ). Then we adjust the value of R 1 for Italy and Romania so that the maximum value of the first wave v 1 * appears at real time T 1 , which in turn was read from the graph [16,17] according to the real data. For the transition-time of R 0 = R 0 ( t ) from R 0 to R 1 and from R 1 to R 2 we assumed a value comparable to 1 γ , but still in accordance to the actual measures taken by authorities. The value of R 1 and γ was determined to adjust the evolution of v = v ( t ) the to the real data. The values for R 1 were checked for values R 1 [ 0.4 ,   0.9 ] , according to [18].
The crucial part of the model is the variation of the value for R 3 (see also [19]). Additionally, we vary the value of t 2 to get its effect on T 2 and on the value of the extent fraction v 2 * v 1 * .
In the paper we answer the following questions:
  • (Q1) How do the publicly announced interventions/restrictions affect parameter R 1 ?
  • (Q2) What is the time delay, Δ T = T 2 T 1 , in the evolution of the spread of the disease for both countries?
  • (Q3) How does parameter R 2 affect the fraction v 2 * v 1 * ?
  • (Q4) What is the effect of the parameter R 2 on the time T 2 ?
  • (Q5) How does the parameter Δ t = t 2 t 1 (the duration of the restrictions) affect the fraction v 2 * v 1 * ?
  • (Q6) How does the parameter Δ t affect the time T 2 ?
  • (Q7) What is the effect of possible re-restrictions on the extent and datum of the second peak?
To answer these questions, we look at the system (3) with the real data output and vary the values of R 1 , R 2 , Δ t for Italy and for Romania. For 1 γ we used/assumed the value of 11–14 days, as reported in news. To obtain numerical solutions to Equations (3) we used Mathematica 10.4. In the sequel model (3) with R 0 = R 0 ( t ) as shown in Figure 3.
  d u d t = R 0 ( t ) u v , d v d t = v ( R 0 ( t ) u 1 ) , d w d t = v .
were adjusted to real data for Italy [16] and Romania [17].
In the time dependent SIR model (4) the corresponding systems of ordinary differential equations are non-linear and therefore it is impossible to obtain the exact solution to a particular initial condition problem. Therefore we use the Mathematica 10.4. command NDSolve , which gives results in terms of interpolating function objects. These interpolating functions are representing the numerical solution to a particular initial condition problem as an approximate function. For the time integration the Runge–Kutta 4th order method was used, while for the interpolation of the solution function the Hermite function was used. We used the following Mathematica code:
In [ 1 ] : = R 0 [ t _ ] : = Piecewise [ { { R 0 ,   0 t N 1 } , { ( R 0 R 1 / N 2 ( t N 1 ) + R 0 , N 1 < t   N 1 + N 2 } , { R 1 , N 1 + N 2 < t   N 1 + N 2 + N 3 } , { ( R 2 R 1 / N 4 ( t ( N 1 + N 2 + N 3 + N 4 ) ) + R 2 , N 1 + N 2   + N 3 < t N 1 + N 2 + N 3 + N 4 } , { R 2 , N 1 + N 2 + N 3 + N 4 < t   N 1 + N 2 + N 3 + 15 N 4 } } ]
In [ 2 ] : = sol 1 = NDSolve [ { v [ t ] = = v [ t ] ( R 0 [ t ] u [ t ] 1 ) , u [ t ] =   = R 0 [ t ]   u [ t ]   v [ t ] , w [ t ] = = v [ t ] , v [ 0 ] = = r 1 , u [ 0 ] = = 1 r 1 , w [ 0 ]   = = 0 } , { v , u , w } , { t , 0 , N 1 + N 2 + N 3 + 15   N 4 } ]
In the model r 1 is the initial value of infected fraction v ( 0 ) = r 1 .
Theoretical Example 1. In this example, we show the sensitivity of the output in model (4) on Δ t , with R 2 = R 0   = 4.0 and R1 = 0.5; i.e., R 0 jumps back to its »original« value. The evolution times t 1 and t 2 under each level are very important. For simplicity, in Figure 4, Figure 5, Figure 6 and Figure 7 we have two slightly different cases of jumps from R 0 and the »boomerang effect« looks completely different. In Figure 5, the second peak is larger than the first one, while in Figure 7, the second peak much smaller. We see that the SIR model with the time-dependent R 0 is very sensitive to the duration of the restrictions (i.e., to the value of Δ t = t 2 t 1 ).
Note that the only difference of R 0 = R 0 ( t ) in Figure 5 and Figure 7 are the time intervals (the levels of higher R 0 H = 4.0 and lower R 0 L = 0.5 is the same).
The time intervals in Figure 4 are as follows: for t [0,3.5] we have R 0 = 4.0 , for t [3.5,5.5] we have the descending line, for t [5.5,10] we have R 1 = 0.5 , for t [10,15] we have the ascending line and for t [15,25] we have again R 2 = 4.0 .
The time intervals in Figure 6 are as follows: for t [0,4] we have R 0 = 4.0 , for t [4,6] we have the descending line, for t [6,10] we have R 1 = 0.5 , for t [10,15] we have the ascending line and for t [15,25] we have again R 2 = 4.0 .
  • For Italy, the parameters of the model (4) are adjusted to real data are found to be
    r 1 = 2.319 · 10 7 ;   N 1 = 1.183 ;   N 2 = 3.630 ;   N 3 = 0.726 ;
    N 4 = 1.6 ; R 0 = 4.28 ; R 1 = 0.4 ; 1 γ = 13   d a y s .
The parameter R 2 was varied in the range of [1.2,3]. Parameter value N 2 + N 3 = Δ t was varied in the range of [4.3,6.3]. Note that in the dimensionless model (4) this corresponds to a range from almost two months to more than three months.
For Italy, the time is in units τ × 13 days, for Romania the time is in units τ × 14 days.
  • For Romania, the parameters of the model (4) adjusted to real data are found to be:
    r 1 = 1.5 · 10 7 ;   N 1 = 1.42 ;   N 2 = 2.60 ;   N 3 = 1.33 ;
    N 4 = 1.6 ; R 0 = 4.80 ; R 1 = 0.83 ; 1 γ = 14   d a y s .
The parameter R 2 was varied in the range of [1.5,2.5]. Parameter N 2 + N 3 = Δ t was varied in the range of [2,8]. Note that in the dimensionless model (4) this corresponds to a range from one to almost four months.

3. Results

Above best-fit-values with R 0 = R 0 ( t ) presented in Figure 8 and Figure 9 for Italy and Romania are varied slightly (i.e., near the best fit values).

3.1. Italy

The overall evolution of the infected fraction in Italy [20] for COVID-19 is presented in the following Figure 10, Figure 11 and Figure 12.
The estimated time of formation of the 2nd peak is 7.8 months after the beginning (in our model 22 February 2020) of the epidemic. According to [16] our model predicts the 2nd peak about the 15th of October 2020, if R 2 = 1.6 . Assuming R 2 = 1.6 , the extent of the 2nd wave is 44.5-times larger than the extent of the first wave:
v 2 * v 1 * = 44.5
The 1st peak is at 1.79‰, while the 2nd peak is at 7.99%, if R 2 = 1.6 . Figure 13 shows the relation between fraction v 2 * v 1 * and R 2 in range of R 2 [ 1.5 , 2.1 ] .
Figure 14 shows that T 2 tends to 5 months, if R 2 exceeds 3 (i.e., the 2nd peak could appear in August). But, if R 2 [ 1.6 , 2.0 ] the 2nd peak might appear between September and November.
Varying R2 and N2 in model (4) we generated the results for extent v 2 * and time T 2 . For fitting to linear regression line we used the Mathematica command LinearModelFit while for fitting to exponential dependence f ( x ) = a + b e c x the command NonlinearModelFit of Mathematica 10.4 was used.
The fraction v 2 * v 1 * depends linearly on R 2 . The (linear) correlation (Pearson) coefficient is found to be 0.999865 . The functorial dependence is found to be v 2 * v 1 * = 97.31036 + 89.0063 · R 2 . In Table 1 the estimates, standard errors, t-Statistics and the corresponding p-value for coefficients a and b in y = a + b x are presented. The p-value stands for the null-hypothesis that the coefficients a , b are equal to zero ( y v 2 * v 1 * and x R 2 ).
The fitted line and the simulated data are shown in Figure 13.
Next, we fitted T 2 in months vs. R 2 as follows T 2 ( R 2 ) = a + b e c · R 2 . We found a = 4.7714696 , b = 242.350437 and c = 2.80811033 . The fitted line and the simulated data are shown in Figure 14.
We fitted the extent fraction v 2 * v 1 * in relation to Δ t and found v 2 * v 1 * = 44.8798615 0.09911836 · Δ t . In Table 2 the estimates, standard errors, t-Statistics and the corresponding p-value for coefficients a and b in y = a + b x are presented. The p-value stands for the null-hypothesis that the coefficients a , b are equal to zero ( y v 2 * v 1 * and x R 2 ).
The fitted line and the simulated data are shown in Figure 15. The correlation coefficient was found to be 0.943522 .
Finally, we fitted the date of the second peak vs. Δ t in terms of τ . The correlation was found to be T 2 = 3.6552199083 + 0.8801832505 · Δ t . The functorial dependence was found to be 0.999797 with the statistical table for the estimates, standard errors, t-Statistics and the corresponding p-value for coefficients a and b in y = a + b x as stated in Table 3. The p-value stands for the null-hypothesis that the coefficient s a , b are equal to zero ( y T 2 and x Δ t ).
The fitted line and the simulated data are shown in Figure 16.

3.2. Romania

The evolution of the infected fraction in Romania [21] for COVID-19 for R 2 = 1.6 is presented in the Figure 17, Figure 18 and Figure 19. In all figures, τ represents time multiplied by 14 days; i.e., 15 τ means 210 days or cca. 7 months (after 1 March).
According to [17] and [22] the estimated time of the 2nd peak for R 2 = 1.6 is 7 months after the beginning (in our model 1 March 2020) of the epidemic; i.e., about 1 October 2020. For R 2 = 2.5 the 2nd peak is estimated to about 25 July. If R 2 = 1.5 the extent of the 2nd wave is about 15 -times larger than the extent of the first wave (the 1st peak is at 3.9 ‰, while the 2nd peak is estimated at 6.1%):
v 2 * v 1 * = 15.64
Figure 20 shows the relation between fraction v 2 * v 1 * and R 2 in range of R 2 [ 1.5 , 2.5 ] .
In Figure 21 we see that the 2nd peak could appear also about 1 August 2020, if R 2 3 . Therefore, most likely, the 2nd peak will appear between the beginning of August 2020 and end of October 2020, if R 2 1.3 .
Varying R2 and N2 in model (4) we generated the results for extent v 2 * and time T 2 . For fitting to linear regression line we used the Mathematica command LinearModelFit , while for fitting to exponential dependence f ( x ) = a + b e c x the command NonlinearModelFit of Mathematica 10.4 was used [15].
The fraction v 2 * v 1 * depends linearly on R 2 . The linear (Pearson) correlation coefficient is found to be 0.999273 . The functorial dependence is found to be v 2 * v 1 * = 48.095079 + 42.868254 · R 2 . In Table 4 the estimates, standard errors, t-Statistics and the corresponding p-Value for coefficients a and b in y = a + b x are presented. The p-Value stands for the null-hypothesis that the coefficients a , b are equal to zero ( y v 2 * v 1 * and x R 2 ).
The fitted line and the simulated data are shown in Figure 20.
Next, we fitted T 2 in months vs. R 2 as follows T 2 ( R 2 ) = a + b e c · R 2 . We found a = 4.409702643 , b = 47.626146393 and c = 1.852615293 . The fitted curve and the simulated data are shown in Figure 21.
We fitted the extent ratio v 2 * v 1 * = a + b · Δ t and found a = 30.490269 and b = 0.303047 . In Table 5 the estimates, standard errors, t-Statistics and the corresponding p-value for coefficients a and b are presented. The p-value stands for the null-hypothesis that the coefficients a , b are equal to zero.
The fitted line and the simulated data are shown in Figure 22. The correlation coefficient was found to be 0.99542 .
Finally, we fitted the date of the second peak vs. Δ t . The correlation was found 0.99981 . The functorial dependence was found to be T 2 = 3.758780098 + 0.59492099 · Δ t . In Table 6 the estimates, standard errors, t-Statistics and the corresponding p-value for coefficients a and b in y = a + b x are presented. The p-value stands for the null-hypothesis that the coefficients a , b are equal to zero ( y T 2 and x Δ t ).
The fitted line and the simulated data are shown in Figure 23. The dependence T 2 = T 2 ( Δ t ) is in units of τ for both variables.
Theoretical Example 2. In this example, we show the sensitivity of the output in the modified model in case of re-restrictions after the beginning of the 2nd peak, as shown in Figure 24. In Figure 25 and Figure 26 we see that the extent of the 2nd peak is always lower but depends on the time t 3 of the reducing of R 2 to R 3 . The evolution for case in Figure 24 and Figure 27 is shown in Figure 25 and Figure 26, respectively.

4. Remarks on Immunity

The coronavirus SARS-CoV-2 has only been circulating in human hosts since December 2019, which means that it is simply impossible to know whether immunity to the disease will last longer than nine months. In the meantime, the results only confirm that COVID-19 patients can maintain the adaptive immunity to SARS-CoV-2 for two weeks post-discharge. Evidence from other coronaviruses (e.g., 2002s SARS 1, or 2012s middle Easter respiratory syndrome) suggests that immunity probably lasts longer than that [23]. The duration of immunity to SARS-CoV-2 by a vaccine (once developed) is estimated to be 6–18 months [24].
However, it is still unclear whether someone with immunity could spread the coronavirus to others while fighting off a second infection. If the immune response were strong enough to crush the virus quickly, the person probably wouldn’t transmit it further. A weaker response that allowed some viral replication might not prevent transmission, especially as people without symptoms are known to transmit the coronavirus [23].
Let us just comment on the evolution in Romania: in the 1st wave, the total number of infected persons was approximately 1% of the population, but, after the 2nd wave the (maximum) total number of infected is estimated to be about 60%. Assuming immunity for approximately 12 months this yields the reduction of the number of susceptible persons to 4 10 N . Possibly the extent of the estimated 3rd wave (again assuming constant R 2 = 1.6 ) will be at most 3% after approximately 20 τ (see green curve in Figure 28). It is possible that due to the continuous (and prompt) reactivation of restrictions already the 2nd wave will not be of such extent, as predicted by the proposed model.
Note that the proposed model must be modified after the 2nd wave, because the regime described in Figure 8 and Figure 9 is no longer valid after the 2nd wave. Most likely the value of parameter R 0 (with a slight random error) will remain above R 0 = 1 . Such a regime is likely to be in operation in Sweden from the beginning (in February 2020) [9,10].
The evolution of COVID-19 in Sweden (i.e., the function of removed persons R = R ( t ) ) is a strictly increasing function (see [25]; real time evolution of the removed persons R = R ( t ) ), while for Italy the same curve runs almost horizontally between May and August (see [26]; real time evolution of the removed persons R = R ( t ) ).
Finally, note that also in Israel [9] the evolution of the COVID-19 disease (until September 2020) is composed of two waves. According to the shape of R = R ( t ) (a large horizontal line between 25 April and 15 June 2020 corresponds to the time period between the end of the 1st and the beginning of the 2nd wave) the evolution in Israel also obeys the proposed model (i.e., the variation of R 0 ). The large difference between the extent of the 1st and 2nd wave in Israel is probably due to larger value of R 2 compared to Romania and Italy.

5. Conclusions

Since at the time of the outbreak of the COVID epidemic in Europe, many Romanian citizens returned to Romania en masse from Italy, Spain, Great Britain, etc. when the first wave of the epidemic started in Romania, the first wave was relatively more pronounced in Romania than in Italy. However, the number of deaths in Romania during the first wave was far from being as high as in Italy. The proposed model answers the questions asked in the introduction and forecasts about the same extend of the second peak (maximum at about 8%) for both countries (in relative terms).
Note that the proposed model cannot answer the question when (if) will the epidemics end—it just rises up the question on the relation between the time and extent of the 1st and 2nd peak. Neither the prevalence of the disease can be considered by SIR models. To explain how R 0 changes, note that in the SIR model the population is assumed to be closed, and N cannot change, neither the rate of recovery γ . Therefore, the real change appears on the parameter β . A simple application of the definition of conditioned probability yields β =   p c , where c is a contact parameter (i.e., the probability of contact between a single infected and a single susceptible) and p is a contagion parameter (i.e., the probability that if a contact occurs, then the susceptible becomes infected). However, p cannot change with restrictions, since it depends on the infectivity of the disease. Therefore, a lockdown-type restriction makes c smaller, since it reduces the contacts. The duration of restrictions (i.e., Δ t ) and the social distance approach as considered in the model is reflected linearly in terms of the fraction of the maximal extent and in terms of the expected datum of the 2nd peak. On the other hand, the relation T 2   =   T 2 ( R 2 ) seems to be exponential. The model predicts the influence of the beginning and duration of the restrictions to the extent and duration of the 2nd wave.
The number of total infected in the 2nd wave seems to be much higher than the number infected in the 1st wave (note that the 1st peak is measured in per mils while the 2nd peak in per cents). A reasonable question is of course also how many waves will occur. That depends, of course, on the duration of the actual immunity and on whether/when we will have a coronavirus vaccine. Similar to the assessment of the duration of immunity, the answer to this question should be compared with cases from real-life. From history we know of cases where the virus has decimated the population, so caution is more than appropriate in combating COVID-19. One of the most affected diseases in Europe was Spanish flu. In 1918, the Spanish flu infected 500 million people with the H1N1 flu virus in four waves, which was one third of the world’s population at that time (and killed 17-50 million people). If there are many waves, and if the immunity to SARS-CoV-2 is as strong as the immunity to SARS-CoV-1, the next waves will not be as strong as the 2nd wave. However, this is left for the future work.
To conclude we give answers to the questions (Q1)–(Q7):
  • Answer to (Q1): In Italy, the final R 2 =0.4, however the average R 2 = 2.017 , with 1 γ = 13 days.
  • In Romania, the final R2 = 0.83 , however the average R 2 = 2.14, with 1 γ = 14 days.
  • Answer to (Q2): The time delay Δ T varies with R 2 and t 2 , as it can be seen in Figure 16 and Figure 23. Varying of R 2 makes sense for the future evolution, while varying of t 2 , once the restrictions were eased, has just a theoretical meaning. For Italy, Δ T varies from 2.5 months (if R 2 = 3) to 7.3 months (if R 2 = 1.5).
  • Answer to (Q3): From Figure 13 and Figure 20 we see that the relation is linear and ascending. For R 2 > 1.2 the extent fraction is higher than 10. From Figure 25 and Figure 26 we see that after repeating the restrictions the extent of the 2nd peak can be much lower (i.e., one third) than predicted by Figure 2 and Figure 3.
  • Answer to (Q4): The relation T 2 = T 2 ( R 2 ) seems to obey the exponential law T 2 ( R 2 ) = a + b e c R 2 . The relation is shown in Figure 14 and Figure 21. The exponential dependence is more considerable in case of Romania (see Figure 21). For higher values of R 2 the period between the peaks, Δ T , is lower. Furthermore, the second peak appears faster and has a greater extent.
  • Answer to (Q5): The longer the restriction period, the smaller the fraction v 2 * v 1 * . From Figure 15 and Figure 22 we see that the relation is linear (descending). This means that prolonging the duration of the restrictions leads to decreasing of the extent of the second wave of infection.
  • Answer to (Q6): The parameters Δ t and T 2 are linearly related (ascending). For longer restriction periods one obtains the 2nd peak later. Also, the appearance of the second wave is of lower extent.
  • Answer to (Q7): The extent of the 2nd wave in case of re-restrictions is always lower and depends on the time of reacting on the start of the 2nd wave. In the middle of September in Italy, a restriction concerning the closure of nightclubs and the obligatory use of face-masks after 18:00 was introduced. This is a concrete example of lowering R 2 to R 3 , as noted in the Theoretical Example 2.

Author Contributions

C.C. and E.B. contributed mostly the statistical analysis, to the literature re-view, conceptual design, J.A. contributed to the statistical interpretation of results, conceptual design, quality check of paper, and interpretation of results, B.F. and M.M. were responsible for mathematical formulation of the model, programming in Mathematica, design of the algorithms and quality check of the entire paper. All authors have read and agreed to the published version of the manuscript.

Funding

The last three authors acknowledge the financial support of the Slovenian Research Agency. In particular, the third author was supported by research core funding by Slovenian Research Agency, Research Program Applied Electromagnetics No. P2-0114, the fourth author was supported by research core funding by Slovenian Research Agency, Research Program Applied Mathematics, Theoretical Physics and Intelligent Systems No. P1-0306 and the fifth author was supported by research core funding by Slovenian Research Agency, Research Program Algebras and Rings No. P1-0288 and by the project by Slovenian Research Agency, Research Program Algebraic Methods and use of differential Equations, No. N1-0063.

Acknowledgments

The authors would like to thank to all the referees of MDPI for their remarks, which considerably improved the paper.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Nomenclature

S ( τ ) number of susceptible persons
I ( τ ) number of infected persons
R ( τ ) number of persons recovered at time τ
Nthe population size
β N the (average) number of susceptible persons an infected person
1 γ the (average) time an infected person takes to recover
R 0 = β N γ the number of new infections caused by an infected person before he/she recovers
S 0 the starting number of susceptible
t = γ τ the modified time
v ( t ) the infected-fraction
R 1 , R 2 , R 3 different values for R 0
( u , v ) invariant plane
u = S N , w = R N non-dimensional variables of SIR model
Δ t = t 2 t 1 the duration of the restrictions
Δ T = T 2 T 1 the time delay between successive waves
r 1 the initial value of infected fraction v ( 0 ) = r 1
a, b, cstatistical coefficients

References

  1. Zhang, X.; Ma, R.; Wang, L. Predicting turning point, duration and attack rate of COVID-19 outbreaks in major Western countries. Chaos Solitons Fractals 2020, 135, 109829. [Google Scholar] [CrossRef] [PubMed]
  2. Fanelli, D.; Piazza, F. Analysis and forecast of COVID-19 spreading in China, Italy and France. Chaos Solitons Fractals 2020, 135, 109761. [Google Scholar] [CrossRef] [PubMed]
  3. Lin, Q.; Zhao, S.; Gao, D.; Lou, Y.; Yang, S.; Musa, S.S.; Wang, M.H.; Cai, Y.; Wang, W.; Yang, L.; et al. A conceptual model for the coronavirus disease 2019 (COVID-19) outbreak in Wuhan, China with individual reaction and governmental action. Int. J. Infect. Dis. 2020, 93, 211–216. [Google Scholar] [CrossRef] [PubMed]
  4. Kucharski, A.J.; Russell, T.W.; Diamond, C.; Liu, Y.; Edmunds, J.; Funk, S.; Eggo, R.M.; Sun, F.; Jit, M.; Munday, J.D.; et al. Early dynamics of transmission and control of COVID-19: A mathematical modelling study. Lancet Infect. Dis. 2020, 20, 553–558. [Google Scholar] [CrossRef] [Green Version]
  5. Rafiq, D.; Suhail, S.A.; Bazaz, M.A. Evaluation and prediction of COVID-19 in India: A case study of worst hit states. Chaos Solitons Fractals 2020, 139, 110014. [Google Scholar] [CrossRef] [PubMed]
  6. Okuonghae, D.; Omame, A. Analysis of a mathematical model for COVID-19 population dynamics in Lagos, Nigeria. Chaos Solitons Fractals 2020, 139, 110032. [Google Scholar] [CrossRef] [PubMed]
  7. Istituto Nazionale di Statistica (ISTAT). Available online: https://www.istat.it/en/archivio/immigrants (accessed on 23 April 2020).
  8. Istituto Nazionale di Statistica (ISTAT). Demographic Indicators, 7 February 2019. Available online: https://www.istat.it/it/files//2019/02/Indicatoridemografici2018_EN.pdf (accessed on 1 May 2020).
  9. European Centre for Disease Prevention and Control. Available online: https://www.ecdc.europa.eu/en/publications-data/download-todays-data-geographic-distribution-covid-19-cases-worldwide (accessed on 5 May 2020).
  10. Available online: https://github.com/owid/covid-19-data/tree/master/public/data (accessed on 13 October 2020).
  11. Transmission of SARS-CoV-2: Implications for Infection Prevention Precautions. Available online: https://www.who.int/news-room/commentaries/detail/transmission-of-sars-cov-2-implications-for-infection-prevention-precautions (accessed on 10 September 2020).
  12. Institutul Naţional de Sănătate Publică. Raport Naţional al Stării de Sănatatea Populaţiei 2017. Available online: http://insp.gov.ro/sites/cnepss/wp-content/uploads/2018/11/SSPR-2017.pdf (accessed on 1 May 2020).
  13. Britton, N.F. Essential Mathematical Biology; Springer Undergraduate Mathematics Series; Springer Ltd.: London, UK, 2003; ISBN 185233536X. [Google Scholar]
  14. Boccaletti, S.; Ditto, W.; Mindlin, G.; Atangana, A. Modeling and forecasting of epidemic spreading: The case of Covid-19 and beyond. Chaos Solitons Fractals 2020, 135, 109794. [Google Scholar] [CrossRef] [PubMed]
  15. Chasnov, J.R. Mathematical Biology, Lecture Notes for MATH 4333; The Hong Kong University of Science and Technology: Hong Kong, China, 2009. [Google Scholar]
  16. Available online: https://covid19.who.int/region/euro/country/it (accessed on 9 June 2020).
  17. Available online: https://covid19.who.int/region/euro/country/ro (accessed on 9 June 2020).
  18. Available online: https://www.weforum.org/agenda/2020/05/covid-19-what-is-the-r-number/ (accessed on 9 June 2020).
  19. Giordano, G.; Blanchini, F.; Bruno, R.; Colaneri, P.; Di Filippo, A.; Di Matteo, A.; Colaneri, M. Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy. Nat. Med. 2020, 26, 855–860. [Google Scholar] [CrossRef] [PubMed]
  20. Coronavirus Disease (Covid-19) Outbreak. Available online: https://www.who.int/emergencies/diseases/novel-coronavirus-2019 (accessed on 1 January 2020).
  21. The Government of Romania through the Authority for the Digitization of Romania. Available online: https://stirioficiale.ro/informatii (accessed on 1 May 2020).
  22. Chen, Y.-C.; Lu, P.-E.; Chang, C.-S.; Liu, T.-H. A Time-dependent SIR model for COVID-19 with Undetectable Infected Persons. arXiv 2003, arXiv:2003.00122v6. Available online: https://arxiv.org/pdf/2003.00122.pdf (accessed on 28 April 2020).
  23. After Recovering from COVID-19, Are You Immune? Available online: https://www.livescience.com/covid-19-immunity.html (accessed on 13 October 2020).
  24. Amanat, F.; Krammer, F. SARS-CoV-2 Vaccines: Status Report. Immunity 2020, 52, 583–589. [Google Scholar] [CrossRef] [PubMed]
  25. Available online: https://www.worldometers.info/coronavirus/country/sweden/ (accessed on 13 October 2020).
  26. Available online: https://www.worldometers.info/coronavirus/country/italy/ (accessed on 13 October 2020).
Figure 1. The comparative evolution of the number of cases of diseases per day of evolution of the COVID-19 infection: (a) Romania, (b) Italy.
Figure 1. The comparative evolution of the number of cases of diseases per day of evolution of the COVID-19 infection: (a) Romania, (b) Italy.
Asi 03 00044 g001
Figure 2. Time dependent R 0 , as following from the assumptions described above. The restrictions/interventions start at time t 1 and end at time t 2 . During the restrictions it is assumed that the value of R 0 decreases, after the restrictions have ended it is assumed that the value of R 0 increases again, but to a slightly lower value than before the restrictions were introduced.
Figure 2. Time dependent R 0 , as following from the assumptions described above. The restrictions/interventions start at time t 1 and end at time t 2 . During the restrictions it is assumed that the value of R 0 decreases, after the restrictions have ended it is assumed that the value of R 0 increases again, but to a slightly lower value than before the restrictions were introduced.
Asi 03 00044 g002
Figure 3. First and second wave of the evolution, as derived from the assumption R 0 = R 0 ( t ) and the susceptible–infected–recovered (SIR) model; Equation (3). The start/end time of the restrictions: t 1 and t 2 are followed by two peaks at time T 1 and T 2 , respectively. The extension of the successive waves are v 1 * and v 2 * respectively.
Figure 3. First and second wave of the evolution, as derived from the assumption R 0 = R 0 ( t ) and the susceptible–infected–recovered (SIR) model; Equation (3). The start/end time of the restrictions: t 1 and t 2 are followed by two peaks at time T 1 and T 2 , respectively. The extension of the successive waves are v 1 * and v 2 * respectively.
Asi 03 00044 g003
Figure 4. Double transition R 0 = R 0 ( t ) with the parameter values R2 = R0 = 4.0 and R1 = 0.5. The start of restrictions is at time t 1 = 3.5 , the end of the restrictions is at time t 2 = 10 .
Figure 4. Double transition R 0 = R 0 ( t ) with the parameter values R2 = R0 = 4.0 and R1 = 0.5. The start of restrictions is at time t 1 = 3.5 , the end of the restrictions is at time t 2 = 10 .
Asi 03 00044 g004
Figure 5. The evolution of the infected fraction v = I N , as follows from the double transition of R 0 = R 0 ( t ) defined in Figure 4. In this case the 2nd wave is much larger than the first one.
Figure 5. The evolution of the infected fraction v = I N , as follows from the double transition of R 0 = R 0 ( t ) defined in Figure 4. In this case the 2nd wave is much larger than the first one.
Asi 03 00044 g005
Figure 6. Double transition R 0 = R 0 ( t ) with the parameter values R 2 = R 0   = 4.0 and R 1 = 0.5 . The start of restrictions is at time t 1 = 4 , the end of the restrictions is at time t 2 = 10 .
Figure 6. Double transition R 0 = R 0 ( t ) with the parameter values R 2 = R 0   = 4.0 and R 1 = 0.5 . The start of restrictions is at time t 1 = 4 , the end of the restrictions is at time t 2 = 10 .
Asi 03 00044 g006
Figure 7. The evolution of the infected fraction v = I N , as follows from the double transition of R 0 = R 0 ( t ) defined in Figure 6. In this case the 2nd wave is much smaller than the first one.
Figure 7. The evolution of the infected fraction v = I N , as follows from the double transition of R 0 = R 0 ( t ) defined in Figure 6. In this case the 2nd wave is much smaller than the first one.
Asi 03 00044 g007
Figure 8. Fitted R 0 ( t ) for Italy/time unit represents 13 days (case A). The timescale starts on 22 February 2020.
Figure 8. Fitted R 0 ( t ) for Italy/time unit represents 13 days (case A). The timescale starts on 22 February 2020.
Asi 03 00044 g008
Figure 9. Fitted R 0 ( t ) for Romania/time unit represents 14 days (case B). The timescale starts on 1 March 2020.
Figure 9. Fitted R 0 ( t ) for Romania/time unit represents 14 days (case B). The timescale starts on 1 March 2020.
Asi 03 00044 g009
Figure 10. Relation between the 1st and 2nd wave in Italy, if R 2 = 1.6 . The overall picture. One time unit corresponds to 13 days. The timescale starts on 15 February 2020.
Figure 10. Relation between the 1st and 2nd wave in Italy, if R 2 = 1.6 . The overall picture. One time unit corresponds to 13 days. The timescale starts on 15 February 2020.
Asi 03 00044 g010
Figure 11. Italy: first wave of magnitude 1.79‰. One time unit corresponds to 13 days. The timescale starts on 22 February 2020.
Figure 11. Italy: first wave of magnitude 1.79‰. One time unit corresponds to 13 days. The timescale starts on 22 February 2020.
Asi 03 00044 g011
Figure 12. Italy: estimated second wave of magnitude 7.99%. One time unit corresponds to 13 days. The timescale starts on 22 February 2020.
Figure 12. Italy: estimated second wave of magnitude 7.99%. One time unit corresponds to 13 days. The timescale starts on 22 February 2020.
Asi 03 00044 g012
Figure 13. Italy: fitted linear dependence v 2 * v 1 * 97.31 + 89.01 · R 2 . One time unit corresponds to 13 days.
Figure 13. Italy: fitted linear dependence v 2 * v 1 * 97.31 + 89.01 · R 2 . One time unit corresponds to 13 days.
Asi 03 00044 g013
Figure 14. Italy: fitted exponential dependence T 2 ( R 2 ) 4.77 + 242.35 e 2.81 · R 2 .
Figure 14. Italy: fitted exponential dependence T 2 ( R 2 ) 4.77 + 242.35 e 2.81 · R 2 .
Asi 03 00044 g014
Figure 15. Italy: fitted linear dependence v 2 * v 1 * 44.88 0.10 · Δ t . One time unit corresponds to 13 days.
Figure 15. Italy: fitted linear dependence v 2 * v 1 * 44.88 0.10 · Δ t . One time unit corresponds to 13 days.
Asi 03 00044 g015
Figure 16. Italy: fitted linear dependence T 2 = 3.66 + 0.88 · Δ t . One time unit corresponds to 13 days.
Figure 16. Italy: fitted linear dependence T 2 = 3.66 + 0.88 · Δ t . One time unit corresponds to 13 days.
Asi 03 00044 g016
Figure 17. Relation between 1st and 2nd wave in Romania, if R 2 = 1.6 . One time unit corresponds to 14 days. The timescale starts on 1 March 2020.
Figure 17. Relation between 1st and 2nd wave in Romania, if R 2 = 1.6 . One time unit corresponds to 14 days. The timescale starts on 1 March 2020.
Asi 03 00044 g017
Figure 18. Romania: fitted first wave of magnitude 3.9 ‰. One time unit corresponds to 14 days. The timescale starts on 1 March 2020.
Figure 18. Romania: fitted first wave of magnitude 3.9 ‰. One time unit corresponds to 14 days. The timescale starts on 1 March 2020.
Asi 03 00044 g018
Figure 19. Romania: estimated second wave of magnitude 7.8%. One time unit corresponds to 14 days. The timescale starts on 1 March 2020.
Figure 19. Romania: estimated second wave of magnitude 7.8%. One time unit corresponds to 14 days. The timescale starts on 1 March 2020.
Asi 03 00044 g019
Figure 20. Romania: extent fraction vs. R 2 fitted to linear dependence v 2 * v 1 * 48.10 + 42.87 · R 2 .
Figure 20. Romania: extent fraction vs. R 2 fitted to linear dependence v 2 * v 1 * 48.10 + 42.87 · R 2 .
Asi 03 00044 g020
Figure 21. Romania: exponential dependence T 2 ( R 2 )   = a + b e c · R 2 . Fitted parameters: a 4.41 , b 47.63 and c 1.85 .
Figure 21. Romania: exponential dependence T 2 ( R 2 )   = a + b e c · R 2 . Fitted parameters: a 4.41 , b 47.63 and c 1.85 .
Asi 03 00044 g021
Figure 22. Romania: fitted linear dependence v 2 * v 1 * 30.50 0.30 · Δ t ( Δ t in units of τ = 14 days).
Figure 22. Romania: fitted linear dependence v 2 * v 1 * 30.50 0.30 · Δ t ( Δ t in units of τ = 14 days).
Asi 03 00044 g022
Figure 23. Romania: fitted linear dependence T 2 3.76 + 0.60 · Δ t ( Δ t in units of τ =14 days).
Figure 23. Romania: fitted linear dependence T 2 3.76 + 0.60 · Δ t ( Δ t in units of τ =14 days).
Asi 03 00044 g023
Figure 24. Modified R 0 = R 0 ( t ) . After some time (approximatively 2 units) R 2 is reducing to R 3 < R 2 (blue line).
Figure 24. Modified R 0 = R 0 ( t ) . After some time (approximatively 2 units) R 2 is reducing to R 3 < R 2 (blue line).
Asi 03 00044 g024
Figure 25. The response to modified R 0 ; according to reducing R 2 to R 3 < R 2 (blue line).
Figure 25. The response to modified R 0 ; according to reducing R 2 to R 3 < R 2 (blue line).
Asi 03 00044 g025
Figure 26. The response to modified R 0 = R 0 ( t ) ; according to reducing R 2 to R 3 < R 2 (blue line).
Figure 26. The response to modified R 0 = R 0 ( t ) ; according to reducing R 2 to R 3 < R 2 (blue line).
Asi 03 00044 g026
Figure 27. Modified R 0 = R 0 ( t ) . After some time (approximatively 4 units) R 2 is reducing to R 3 < R 2 (blue line).
Figure 27. Modified R 0 = R 0 ( t ) . After some time (approximatively 4 units) R 2 is reducing to R 3 < R 2 (blue line).
Asi 03 00044 g027
Figure 28. Romania: the estimated 3rd wave, assuming the jump after the 2nd wave from R 1 to R 2 remains the same (i.e., from 0.83 to 1.6). Compare this result with the extent of the 2nd wave in Figure 19.
Figure 28. Romania: the estimated 3rd wave, assuming the jump after the 2nd wave from R 1 to R 2 remains the same (i.e., from 0.83 to 1.6). Compare this result with the extent of the 2nd wave in Figure 19.
Asi 03 00044 g028
Table 1. Statistics table for coefficients.
Table 1. Statistics table for coefficients.
EstimateStandard Errort-Statisticp-Value
1−97.31364.16505−23.36432.67331 × 10−6
x89.00632.0659943.08171.27156 × 10−7
Table 2. Statistics table for coefficients.
Table 2. Statistics table for coefficients.
EstimateStandard Errort-Statisticp-Value
144.87990.0871619514.9035.24398 × 10−13
x−0.09911840.015565−6.368010.00141214
Table 3. Statistics table for coefficients.
Table 3. Statistics table for coefficients.
EstimateStandard Errort-Statisticp-Value
13.655220.044430382.26875.02856 × 10−9
x0.8801830.00793419110.9351.12869 × 10−9
Table 4. Statistics table for coefficients.
Table 4. Statistics table for coefficients.
EstimateStandard Errort-Statisticp-Value
1−48.09511.27627−37.68412.33102 × 10−8
x42.86830.66791564.18229.61961 × 10−10
Table 5. Statistics table for coefficients.
Table 5. Statistics table for coefficients.
EstimateStandard Errort-Statisticp-Value
130.49030.0560693543.7962.61013 × 10−15
x−0.3030470.011882−25.50482.39391 × 10−7
Table 6. Statistics table for coefficients.
Table 6. Statistics table for coefficients.
EstimateStandard Errort-Statisticp-Value
13.758780.0223477168.1962.9797 × 10−12
x0.5949210.00473583125.6211.7159 × 10−11
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chiruţă, C.; Bulgariu, E.; Avsec, J.; Ferčec, B.; Mencinger, M. Comparison of the Evolution of the COVID-19 Disease between Romania and Italy. Appl. Syst. Innov. 2020, 3, 44. https://0-doi-org.brum.beds.ac.uk/10.3390/asi3040044

AMA Style

Chiruţă C, Bulgariu E, Avsec J, Ferčec B, Mencinger M. Comparison of the Evolution of the COVID-19 Disease between Romania and Italy. Applied System Innovation. 2020; 3(4):44. https://0-doi-org.brum.beds.ac.uk/10.3390/asi3040044

Chicago/Turabian Style

Chiruţă, Ciprian, Emilian Bulgariu, Jurij Avsec, Brigita Ferčec, and Matej Mencinger. 2020. "Comparison of the Evolution of the COVID-19 Disease between Romania and Italy" Applied System Innovation 3, no. 4: 44. https://0-doi-org.brum.beds.ac.uk/10.3390/asi3040044

Article Metrics

Back to TopTop