key: cord-1047545-catapr59 authors: Baruah, H. K. title: The COVID-19 Spread Patterns in Italy and India: A Comparison of the Current Situations date: 2020-06-23 journal: nan DOI: 10.1101/2020.06.21.20136630 sha: 8e4763e762dbd4a0090ae32791b78a516e899ddd doc_id: 1047545 cord_uid: catapr59 Epidemiological mathematical models and time series models can be used to forecast about the spread of an infectious disease. In this article, without using such models, we are going to show how exactly the pattern evolves day by day once a pattern is seen to be approximately followed by the data. Although in Italy as well as in India the novel corona virus appeared on the same day, in Italy the spread is nearly logarithmic by now and in India it is nearly exponential even now. Various epidemiological models to study the spread of infectious diseases are available in the literature. Meyers [1] studied the use of the Susceptible-Infected-Recovered (SIR) model for infectious diseases affecting a small region. Kamp [2] concluded that the transmission network strongly influences speed and range of spread of epidemics. Liu et. al. [3] studied theoretically the Susceptible-Infectious-Susceptible (SIS) pandemic model using time delays in the thresholds. Majeed and Shawka [4] studied mathematically the SIS model of epidemic growth. Wu et. al. [5] used the Susceptible-Exposed-Infectious-Recovered (SEIR) model in a simulation study of the COVID-19 spread. Recently Grant [6] has shown that SEIR model does not work properly in the case of the spread of COVID-19. There are certain other models too such as the Susceptible-Infectious-Recovered-Dead (SIRD) model [7] . Indeed, quite a few more epidemiological models are in use to study the spread of infectious diseases. It has been seen that some models do forecast about the spread very nicely while some others have not been very successful. Time series models using the auto-regressive integrated moving average (ARIMA) method have also been used successfully by a few authors for forecasting the COVID-19 spread in India. Poonia and Azad [8] and Azad and Poonia [9] studied the situation in two phases using the ARIMA method. Their forecasts were close to the actual values observed later. Basu [10] studied time dependent spread of the virus in India using his own model. As per his forecasts expressed in the form of a graph, the total number of cases in India should have crossed 200,000 in the beginning of June, and his forecast has been found to be true. In this work, we have two objectives. First, we shall compare the current spread patterns in Italy and India. From the graphs of the current data on total number of cases that includes active cases, recovered cases and deaths, some mathematical pattern would be apparent. We would find out the current approximate patterns of the spread in these two countries. Secondly, we would like to show how the patterns evolve day by day. Hypothesizing about a particular mathematical model as the underlying spread pattern is one thing, while studying the changes day by day looking into the recent data is quite another. We would show how the logarithmic function is being followed by the total number of cases in Italy, and how in India it is following a nearly exponential function, while the patterns are changing slowly and steadily. We would like to state that while trying to find a mathematical model regarding the spread pattern of a pandemic, the epidemiological models assume that in the first phase of the spread the mathematical pattern is exponential. In fact, this assumption is a very simplified one. Initially, the pattern would be nonlinear but assuming it to be exponential right from the beginning is not actually valid. Indeed, it can in reality be only approximately exponential, because if the pattern is exponential, it would only mean that the pattern would remain so until everything is finished. When the pandemic happens to retard, it is said that the curve is getting flattened. It is important to examine following what pattern the growth starts retarding. In any phase, we cannot say that a particular growth pattern is being strictly followed. When the pandemic would finally come to an end, it would mean that the pattern is a straight line parallel to the time axis. Before that happens, the pattern may be approximately logarithmic. Indeed a logarithmic growth is also a growth, however slow it is. We are in this article going to discuss and compare the current growth patterns of COVID-19 in Italy and India. We would show that in Italy the current growth is approximately logarithmic while in India it is still growing approximately exponentially, although as per data published by Worldomters.info [11], on February 15, 2020, both in Italy and in India the total numbers of cases were equal to 3. where m and c are constants. If we need to fit such a curve from some observed data, we have to estimate the parameters from the data. That would need a method such as the method of least squares to estimate the two parameters. However, in the case that we are currently going to discuss, we would assume that from a certain value of x the concerned curve is logarithmic. Therefore taking that starting value of x as 1, if we want to verify whether the curve is approximately logarithmic, then the value of the constant c is already known to us, and therefore in that kind of a situation we would need to estimate only the value of m. If we now see that the estimated values of m found from . CC-BY-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted June 23, 2020. Although there are two parameters involved, when we proceed to apply this model in the exponential phase of a pandemic, we can take some particular date as the base to describe a few dates thereafter with reference to the base date. Accordingly, the value of the constant a would be available to us already and we would have to find an estimate of the parameter b. When we would observe that the estimated values of b are very nearly constant, we would be able to say that the pattern is approximately exponential. As per the Worldometers.info data [11], on February 15 there were 3 cases reported in Italy. From the graphs it can be seen that by the middle of March the spread pattern became highly nonlinear. By the end of April, the process of curve flattening had started. In other words, the pattern changed by the end of April from an approximately exponential form to a logarithmic form. We are explaining the matters in a rather simple way. But it can actually be checked whether it was really so. In Table- . CC-BY-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted June 23, 2020. . CC-BY-ND 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted June 23, 2020. . https://doi.org/10.1101/2020.06.21.20136630 doi: medRxiv preprint Contact network epidemiology bond percolation applied to infectious disease prediction and control Demographic and behavioural changes during epidemics SIS model on homogeneous networks with threshold type delayed contact reduction The dynamics of an eco-epidemiological model with (SI), (SIS) epidemic disease in prey Nowcasting and forecasting the potential domestic and international spread of the 2019 nCoV outbreak originating in Wuhan, China: A modeling study Dynamics of COVID-19 epidemics: SEIR models underestimate peak infection rates and overestimate epidemic duration Data-based analysis, modeling and forecasting of the COVID-19 outbreak Short term forecasts of COVID-19 spread across Indian States until Short term forecasts of COVID-19 spread across Indian States until 29 May, 2020, under the worst case scenario fits the data of spread in India [12] approximately.To estimate the value of the parameter b at some point of time we would need data about the total number of cases for a few days prior to that. If the values of , the first order differences of z are seen to be nearly constant, then we can say that the pattern is nearly exponential. It was observed in [12] that the values of ∆ ‫ݖ‬ have been following a reducing trend with some irregularities in between which is inherent in the case of a time series of this type of a pandemic. In Table- This shows that at least up to June 17 the pattern of spread in India is nearly exponential. Our objective was not to forecast, but to show that the nearly exponential pattern is still continuing in India. Indeed, our method is valid for short term forecasting only. As has been mentioned earlier, the values of ∆ ‫ݖ‬ were seen to be in a reducing trend. So this method is unsuitable for long term forecasting. Although the pandemic had started in Italy and India on the same date, February 15, 2020, the situations soon became very different in the two countries. In Italy, the nearly exponential pattern could be observed before the middle of April, and currently it is increasing logarithmically. In India however, in April the spread was just highly nonlinear and currently it is nearly exponential. This shows that a study regarding the total number of cases in the world as a whole cannot follow one single mathematical model, because whereas in India the spread pattern is continuing to be nearly exponential, in Italy it is the inverse function -the logarithmic functionbeing followed by the data.