The question of how and when to open up the economy as Covid-19 rates drop is fraught with great risk on both sides. How the question is answered is likely the most critical public policy decision in the last few decades. Unfortunately, most of the arguments made so far have been based more on philosophy than science. Reliable data has been sparse, but modern technology provides opportunities to make quantitative arguments.
This paper attempts to find relationships between Covid-19 infection rates in the United States and mobility data collected from mobile devices. Google has recently made this Mobility Data publically available for use in research on the Virus. The analysis demonstrates that Google Mobility Data is a reasonable proxy for social interaction that correlates significantly with infection rates.
The Google reports utilize aggregated, anonymized global data from mobile devices to quantify geographic movement trends over time across 6 area categories: retail and recreation, groceries and pharmacies, parks, transit stations, workplaces, and residential areas (1). The numbers are percentages that represent changes above or below the long term trend. The data exist for 131 countries and regions, but I am using only data for the United States in order to compare with relatively consistent Covid-19 epidemiological data. Google’s definitions of the area categories are in Table 1.
Table 1. Mobility area category definitions
The U.S. aggregates since February 15 are shown below. A couple of interesting things to note: Grocery and Pharmacy spiked up in early March, as people stocked up for the lockdowns and “Social Distancing”. Parks and Retail/recreation did also though to a lesser extent, suggesting people wanted to carry out these activities before lockdowns were put in place. Workplaces and Residential are clearly inversely correlated, as workplaces shut down people spent more time travelling near the home.
Figure 1. U.S. aggregate mobility by date since Feb. 15 for 6 different area categories.
The New York Times has published State and County level data to github (2). Cases and Deaths are cumulative by Date, going back to Washington State on 1/21/2020. Figure 2 shows cumulative cases for 4 counties, Westchester (NY), Los Angeles (CA), Dallas (TX), and Snohomish (WA). Note that because the cases are cumulative, no new cases are being added when the slope becomes horizontal. Snohomish and Westchester are closer to this than Los Angeles and Dallas, which experienced later onsets of the disease.
The data represent verified cases only. Some recent antibody studies in Germany, Norway, and The United States suggest that as many as 20% of certain populations have already been infected by the virus (3). In that light, the numbers being used here are almost certainly a significant underrepresentation, but they are useful for two reasons:
Death counts are likely far less ambiguous than case counts, and it is possible to do this analysis with them, but the data for deaths is also far more sparse and more truncated, as it is usually 1-2 weeks from diagnosis to mortality. This leads to more numerical problems in regressing the data.
Figure 2. Cumulative Covid-19 cases in 4 representative U.S. counties.
In order to tie the Mobility data to outcomes, we need robust metrics to represent each. A simple dependent variable is simply the percent increase in cases over a specific time period. This is unstable in the early days of the viral spread, when case counts are low in a specific county, but can be regularized by weighting the regression on the number of cases. The choice of a “lookahead window” is somewhat subjective, you need one long enough to capture any changes influenced by mobility, but if it is too long you truncate your data. According to the CDC, people who get symptoms nearly always do so in the first 2-14 days (4), with the 97.5% experiencing symptoms in the first 11.5 days (6), so a 12 day lookahead is probably adequate to compute the percent increase.
Defining the Independent variable is also somewhat subjective. I chose to look at Mobility for the 12 days leading up to the lookahead, but filter it with a 12 period Gaussian (mean = 3, sd = 2.0) (Figure 3). This gives the greatest weight to mobility 7-9 days before the lookahead, and slowly deprecates the effects to nearly zero a couple days before the window.
Combining the datasets above produced 47,847 rows of data, of which 20,609 were removed because of missing mobility values. To this data I added several time-independent covariates from the U.S. Census data (5) which are sometimes associated with variance in epidemiology:
Lastly, I added an independent variable measuring the previous 5 days viral growth rate. This allows the model to make more accurate projections of the growth rates 12 days into the future.
The regression results are shown in Table 2 below. The model has an R Squared of 0.596, meaning that most of the results are explained by these covariates, although their individual contributions vary significantly.
All of the covariates except for “PctAsian” are significant beyond the 99% confidence level. It is clear that the most predictive is “PreviousFiveDaysPctChangeCases”, which just means the future slope of the curve is related to the current slope for each county.
By changing one variable at a time while holding the others constant, we get an estimate of the influence of the time dependent covariates (Table 3) and the time independent ones (Table 4).
Table 3. Time dependent covariates and their predicted effects on infection rates.
Among the mobility variables, the strongest predictor of increase in infection rate is mobility around the workplace, followed closely by mobility around retail and recreation areas. Curiously, Residential mobility was third, suggesting that lockdowns and “sheltering in place” measures are not as effective as suggested, or are at least are being sabotaged by some amount of interaction with housemates or friends/neighbors. This may explain why Sweden, which did not enforce strict lockdowns, has not had significantly higher rates of infection than other European countries. The model also suggests that greater mobility in the areas of grocery/pharmacy and parks/recreation would not increase infection rates. The one thing that Retail/Recreation (which includes bars, restaurants, concerts, etc.) and Workplaces have in common is close social interaction, which Parks and Grocery stores have less of. This suggests it may be more common to get the virus from respiration rather than touching it.
Most of the time-independent factors seem to have very little influence on rates of infection. The exceptions are Latitude, which might suggest warmer weather has a small effect, as does persons per household (this is not surprising), and the percentage of foreign born in the county (possibly due to more visitors from their native countries). Race does not seem to have a large effect, nor does income. It is widely known that those over 65 are more at risk of death from Covid-19, but as far as infection rates goes it appears that having a large percentage of seniors in the county is a slight deterrent, possibly because they take the social distancing guidelines more seriously.
Table 4. Time independent covariates from Census data and their predicted effects on infection rates.
Table 5 contains a county-by-county breakdown of weighted average mobility trends and the projected changes in cumulative infected rates for the Baseline scenario (current status quo), Scenario 1 (returning to 50% of historical mobility), and Scenario 2 (returning to 100% of historical mobility). The most populous 30 counties in the U.S. are shown. A change of 200% in infection rate represents a doubling of cumulative cases over the 12 day lookahead period. Because 2 weeks is roughly the time it takes for an infected patient to either die or recover, a 200% growth rate is roughly keeping a constant rate of infection.
Table 5. Mobility and predicted 12 day infection growth rates (last 3 columns) as of May 1, 2020. The Baseline projections are for 12 days in the future with current mobility and can be compared with Scenario 1 (return 50% to long term mobility) and Scenario 2 (returning 100% to long term mobility).
It should be noted that these projections are based on pre-Covid19 norms of social contact, and do not take into account mitigation like social distancing. For example, it is probably possible to return to historical norms in the workplace without dramatically increasing infection rates if social distancing is used and large meetings are avoided.
Because Mobility can be a proxy for social interaction, it is clearly a significant factor in the transmission of Covid-19. Using Google’s mobility data allows us to see the relationships between mobility in different geographical areas and their corresponding increase in infection rates. Regressing the data suggests that it is possible to achieve previous levels of mobility but doing so must be undertaken with caution and mitigation, especially in the workplace and in retail/entertainment venues. Future work can utilize the Global dataset in order to see correlations by country.
Comment
Thank you for doing this work and for sharing it!
I have one question about "The model also suggests that greater mobility in the areas of grocery/pharmacy and parks/recreation would not increase infection rates. ": do you refer to lower correlation between these data and the covid-19 cases or do you actually mean a cause-effect?
Kind Regards,
Hugo
Hi Sana,
The trends are definitely upward because this is a cumulative rate of infection. If you plot the new daily cases (cases[n] - cases[n-1]) you will see peaks for most counties. As far as modeling goes though you can still measure the slope, and the differences in the slope vs. time, which is what I related to the mobility (with a 12-18 day time lag).
I plotted the infection rate and it seems to have a pretty steady upward trend. How did you manage to see the impact mobility has on the infection rate if the trend shows no change across a lot of days?
Sorry it took me so long to get back with you. I emailed the data I regressed on. The best performing model I found to be a RandomForest, closely followed by Light Gradient Boosted Trees. A couple of things to keep in mind, here are the features I used:
pct_chg_cases ~ retail_and_recreation_percent_change_from_baseline_score
+ grocery_and_pharmacy_percent_change_from_baseline_score
+ parks_percent_change_from_baseline_score
+ transit_stations_percent_change_from_baseline_score
+ workplaces_percent_change_from_baseline_score
+ residential_percent_change_from_baseline_score
I weighted the regression by "current_cases" because the rows with very few cases (small counties early in the pandemic) tend to have very high variance. You can experiment with using an exponent other than one to improve performance. Also, because there are many rows for each county (one per day) many of the rows look very similar and it is possible to get target leakage. That is not a problem with something like linear regression, but with a tree-based method which has many degrees of freedom, it is definitely a problem. I used 5-fold cross validation and grouped all rows for a given county in the same fold to prevent any leakage. I originally compiled this data about 3 weeks ago, the data sources have been updated since then, it would be great to update the regression also.
Hi Paul I don't know how much the datasets are secret that people publish their datasets on the GitHub.
Hi Paul,
Have you performed a polynomial linear regression or just a basic one? Im not sure but wouldn't a polynomial one fare better in this case?
Im confused as to how exactly you constructed the Gaussian filter. Would you mind giving me more details on it. Also I would really appreciate it if you could also provide me the manipulated data after you applied the Gaussian filter, if it's not too much trouble. My email: [email protected]
Hi Paul,
My concern was also linear regression. I. really appreciate that if you give me the final data ( manipulated data) to play with.
Also explaining the Gaussian filtering. Here is my email
Thanks for your suggestions Patrick - I like the suggestion about the control set although what I normally do is regress on 10 datasets where I randomly mix the dependent variables, and in this case I got no better than 0.07 RSquared.
The limitations of Google's data are spelled out on their URL. I am not sure about the accuracy beyond that, but when trying to glean information about Coronavirus infection rates, the question has to be asked, compared to what? It is very difficult to find anything beyond anecdotal data.
The choice of linear regression has to do with what I was looking for. I assumed when it came to mobility around certain potential contact areas, there was a proportional relationship. I was more interested in finding which factors were the most robust predictors than simply fitting a tree based model to every inflection of the data, which could be deceptive where the data is sparse. That said, I did build GBT and RF models with better fits, but similar relationships between the variables.
(county level; not state level data; just re-read)
I suppose I am quite a bit more cautious about the data sources. Even Google's tracking. But a few questions/comments:
Is this OLS / linear regression? Probably not the best way.
But a valiant effort at data integration, etc. State level time series for 8 weeks.
A few suggestions:
1. Put in dummy variables for each state, perhaps based on their policy reactions (if any)?
2 Grab the CDC weekly mortality data from prior years. That would be an interesting control.
© 2020 Data Science Central ® Powered by
Badges | Report an Issue | Privacy Policy | Terms of Service
Upcoming DSC Webinar
Most Popular Content on DSC
To not miss this type of content in the future, subscribe to our newsletter.
Other popular resources
Archives: 2008-2014 | 2015-2016 | 2017-2019 | Book 1 | Book 2 | More
Upcoming DSC Webinar
Most popular articles
You need to be a member of Data Science Central to add comments!
Join Data Science Central