Evaluating the policy of closing bars and restaurants in Cataluña and its effects on mobility and COVID19 incidence

The world has gone through unprecedented changes since the global pandemic hit. During the early phase of the pandemic, the absence of known drugs or pharmaceutical treatments forced governments to introduce different policies in order to help reduce contagion rates and manage the economic consequences of the pandemic. This paper analyses the causal impact on mobility and COVID19 incidence from policy makers in Cataluña, Spain. We use anonymized phone-based mobility data together with reported incidence and apply a series of causal impact models frequently used in econometrics and policy evaluation in order to measure the policies impact. We analyse the case of Cataluña and the public policy decision of closing all bars and restaurants down for a 5 week period between 2020-16-10 and 2020-23-11. We find that this decision led to a significant reduction in mobility. It not only led to reductions in mobility but from a behavioural economics standpoint, we highlight how people responded to the policy decision. Moreover, the policy of closing bars and restaurants slowed the incidence rate of COVID19 after a time lag has been taken into account. These findings are significant since governments worldwide want to restrict movements of people in order to slow down COVID19 incidence without infringing on their rights directly.

New cases per day Figure 2. New cases per day for each CCAA: number of COVID19 cases for each autonomous community (CCAA) in Spain. Madrid experienced its second peak before the rest of Spain, with Cataluña (as the bold line) experiencing its second peak just after. The policy was a direct response to control the outbreak of the second peak in Cataluña. On the x-axis, ene refers to Enero (January).

Literature review
This section discusses relevant literature. We first break the section down to related literature using interrupted time series for policy decisions. Secondly, we link mobility with public policy decision making and finally we discuss how mobility has been linked to COVID19 incidence.
Interrupted time series. Interrupted time series models are being increasingly applied when analysing the impact of public health interventions. Reference 16 applied Box-Tiao intervention analysis to analyse the effects of the introduction of U.S. legislation requiring the use of mandatory seat belts across 8 states between 1976 and 1986. Reference 17 applied interrupted time series models to study the association between the introduction of helmet legislation and admissions to hospital for cycling-related head injuries between 1994 and 2003 in Canada. Reference 18 used interrupted time series to assess the effect of U.K. legislation of reducing package sizes of paracetamol on deaths from paracetamol poisoning between 1998 and 2009. Reference 19 used interrupted time series to study the effect of the introduction of 20 mph (32 km an hour) traffic speed zones on-road collisions, injuries and fatalities in London between 1986 and 2006. Reference 20 used interrupted time series to analyse the association between the 2008 financial crisis and suicide rates in Spain, using data between 2005 and 2010. Interrupted time series models have long be used to analyse public policy and time-series related events in economics and provide a suitable methodology for analysing policy events related to COVID19.
Mobility and public policy. Reference 21 used Google Community Mobility Data to analyses mobility changes from a baseline mobility level (pre-pandemic). They found three distinct patterns of societal reaction to social restrictions. In Australia (which implemented a near-complete lockdown) people did not go to their workplaces and stayed at home. In Sweden (which implemented relaxed lockdown with preserved workplace activity) the change in workplace mobility was smaller suggesting preservation of workplace activity. South Korea (which implemented minimal lockdown and preservation of workplace activity) changes in workplace mobility was even smaller. These findings suggest that populations respond directly to governmental interventions and highlights the important role of government decisions on day-to-day lives of its populations. Reference 22 analysed public policy effects on driving, transit and walking mobility behaviour, finding that they dropped to lower levels in Canada than the U.S. during March 2020 and show strong evidence that policy affects mobility behaviour. Reference 23 used a difference-in-difference approach and mobility data to study the causal impact of policies finding that statewide stay-at-home orders had the strongest causal impact on reducing social interactions and affects both mobility patterns and subsequent infection rates. Reference 24 quantifies the effect of U.S. state reopening policies on daily mobility levels, finding that four days after reopening mobility increased between 6% to 8%. Public policy decisions therefore can have a direct effect on the behavioural patterns of inhabitants mobility in regions directly affected by the policy.
Mobility and its relation to COVID19 incidence. Reference 25 show that mobility patterns are strongly correlated with decreased COVID-19 case growth rates in the USA. They collected data from 1st January 2020 to 20th April 2020. They define a Mobility Ratio (MR) to quantify the change in mobility patterns when compared to a baseline day before the pandemic when travel patterns were stable. They use this as a proxy for social distancing such that when an individual makes fewer trips, they interact less. They link mobility data with data on cases and construct a COVID-19 Growth Ratio (GR) in order to capture the complex and time-dependent dynamics between mobility and cases. They show a statistically significant correlation between their socially distancing metric and reductions in COVID-19 growth rates, showing that the effect of social distancing on case growth is not likely to be noticed for at least 9-12 days after implementation. Moreover 26 , extended their own research to include more regions and expanded the time horizon from 16th March 2020 to 16th September 2020. They found that the linear association between mobility and case growth rates previously observed is absent after April 2020 and that mobility has a less significant role in the transmission of COVID-19 than other adopted behavioural changes and NPIs such as wearing face masks, hand-washing, maintaining physical distance, avoiding large gatherings and closing schools. The strong association revealed in March-April is related to the adoption of NPIs in parallel and after an introduction of varying policies and changes in individuals, mobility behaviours confound the role of mobility. They conclude that using mobility data alone is likely to result in inaccurate models and forecasts and that there are more critical factors than mobility for controlling COVID-19. Reference 27 also identified a strong correlation between decreased mobility and reduced COVID-19 case growth during the period 27th March 2020 and 20th April 2020. They found that when they extended their time-horizon between 21st April 2020 and 24th May 2020 and later 25th May 2020 to 22nd July 2020, there was only a weak correlation between daily distance difference and case growth. They find that mobile phone data only captures a small component of the behaviours associated with social distancing and reduced case growth rates and that other NPIs such as wearing masks, maintaining distance are likely to be more important than mobility alone. Reference 28 used mobile phone data and modelled the relationship between mobility inflow and infections across counties in the U.S. between March 1 and June 9 2020. They found that travel between counties decreased 35% after the U.S. entered a lockdown but recovered rapidly once the lockdown started to ease. Using a system of equations they find a strong and positive relationship between mobility and the number of infections across counties with an average coefficient of 0.243, that is, a 10% increase in mobility is associated with a 2.34% increase in infections a week later. Reference 29 analysed mobility data in China to track population outflows from Wuhan between January 1st and January 24th 2020 and linked it with COVID-19 infection counts by location. They find a strong correlation between total population flow and the number of infections across regions. Reference 30 construct network maps of hourly movements of people to and from non-residential locations and apply a Scientific Reports | (2022) 12:9132 | https://doi.org/10.1038/s41598-022-11531-y www.nature.com/scientificreports/ meta-population SEIR model with susceptible (S), exposed (E), infectious (F) and recovered (R) compartments which tracks the trajectories of infections. Reference 31 used aggregated mobile phone data to build a SEIR model for the city Shenzhen, China. They simulate how the spread of COVID-19 changes when the type and magnitude of mobility restrictions vary. They found that reducing mobility by 20% delayed the epidemic peak by around 2 weeks with the peak incidence decreasing to 33%, a 40% reduction in mobility was associated with a delay of 4 weeks and reduced the peak by 66% and a 60% reduction was associated with a delayed peak of 14 weeks, decreasing the magnitude by 91%.

Methods
Mobility data records. Mobility data records come from a study conducted by the Ministerio de Transportes, Movilidad y Agenda Urbana (MITMA) https:// www. mitma. gob. es/ minis terio/ covid-19/ evolu cion-movil idad-big-data. The study collects data on mobility and distribution of the population in Spain from 13 million anonymised mobile-phone lines provided by a single mobile operator whose subscribers are evenly distributed. More specifically, the data is reported on a geographical layer composed of 2850 mobility zones across the whole of Spain where each mobility zone corresponds to a district or group of districts in densely populated areas, and municipalities or groups of municipalities in regions with low population density (see Fig. 1). The unit to measure mobility is the trip and the data contains the number of trips between and within mobility zones reported on an hourly basis. The start of a trip event is defined as when a user moves more than 500 m and the end of the trip is defined as when that user remains in an antenna coverage area for more than 20 min. Thus, a person moving from region A to region B and remains in region B for longer than 20 min is defined as a trip between regions. Moreover, a person who remains in region A and is connected to antenna A i then moves to antenna A j is defined as a trip within region. We collected mobility data from 01-09-2020 to 30-11-2020 for all 17 autonomous communities in Spain using the Flow-Maps systems. https:// github. com/ bsc-flowm aps/ Note: due to data handling issues, País Vasco and Castilla-La Mancha were omitted. Figure 1 shows the different mobility zones for all of Spain, the colours do not hold meaning and are simply coloured for easier illustration of the different MITMA areas. After retrieving the raw mobility data reported by the hour, we aggregate it by summing the total trips between and within mobility zones for each day, to obtain mobility reported on a daily basis. Using the daily mobility data we construct origin-destination (OD) matrices for each day t. Table 1 represent an example of OD matrix where rows correspond to the origin zone and the columns correspond to the destination zone and each entry corresponds to the total daily trips from origin to destination.
Following the example of Table 1, trips from A o to B d are associated with trips leaving MITMA zone A o and going to MITMA zone B d and are the outgoing trips from MITMA zone A o . Incoming trips are defined as the opposite and thus trips which began at B d and went to A o . The trips going from A o to A d are internal trips (the diagonal of the matrix) or the number of trips within MITMA zone A . Furthermore, for each mobility zone on a given date we define the following four mobility indexes: (i) incoming mobility; (ii) outgoing mobility; (iii) internal mobility; and (iv) the total mobility. More formally, the incoming data takes the column sums of the OD matrix minus the diagonal, the outgoing data takes the row sums of the OD matrix minus the diagonal, the internal is just the diagonal and the total is incoming + outgoing + internal.
Linear regression. We first run a linear regression model in order to see the effects on weekday and weekend mobility. We expect to see that the closure of bars and restaurants will have a greater effect on the weekend than on the weekdays since mobility is largely unaffected by people commuting to and from work but is affected by peoples decision to go and socialise on the weekends. We estimate the following for policy status j, at time t: where, outcome jt is mobility for MITMA region j at time t, β 0 is the intercept of the existing level at point 0, β 1 gives the existing trend in mobility before the policy, level takes on a value of 0 before the policy and 1 after the policy with its β 2 capturing the impact of the policy, β 3 captures the change in trend.

Difference-in-differences (DiD).
We use a Difference-in-difference technique to infer the causal impact of a policy. That is, with two groups and two periods, the DiD estimator is defined as the difference in average outcome in the treatment group before and after treatment minus the difference in average outcome in the control group before and after treatment. Here the outcome variable is a daily time-series of mobility data and therefore we are considering mobility before the policy intervention and mobility after the policy intervention. Consider  Table 2) in which Cataluña corresponds to the treatment group, denoted as T and other CCAAs, for instance, Madrid corresponds to the control group, denoted as C. Therefore, ȳ T 0 corresponds to the average mobility for a given mobility zone in Cataluña before the policy, ȳ T 1 the average mobility in a given mobility zone in Cataluña after the policy and ȳ C 1 and ȳ T 0 corresponds to the average mobility in a given mobility zone in a CCAA (excluding Cataluña) before and after the policy. Thus, the causal impact of the policy can be given as the difference between these two differences. It is important to note that the data is normalised by; x , bounded by [0, 1].
Bayesian structural time-series. We use Bayesian structural time series to infer the causal impact of a policy by explicitly modelling the counter-factual observed before and after an intervention (see 32,33 ). The model assumes that the outcome time series can be explained in terms of a set of control time series that were themselves not affected by the intervention. Furthermore, the relation between the treated series and control series is assumed to be stable during the post-intervention period. This allows us to generalise the results obtained using the difference-in-difference approach. The model is first estimated using pre-intervention data, and then an intervention period occurs and the model then tries to predict the post-intervention period. The difference between the prediction and the observed data can be thought of as the causal impact of the policy. For the preintervention period, we take mobility before the policy and for the post-intervention period, we take mobility after the policy. Consider the following state-space model for time-series data.
In which, ǫ t N(0, σ 2 t ) and η t N(0, Q t ) are independent of all other unknowns. The observation equation (3) links the observed data y t to a latent d-dimensional state vector α t and the state equation (4) governs the evolution of the state vector alpha t through time. y t is a scalar, Z t is a d-dimensional output vector, T t is a d × d transition matrix, R t is a d × q control matrix, ǫ t is a scalar observation error with noise variance σ t and η t is a q-dimensional system error with a q × q state-diffusion matrix Q t , where qd. More details can be found in 33 .

Results
This section discusses the results in which we assess the effect of NPIs on mobility. We compare the mobility in Cataluña with that of other regions in Spain where other policies or no policies have been applied. We apply an OLS regression, difference-in-difference and Bayesian structural time series model in order to try and capture the causal effect of this policy intervention. We first report some analysis of the mobility data and OD matrices for Cataluña and Madrid, two of the most populous regions in Spain and where one region introduced the policy whereas the other did not. We then report the main results of the paper from the different econometric models. Finally, we try to link the policy to reductions in COVID19 growth rates.

Analysis of mobility.
We first report summary statistics for Cataluña and Madrid for each mobility type, Table 3. We can observe that the average mobility levels fall for each of the mobility types for Cataluña when comparing before and after the policy intervention. We additionally observe that the mobility levels increased from before and after the policy intervention in Cataluña for Madrid. However, this increase has no relation to the policy in Cataluña. Figure 2 shows that the number of cases were falling in Madrid before Cataluña introduced the policy and thus mobility began to increase when the risk of contagion was falling in Madrid.(Since no one single region can act as a perfect control group for Cataluña due to varying differences between the policies, socio-demographic, population etc. we use all Autonomous Communities in our analysis and control groups). Figure 3 shows a time series plot of the total number of trips for Cataluña and Madrid. The policy came into effect right in the centre of the time series and lasted until the end of the time series (shaded as yellow) thus we see a shift down in the number of trips for Cataluña after the introduction of the policy and it remained down for the duration of the policy. Contrasted with Madrid in which we see no shift downward, which is, as we would expect since Madrid did not introduce the same policy as Cataluña. Figure 4 shows the mobility matrices for Cataluña and Madrid for each day for incoming trips into MITMA regions.(MITMA: Ministry of Transportation, Mobility and Urban Agenda). The darker colours represent higher mobility and the lighter regions lower mobility. The plot has been arranged to show MITMA regions with overall Table 2. Diff-in-diff: difference-in-difference estimates.

Pre
Post Post-pre difference Scientific Reports | (2022) 12:9132 | https://doi.org/10.1038/s41598-022-11531-y www.nature.com/scientificreports/ higher mobility on the right and MITMA regions with overall lower mobility on the left. There are distinct periods of higher and lower mobility, the higher mobility corresponds to weekdays and thus this mobility incorporates people travelling to and from work and the lighter periods correspond to weekends where people travel less. Additionally, there are other days of lower mobility. At point (A) in Fig. 4 there are 3 days of lower mobility, 2 of these days are a weekend and the 3rd day is a public holiday in Cataluña falling on a Friday, 2020-09-11 (Fiesta Nacional de Cataluña). The same region in the Madrid plot does not show this reduced mobility since the public holiday is only specific to Cataluña and thus not Madrid, therefore people in Madrid went to work on this day. There are other public holidays in this data, specific to Madrid which fell on Monday, 2020-11-01 (Fiesta de Todos los Santos) and Monday, 2020-11-09 (Fiesta de la Almudena) which are shown in the mobility matrix data as (C) and (D) respectively. Finally, there is a public holiday which is celebrated across the whole of Spain, on the 2020-10-12 (Fiesta Nacional de España) denoted as point (B). Therefore, mobility for this day was reduced in both Cataluña and Madrid along with other Autonomous Communities. The next sections report the main findings of this study.
Linear regression. We first applied a series of linear regression models. Firstly, using all MITMA zones in Cataluña and then secondly, we aggregate the data in order to analyse the effect for the whole of Cataluña.   We see a reduction in activities which imply the movement of people. We expect to see less of a drop during the weekdays than on weekends since people still need to go to work. However, the 3.9% drop during the weekdays may suggest that people are socialising less in bars and restaurants after work and are going straight home. After the policy impact, there appears to be a sustained drop in daily mobility of 6.9% after controlling for weekend effects and therefore the policy caused a shock to mobility and then levelled out, indicated by the trend of 0.1% afterwards. Therefore, using a regression model, we have aimed to quantify the drop in mobility as seen in the Cataluña panel of Fig. 3. Linear regression aggregated MITMA zones (Cataluña). In order to visualise and illustrate the previous regression results in a more intuitive way, We aggregate the MITMA regions for all of Cataluña into daily mobility totals and thus we have a single observation for each day for Cataluña. Figure 5 plots the regression fitted values for the incoming mobility type before and after the policy with and without a weekday control variable. We note that the level variable of interest for Panel (A) is not statistically significant but the level variable in Panel (B) is   Table 4. OLS regression results: the Weekend regression is only run on the data points on the weekend whereas the Weekday regression is only run on the data points on the weekdays. We introduce a weekday control, measuring both the weekends and weekdays mobility. * p<0.1; * * p<0.05; * * * p<0.01. Difference-in-differences. We next apply a difference-in-difference model in order to additionally quantify the effect on mobility after the policy of closing the bars and restaurants in Cataluña was introduced. A more formal description of the model is left to "Difference-in-differences (DiD)". Firstly, we use Madrid as a control group in order to visualise the difference-in-difference model since Madrid did not introduce the same policy of closing the bars and restaurants as in Cataluña, reported in Fig. 6. Finally, we present the results using each autonomous community.. We applied the difference-in-difference model to the four different mobility types incoming, internal, outgoing and total and extended the control groups to each autonomous community in Spain. (Note: due to data handling issues, País Vasco and Castilla-La Mancha were omited.) Table 5 reports the difference-in-difference estimators for each mobility type and autonomous community. Interestingly, the internal mobility does not show statistically significant results for all but three CCAAs (i.e. the cells in the table are not coloured). Recall, that internal mobility corresponds to the movements of people within a given MITMA zone. One economic interpretation could be that since bars and restaurants were closed people chose not to travel as much to different districts to socialise. For internal mobility, people still went about their daily business closer to their home, i.e. continue to shop at the supermarket and go to the pharmacy etc. but they had less incentive to travel to different districts (or MITMA regions) to meet with friends and family (which may be the reason for the non-statistically significant results in the internal mobility column). This result can be further seen in Fig. 12 in the Appendix which plots a 14-day rolling moving average for 4 CCAA's, the internal plot shows relatively flat lines for mobility, with Cataluña's internal mobility only dropping slightly, whereas incoming and outgoing showed a much steeper drop for Cataluña relative to the other CCAA's. Therefore, people still moved within their MITMA regions in order to go about their daily lives but after the policy people stopped migrating to other regions and therefore longer distance mobility habits changed.   The difference can be seen as the reduction in mobility. Panel (B) also shows the fitted values before and after the policy. The darker green line represents the fitted values before the policy and the transparent green line after the policy represents the fitted values had the policy not been introduced. The orange line represents the fitted values on the data points after the policy was introduced and the difference can be seen as the reduction in mobility due to the policy. In both panels, the blue points represent the observed values of the mobility levels on a daily basis. That is, the data started on Tuesday, September 1st 2020 and thus the first four points correspond to mobility levels from Tuesday 1st to Friday 4th. We then see two points significantly drop off on the weekend. The following week, we observe three significantly lower points. One of the data-points being Friday 11th September, also shown in Fig. 4 previously (2020-09-11, Fiesta Nacional de Cataluña) and the other two being the weekend reduction in mobility. On the x-axis, dic refers to Diciembre (December).  Figure 7 shows the Bayesian structural time-series model in which the model uses Madrid's mobility data in order to build a predictive model for Cataluña's mobility. A more formal description of the model is left to "Linear regression". The model is able to build a very predictive model during the training phase. The large spike downwards in the pointwise panel, at the date 2020-09-11 in the training phase, corresponds to La Diada Nacional de Catalunya a public holiday specific only to Cataluña and thus affects Cataluña mobility and not Madrid (Fig. 4 shows that on public holidays mobility is reduced.). Aside from this single day, the differences between the predictive model and the observed data lie in and around zero before the policy is implemented. After the introduction of the policy, these differences become negative except on two dates 2020-11-02 and 2020-11-09 which correspond to Día de todos los Santos and Virgen de la Almudena two public holidays specific to Madrid and not Cataluña. We find that in absolute terms during the post-policy period, the average mobility in Cataluña was 12.39 million trips, in the absence of the policy the model expected an average of 14.19 million trips with a 95% confidence interval of (13.79, 14.58) million. The causal effect is 12.39 − 14.19 = −1.8 million reduction in the number of trips in Cataluña with a 95% confidence interval (− 2.19, − 1.39) million. In relative terms, mobility reduced by − 13% with a 95% confidence interval (− 15%, − 10%) which suggests that this causal effect in the reduction in mobility is statistically significant with a posterior Bayesian one-sided tail-area probability of 0.001. This finding is consistent with the linear regression and difference-in-difference models presented earlier. Figure 7 reports the model for Cataluña and Madrid, Fig. 13 in the Appendix reports a case where the model did not work using Andalucia in place of Madrid as a control group. Finally, Table 6 shows the results for all CCAAs for the incoming mobility type, the internal, outgoing and total are left to the appendix. Tables 7, 8 and 9 reports the tables for the Internal, Outgoing and Total mobility respectively. Mobility and COVID19. Thus far, this paper has shown that the policy has had a clear effect on the reduction of mobility in Cataluña. In this section we aim to assess the policies effect on daily COVID19 incidence, however, some comments are required.

Madrid
Quantifying the policies impact on reductions in the number of new cases is more difficult than quantifying its impact on mobility. That is, mobility levels over time are roughly constant and predictable which allows us to build a strong counter-factual and thus when a policy shock has been implemented we can use this counterfactual in order to measure the real causal impact of that policy. Using new COVID19 case data presents more difficult problems since outbreaks are unpredictable and different regions experience fluctuations and peaks in COVID19 cases at different times, this renders the models used in this paper obsolete when applied to COVID19 case data. Consider Fig. 2 which plots the number of new cases over time since the start of the pandemic. There are three regions with high peaks, Madrid experienced its peak before that of Cataluña and would therefore not make a good counter-factual after the policy. Andalucía and Cataluña peaked during the same period however, both CCAAs introduced strict countermeasures to reduce the number of cases and both their number of new cases fell at the same time and therefore Andalucía would not make a suitable counter-factual either.  where C i is the number of cases for an autonomous community i at time t. The GR-ratio aims to smooth out the case incidence data by using a 14 and 21 day rolling average. From here the ratio between the two smoothed average incidence data is taken to try to capture the growth in the number of cases over time. Figure 8 shows the growth rate ratio for the number of cases over time for all autonomous communities. The growth rate began to increase in the month leading up to when regional governments introduced different policies. The growth rate in Cataluña began to decrease after an initial 2-3 week lag once the policy came into effect. Figure 14 shows that the growth rate in the number of cases began to increase in December, before Christmas and when the policy was not enforced. Figure 9 shows the same calculation as in Eq. (5) but applied to the mobility time series. The growth rate ratio for mobility began to increase from September for many CCAAs, this could be attributed to people returning from their second homes after the summer break.
We next anticipate that there is a lag between the movement of people and becoming infected with COVID19, due to the lag in symptoms, reporting etc. We firstly normalise the mobility data, fixed to the first week of October as follows in Eq. (6). The first week of October is the first week in our sample data. Moreover, the data is normalised such that, each Monday is normalised to the first Monday of October, each Tuesday is normalised to the first Tuesday of October and so on. Table 5. Diff-in-diff estimates for mobility type: the mobility types incoming and outgoing appear correlated with all of the coefficients being similar to each other. Since these two mobility types are correlated, the total is also somewhat correlated. The CCAAs whose coefficients were statistically significant at the 0.1% level ranged between − 9 and − 17.5% and between − 8.9 and − 17.4% for incoming and outgoing respectively suggesting that the policy reduced mobility anywhere between − 9 and − 17.5% depending on the control group used. Total mobility at the same significance level ranged from − 8.6 to − 17.1%. These findings indicate that the policy had a real and direct effect on reducing the movement of people across MITMA regions. where t o corresponds to the mobility data for each corresponding weekday at the beginning of October. Next, we compute the Pearson correlation coefficients with a 95% confidence interval for lags 1:40 as shown in Fig. 10. We finally take the optimal lag of 21 and plot the scatter-plot between the growth rate ratio defined in Eq. (5) and the normalised mobility defined in Eq. (6). The results can be seen in Fig. 11. There is an association between increased mobility and an increase in the growth rate ratio.

Conclusion
This paper quantifies the impact of mobility on the closing of bars and restaurants in Cataluña using a number of causal inference models. Overall we find that this policy reduced mobility in Cataluña, but it not only reduced mobility, it caused people to change the way they behave and respond to the policy, with trips across regions being affected more than trips within regions. This finding is significant since people did not substitute meeting friends and family at bars and restaurants with meeting friends and family in other locations, people simply stayed within their own MITMA region, reducing their overall mobility and changed the way in which people went about their daily business. The findings of this paper are 3-fold. (a) we quantify the impact of the policy of closing bars and restaurants down during a 5 week period on mobility. (b) we find that people changed their  Correlations between GR Ratio and Normalised Mobility Figure 10. Normalised correlations: relationship between the growth rate ratio and mobility-normalised. The optimal lag occurs around 21 days. This may account for lags in the time people notice symptoms of COVID19 after some days and lags in the reporting of cases for each regional healthcare systems. Scatter−plot between the GR Ratio and Normalised Mobility Figure 11. Scatter-plot between the growth rate ratio in the number of cases and the normalised mobility: there appears to be some relationship between an increase in mobility and an increase in the growth rate ratio for the optimal lag of 21 days. www.nature.com/scientificreports/ When the GR-ratio is not correlated with mobility. Figure 8 of the manuscript shows that a number of Autonomous Communities saw a reduction in the number of reported cases. The drop-off in the number of cases can also be seen in Fig. 1 Fig. 7 in which we report a case where the control was not suitable for Cataluña. The pairwise differences hover around zero and the cumulative mobility returns to zero suggesting that causal inference cannot be inferred from this control group. As Table 6 shows, Andalucia is not statistically significant and the sign of the relative effect is positive, not negative as one would expect. On the x-axis, dic refers to Diciembre (December).
Scientific Reports | (2022) 12:9132 | https://doi.org/10.1038/s41598-022-11531-y www.nature.com/scientificreports/ aimed to explain the relationship between the fall (rise) in the number of cases and a low (high) level of mobility. There does appear to be a correlation between the two variables for the case of Cataluña.
In order to illustrate this point we ran the analysis again using the Community of Valencia in place of Cataluña. Figures 15 and 16 in this document shows the same analysis as Figs. 10 and 11 in the manuscript. For the Community of Valencia we can see that the relationship between the GR Ratio and Mobility are not statistically significant for all but two lagged periods. The Community of Valencia did not impose any policies during this time period which affected mobility levels. Figure 16 also shows no relationship between mobility and the GR ratio.
As previously discussed, not all regions introduced a policy. For example, Madrid did not introduce any policies at the time Cataluña introduced their policy of closing down leisure activities. Therefore, Madrid can be considered a suitable control group in terms of not having a policy Madrid also has other favorable characteristics   www.nature.com/scientificreports/ to make it a suitable control group, such as, population levels, economic stability etc.). Moreover, other regions such as Andalucia introduced a policy to try to mitigate the spread of COVID19 in the region. This region turns out to be a poor control group in the Bayesian Structural Time-Series model given in Fig. 13. Moreover, this region is not statistically significant in Table 5 (DIff-in-Diff model) and Table 6 (Bayesian Structural Time-Series). Correlations between GR Ratio and Normalised Mobility