Association between mobility, non-pharmaceutical interventions, and COVID-19 transmission in Ghana: a modelling study using mobile phone data

Background: Governments around the world have implemented non-pharmaceutical interventions to limit the transmission of COVID-19. While lockdowns and physical distancing have proven effective for reducing COVID-19 transmission, there is still limited understanding of how NPI measures are reflected in indicators of human mobility. Further, there is a lack of understanding about how findings from high-income settings correspond to low and middle-income contexts. Methods: In this study, we assess the relationship between indicators of human mobility, NPIs, and estimates of Rt, a real-time measure of the intensity of COVID-19 transmission. We construct a multilevel generalised linear mixed model, combining local disease surveillance data from subnational districts of Ghana with the timing of NPIs and indicators of human mobility from Google and Vodafone Ghana. Findings: We observe a relationship between reductions in human mobility and decreases in Rt during the early stages of the COVID-19 epidemic in Ghana. We find that the strength of this relationship varies through time, decreasing after the most stringent period of interventions in the early epidemic. Interpretation: Our findings demonstrate how the association of NPI and mobility indicators with COVID-19 transmission may vary through time. Further, we demonstrate the utility of combining local disease surveillance data with large scale human mobility data to augment existing surveillance capacity and monitor the impact of NPI policies.


Research in Context
Evidence before this study We searched PubMed and preprint archives for articles published in English that contained information about the COVID-19 pandemic published up to Nov 1, 2021, using the search terms "coronavirus", "CoV", "COVID-19", "mobility", "movement", and "flow".
The data thus far suggests that NPI measures including physical distancing, reduction of travel, and use of personal protective equipment have been demonstrated to reduce COVID- 19 transmission. Much of the existing research focuses on comparisons of NPI stringency with COVID-19 transmission among different high-income countries, or on high-income countries, leaving critical questions about the applicability of these findings to low-and middle-income settings.

Added value of this study
We used a detailed COVID-19 surveillance dataset from Ghana, and unique high resolution spatial data on human mobility from Vodafone Ghana as well as Google smartphone GPS location data. We show how human mobility and NPI stringency were associated with changes in the effective reproduction number (R t ). We further demonstrate how this association was strongest in the early COVID-19 outbreak in Ghana, decreasing after the relaxation of national restrictions.

Implications of all the available evidence
The change in association between human mobility, NPI stringency, and R t may reflect a Introduction Nations around the world introduced a range of non-pharmaceutical interventions (NPIs) to limit the spread of COVID-19 in the early phases of the epidemic [1]. These NPIs have been diverse, and have included the use of personal protective measures, environmental measures, physical distancing, restricting movement, and limiting the gathering of people.
NPIs have been implemented at different times in relation to the progression of local and national disease outbreaks, with some put in place before transmission was established, and others reactive to rises in cases. NPI measures have also overlapped one another in the timing of their application [1,2]. Previous research has attempted to quantify the relative effectiveness on COVID-19 transmission of different NPIs [3][4][5][6], but modelling the impact of . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.  [7].
One approach used by researchers and policymakers to measure the impact of NPIs during the COVID-19 pandemic has been to observe changes in measurements of human behaviour under individual interventions or under a combination of interventions [8][9][10][11][12].
Perhaps the most common way to quantify varying patterns of human behaviour is the use of human mobility datasets, which measure the locations of individuals using GPS or Call Detail Records (CDRs) [13,14]. These mobility datasets have been made available by a variety of network service and mobile application providers [15][16][17]. Mobility data has been used widely during the COVID-19 pandemic to predict the introduction of COVID-19 cases, and to monitor and estimate adherence to NPIs including travel restrictions [8,9,[18][19][20], but questions remain about how patterns of mobility and NPI stringency relate to transmission in LMIC settings.
Previous research has been conducted in Africa on the implications of mobility patterns for transmission of infections other than COVID-19 [21,22] and during the COVID-19 epidemic, analysis of movement patterns in Ghana has been conducted to inform policy makers about the volume of reductions coinciding with lockdown interventions in Accra and Kumasi [23].
These indicators may be used as a proxy for social contact [13] and therefore, for potential COVID-19 transmission, although the "link" between movement and disease transmission may decrease due to greater adherence to social distancing or personal protective equipment guidelines [24].There remain questions about how mobility indicators can be used to estimate COVID-19 transmission and how these indicators reflect behavioural responses to NPI measures, particularly in an LMIC context. In this paper, we combine human mobility and NPI data to estimate their relationship over time to the progression of the COVID-19 epidemic in Ghana.
Here, we used surveillance data in a sample of 27 districts collected by the Ghana Health Service on PCR confirmed COVID-19 patients at the district level (administrative level 2, 261 total districts) in 11 of the 16 regions of Ghana between March and December 2020 to produce individual estimates of R t , the real-time reproduction number, for individual districts.
R t is a time-varying parameter describing the average number of infections derived from a . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted November 17, 2021. ; https://doi.org/10.1101/2021.11.01.21265660 doi: medRxiv preprint single infection and indicates whether an epidemic is growing (R t > 1) or decreasing (R t < 1).
We combined estimates of R t with subnational human movement data from Vodafone Ghana and Google. These data measure the volume of movement activity (a proxy for social contact), in each district. We then modelled the relationship between human mobility indicators and R t using a multilevel generalised linear mixed model to assess whether changes in mobility and NPI stringency were associated with changes in R t in during the COVID-19 epidemic in Ghana.

Study Setting
The first cases of COVID-19 in Ghana were reported on 12th March 2020 [25]. These cases were detected in Accra, the capital city of Ghana and were imported via international travel [25]. Following the announcement of the first COVID-19 cases, the Ghanaian government announced the suspension of international travel and the closure of land borders to reduce the risk of further introduction [26]. Domestic case numbers grew in March and April 2020, leading to the closure of universities and high schools and the announcement of a partial lockdown in the Ashanti and Greater Accra regions, the two most populous regions of Ghana [26]. This lockdown introduced a stay at home order except for essential travel including shopping, healthcare, and use of public toilets. Almost all COVID-19 NPI restrictions were lifted by July, although restrictions on international travel and mandated use of facemasks remained in place until September 2020.

COVID-19 surveillance data
We used line list data recording lab-confirmed COVID-19 cases collected by Ghana Health Services each day between 12th March and 1st September 2020. These data were collected in 11 of the 16 regions in Ghana (administrative level 1: excluding "Ahafo", "Bono", "Upper West", "Volta", and "Western North" districts due to limited detail of data in these districts).
Patient-level records were referenced to a standard spatial reference provided by the Ghana Statistical Service using patients' reported district of residence. Using the date of case confirmation and the district of residence, we aggregated individual records into daily case counts of confirmed COVID-19 cases per district (Supplemental Figure 1). Through visual inspection, we replaced three outliers in reporting in two districts with a linear interpolation between the preceding and following records. For these outliers, the number of cases . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted November 17, 2021. ; https://doi.org/10.1101/2021.11.01.21265660 doi: medRxiv preprint reported in a district clearly exceeded the overall trend of case reporting (reported cases greater than 5x higher than all previous reports) (Supplemental Table 1). We assumed that these records reflected "late reporting" with samples collected on multiple days reported on the same date. This assumption potentially underestimates the total number of COVID-19 cases in these two districts, but it is not possible to approximate when the cases reported in these intervals may have been originally tested.

Defining stringency indices for NPIs
Data on the dates of NPIs implemented in Ghana were provided by Ghana Health Service detailing the starting dates of public health interventions. We used this to define the start dates of interventions and augmented it with available news sources and government press releases to create a dataset of the start and end dates of nine intervention measures. Using these intervention data, we constructed a stringency index to measure the stringency of COVID-19 interventions through time, defined daily as the number of active interventions divided by the total number of interventions. This stringency index assigns a uniform level of stringency to each intervention measure and records the length of time that the measure was implemented. We also used the OxCGRT stringency index, calculated from a global database of NPIs, which is used to construct a stringency index based on a taxonomy of government interventions [1] (Supplemental Figure 2). We used the most recent version of the stringency index (as of June 2021), rather than the OxCGRT "legacy stringency index." OxCGRT data also records nine interventions resulting in a change in the stringency index in the study period. We extracted the date of maximal intervention from both indices to compare both stringency indicators, defined as the first date with the highest stringency for each index. Although both sources of intervention data reported interventions at a national level, intervention measures were introduced at different spatial scales in Ghana. School closures (including different educational tracks) and mask mandates, for example, were implemented nationally, while partial lockdown measures were introduced only in Ashanti and Greater Accra regions.

R t estimation
To ensure data coverage during the early stages of the COVID-19 outbreak in Ghana, we limited R t estimation to 27 districts with reported cases before March 30th, 2020 (the beginning of the partial lockdown in Ashanti and Greater Accra), and at least 100 reported cases during the entire study period (12th March to 1st September 2020). We chose this threshold to select districts with surveillance resources capable of detecting cases in the . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted November 17, 2021. ; https://doi.org/10.1101/2021.11.01.21265660 doi: medRxiv preprint early epidemic, due to uncertainty about whether case reporting followed the path of the epidemic or the availability of testing in districts with irregular case reporting. R t estimates were calculated at the district level using the EpiNow2 R package (1.3.2) using MCMC, as implemented in Stan [27,28], based on weekly reported cases [29]. Expected daily cases were estimated using the renewal equation to weight prior expected cases multiplied by the estimated R t . Variation in R t over time was modelled using a mean intercept and an approximate Gaussian process with a 3/2 matern kernel on the log scale [30,31]. Unlike in Abbott et. al. [30] we modelled R t explicitly with the gaussian process and not as a first order difference. This has little impact on retrospective R t estimates and substantially reduces the computational overhead. We used a generation time modelled as a gamma distribution with mean: 3.6 (standard deviation of mean: 0.7), standard deviation 3.1 (standard deviation of standard deviation: 0.76) and maximum: 15 [5,32]. We assumed a negative binomial observation model for reported cases with a day of the week effect modelled as a simplex allowing us to model weekly reported cases without manual specification. R t estimates did not include an estimate of reporting delays as lags were estimated in subsequent analyses.
Therefore, estimates of R t on a given date vary as a result of the reported cases on that date. Inference was performed across 4 chains for 2000 samples with a burn-in period of 250 samples. Convergence was diagnosed using the R hat diagnostic [28].

Defining a Mobility Indicator from Vodafone Data
We used Vodafone Ghana Call Detail Records (CDRs) aggregated by the Flowminder Foundation prior to data sharing [33]. CDRs record mobile phone connections to the cellular towers routing a call or SMS. CDRs are used by mobile network providers for billing purposes. CDRs provide the location of the mobile phone and SIM card based on the location of the cellular tower routing the signal, most often the nearest one. The precision of estimated mobile phone locations depends on the density of cellular towers in an area and signal coverage. It can reach up to 3 km in average coverage conditions and up to 8 km in good coverage conditions. Individual mobility can be estimated by recording a series of mobile phone connections to cellular towers over time. We used CDR data aggregated into an origin-destination matrix, based on the locations recorded within 24 hours for individual mobile phones. The data were censored to remove daily counts of 10 or fewer subscribers recorded for an origin-destination pair, in order to reduce the risk of statistical disclosure of personally identifiable information.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted November 17, 2021. ; https://doi.org/10.1101/2021.11.01.21265660 doi: medRxiv preprint We generated a normalised index of movement outside of individual districts relative to baseline values from the origin-destination matrix. We used two metrics to calculate this normalised measure of mobility: (1) trips between districts: the daily number of subscribers travelling between pairs of districts, and (2) total subscribers per district: the total number of unique mobile phone subscribers recorded in a district on each day. Because of inconsistencies between the spatial references used for the mobility data and the case data, we aggregated mobility indicators for Accra, Tema, and Kumasi Metropolitan Areas by removing trips between aggregated districts and calculating the sum of subscribers for these districts. The mobility data also contained 10 missing dates (4.1%) and we performed linear interpolation for each district for both metrics (trips and subscriber counts) on these dates (Supplemental Figure 3). This resulted in a normalised mobility indicator (Supplemental Figure 4).

Google Mobility Indicator
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted November 17, 2021. ; https://doi.org/10.1101/2021.11.01.21265660 doi: medRxiv preprint We used mobility data from Google as a second measure of human movement [16]. This data records the GPS location of individuals actively using Google services who have chosen to share their location data with Google. The data is provided as a measure of changes in activity relative to a baseline in different settings (Residential, Grocery & Pharmacy, Retail & Recreation, Transit Stations, Workplaces). The dataset documentation recommends consideration of the specifics of mobility in different settings. We chose to use only the mobility indicator from the "Residential" setting because we considered this setting to be the most clearly defined setting in the context of Ghana, and because of the relatively lower variance of this indicator (Supplemental Figure 5). We calculated the inverse of the percent change in residential mobility. This percentage is relative to a baseline period between 3rd January and 6th February, 2020 which is defined by Google prior to data sharing.
Google mobility data is not referenced to known administrative areas but rather to custom boundary polygons, which do not closely align with administrative districts in Ghana. To combine this mobility data with the other data sources used in this study, we manually digitized (traced) these features to create a spatial representation of the coverage area of each metric. To align Google mobility data with our spatial reference, we assigned Google mobility data to those districts with greater than 50% overlap with the administrative areas defined by Google for Accra and Kumasi. This restricted the coverage of Google mobility data to central districts in the Accra and Kumasi Metropolitan Areas (Supplemental Figure   6).

Statistical Analysis
We assessed the association between NPIs, mobility and median R t while adjusting for public holidays using a two-level multilevel generalised linear mixed model (using a Gaussian observation model), with random intercepts for individual districts to account for local variation. We used data from 12th March 2020 to 1st September 2020, during the period of interventions in Ghana, and before the detection of the Alpha or Beta variants [34].
Level 1 of the multi-level model included the relationship between mobility indicators and the inverse NPI stringency with district-level random effects, x to account for correlation between mobility and NPI stringency:

Level 1: Mobility a ~ x a + NPI stringency
Where a indicates individual districts.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

+ Holidays , i -j
Where x is a random effect, a indicates individual districts, j is a lag between 1 and 30 days, and i is the original date of data collection. We trained models for each time period and lag values to determine the optimal lag between R t and mobility. We assessed the different models by comparing the marginal R 2 , which represents the contribution of fixed effects only, for different lag values and time periods. We also calculated the Median Absolute Error of each model: For i in 1...N values (either predicted or observed) where N is the total number of observations.
Coefficients of the Level 2 model estimate R t given inverse NPI stringency, holiday events, and the residuals of the Level 1 model (which can be interpreted as "mobility not explained by NPI stringency"). The use of inverse NPI stringency means that positive coefficients can be interpreted similarly for NPI stringency and residuals of the Level 1 model. For example, a positive coefficient indicates that R t will increase as NPI stringency decreases. Independent variables were centred and scaled for all models to allow for comparison between model coefficients.
Holiday periods included Easter, Eid al-Fitr, Eid al-Adha, and National holidays. We used the custom NPI stringency index and performed a sensitivity analysis using the OxCGRT index (Supplemental Section 2). We used the Vodafone mobility index in the main model since it is available in more districts, and performed a sensitivity analysis using the Google mobility index (Supplemental Section 3).
To determine if the association between mobility, NPIs, and R t at different points of the epidemic is time-varying, we repeated model training for varying-length periods from 12th March to t for t in 19th March to 1st September (7-137 days). To understand the influence of varying time periods on model training, we also conducted a sensitivity analysis training the model in rolling fixed-length periods of 30, 60, and 90 days (Supplemental Section 4).
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

COVID-19 epidemic in Ghana
Broadly, the first wave of the national COVID-19 epidemic in Ghana was characterised by an early increase in cases in March and April 2020, followed by a decline in cases over the summer and a resurgence in June and July 2020 (Figure 1). Using aggregated surveillance data for 27 districts included in the estimation of R t , we observed variations in the progression of local epidemics in individual districts (Supplemental Figure 1). Patterns in each district varied, with case reports ranging from 1 to 250 cases per day, with districts reporting cases on average in 89 of 173 days. Ghana introduced a series of NPIs in response to the growing number of COVID-19 cases in March 2020 ( Figure 1). On 1st April 2020, a partial lockdown was introduced in the Ashanti and Greater Accra regions requiring individuals to remain at home except for essential errands (shopping, healthcare, use of public toilets). The restrictions also prohibited inter-city travel except for essential services.
Lockdown restrictions remained in place until 28th April 2020.

Figure 1: Confirmed COVID-19 cases and Non-pharmaceutical interventions. The total number of confirmed cases of COVID-19 in districts included in this study. The timeline of different nonpharmaceutical interventions are indicated with dashed lines.
Both the custom and OxCGRT stringency indices reflect similar patterns in Ghana: a peak stringency coinciding with the introduction of the partial lockdown in Ashanti and Greater Accra regions, and a following reduction beginning in July 2020 (Supplemental Figure 2).
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Figure 2. Mobility indicators in Accra and Kumasi Metropolitan Areas. a) A comparison of the Vodafone and Google mobility indicators in Accra and Kumasi Metropolitan Areas. b) The correlation between Vodafone and (inverse) Google mobility indicators. Blue dashed line indicates the best fit line. This shows a strong correlation between both mobility indicators across the study period.
Note that these data are collected from two different sources (Google: GPS, Vodafone: CDRs) and describe different aspects of mobility (Google: activity in "residential" areas, Vodafone: travel between administrative districts).

District-level estimates of R t
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted November 17, 2021. ; https://doi.org/10.1101/2021.11.01.21265660 doi: medRxiv preprint After the first detection of cases on March 12th, reported cases and the number of districts reporting cases grew until the announcement of the partial lockdown (Figure 3a-c). We found R t above 1 (indicating a growing epidemic) after the first reported COVID-19 cases in individual districts and a subsequent decline coinciding with the period of maximal interventions (Figure 3b-c). This was followed by an increase in R t during the summer of 2020. While district-specific epidemics followed a broad trend, transmission in individual districts was characterised by varying patterns of epidemic progression (Supplemental Figure 8). We compared R t estimates one week before and after the announcement of a partial lockdown in Ashanti and Greater Accra regions (the date of maximum intervention), finding that between 25th March and 8th April, R t decreased in 16 of the 27 districts (7 missing).

Association between NPI stringency, mobility, and R t
We found an optimal lag of mobility, NPI stringency, and holiday dates of 22 days associated with R t , measured by the maximum marginal R 2 of the multilevel model training across all periods. Because R t estimates did not include estimated delays from infection to reporting, this lag reflects the delay between mobility and NPI, infection, case detection and reporting.
Sensitivity analysis using OxCGRT stringency index and Vodafone mobility indicator, as well as Google mobility indicator and Custom stringency index identified optimal lags of 21 and 19 days, respectively. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted November 17, 2021. ; https://doi.org/10.1101/2021.11.01.21265660 doi: medRxiv preprint We identified correlation between mobility and inverse NPI stringency using the Level 1 model (Supplemental Table 2). We found that the Level 2 model explained a greater amount of variance in R t and that the strength of association between R t and NPI stringency was highest during the early epidemic (Figure 4a, Supplemental Figure 9). Across all training periods (173) and lag values (0-30), the maximum marginal R 2 was 0.51 (conditional R 2 : 0.64) using data between March 12th and May 12th (Figure 4c). We observed higher absolute error in the beginning of case reporting in specific districts (Supplemental Figure   10). The model identified strong evidence of a positive association between R t and both NPIs and Residual Mobility (Table 1). Positive coefficients indicate an association between NPI stringency and Residual Mobility where R t increases as NPI stringency decreases and mobility increases. Note that a positive coefficient for NPI stringency results from the use of inverse stringency in the model. We did not find evidence of association between R t and holidays (Table 1). For this model, we found 15 district-specific random effects distinguishable from 0 (55.6%) (Supplemental Figure 11).

Discussion
We found that R t was associated with human mobility and NPI stringency in the early stages of the COVID-19 pandemic in Ghana and that this association decreased through time. We also identified a positive association between residual mobility (mobility not explained by NPIs) and R t . In our sensitivity analyses, we found similar optimal lags for both the Custom and OxCGRT stringency indices. These lags were greater than those reported in other studies, (Badr et. al. for example, identified an optimal lag of 14 days between mobility and COVID-19 transmission in the USA). In sensitivity analyses included in Supplemental Section 2, using Vodafone and OxCGRT data, we found that the OxCGRT stringency index explained a greater amount of variance in R t in the end of the study period, which may be due to higher stringency value of the OxCGRT index in this period. We also performed sensitivity analyses using mobility indicators from two providers, Google and Vodafone Mobility. The declining relationship may indicate a disconnection between mobility and R t as the effect of mobility was mediated by other behavioural changes. It is also notable that we did not observe a decrease in mobility preceding the end of the second wave of the pandemic in July. This may indicate that mobility data is most useful in the beginning of the pandemic, when mobility patterns reflect behavioural changes relevant to disease transmission.
These findings are in line with those from high-income countries which find associations between decreases in human movement and a reduction in R t [12,20,35,36]  We used an R t estimation method that supports uncertain generation times via a Bayesian prior with mean 3.6 days (standard deviation of mean 0.7 days) for calculating R t . The use of longer generation times will lead to greater variance in estimates of R t . This could translate into larger effect sizes (positive or negative) in the statistical model.
Analysing only districts with case counts which were available early in the COVID-19 epidemic in Ghana may bias our estimates towards urban populations or populations with greater disease surveillance resources. It is also not possible to determine whether the timing of the first reported cases of COVID-19 in individual districts is related to the progression of local epidemics or to the first availability of PCR testing resources in each district. The mobility indicators used in this study rely on the aggregated locations of subscribers to mobile networks (Vodafone) and users of internet services (Google). The volume and reporting of these locations may be influenced by varying patterns of mobile device usage. The demographics of users of either service may also be different from the demographic of the population of Ghana, particularly for Google data which relies on data collected from internet-connected smartphones [37]. Additionally, we used national, not district-specific indices of NPI stringency. Neither index includes intervention measures . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted November 17, 2021. ; which may have been implemented in local districts but which are not recorded at a national level.
In this study, we identified evidence of positive associations between mobility, NPI stringency, and R t and show how the strength of this relationship changed through time. We found that mobility and NPI stringency was able to explain variance in R t during the early epidemic but this pattern declined as the epidemic progressed. This decline may reflect a disconnection between disease transmission and behavioural changes measured by mobility and NPI indicators. For policymakers and public health decision makers responding to the COVID-19 pandemic, our findings demonstrate that mobility and NPIs were effective for estimating disease transmission during the early epidemic, but that subsequent outbreaks may be more related to factors that are not captured in these data.

Data Availability
Data used in this study included individual Line List data shared with authors by the Ghana Health Service. Use of this data was approved by the LSHTM Research Committee (Ref: 22477) and the Noguchi Memorial Institute of Medical Research (Ref: 048/20-21). Vodafone mobility data used in this study is proprietary data shared by Vodafone Ghana in partnership with the Flowminder Foundation and Ghana Statistical Service. This mobility data is available to researchers by application. Google mobility data used in this study is available in the public domain. Downloads of this data can be found in the references. Code used in this study is available from: https://github.com/hamishgibbs/ghana_rt_npi_mobility.
The following funding sources are acknowledged as providing funding for the named authors. This research was partly funded by the Bill & Melinda Gates Foundation (INV-003174: YL). EDCTP2 (RIA2020EF-2983-CSIGN: HPG, OM, RME, MM). This project has received funding from the European Union's Horizon 2020 research and innovation programme -project EpiPose (101003688: YL). HDR UK (MR/S003975/1: RME). This research was partly funded by the National Institute for Health Research (NIHR) using UK aid from the UK Government to support global health research. The views expressed in this publication are those of the author(s) and not necessarily those of the NIHR or the UK Department of Health and Social Care (16/137/109: YL; NIHR200908: RME). UK DHSC/UK Aid/NIHR (PR-OD-1017-20001: HPG). UK MRC (MC_PC_19065: RME, YL). Wellcome Trust (210758/Z/18/Z: SA).