Retrospective study of the first wave of COVID-19 in Spain: analysis of counterfactual scenarios
================================================================================================

* Benjamin Steinegger
* Clara Granell
* Giacomo Rapisardi
* Sergio Gómez
* Joan T. Matamalas
* David Soriano-Paños
* Jesús Gómez-Gardeñes
* Alex Arenas

## Abstract

One of the most important questions on the COVID-19 pandemic is ascertaining the correct timing to introduce non-pharmaceutical interventions (NPIs), based mainly on mobility restrictions, to control the rising of the daily incidence in a specific territory. Here, we make a retrospective analysis of the first wave of the epidemic in Spain and provide a set of useful insights to optimize actions in the near future. We have reconstructed the exposure times, from infection to detectability, to correctly estimate the reproduction number Rt. This enables us to analyze counterfactual scenarios to understand the impact of earlier or later responses, decoupling containment measures from natural immunity. Our results quantify the differences in the number of fatalities for earlier and later responses to the epidemic in Spain.

**Teaser** *“We propose a backward analysis of pandemic incidence in a region to determine the correct timing of authorities’ non-pharmaceutical interventions to fight COVID-19”*

Spain was among the strongest hit countries globally during the first wave with officially 28.000 fatalities caused by COVID-19, in excess deaths even 40% more, and an attack rate between 4% and 5% (1). Surprisingly, the Spanish authorities implemented very stringent containment measures, which contrasted with the high fatality and attack rates. The national government implemented the first lockdown on March 15, which was reinforced on March 29 until April 12. The confinement was then gradually lifted from May 2 onward.

In short, the non-pharmaceutical interventions NPIs issued by the authorities, all together with physical distancing, face masks, and hands hygiene adopted by the population, were sufficient to mitigate the daily infections. Therefore, a natural retrospective question arises: how to identify the impact of the various NPIs implemented at the timing they were implemented. The answer is essential to design the future mitigation strategies of COVID-19, as it has been done for China or Italy (2,3).

To assess the impact of NPIs, one can either directly analyze the epidemiological data such as case numbers, hospitalizations, and fatalities (4) or perform model-based inference (5, 6) to extract cornerstone conclusions. Here, we blend both approaches. Analyzing the epidemiological data allows us to propose a hypothesized functional form for Rt during Spain’s first wave. Then, we combine this information with a model-based approach to infer the behavior of Rt. With this analysis, we can explore counterfactual scenarios.

## Results

One of the most pressing difficulties when analyzing an epidemic’s impact retrospectively is estimating the unobserved past incidence of the infection (7). Estimating the past incidence of COVID-19 relies on the accurate determination of the period’s distribution between the exposure to the virus to affliction (or detectable symptomatic infection).

In Spain, the ISCIII provides a time series with symptom onset dates for reported cases (8). Using current estimates of the incubation period, we can infer the daily infections using a non-parametric back-projection method based on a maximum likelihood estimation (7) (see Materials and methods).

In Fig. 1A we show the results of the estimated infections for the Comunidades Autónomas (CCAA) in Spain, as well as for the whole country, together with the corresponding official data for symptoms onset and reporting date. The substantial delays between infection and report are evident for all regions, as corroborated in Fig. 1B showing the day when the infections and reported cases peaked, andthe difference between them. On average, infections and reported cases peaked with a 15 days difference (sd: 3). We find the maximal difference of 20 days for Catalonia, together with Castilla La-Mancha.

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F1)

Figure 1: Reconstructed infections vs. reported notifications in Spain
(**A**) Time series for reported cases in blue and symptom onset in yellow for Spain and the CCAA. Solid lines indicate a centered, 7 day rolling average. In red we show the reconstructed exposure times. The shaded area in light gray indicates the lockdowns 1 and 3. In dark gray, we indicate the lockdown 2 where, in addition, all non essential economic activity was shut down. (**B**) Date on which the reconstructed exposure (red) times and reported cases (blue) peaked. Additionally, we explicitly indicate the delay between the two peaks. Shaded areas and bars around lines indicate 95% credible intervals in the entire manuscript.

Then, we estimated Rt’s evolution in Fig. 2A using EpiEstim (9). In all CCAA, we observe a substantial reduction of Rt before the first confinement. Close to this first confinement, Rt is for the first time below one in almost all CCAA, and stayed mainly below one until the end of the confinement. Fig. 2B shows the day on which Rt was first below one in the CCAA. From March 13 to 21, Rt became smaller than 1 in all CCAA. On average, Rt crosses 1 on March 16 with a standard deviation of 2,1 days, which is consistent with other estimations (14, 15). The Community of Madrid, was the first CCAA to push Rt below 1 as a consequence of the implemented containment measures on March 9, before any other CCAA. Similar results are found for other CCAA (see SM for detailed information about the evolution across CCAA).

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F2)

Figure 2: Estimation of Rt from reconstructed infections
(**A**) Estimation of Rtfor Spain and the CCAA from the reconstructed exposure times. The shaded area in light gray indicates the lockdowns 1 and 3. In dark gray, we indicate the lockdown 2 when all non-essential economic activity was shut down. (**B**) Chronological order for the date when Rt dropped first below 1 in Spain and the CCAA.

After obtaining the retrospective values of Rt, and given that the NPIs impact starts with the behavior of people, we wonder when the epidemic response was initiated. The epidemic response should result in a breakpoint (17) in Rt after which we observe a substantial reduction of Rt. To determine this breakpoint, R1, we separate (see Materials and Methods) Rt’s evolution into three linear segments: the “free” spreading phase before restrictions, the initiation of the epidemic response, and the lockdown. The result is shown in Fig. 3A. The second breakpoint pretty much coincides with the implementation of the first lockdown by the authorities in Spain. More interestingly, the reproduction number starts already to decrease between March 5 and 6 (first breakpoint). The early decrease precedes the introduction of any containment measures, also on a regional level, as well as the mobility reduction, which started between March 9 and 10 (see Fig. S1).

![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F3)

Figure 3: Early decrease of Rt and saturation of testing capacity
(**A**) Estimated evolution of Rtin Spain. Yellow dashed lines and crosses indicate the identified segments and breakpoints, respectively, in the time series. We explicitly fixed three segments to separate the “free” spreading phase, the epidemic response, and the lockdown. (**B**) Time evolution of the per capita incidence (top), number of tests (middle) and test positivity (bottom) in Aragón (18) and Catalunya (19). While cases rise exponentially, the number of tests only increases linearly. Since test positivity does not significantly change, this strongly indicates saturation in test capacity. (**C**) Fit between the model and the daily fatalities.

However, this early decrease is found based on the reported case data. Therefore, the estimation of Rt is subjected to changes in the testing capacity and reporting criteria. The preceding decrease of Rt concerning mobility and containment measures could be due to a saturation of the test capacity, as illustrated in Fig. 3B for the CCAA Catalonia and Aragon. To corroborate this hypothesis, we compare the estimation with the most reliable epidemiological data in Spain, daily fatalities. In particular, we use our estimation to inform a simple, but robust, model (Methods Materials) about Rt’s general form, and adjust it to the daily fatalities. The first breakpoint, BP, which indicates the start of the epidemic response, as well as the initial number of infected individuals, I0, are left as free parameters. To reduce the number of free parameters, we fix the second breakpoint on March 15, one of two break points in the time series for Rt, Fig. 3A. Additionally, we assume a constant value R2 of Rt during the lockdown.

Figure 3C shows the fit of the model to the time series of daily fatalities. The parameters were calibrated as: I0 = 159(CI: 80-276), R1 = 4.75(CI: 4.3-5.3), R2 = 0.713(CI: 0.704 – 0.721), BP = 2020-03-12(CI: 2020-03-10-2020-03-13). Figure S2 presents the summary statistics. The corresponding evolution of Rt is shown in Fig.4A. The estimations of R1 and R2 lead to a doubling time and half-life time of 1.84(CI: 1.69– 1.99) days and 11.2(CI: 10.8– 11.6) days, respectively (27).

![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F4.medium.gif)

[Figure 4:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F4)

Figure 4: Underdetection and counterfactual scenarios
(**A**) Estimated Rt from the reconstructed exposure times together with the piecewise linear Rt, which we inferred through a fit to the daily fatalities. (**B**) Attack rate (top) and total number of fatalities (bottom) as a function of the shift in number of days of the epidemic response. (**C**) Estimated ascertainment bias according to the ratio between the daily infections from the model and the reconstructed exposures from the reported cases.

In Fig. 4A, we compare the Rt estimated from official incidence data, with the one fitted to reported daily fatalities. The first observation is that the actual Rt seems to have been substantially higher than what can be estimated from the reported incidence data. Additionally, we observe a reduction in Rt between March 10 and 13 substantially later than inferred from the reported data. Importantly, this indicates a decrease in Rt, preceding the implementation of the national lockdown. This seems to be the fingerprint of a spillover effect (16) from the earlier introduction of local restrictions as in Madrid, for example. These earlier restrictions raised the individual awareness of the population as indicated by an increase of COVID-19 related searches in Google from March 9 onward (see Fig. S3).

Figure 4C shows the estimated ascertainment bias according to the ratio between the daily infections from the model and the exposure times. Our results indicate that, initially, ascertainment was only around 10%. Additionally, ascertainment continuously deteriorated towards the first confinement, when it almost reached 5%. This decreasing ascertainment indicates that the ramping up of the test capacity could not keep up with the propagation of the epidemic.

Finally, the inferred form of Rt serves us to build counterfactual scenarios. We evaluate the accumulated fatalities, as well as the attack rate, if the epidemic response (including confinement, an earlier adoption of social distancing and the specific measures taken by the CCAA) had occurred earlier or later by shifting the function of Rt.

Figure 4B presents the total attack rate and accumulated fatalities if the entire epidemic response had occurred some days earlier/later. For a 7 days earlier response, we find an attack rate of 1.2% (CI: 1.1– 1.4%) and a total number of fatalities of 4,782 (CI: 4,146-5,503). In contrast, supposing the epidemic response 7 days later, we find an attack rate of 30.8% (CI: 28.2-34.0%) and a total number of fatalities of 120,353 (CI: 110,165-132,545). Due to the high Rt in Spain, which translates into a minimal doubling time, any delay in the epidemic response would have had drastic consequences.

## Discussion

Our results indicate the existence of a substantial delay between exposure and notification in Spain. An average delay of 15 days proves how late individuals with symptoms were tested, the substantial test turnaround time, and the significant delay to be notified. Such pronounced delay hinders an actual evaluation of the epidemiological situation and therefore impedes effective management of the epidemic response. For illustration, as the first confinement was implemented, our results indicate that the daily infections had already peaked. Similarly, on March 28, all non-essential economic activity was shut down despite still rising case numbers. However, we find that daily infections had already peaked more than 10 days before in most CCAA.

Besides the homogeneous, impactful containment measures implemented on a national level, we found that the earlier introduction of containment measures on a regional level led to differences when the daily case numbers first stopped growing for upto a week. In general, the CCAA,which first pushed Rt below 1,coincides with the ones that implemented containment measures the earliest. The most prominent example is the Community of Madrid. Despite the regional heterogeneity, the first confinement on March 15 fell into the range when Rt became smaller than 1.

By fitting a hypothesized functional form of Rt to reproduce the evolution of daily fatalities, we found that Rt started decreasing between March 10 and 13. The first containment measures taken from March 9 onwards and the corresponding mobility reduction which started around March 10, fall into this range. This reduction of Rt, before the first confinement, highlights the role of individual awareness of the population, a factor whose consideration is crucial (16). Nonetheless, the observation based on reconstructed infections contrasts with the reported data, which indicates a constant decrease in Rt from March 5 onwards. Available data from Aragón and Catalunya indicates that this earlier decrease can be attributed to saturation in the testing capacity. Similarly, the higher Rt value inferred from the model suggests that the expansion of test capacity could not keep up with the propagation speed of the epidemic in Spain. Accordingly, constantly fewer cases were detected towards lockdown before eventually only capturing approximately 5% of cases. In contrast, Althaus et al. find the opposite for Switzerland, where an increasing percentage of cases was detected towards lockdown (15). A decreasing detection rate, a substantial delay between exposure and notification, together with a doubling time of about 1.8 days, hampered the evaluation of Spain’s epidemiological situation and may have led to an underestimation regarding the severity of the situation.

Our counterfactual scenarios indicate that only a minimal delay in the epidemic response leads to a substantially different outcome. For an epidemic response 7 days earlier, we find approximately 5,000 fatalities and an attack rate of 1%, while a 7 days later response results in 120,000 fatalities and an attack rate of 30%. In this sense, the epidemic response and NPIs prevented many fatalities. Nevertheless, an earlier introduction of these measures would have substantially reduced COVID-19 related mortality.

## Data Availability

Only public data was used.

## Funding

A.A., B.S., and S.G. acknowledge financial support from Spanish MINECO (Grant No. PGC2018-094754-B-C21), Generalitat de Catalunya (Grant No. 2017SGR-896), and Universitat Rovira i Virgili (Grant No. 2019PFR-URV-B2-41). A.A. also acknowledges support from Generalitat de Catalunya (PDAD14/20/00001), ICREA Academia, and the James S. McDonnell Foundation (Grant No. 220020325). J.G.-G. and D.S.-P. acknowledge financial support from MINECO and FEDER funds (Project No. FIS2017-87519-P) and from the Departamento de Industria e Innovación del Gobierno de Aragón y Fondo Social Europeo (FENOL group E-19). C.G. acknowledges financial support from Beatriz de Galindo (Ministerio de Ciencia, Innovación y Universidades). B.S. acknowledges financial support from the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie Grant Agreement No. 713679 and from the Universitat Rovira i Virgili (URV).

## Author contributions

A.A., S.G. and J.G.-G. are joint senior authors. B.S. and A.A. designed the project. B.S., J.T.M., G.R. and D.S.-P. collected data. B.S., J.T.M., D.S.-P., G.R. and C.G. analyzed the data. B.S., J.T.M., D.S.-P., G.R., S.G, C.G. J.G.-G. and A.A. interpreted the results. B.S., J.T.M., D.S.-P., G.R., S.G, C.G. J.G.-G. and A.A. wrote the manuscript.

## Competing interests

Authors declare no competing interests.

## Data and materials availability

All data is available in the main text or the supplementary materials.

## Supplementary Material for

### Materials and Methods

#### Reconstruction of exposure times

To reconstruct the exposure times in the different CCAA, we use the function backprojNP (*17*) from the surveillance package (*18, 19*) in R. The method, initially proposed by Becker et al. (*7*), infers the expected number of exposures given the probability mass function of the incubation period through a maximum likelihood approach.

Credible intervals are calculated based on a bootstrap procedure (*20*). We fixed the smoothing factor *k* = 6, which corresponds to a centered rolling average of 7 days. The bootstrap procedure made use of 1000 samples (*B* = 1000). The incubation period distribution was fixed as a gamma distribution with mean 5.2 days and standard deviation 2.8 days (*21*). The time series with the symptom onset is provided by the Centro Nacional de Epidemilogía (*22*).

#### Estimation of *R**t*

From the median incidence, obtained by the reconstruction of the exposure times, we estimate the evolution of *R**t* using the EpiEstim package (*9, 23*). Our interest lies in the impact of the NPI. Accordingly, we estimate the instantaneous *R**t* (*24*) rather than the case *R**t* (*25*), since the former reflects changes in the transmission dynamics more abruptly.

For the infectivity function, we choose the generation time estimated by Ganyani et al. (*12*). To be more precise, we assume a generation time following a gamma distribution with mean 5.2 (CI : 3.78 – 6.78) days and standard deviation 1.72 (CI : 0.91 −3.93) days. This generation time distribution corresponds to the estimation by Ganynai et al. for Singapore while assuming the same incubation period distribution (*21*) as we do in the reconstruction of the exposure times. We assume a standard deviation for the mean and the standard deviation of the generation time of 1.0 days and 1.2 days, respectively. However, we bound the values for the mean and standard deviation by the estimations of Ganyani et al..

We fixed a centered rolling average of 7 days. To bootstrap the credible intervals, we took 100 samples of the generation time distribution and considered 100 posteriors for each of these samples (n1 = 100, n2 = 100).

#### Identifying segments of *R**t*

To identify the linear segments we use the R package segmented (*13,26*). The method proposed by Muggeo implements a maximum likelihood approach using linear predictors. The credible intervals are obtained through bootstrap. As previously pointed out, we assume three segments of *R**t*. An initial constant value R1 as the disease was spreading freely in Spain as well as two linear evolving parts corresponding to the decrease of *R**t* towards lockdown and the constant decrease observed during the lockdown.

#### Adjusting the model to the daily fatalities

To evaluate the validity of our estimation of *R**t*, we fit a piece-wise linear *R**t* to the epidemiological data in Spain. We then compare this fit with our estimation of *R**t* from the exposure times. Unfortunately, the only reliable time series available are the daily fatalities (*27*). Both, hospitalizations and ICU occupation, were subject to various changes of reporting criteria from March until May 2020. As an example, CCAA switched between reporting the accumulated number of admissions and the current number of ICUs occupied without announcement. For this reasons, we refrain from using hospitalization data and focus on the daily fatalities.

We assume *R**t* to have three segments as in the previous analysis. Motivated by the identified segments and to reduce the number of parameters to fit, we fix the second breakpoint on March 15, the day the lockdown was implemented. It seems reasonable that, as the first lockdown was implemented, *R**t* reached a relatively stable value. In contrast, the first breakpoint BP, which initiates the decrease of *R**t* and thus the epidemic response, we leave as a free parameter. In this way, we are able to evaluate to which extent the decrease of *R**t* before lockdown actually took place. Additionally, we assume a constant value R1 and R2 of *R**t* before and after the second breakpoint.

We then inform an epidemic model with the assumed form of *R**t*. Given *R**t*, the generation time distribution *w*(*t*) and the size of the population *N*, the daily incidence *I**t* on day *t* evolves as (*9, 24*): ![Formula][1]</img>  

In the above equation, we decoupled the effect of containment measures and a reduction of *R**t* due to natural immunity. This separation serves us later to build the counterfactual scenarios. We use the same generation time distribution as for estimating *R**t* (*12*). From the daily incidence, we propagate the symptom onset as well as the daily fatalities through convolution (*28*). Given the incubation time distribution *P* (*t*) as well as the time from symptom onset to decease *D*(*t*), the daily number of individuals developing symptoms *S**t* and the daily fatalities *F**t* evolve as: ![Formula][2]</img>  ![Formula][3]</img>  

We fix the incubation time period distribution as for the estimation of the exposure times (*21*). We set the IFR as 0.83% according to the nationwide seroprevalence study (*29*). The average time from symptom onset to decease is fixed as 11 days according to the information by the public health authorities (*30*). Since no information is available on the type of distribution which fits best the data from symptom onset to decease, we make use of a gamma distribution (*5, 31*). We fix the shape factor (0.63) to match the interquartile range of the empirical distribution of 7 and 17 (*30*). In the Supplementary Text we perform a sensitivity analysis by considering a log-normal distribution.

Finally, in order to fit the model to the data, we make use of a MCMC approach. To be more precise, we use the DREAMzs sampler (*32, 33*) from the BayesianTools package (*34*). We use all flat priors, namely I0 ∈ [1, 1000], R1 ∈ [3, 8], R2 ∈ [0, 1] and BP between March 1st and 15. The number of iterations is fixed to 3 ·106 with three chains. We fix the snooker update probability as 0.5. We use 20% of the iterations for adaption and discard the first 30% of the iterations. The posteriors and their correlation are shown in Fig. S2, whereas the actual fit is presented in Fig. 3C.

#### Counterfactual scenarios

To analyze what would have happened if the epidemic response had occurred earlier or later, we make use of the previous fit of the piece-wise linear *R**t*. By shifting the breaking points, we can analyze straightforwardly the impact of an earlier or later response since we decoupled the containment from the effect of natural immunity. This approach is strongly inspired by Althaus et al. (*11*). We make all the simulations sampling 5000 parameter sets from the calibration’s posterior to bootstrap the credible intervals.

## Supplementary Text

### Heterogeneity in the epidemic response between CCAA

In the main text, Fig. 2B shows the precise day on which *R**t* was first below one in the CCAA. Here, we describe in more detail how this day correlates with the measures implemented in the CCAA. The Community of Madrid was the first CCAA to push *R**t* number below one. Madrid was the strongest hit CCAA early on and implemented first containment measures on March 9 before all other CCAA. Similarly, other CCAA, which first pushed *R**t* below one, were among the first to introduce containment measures such as Asturias (*35*), Valencian Community (*36*) or Basque Country (*37*). Interestingly, the last CCAA to push *R**t* below one, are in geographical proximity to Madrid such as Aragon, Castile-La Mancha, Castile and León, and Extremadura. This is consistent with a study by Mazzoli et al. that indicates how the local incidence strongly correlated with mobility in the early phase from and to Madrid, due to multi-seeding (*38*). Furthermore, Castile-La Mancha, Castile and León, and Extremadura have the highest percentages in Spain of their students attending university in the Community of Madrid; 27%, 14%, and 12%, respectively (*39*). Even though it is speculative, it seems reasonable that, after the closure of the universities, the outflux of students from Madrid to their hometowns added various seeds with mostly mild symptoms, accelerating the evolution of the epidemic in these regions.

### Cross-correlation between mobility and the regional *R**t*

The SARS-CoV-2 epidemic led to a substantial reduction in mobility. Figure S4 presents this data from the Google mobility report. We note a strong decrease of all types of mobility before, and shortly after the lockdown. Afterwards, the mobility stays very stable, but an additional decrease can be seen coinciding with the second lockdown as all-non essential mobility was closed. This additional reduction of mobility during the closure of all non-essential is better illustrated in Fig. S5, which shows the average mobility during the three subsequent lockdowns. We also note that the third lockdown shows a slightly stronger mobility reduction compared to the first one despite the same restrictions being applied.

Comparing the mobility reduction with the decrease in *R**t*, we find an elevated cross-correlation. Figure S6 presents the cross-correlation between the mobility types and *R**t* for the CCAA. A detailed distribution for the maximal cross-correlation is shown in Fig. S7. Additionally, with very little exceptions, the maximum cross-correlation is reached without any lag between the two time series. The exception here is the mobility related to grocery shops and pharmacies. This negative lag is due to the strong increase in visits of grocery stores and pharmacies due to fear of supply shortages as shown in Fig. S4.

### Detailed analysis of segments in the national *R**t*

In the main text we separated *R**t* into three linear segments: the “free” spreading phase before restrictions, the initiation of the epidemic response, and the lockdown. Here, we present a more detailed analysis regarding the epidemic response and the lockdown. Therefore, we cut the time series of *R**t* at the turning point as the decrease is initiated in the three segment analysis. Subsequently, we analyze the different segments observed until the end of lockdown. We choose the number of segments according to the Bayesian information criterion (*40*) and corroborated with the Akaike information criterion (*41*) by iterating from 1 to 15 segments.

The identified segments are shown in Fig. S8 (bottom) together with the evolution of mobility according to the COVID-19 Community Mobility report from Google (*42*) (top). Interestingly, between the first and second breakpoint, we observe a substantial reduction of *R**t* before mobility starts to decrease. The second breakpoint then coincides with the initiation of mobility reduction. The third breakpoint matches the implementation of the first lockdown on March 15. At breakpoint 4, *R**t* starts a plateau together with the mobility level. The breakpoint 6 then corresponds to the shutdown of all non-essential activity, the second lockdown. This also coincides with a further reduction in overall mobility as shown in Fig. S5.

The lift of the second lockdown and the associated mobility increase, seem not to have caused any increase in *R**t*. It is probable that the good weather conditions, as well as an increased individual awareness, raised by the elevated death toll and reflected in the reduced mobility (see Fig. S5), have contributed to a continuous decrease in *R**t*. Finally, at breakpoint 9 the increase in *R**t* coincides with the lift of the lockdown and corresponding increase in mobility. Various identified breakpoints cannot be associated with any nation-wide introduced measures put in place by the national authorities or changes in mobility. They may stem from measures introduced in the different CCAA.

### Sensitivity Analysis

For adjusting the model to the daily fatalities we assumed a gamma distribution. The mean of the distribution was fixed to 11 days according to the information provided by the health authorities (*30*). Furthermore, we adjusted the shape factor to match the interquartile range (7 and 17 days) since this is, besides the mean, the only information we have at our disposal of the empirical distribution. Even though gamma distributions have been shown to describe particularly well the time between symptom onset and disease (*31*), we performed the same analysis making use of a log-normal distribution. We assumed the same priors as in the case of the gamma distribution.

Adjusting the log-normal distribution to the interquartile range, we found a logarithmic standard deviation of 0.54. The log-normal distribution is compared to the gamma distribution in Fig. S9. The parameters we find adjusting to the daily fatalities are: I0 = 261 (CI: 139 – 434), R1 = 4.5 (CI: 4.0 −5.0), R2 = 0.732 (CI: 0.724 −0.740), BP = 2020-03-11 (CI: 2020-03-09 – 2020-03-12). The estimations of R1 and R2 lead to a doubling time and half-life time of 1.94 (CI: 1.75 – 2.11) days and 12.1 (CI: 11.7−12.5) days, respectively (*12*). The results are not substantially altered with respect to a gamma distribution. The smaller *R**t* before the epidemic leads to a higher number of initially infected individuals. Furthermore, the shape of the log-normal distribution causes *R**t* to decrease a day earlier. The posteriors and their correlations are shown in Fig. S10.

Figure S11 shows the obtained fit between the model and the data. Figure S12 shows the inferred form of *R**t*, while in Fig. S13 we present the evolution of under detection. The figures show very similar results as in the main text. Regarding the counterfactual scenarios, we find 5, 274 (CI: 4, 471 – 6, 005) fatalities and an attack rate of 1.4% (CI: 1.1− 1.5%) for a 7 day earlier epidemic response. Considering a 7 day later epidemic response, the model shows 112, 416 (CI: 102, 537 – 125, 562) fatalities and an attack rate of 28.8% (CI: 26.2 −32.%). The complete results are shown in Fig. S14.

![Figure S1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F5.medium.gif)

[Figure S1:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F5)

Figure S1: 
First breakpoint in *R**t* and workplace related mobility for Spain and all CCAA. Points indicate the mean, whereas horizontal bars indicate 95% credible intervals. We explicitly separated the time series into three segments, and identified them with the segmented package (*26*). Similarly as in the main text, we assumed the first segment of the time series was constant, and the others linear. To discard initial fluctuations in *R**t* due to a small incidence, we did not consider values before March 1. Additionally, we cut the series on March 28 to exclude the additional decrease in mobility due to lockdown 2. Overall, *R**t* decreases substantially earlier than mobility.

![Figure S2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F6.medium.gif)

[Figure S2:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F6)

Figure S2: 
Posterior distribution of the model parameters and their correlations, assuming the time between symptom onset to death follows a gamma distribution. The range of the flat priors was larger than the posteriors in all cases. For the variable BP, the integer values refer to the dates in March 2020.

![Figure S3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F7.medium.gif)

[Figure S3:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F7)

Figure S3: 
Frequency of different COVID-19 related search queries according to Google Trends. Note that Google Trends normalizes the frequency of queries between 0 and 100 during the considered period. All key words see a sharp increase from March 9 onward. On the same day, the Community of Madrid (*43*) closed all educational centres and the Basque Country (*37*) announced their closure for the next day. These closures correspond to the first NPIs introduced in Spain.

![Figure S4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F8.medium.gif)

[Figure S4:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F8)

Figure S4: 
Visualization of the mobility data from the COVID-19 Community Mobility report from Google (*42*). Dots indicate the data points, while solid lines represent a centered, 7 day rolling average. The shaded areas in light gray indicate the lockdowns 1 and 3. In dark gray, we indicate the lockdown 2.

![Figure S5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F9.medium.gif)

[Figure S5:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F9)

Figure S5: 
Average mobility according to the COVID-19 Community Mobility report from Google (*42*) during the three lockdown stages in Spain. Mobility is substantially lower during the second lockdown since all non-essential economic activity was shutdown. Additionally, mobility was lower during lockdown 3 when compared to lockdown 1, even though restrictions were the same. This indicates an increased awareness among the population, probably due to the huge death toll to this point.

![Figure S6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F10.medium.gif)

[Figure S6:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F10)

Figure S6: 
Lag and value of maximal cross-correlation between different mobility types, according to the COVID-19 Community Mobility report from Google (*42*) and *R**t*. Negative lag indicates that *R**t* precedes mobility.

![Figure S7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F11.medium.gif)

[Figure S7:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F11)

Figure S7: 
Maximal cross-correlation among CCAA between different mobility types, according to the COVID-19 Community Mobility report from Google (*42*) and *R**t*. For the Residential mobility category we take the absolute value to facilitate visualization.

![Figure S8:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F12.medium.gif)

[Figure S8:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F12)

Figure S8: 
(Top) Average national mobility according to the COVID-19 Community Mobility report from Google (*42*). Dots indicate the data points and solid lines represent a centered, 7 day rolling average. (Bottom) *R**t* and the identified segments. Purple dots are the estimated values of *R**t* while bars indicate credible intervals. Yellow crosses refer to the breakpoints, and green dashed rectangles indicate their uncertainty range. The shaded area in light gray indicates lockdowns 1 and 3. In dark gray, we indicate lockdown 2.

![Figure S9:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F13.medium.gif)

[Figure S9:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F13)

Figure S9: 
Gamma and log-normal distributions used for the time between symptom onset and decease. The average of both distributions is 11 days. The gamma distribution has a shape factor of 0.63, whereas the logarithmic standard deviation of the log-normal distribution is 0.54.

![Figure S10:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F14.medium.gif)

[Figure S10:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F14)

Figure S10: 
Posterior distribution for the model parameters and their correlations, assuming the time between symptom onset to death follows a log-normal distribution. The range of the flat priors was larger than the posteriors in all cases. For the variable BP, the integer values refer to the dates in March 2020.

![Figure S11:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F15.medium.gif)

[Figure S11:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F15)

Figure S11: 
The daily fatalities in Spain from March 8 to May 8 compared to the daily fatalities for the sensitivity analysis. In the sensitivity analysis, we used a log-normal distribution describing the time between from symptom onset to decease. The shaded area represents 95% credible intervals.

![Figure S12:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F16.medium.gif)

[Figure S12:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F16)

Figure S12: 
*R**t* estimated from the reconstructed exposure times vs. the model inferred as part of the sensitivity analysis. This is equivalent to Fig. 4A in the main text when using the log-normal distribution.

![Figure S13:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F17.medium.gif)

[Figure S13:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F17)

Figure S13: 
Ratio between the reconstructed exposure times and the model inferred infections as part of the sensitivity analysis. This is equivalent to Fig. 4B in the main text when using the log-normal distribution.

![Figure S14:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.16.21251832/F18.medium.gif)

[Figure S14:](http://medrxiv.org/content/early/2021/02/18/2021.02.16.21251832/F18)

Figure S14: 
Counterfactual scenarios for the attack rate as well as the total number of fatalities. This is equivalent to Fig. 4C in the main text when using the log-normal distribution.

## Acknowledgments

The authors acknowledge Prof. Dan Larremore and Prof. Miguel Hernán for their helpful comments on the manuscript.

*   Received February 16, 2021.
*   Revision received February 16, 2021.
*   Accepted February 18, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Prevalence of SARS-CoV-2 in Spain (ENE-COVID): a nationwide, population-based seroepidemiological study. The Lancet, 396(10250): 535–544, 2020.
    
    

2.  2. F. Li,  Y. Y. Li,  M.J. Liu,  L.Q. Fang,  N.E. Dean,  G.W.K. Wong,  X.B. Yang,  I. Longini,  M.E. Halloran,  H.J. Wang,  P.L. Liu,  Y.H. Pang,  Y.Q. Yan,  S. Liu,  W. Xia,  X.X. Lu,  Y. Yang,  S.Q. Xu. Household transmission of SARS-CoV-2 and risk factors for susceptibility and infectivity in Wuhan: a retrospective observational study. Lancet Infectious Diseases, 2021.
    
    

3.  3. V. Marziano,  G. Guzzetta,  B.M. Rondinone,  F. Boccuni,  F. Riccardo,  Antonino Bella,  Piero Poletti,  Filippo Trentini,  Patrizio Pezzotti,  Silvio Brusaferro,  Giovanni Rezza,  Sergio Iavicoli,  Marco Ajelli,  Stefano Merler. Retrospective analysis of the Italian exit strategy from COVID-19 lockdown. Proceedings of the National Academy of Sciences USA 118(4):e2019617118, 2021.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxNzoiMTE4LzQvZTIwMTk2MTcxMTgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMi8xOC8yMDIxLjAyLjE2LjIxMjUxODMyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

4.  4. N. Haug,  L. Geyrhofer,  A. Londei,  E. Dervic,  A. Desvars-Larrive, V. Loreto, B. Pinior,  S. Thurner, P. Klimek. Ranking the effectiveness of worldwide COVID-19 government interventions. Nature Human Behaviour, 4:1303–1312, 2020.
    
    

5.  5. S. Flaxman,  S. Mishra,  A. Gandy,  H. Juliette  T. Unwin,  T.A. Mellan,  H. Coupland,  C. Whittaker,  H. Zhu,  T. Berah,  J.W. Eaton,  M. Monod, Imperial College COVID-19 Response Team,  A.C. Ghani,  C.A. Donnelly,  S. Riley,  M.A.C. Vollmer,  N.M. Ferguson, L.C. Okell,  S. Bhatt. Estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe. Nature, 584(7820):257– 261, 2020.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.16.21251832.atom) 

6.  6. J. Dehning,  J. Zierenberg,  F. P. Spitzner,  M. Wibral,  J. P. Neto,  M. Wilczek,  V. Priesemann.Inferring change points in the spread of COVID-19 reveals the effectiveness of interventions.Science, 369(6500):eabb9789, 2020.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE3OiIzNjkvNjUwMC9lYWJiOTc4OSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAyLzE4LzIwMjEuMDIuMTYuMjEyNTE4MzIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

7.  7. N.G. Becker,  L.F. Watson,  J.B. Carlin. A method of non-parametric back-projection and its application to AIDS data. Statistics in Medicine, 10(10):1527–1542,1991.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/sim.4780101005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1947509&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.16.21251832.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1991GH68700004&link_type=ISI) 

8.  8.Centro Nacional de Epidemiología. COVID-19 en España, 2020.[https://cnecovid.isciii.es/covid19/](https://cnecovid.isciii.es/covid19/)
    
    

9.  9. A. Cori,  N.M. Ferguson,  C. Fraser,  S. Cauchemez. A new framework and software to estimate time-varying reproduction numbers during epidemics. American Journal of Epidemiology, 178(9):1505–1512, 2013.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwt133&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24043437&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.16.21251832.atom) 

10. 10. A. Arenas,  W. Cota, J.Gómez-Gardeñes, S. Gómez, C. Granell,  J.T. Matamalas,  D. Soriano-Paños,  B. Steinegger. Modeling the spatiotemporal epidemic spreading of COVID-19 and the impact of mobility and social distancing interventions. Physical Review X, 10:041055, 2020.
    
    

11. 11. C.L. Althaus,  D. Probst,  A. Hauser,  J. Riou. Time is of the essence: containment of the SARS-CoV-2 epidemic in Switzerland from February to May 2020. medRxiv, 2020.07.21.20158014, 2020.
    
    

12. 12. T. Ganyani,  C. Kremer,  D. Chen,  A. Torneri,  C. Faes,  J. Wallinga,  C. Hens. Estimating the generation interval for coronavirus disease (COVID-19) based on symptom onset data, march 2020. Eurosurveillance, 25(17): 2000257, 2020.
    
    

13. 13. V.M.R. Muggeo. Estimating regression models with unknown break-points. Statistics in Medicine, 22: 3055–3071, 2003.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/sim.1545&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12973787&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.16.21251832.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000185451700006&link_type=ISI) 

14. 14.Departamento de Sanidad. Gobierno de Aragón. Datacovid: [https://datacovid.salud.aragon.es/covid/](https://datacovid.salud.aragon.es/covid/)
    
    

15. 15.Departament de Salut de la Generalitat de Catalunya. Dades COVID: [https://dadescovid.cat/](https://dadescovid.cat/)
    
    

16. 16. D. Holtz,  M. Zhao,  S.G. Benzell,  C.Y. Cao,  M.A. Rahimian,  J. Yang,  J. Allen,  A. Collis,  A. Moehring,  T. Sowrirajan,  D. Ghosh,  Y. Zhang,  P.S. Dhillon,  C. Nicolaides,  D. Eckles,  S. Aral, Interdependence and the cost of uncoordinated responses to COVID-19, Proceedings of the National Academy of Sciences 117, 33, 19837--19843 (2020)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzMzLzE5ODM3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDIvMTgvMjAyMS4wMi4xNi4yMTI1MTgzMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

## References

1.  1. M. Pollán,  B. Pérez-Gómez,  R. Pastor-Barriuso,  J. Oteo,  M. A. Hernán,  M. Pérez-Olmeda,  J. L. Sanmartín,  A. Fernández-García,  I. Cruz,  N. Fernández de Larrea,  M. Molina,  F. Rodríguez-Cabrera,  M. Martín,  P. Merino-Amador,  J. León Paniagua,  J. F. Muñoz-Montalvo,  F. Blanco,  R. Yotti,  R. Gutiérrez Fernández,  S. Mezcua Navarro,  J. F. Muñoz-Montalvo,  M. Salinero Hernández,  J. L. Sanmartín,  M. Cuenca-Estrella,  J. León Paniagua,  P. Fernández-Navarro,  A. Avellón,  G. Fedele,  J. Oteo Iglesias,  M. T. Pérez Olmeda,  M. E. Fernandez Martinez,  F. D. Rodríguez-Cabrera,  M. A. Hernán,  S. Padrones Fernández,  J. Rumbao Aguirre,  J. Navarro Marí,  B. Palop Borrás,  A. B. Pérez Jiménez,  M. Rodríguez-Iglesias,  A. M. Calvo Gascón,  M. L. Lou Alcaine,  I. Donate Suárez,  O. Suárez Á lvarez,  M. Rodríguez Pérez,  M. Cases Sanchís,  C. J. Villafáfila Gomila,  L. Carbo Saladrigas,  A. Hurtado Fernández,  A. Oliver,  E. Castro Feliciano,  M. N. González Quintana,  J. Bar-rasa Fernández,  M. A. Hernández Betancor, M. Hernández Febles,  L. Martín Martín,  L.-M. López López,  T. Ugarte Miota,  I. De Benito Población,  M. S. Celada Pérez,  M. N. Vallés Fernández,  T. MatéEnríquez,  M. Villa Arranz,  M. Domínguez-Gil González,  I. Fernández-Natal,  G. Megías Lobón,  J. L. Muñoz Bellido,  P. Ciruela,  A. Mas i Casals,  M. DoladéBotías,  M. A. Marcos Maeso,  D. Pérez del Campo,  A. Félix de Castro,  R. Limón Ramírez,  M. F. Elías Retamosa,  M. Rubio González,  M. S. Blanco Lobeiras,  A. Fuentes Losada,  A. Aguilera,  G. Bou,  Y. Caro,  N. Marauri,  L. M. Soria Blanco,  I. del Cura González,  M. Hernández Pascual,  R. Alonso Fernández,  N. Cabrera Castro,  A. Tomás Lizcano,  C. Ramírez Almagro,  M. Segovia Hernández,  N. Ascunce Elizaga,  M. Ederra Sanz,  C. Ezpeleta Baquedano,  A. Bustinduy Bascaran,  S. Iglesias Tamayo,  L. Elorduy Otazua,  R. Benarroch Benarroch,  J. Lopera Flores,  A. Vázquez de la Villa, Prevalence of SARS-CoV-2 in Spain (ENE-COVID): a nationwide, population-based seroepidemiological study, The Lancet 396, 535 (2020).
    
    

2.  2. F. Li,  Y.-Y. Li,  M.-J. Liu,  L.-Q. Fang,  N. E. Dean,  G. W. K. Wong,  X.-B. Yang,  I. Longini,  M. E. Halloran,  H.-J. Wang,  P.-L. Liu,  Y.-H. Pang,  Y.-Q. Yan,  S. Liu,  W. Xia,  X.-X. Lu,  Q. Liu,  Y. Yang,  S.-Q. Xu, Household transmission of SARS-CoV-2 and risk factors for susceptibility and infectivity in Wuhan: a retrospective observational study, The Lancet Infectious Diseases (2021).
    
    

3.  3. V. Marziano,  G. Guzzetta,  B. M. Rondinone,  F. Boccuni,  F. Riccardo,  A. Bella,  P. Poletti,  F. Trentini,  P. Pezzotti,  S. Brusaferro,  G. Rezza,  S. Iavicoli,  M. Ajelli,  S. Merler, Retro-spective analysis of the italian exit strategy from COVID-19 lockdown, Proceedings of the National Academy of Sciences (USA) 118 (2021).
    
    

4.  4. N. Haug,  L. Geyrhofer,  A. Londei,  E. Dervic,  A. Desvars-Larrive,  V. Loreto,  B. Pinior,  S. Thurner,  P. Klimek, Ranking the effectiveness of worldwide COVID-19 government interventions, Nature Human Behaviour (2020).
    
    

5.  5. S. Flaxman,  S. Mishra,  A. Gandy,  H. J. T. Unwin,  T. A. Mellan,  H. Coupland,  C. Whittaker,  H. Zhu,  T. Berah,  J. W. Eaton,  M. Monod,  P. N. Perez-Guzman,  N. Schmit,  L. Cilloni,  K. E. Ainslie,  M. Baguelin,  A. Boonyasiri,  O. Boyd,  L. Cattarino,  L. V. Cooper,  Z. Cu-cunubá,  G. Cuomo-Dannenburg,  A. Dighe,  B. Djaafara,  I. Dorigatti,  S. L. van Elsland,  R. G. FitzJohn,  K. A. Gaythorpe,  L. Geidelberg,  N. C. Grassly,  W. D. Green,  T. Hallett,  A. Hamlet,  W. Hinsley,  B. Jeffrey,  E. Knock,  D. J. Laydon,  G. Nedjati-Gilani,  P. Nouvel-let,  K. V. Parag,  I. Siveroni,  H. A. Thompson,  R. Verity,  E. Volz,  C. E. Walters,  H. Wang,  Y. Wang,  O. J. Watson,  P. Winskill,  X. Xi,  P. G. Walker,  A. C. Ghani,  C. A. Donnelly,  S. Riley,  M. A. Vollmer,  N. M. Ferguson,  L. C. Okell,  S. Bhatt, Estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe, Nature 584, 257 (2020).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.16.21251832.atom) 

6.  6. J. Dehning,  J. Zierenberg,  F. P. Spitzner,  M. Wibral,  J. P. Neto,  M. Wilczek,  V. Priesemann, Inferring change points in the spread of COVID-19 reveals the effectiveness of interventions, Science 369 (2020).
    
    

7.  7. N. G. Becker,  L. F. Watson,  J. B. Carlin, A method of non-parametric back-projection and its application to aids data, Statistics in Medicine 10, 1527 (1991).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/sim.4780101005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1947509&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.16.21251832.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1991GH68700004&link_type=ISI) 

8.  8.Centro Nacional de Epidemiología, COVID-19 en España (2020).
    
    

9.  9. A. Cori,  N. M. Ferguson,  C. Fraser,  S. Cauchemez, A new framework and software to estimate time-varying reproduction numbers during epidemics, American Journal of Epi-demiology 178, 1505 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwt133&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24043437&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.16.21251832.atom) 

10. 10. A. Arenas,  W. Cota,  J. Gómez-Gardeñes,  S. Gómez,  C. Granell,  J. T. Matamalas,  D. Soriano-Paños,  B. Steinegger, Modeling the spatiotemporal epidemic spreading of COVID-19 and the impact of mobility and social distancing interventions, Physical Review X 10, 041055 (2020).
    
    

11. 11. C. L. Althaus,  D. Probst,  A. Hauser,  J. Riou, Time is of the essence: containment of the SARS-CoV-2 epidemic in Switzerland from February to May 2020, medRxiv p. 2020.07.21.20158014 (2020).
    
    

12. 12. T. Ganyani,  C. Kremer,  D. Chen,  A. Torneri,  C. Faes,  J. Wallinga,  N. Hens, Estimating the generation interval for coronavirus disease (COVID-19) based on symptom onset data, March 2020, Eurosurveillance 25 (2020).
    
    

13. 13. V. M. Muggeo, Estimating regression models with unknown break-points, Statistics in Medicine 22, 3055 (2003).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/sim.1545&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12973787&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.16.21251832.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000185451700006&link_type=ISI) 

14. 14.G. d. A. Departamento de Sanidad, Datacovid: [https://datacovid.salud.aragon.es/covid/](https://datacovid.salud.aragon.es/covid/).
    
    

15. 15.D. de Salut de la Generalitat de Catalunya, Dades covid: [https://dadescovid.cat/](https://dadescovid.cat/).
    
    

16. 16. D. Holtz,  M. Zhao,  S. G. Benzell,  C. Y. Cao,  M. A. Rahimian,  J. Yang,  J. Allen,  A. Collis,  A. Moehring,  T. Sowrirajan,  D. Ghosh,  Y. Zhang,  P. S. Dhillon,  C. Nicolaides,  D. Eckles,  S. Aral, Interdependence and the cost of uncoordinated responses to covid-19, Proceedings of the National Academy of Sciences 117, 19837 (2020).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzMzLzE5ODM3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDIvMTgvMjAyMS4wMi4xNi4yMTI1MTgzMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

17. 17. D. S. B. Michael Höhle, Non-parametric back-projection of incidence cases to exposure time (2019).
    
    

18. 18. M. Salmon,  D. Schumacher,  M. Höhle, Monitoring count time series in R: Aberration detection in public health surveillance, Journal of Statistical Software 70, 1 (2016).
    
    

19. 19. S. Meyer,  L. Held,  M. Höhle, Spatio-temporal analysis of epidemic phenomena using the R package surveillance, Journal of Statistical Software 77, 1 (2017).
    
    

20. 20. P. S. Yip,  K. F. Lam,  Y. Xu,  P. H. Chau,  J. Xu,  W. Chang,  Y. Peng,  Z. Liu,  X. Xie,  H. Y. Lau, Reconstruction of the infection curve for SARS epidemic in Beijing, China using a back-projection method, Communications in Statistics: Simulation and Computation 37, 425 (2008).
    
    

21. 21. J. Zhang,  M. Litvinova,  W. Wang,  Y. Wang,  X. Deng,  X. Chen,  M. Li,  W. Zheng,  L. Yi,  X. Chen,  Q. Wu,  Y. Liang,  X. Wang,  J. Yang,  K. Sun,  I. M. Longini,  M. E. Halloran,  P. Wu,  B. J. Cowling,  S. Merler,  C. Viboud,  A. Vespignani,  M. Ajelli,  H. Yu, Evolving epidemiology and transmission dynamics of coronavirus disease 2019 outside Hubei province, China: a descriptive and modelling study, The Lancet Infectious Diseases 20, 793 (2020).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.16.21251832.atom) 

22. 22.Centro Nacional de Epidemilogía, COVID-19 en España (2020).
    
    

23. 23. A. Cori, EpiEstim: Estimate time varying reproduction numbers from epidemic curves (2020). R package version 2.2-3.
    
    

24. 24. C. Fraser, Estimating individual and household reproduction numbers in an emerging epidemic, PLoS ONE 2 (2007).
    
    

25. 25. J. Wallinga,  P. Teunis, Different epidemic curves for severe acute respiratory syndrome reveal similar impacts of control measures, American Journal of Epidemiology 160, 509 (2004).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwh255&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15353409&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.16.21251832.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000223938000001&link_type=ISI) 

26. 26. V. M. Muggeo, segmented: an R package to fit regression models with broken-line relationships, R News 8, 20 (2008).
    
    

27. 27. P. De Salazar,  F. Lu,  J. Hay,  D. Gómez-Barroso,  P. Fernández-Navarro,  E. Martínez,  J. Astray-Mochales,  R. Amillategui,  A. García-Fulgueiras,  M. Chirlaque,  A. Sánchez-Migallón,  A. Larrauri,  M. Sierra,  M. Lipsitch,  F. Simón,  M. Santillana,  M. Hernán, Near real-time surveillance of the sars-cov-2 epidemic with incomplete data (2021).
    
    

28. 28. G. N. Wong,  Z. J. Weiner,  A. V. Tkachenko,  A. Elbanna,  S. Maslov,  N. Goldenfeld, Modeling COVID-19 dynamics in Illinois under nonpharmaceutical interventions, Physical Review X 10, 041033 (2020).
    
    

29. 29. R. Pastor-Barriuso,  B. Pérez-Gómez,  M. A. Hernán,  M. Pérez-Olmeda,  R. Yotti,  J. Oteo-Iglesias,  J. L. Sanmartín,  I. León-Gómez,  A. Fernández-García,  P. Fernández-Navarro,  I. Cruz,  M. Martín,  C. Delgado-Sanz,  N. Fernández de Larrea,  J. León Paniagua,  J. F. Muñoz-Montalvo,  F. Blanco,  A. Larrauri,  M. Pollán, on behalf of the ENE-COVID Study Group, Infection fatality risk for SARS-CoV-2 in community dwelling population of spain: nationwide seroepidemiological study, BMJ 371 (2020).
    
    

30. 30.Instituto de Salud Carlos III, COVID-19. Informes previos.
    
    

31. 31. I. Hawryluk,  T. A. Mellan,  H. Hoeltgebaum,  S. Mishra,  R. P. Schnekenberg,  C. Whittaker,  H. Zhu,  A. Gandy,  C. A. Donnelly,  S. Flaxman,  S. Bhatt, Inference of COVID-19 epidemiological distributions from Brazilian hospital data: Inference of COVID-19 epidemiological distributions from Brazilian hospital data, Journal of the Royal Society Interface 17 (2020).
    
    

32. 32. C. J. Ter Braak,  J. A. Vrugt, Differential Evolution Markov Chain with snooker updater and fewer chains, Statistics and Computing 18, 435 (2008).
    
    

33. 33. J. A. Vrugt,  C. J. Ter Braak,  C. G. Diks,  B. A. Robinson,  J. M. Hyman,  D. Higdon, Accelerating Markov chain Monte Carlo simulation by differential evolution with self-adaptive randomized subspace sampling, International Journal of Nonlinear Sciences and Numerical Simulation 10, 273 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1515/Ijnsns.2009.10.3.273&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000264593800001&link_type=ISI) 

34. 34. F. Hartig,  F. Minunno,  S. Paul, BayesianTools: General-purpose MCMC and SMC samplers and tools for Bayesian statistics (2019). R package version 0.1.7.
    
    

35. 35.La Nueva España, El Principado decreta el cierre de todos los centros educativos asturianos ante la recomendación del Gobierno (2020).
    
    

36. 36.La Vanguardia, Qué impacto económico tendrá la cancelación de las Fallas de València? (2020).
    
    

37. 37.El Confidencial, El gobierno vasco decreta el cierre de todos los colegios de Vitoria 15 días por el COVID-19 (2020).
    
    

38. 38. M. Mazzoli,  D. Mateo,  A. Hernando,  S. Meloni,  J. J. Ramasco, Effects of mobility and multi-seeding on the propagation of the COVID-19 in Spain, medRxiv (2020).
    
    

39. 39.Ministerio de Ciencia, Innovación y Universidades, Datos y cifras del sistema universitario espanol. publicacion 2018-2019 (2020).
    
    

40. 40. G. Schwarz, Estimating the dimension of a model, The Annals of Statistics 6, 461 (1978).
    
    

41. 41. P. V. Bertrand,  Y. Sakamoto,  M. Ishiguro,  G. Kitagawa, Akaike Information Criterion statistics., Journal of the Royal Statistical Society. Series A (Statistics in Society) 151, 567 (1988).
    
    

42. 42.Google, COVID-19 Community Mobility Reports (2020).
    
    

43. 43. E. Páis, Madrid deja sin clases a 1,5 millones de alumnos y cancelará operaciones y citas para combatir el coronavirus (2020).

 [1]: /embed/graphic-5.gif
 [2]: /embed/graphic-6.gif
 [3]: /embed/graphic-7.gif