## Abstract

The first case of COVID-19 was confirmed in Israel on February 21, 2020. Within approximately 30 days, the total number of confirmed cases climbed up to 1,000, accompanied by a doubling period of less than 3 days. About one week later, after this number exceeded 4, 000 cases, and following some extreme lockdown measures taken by the Israeli government, the daily infection rate started a sharp decrease from the peak value of 1,131 down to slightly more than 100 new confirmed cases on April 30. Motivated by this encouraging data, similar to the trends observed in many other countries, along with the growing economic pressures, the Israeli government has quickly lifted most of its emergency regulations. Throughout May, the daily number of new cases stayed at a very low level of 20 – 40 until at the end of May it started a steady increase, exceeding 1, 000 by the end of June and 2,000 on July 22. As suggested by some experts and popular media, this disturbing trend may be even a part of a “second wave”. This article attempts to analyze the currently available data on Israel, compared to three European countries (Greece, Italy, and Sweden), in order to understand the local dynamics of COVID-19, assess the effect of the implemented intervention measures, and discuss some plausible scenarios for the foreseeable future.

## Introduction

The first case of COVID-19 was confirmed in Israel on February 21, 2020. On March 9, a two-week self-isolation was imposed on all people coming from abroad and on March 12, schools and universities were closed by the government order, partially switching to distant teaching. Then, on March 19, when the number of daily confirmed cases exceeded 100, all non-essential businesses were ordered to close, the employees were required to increase social distancing of their workers and, if possible, allow them to work from home, whereas people’s movement outside their homes was restricted significantly. Still, by March 23, the number of confirmed cases climbed up to 1, 000, accompanied by a doubling period of less than 3 days. A continuous decrease in the growth rate has started only around March 29, after the total number of confirmed cases exceeded 4,000.

With the peak unemployment rate of 24%, mostly due to people placed on unpaid leave as a result of the pandemic, the Israeli government decided to partially lift the lockdown, which has been maintained throughout the Passover holiday (April 7-16). The “exit strategy” steps under consideration included increasing the allowed percentage of workers at workplaces, opening some shops, partially resuming public transportation, gradually getting children back to school, and selectively releasing people from the movement restrictions based on the level of infection in their neighborhood and their risk group. A list of detailed measures for easing the current restrictions was approved by the Israeli government on May 04, 2020. By April 30, the doubling period went up to several weeks with only 222 death outcomes out of 15, 946 confirmed cases. However, at the end of May, a few weeks after the government lifted most of its restrictions, the daily number of new cases started a new and prolonged increase, approaching 2, 000 around July 22.

Table 1 attempts to compare COVID-19 evolution in Israel and three European countries: Sweden, Greece, and Italy. Each of these countries has some common characteristics with Israel along with some prominent differences. In terms of population size, Sweden and Greece are similar to Israel, whereas Italy is several times bigger. On the other hand, the median age in Israel is much lower than in the other three countries, indicating a higher percentage of young population. The weather conditions in Israel are quite similar to Greece, slightly warmer than in Northern Italy, and radically different from Sweden. Also, Sweden, as a typical Scandinavian country, is normally characterised by a greater social distancing than such Mediterranean countries as Italy, Greece, and Israel.

The first multiple COVID-19 cases were detected in Israel, Sweden, and Greece around February 26, about three weeks later than in Italy. In all four countries, the number of confirmed cases has reached 100 within 1.5 – 2.5 weeks, with a first death case reported 2 - 3 weeks from the apparent start of the epidemic. However, none of the above mentioned characteristics of each country (population size, climate, social distancing culture, age distribution, etc.) can explain such extreme differences in COVID-related mortality rates: 11 times more in Italy and Sweden than in Israel, and 2.5 times higher in Israel than in Greece.

The ongoing discussions of alternative “coronavirus strategies” taken by different governments brings up the following questions:

What is the basic reproduction number

*R*0 in Israel vs. other countries and how is it evolving over time?Which measures, if any, have caused the steep decrease in the daily number of new confirmed cases in Israel, which has been observed during April–May 2020, and what can be done to reverse the current growth in the infection rate ?

What is the true Infection Rate (IR) in Israel, i.e. what is the actual fraction of infected people including those who were not tested for coronavirus during their detectable period ?

What are the Case Fatality Rate (CFR) and the Infection Fatality Rate (IFR), defined as the ratio of COVID-19 attributed deaths to the number of confirmed cases and the total number of infection cases, respectively? Are they changing over time?

Finally, the most significant question: What is the expected evolution of the epidemic in Israel as a function of alternative intervention strategies?

The subsequent sections will focus on each one of these crucial questions. Finally, we will discuss the implications of our data analysis on the prospective decisions to be taken by the government in Israel and other countries.

## Materials and methods

### Data Sources

For daily case count data in Israel and other countries, we rely upon the Novel Coronavirus (COVID-19) Cases Data, which is part of Humanitarian Data Exchange (HDX) platform [1]. The data is compiled since 22 January 2020 by the Johns Hopkins University Center for Systems Science and Engineering (JHU CCSE) from various sources such as the World Health Organization (WHO) and European Centre for Disease Prevention and Control (ECDC). The available fields include Province/State, Country/Region, Last Update, Confirmed, Suspected, Recovered, and Deaths.

The regional mortality data in Italy was downloaded from the website of the Italian Department of Civil Protection [2] on July 23, 2020. The available fields for each date and region include (translated from Italian): hospitalized with symptoms, intensive care, total hospitalized, isolated at home, total positive, total positive variation, new positive, discharged recovered, deceased, suspected cases, screening cases, total cases, swabs, cases tested.

The coronavirus geographical data in Israel was downloaded from the COVID-19 Data Repository of the Israeli Ministry of Health [3] on July 23, 2020. The available fields included: town_code, agas_code, date accumulated_tested, new_tested_on_date, accumulated_cases, new_cases_on_date, accumulated_recoveries, new_recoveries_on_date, accumulated_hospitalized, new_hospitalized_on_date, accumulated_deaths, new_deaths_on_date, town.

### Modeling the Infection Dynamics

Following [6], we assume that the COVID-19 outbreak can be represented by SIR dynamics, which assumes that at any given point of time, each individual in the population belongs to one of the following three states: (I)nfected, (S)usceptible to infection, or (R)emoved from the transmission process. However, the original SIR model and its variations, like SEIR (susceptible-exposed-infected-recovered) [5], are continuous time models based on a set of differential equations. In contrast, we represent the COVID-19 dynamics by a discrete time model, which is more appropriate for the time series data of daily case counts. Fixing the total size of the population under a discrete time SIR model implies that the sum of the daily changes in the amount of individuals belonging to each one of the above three compartments should be equal to zero:
where *I _{t}, S_{t}*, and

*H*stand for the changes in the total amount of infected, susceptible, and healed (’’recovered”) people, respectively, between the days

_{t}*t*− 1 and

*t*. Our discrete time model for the temporal evolution of the number of infected individuals builds upon the basic reproduction number

*R*

_{0}, which represents the average number of secondary infections an infected person will cause in a completely susceptible population before he or she is effectively removed from the population as a result of recovery, hospitalization, quarantine, etc. [6]. At the time of a pandemic, the value of

*R*

_{0}may be temporarily reduced by decreasing the amount of social interaction (a measure known as “social distancing”). Additional parameters used by our model include the average duration of the incubation period

*L*between the exposure and the onset of clinical symptoms and the generation time [

_{I}*L*−

_{I}*max_T*;

_{I}*L*] between the primary case and the secondary case, where

_{I}*max_T*represents the maximum number of days from the start of the infectious period to the onset of the symptoms. We assume that symptomatic individuals are not infectious anymore as they are immediately isolated from the susceptible population and required to take a COVID-19 test. We also assume that asymptomatic people are never referred to a test and that the social behavior of infected individuals does not change during their infectious period. Given the initial size of the susceptible population

_{I}*N*, we can define a discrete time model of infection dynamics using Eq 2.

Eq 4 shows the condition that the number of infected people continues to grow exponentially. From there, using Eq 5 we can extract the herd immunity threshold, which is the minimum cumulative amount of infected individuals that should stop an epidemic under a fixed value of *R*_{0}.

Given a time series of estimated daily infections *I _{i}*, the incubation period

*L*, the generation time range [

_{I}*L*−

_{I}*max_T*;

_{I}*L*], and the population size

_{I}*N*, we can use a stochastic optimization method to find the value of

*R*

_{0}, which should reconstruct that series with a minimal average absolute error. The error can be further reduced by splitting the time series into several segments and estimating

*R*

_{0}for each segment separately.

To monitor the evolution of *R*_{0} on daily basis, we can take several simplifying assumptions. First, we assume that as long as only a small proportion of a country’s population is infected, the number of susceptible individuals *S _{t}* remains close to the total population size

*N*. We also assume the generation time to be fixed to the midpoint day of the interval [

*L*−

_{I}*max*_

*T*;

_{I}*L*]. Applying these assumptions to Eq 2, we can approximate the expected number of new infected cases on day

_{I}*t*as:

Eq (6) is a linear autoregressive process of order *L _{I} − max_T_{I}*/2 + 1 and thus we can use the least squares method to estimate the value of

*R*0 from the daily counts of new infected cases. Since

*R*

_{0}may be affected by multiple factors, such as the level of social distancing, we can monitor its change over time by calculating the moving slope of the autoregression equation (6) over a sliding window of

*n*days.

To project the epidemic dynamics into the future, we just need to choose the expected value of *R*_{0}, initialize the daily amounts of infected people on the last *L _{I}* days before the beginning of the projection period, and start calculating the succeeding values of

*I*recursively using Eq (2). A similar simulation paradigm was implemented in [8] to evaluate several transmission and intervention scenarios for COVID-19 in the United States over the next five years.

_{t}As long as there is no daily testing of the entire population, the exact values of *I _{t}* remain unknown. However, in the next sub-section, we explain how one can estimate the actual infection dynamics from the currently available epidemiological data.

### Estimating the Actual Infection Dynamics from Available Data

We assume that symptomatic individuals take a COVID-19 test on the day of the symptoms onset (the end point of the incubation period), but it takes r days on average to report a positive test result as a “confirmed case”. Given the average incubation period *L _{I}* and the average reporting delay

*r*, the time series of the daily number of new infection cases

*I*is related to the time series of the daily number of confirmed cases

_{t}*C*by Eq (7). where

_{t}*p*is the reporting rate on day t, which is equal to 1 if there is a full and timely detection of all positive cases, based on daily testing of the entire population, and greater than 1 when some under-reporting takes place.

_{t}While the amount of new confirmed COVID-19 cases in various countries, regions and even cities is published on daily basis, its level of under-reporting is unknown and potentially unstable due to inconsistent testing policy, varying reporting time, and other factors and thus we cannot use it as a trustworthy indicator of the epidemic dynamics. On the other hand, the daily number of virus-related deaths *D _{t}* is expected to be much more reliable than the testing results and it should not suffer from any significant reporting delay. Assuming a fixed Infection Fatality Rate

*IFR*, we can estimate the daily amount of infection cases on day

*t*as a function of the amount of death outcomes on day

*t*+

*L*+

_{I}*d*, where d is the average number of days between the onset of the virus symptoms and the patient death (see Eq 8). Hence, disregarding the true value of

*IFR*, we can explore the dynamics of the daily death rate

*D*as a time-shifted substitute for the true infection rate

_{t}*IR*(Eq 9). Given the average reporting delay

_{t}*r*, we can also estimate the Case Fatality Rate (CFR), which is defined in epidemiology as the proportion of people who die from a specified disease among all individuals diagnosed with the disease over a certain period of time [7]. Since the CFR numerator is restricted to deaths among people included in the denominator, we calculate it as the average ratio between the amount of death outcomes on day

*t*+

*d*−

*r*and the number of confirmed cases on day

*t*(Eq 10).

## Results

### The Epidemic Dynamics

Our analysis of COVID-19 evolution is based on the following timing parameters provided in [9,10]:

Mean incubation period

*L*= 5.5 days_{I}Minimum time from the start of the infectious period to the onset of the symptoms

*min*_*T*= 1 day_{I}Maximum time from the start of the infectious period to the onset of the symptoms

*max_T*= 3 days_{I}Mean generation period (average number of days it takes a patient to start infecting others)

*L*/2 = 4 days_{I}− max_T_{I}Mean time between the onset of the virus symptoms and the patient death

*d*=11 daysMean time from infection to death

*L*= 16.5 days_{I}+ dMean time from symptoms to detection

*r*= 4.5 daysMean time from infection to detection

*L*+_{I}*r*= 10 daysMean time from detection to death

*d*−*r*= 6.5 days

Fig 1 shows the overall evolution of the reproduction number *R*_{0} in Israel between March 27 and July 5, 2020 along with the 14-day moving average of the daily number of death cases between March 22 and July 18, 2020. On each day, *R*_{0} was estimated as a moving slope of a linear autoregressive model over a 14-day sliding window of daily death cases^{1} using an autoregression lag of 4 days, which is equal to the mean generation interval, and assuming that the number of deaths on day d represents the number of new infections on day *d* − 16. The chart also shows the main intervention measures taken by the Israeli government during the same period. Since we did not have enough mortality data to estimate the infection rate before March 27, we could not evaluate the direct effect of the school closing, which occurred around March 12, as well as business activity and movement restrictions imposed on March 19. Further measures, such as mandating masks in public (April 1) and Passover lockdown (April 7) were actually introduced when the value of *R*_{0} has already decreased below the critical level of 1.0, a trend which continued throughout the first half of May despite a fast relaxation of most restrictions enforced in March and April. As *R*_{0} went below the value of 1.0, the doubling rate (number of days required to double the number of infected people) climbed highly above the 10 days threshold, which was defined in the Government decision from May 04, 2020 as one of the criteria for rolling back the “exit strategy” steps.

After reaching an extremely low level of 0.3 around May 17 and shortly after reopening of all schools in the country (which stayed open until the end of the school year on June 20), *R*_{0} started to climb up again in a continuous trend, which included four distinct peaks on May 23 (1.1), June 12 (1.8), June 22 (1.3), and July 2 (1.4). The moving average of the daily mortality rate clearly followed the evolution of *R*_{0} with a lag of several weeks, peaking at 9.0 per day on April 19, going down to 0.6 at the end of June, and rising up to 5.0 in the second half of July. However, Fig 1 does not reveal a consistent effect of any specific measures on the actual infection rate represented by the reproduction number.

Fig 2 shows the local evolution of the reproduction number *R*_{0} in four most populous Israeli cities (Jerusalem, Tel-Aviv-Yaffo, Haifa, and Rishon-Le-Zion) along with the city of Bnei-Brak, which, for various reasons, has accumulated the second largest number of confirmed COVID-19 cases in Israel despite being only no. 9 in population size. Due to relatively low daily mortality rates in each city, R_{0} was calculated over a 14-day sliding window of daily *confirmed* cases, smoothed with a five-day moving average. It is noteworthy that between the end of March and the end of May, the reproduction numbers in all five cities followed more or less the same decreasing trend. However, the months of June and July were characterised by distinct local peaks (brief but significant outbreaks), which can be easily related to fluctuations in *R*_{0} at the country level that appear in Fig 1: May 23 - an outbreak in Jerusalem and Tel-Aviv, June 12 - an outbreak in Haifa and Bnei-Brak, June 22 - another outbreak in Jerusalem, and July 2 - an outbreak in Rishon-Le-Zion.

Fig 3 shows the overall evolution of the reproduction number *R*_{0} in Greece, along with the 14-day moving average of the daily number of death cases, for the period between March 17 and June 22, 2020. Similar to Israel, the reproduction number in Greece quickly went from 1.2 down to 0.7, in parallel to domestic travel restrictions and various other lockdown measures. Despite several brief peaks in the infection rate, the overall decreasing trend continued throughout April, May, and June (when thelockdown restrictions were gradually removed), with *R*_{0} staying well below the value of 1.0 most of the time and less than one death case per day on average.

As shown by Fig 4, the evolution of COVID-19 in Italy followed a pattern similar to Greece but on a completely different scale. The initial average value of the reproduction number *R*_{0} in Italy was 2.5, about two times higher than in Israel and Greece, and it took it slightly more than one month to go down below 1.0, whereas the average daily number of deaths approached 800 in the beginning of April. Unlike in Israel, the lifting of lockdown restrictions in Italy was not followed by a significant increase in the value of *R*_{0}, which stayed below 1.0 from mid-March until the end of July. At the same time, the daily mortality rate in Italy faced a steady descent from 700 – 800 in March to nearly 10 in July. Fig 5 indicates that COVID-19 dynamics followed a very similar trend in all six Italian regions that had the highest number of death cases.

Fig 6 shows the 14-day moving average of *R*_{0} in Sweden from March 17 to July 19, 2020. The chart also shows the main intervention measures taken by the Swedish government during the same period. On March 17, the mean reproduction number in Sweden was at its peak value of about 2.0, much higher than its maximum values in Israel and Greece but lower than in Italy. Within the next two weeks it went down below 1.0, while the Swedish Government introduced various restrictions on social distancing rather than enforcing a total lockdown like in Israel, Italy, Greece, and many other countries. During the following months, the Swedish reproduction number fluctuated between 0.6 and 0.8, eventually converging to 0.5 in the second half of June, a much lower value than in Israel during the same period (0.8 − 1.2). While it is evident that Sweden did a better job than Israel in terms of “smoothing the curve”, it experienced a relatively high mortality rate, with a peak number of about 100 cases per day at the end of April, which eventually went down below 15 in the second half of May.

Despite the differences in the demographic and cultural characteristics of the four countries, along with different intervention policies taken by their governments, all of them experienced quite similar dynamics of COVID-19. The initial reproduction number varied between 1.2 and 2.5 and within the following few weeks, it rapidly decreased below 1.0, disregarding the specific social distancing measures introduced during that time. Subsequently, *R*_{0} stayed below 1.0 for at least 1-2 months, without being significantly affected by various lockdown relief steps. A repeated increase in *R*_{0} observed in Israel at a later stage requires further investigation, which may have implications for other countries as well.

The true Infection Rate (IR) of COVID-19, like any other epidemic, is a non-decreasing function of time. Hopefully, it will be estimated in the near future using massive serological tests but its true past values will remain unknown forever. However, we believe, like the authors of [11], that the Test Positivity Rate (TPR) during the first few weeks of the epidemic, calculated as the fraction of positive to total tests, should be an upper bound on the true IR. Fig 7 shows the 7-day moving averages for TPR and for the daily amount of COVID-19 tests performed in Israel between March 11 and April 25. Since the mean TPR reached its peak value of 10% around March 25, when the infection rate was also at its maximum, we can conclude that the true IR in Israeli population also did not exceed 10% at that time.

The results of initial 1,700 serological tests performed recently in Israel [12] suggest that the true infection rate in Israel is at least 10 times higher than the number of confirmed cases. Assuming that this reporting ratio is preserved and having 56, 000 confirmed cases up to July 22, we may estimate that by that date about 6.5% of the country’s population of nearly 8.9 million have been infected with the virus. Using Eq 5 and the initial reproduction number of *R*_{0} = 1.136, which was observed at the beginning of the epidemic, we can predict that the daily number of new cases will start a steady decline after the virus prevalence will exceed 12%.

A similar antibody test conducted between April 23 and June 3 on nearly 10, 000 residents of an Italian city of Bergamo has shown about 57% of positive results [13]. Considering the total amount of 13, 600 cases confirmed in the city of 122k people around the same time, the estimated reporting ratio in Bergamo stands at about 5. Using Eq 5 and the initial reproduction number of *R*_{0} = 3.014, which was observed at the beginning of the epidemic in Lombardia Region, we can predict that the daily number of new cases in Bergamo will start a steady decline after the virus prevalence will exceed 67%, just slightly higher than 61% based on 14, 872 cases confirmed until July 22.

Figs 8, 9, 10, and 11 show the calculation of the Case Fatality Ratio (CFR) in Israel, Greece, Italy, and Sweden, respectively, as a slope of a linear regression model between the cumulative amount of confirmed cases and the cumulative number of death outcomes 7 days later. However, in all four countries, we can observe a certain decrease in the line slope, i.e., in CFR, about 2-3 months from the start of the epidemic: from 1.8% to 0.5% in Israel, from 5.9% to 1.2% in Greece, from 16% to 13% in Italy, and from 14% to 2.9% in Sweden. The high values of R-Square in all regression models indicate that most COVID-19 deaths are related to previously confirmed cases. This result supports the hypothesis that most symptomatic patients are detected by the lab tests. However, we do not have and may never have enough data to test this claim in either country. We also have no plausible explanation for the exceptionally high CFR values in Italy compared to Israel, Greece, and Sweden and to the world average of 3.6% in general. The reasons for the significant gaps between mortality rates in different countries and regions deserve a separate study as they may include differences in lab testing strategies, clinical treatment practices, outcome classification policies, coronavirus strains, etc.

### Projecting into the Future

Based on the infection data estimated up to July 22, 2020 and assuming the reporting rate of *p* = 10 (according to the initial serology test conducted in Israel), Figs 12, 13, and 14 simulate the projected number of daily confirmed, cumulative confirmed, and cumulative death cases, respectively, for several possible values of the reproduction number *R*0. Our simulation covers the period from the end of July to the end of December 2020, when we may be close to the peak of a new flu season in Israel. These projected numbers allow us to forecast three important parameters: maximum expected number of critical patients, total amount of death outcomes, and the start of a decline in the amount of new daily infections as a result of crossing the herd immunity threshold (calculated using Eq 5). Following [10], we assume that the average hospital stay of a critical patient is equal to 20 days and that the average proportion of critical patients in Israel is about 1% of confirmed cases. Consequently, the maximum number of available critical care beds should be equal to 20% of the maximum number of daily infections. The number of critical COVID-19 patients at the end of July 2020 is very close to the above estimator based on the moving average of new daily infections detected during the last two weeks of July. The total number of death outcomes is estimated as a product of the current Case Fatality Ratio (0.5%) by the cumulative number of cases confirmed 7 days earlier.

Here is a brief discussion of each simulated scenario:

*R*0 = 1.00. This value means that the daily infection rate has started a slow decline towards less than 100 cases by the end of October and less than 10 cases by the end of December, resulting in about 300 additional deaths. This scenario will require about 360 critical care beds at its peak, slightly higher than the actual number occupied at the end of July and much lower than the maximum capacity of 1,000 beds.*R*0 = 1.10. This more severe scenario will increase the daily infection rate up to 1, 900 cases in the beginning of August, until, according to Eq 5, the amount of infected people in the Israeli population will reach the threshold of (1 − 1/1.1) * 100 = 9.1%, equivalent to about 79,000 confirmed cases. Then the infection rate will start a decrease down to 100 per day by the beginning of December. At its peak, this scenario will require close to 400 critical care beds. About 700 more people are expected to die from the virus until the end of 2020.*R*0 = 1.20. This may be the most realistic scenario, since as shown on Fig 1, similar values of the reproduction number were observed in Israel both at the beginning of the first wave, when virtually the entire population was susceptible, and during the first outbreak after the exit from the lockdown, when the population immunity was still very small. The peak is expected at the end of August, when the virus prevalence will reach (1 − 1/1.2) * 100 = 16.7%, equivalent to about 144, 000 confirmed cases. The corresponding daily infection rate of 3, 100 cases will require about 620 critical care beds, still well below the current capacity of the Israeli hospitals (1,000 beds). Then the daily infection rate should go down to less than 100 cases by the beginning of December. Under this scenario, about 1; 100 more people are expected to die from the virus until the end of 2020.

## Conclusion

The future evolution of COVID-19 in Israel and elsewhere is hardly predictable. Until a vaccine becomes available, which may take at least several months from now, the governments will have to strike the tough balance between health and economic issues. As shown by Eq 4, a decrease in the basic reproduction number *R*_{0}, e.g., as a result of a strict lockdown, should reduce the herd immunity threshold in terms of the total amount of infected people. However, the economic and psychological price of maintaining extreme social distancing measures until a massive vaccination program can be launched, may be too prohibitive. On the other hand, as shown by recent experience in Israel and some other countries, lifting a majority of restrictions creates a potential for multiple local outbreaks in places where *R*_{0} exceeds 1/*fraction of susceptible population*. In the absence of an effective vaccine, the preparedness for such “second wave” phenomena is absolutely necessary.

A further progress in understanding the current pandemic will be possible with a release of detailed clinical records of COVID-19 patients in Israel and other countries tothe research community. Analyzing these records with state-of—the-art statistical and machine learning algorithms may reveal answers to many important epidemiological questions such as an accurate and early detection of high-risk patients, identification of the most infectious persons (”superspreaders”) and locations, characterization of the most common infection pathways, etc. Many thousands of human lives worldwide are at risk and, as we all know, “data saves lives”.

## Data Availability

The daily case count data in Israel and other countries is available from the Humanitarian Data Exchange (HDX) platform compiled by the Johns Hopkins University Center for Systems Science and Engineering (JHU CCSE). The coronavirus testing data in Israel is available from the COVID-19 Data Repository of the Israeli Ministry of Health.

## Acknowledgments

I thank the Israeli Ministry of Health, the Gertner Institute, and the Johns Hopkins University Center for Systems Science and Engineering for publishing the COVID-19 data used in this research. I am also grateful to Michal Feldman, Sivan Frenkel, Neil Gandal, Isidore Last, Ady Pauzner, Avraham Tabbach, Matan Yonas (Tel-Aviv University), Yoram Louzoun (Bar-Ilan University), and Lev Muchnik (Hebrew University of Jerusalem) for the fruitful discussions of some of the topics presented in this paper.

## Footnotes

The analysis now includes four countries: Israel, Greece, Italy, and Sweden. All results were updated with the data available by July 2020.

↵1 Smoothed with a three-day moving average