SCoVMod – a spatially explicit mobility and deprivation adjusted model of first wave COVID-19 transmission dynamics

Background: Mobility restrictions prevent the spread of infections to disease-free areas, and early in the coronavirus disease 2019 (COVID-19) pandemic, most countries imposed severe restrictions on mobility as soon as it was clear that containment of local outbreaks was insufficient to control spread. These restrictions have adverse impacts on the economy and other aspects of human health, and it is important to quantify their impact for evaluating their future value. Methods: Here we develop Scotland Coronavirus transmission Model (SCoVMod), a model for COVID-19 in Scotland, which presents unusual challenges because of its diverse geography and population conditions. Our fitted model captures spatio-temporal patterns of mortality in the first phase of the epidemic to a fine geographical scale. Results: We find that lockdown restrictions reduced transmission rates down to an estimated 12\% of its pre-lockdown rate. We show that, while the timing of COVID-19 restrictions influences the role of the transmission rate on the number of COVID-related deaths, early reduction in long distance movements does not. However, poor health associated with deprivation has a considerable association with mortality; the Council Area (CA) with the greatest health-related deprivation was found to have a mortality rate 2.45 times greater than the CA with the lowest health-related deprivation considering all deaths occurring outside of carehomes. Conclusions: We find that in even an early epidemic with poor case ascertainment, a useful spatially explicit model can be fit with meaningful parameters based on the spatio-temporal distribution of death counts. Our simple approach is useful to strategically examine trade-offs between travel related restrictions and physical distancing, and the effect of deprivation-related factors on outcomes.


Introduction
The coronavirus disease 2019  pandemic resulted in the introduction of severe non-pharmaceutical interventions (NPIs) to control its spread in most countries around the world. COVID-19 was introduced into Scotland no later than February 2020. Following a series of lesser restrictions and recommendations, on 23rd March a lockdown with widespread NPIs was put in place across the UK, which resulted in cessation of all non-essential activities, many businesses closed and most individuals except key workers were restricted to only short distance trips from their homes, initially with a maximum frequency of only once per day. These measures reduced the transmission of COVID-19 substantially, so that by the first week of April, the average infected individual in Scotland was estimated to be infecting fewer than one other individual each (i.e. causing the reproduction number R to fall below one).
While we therefore know that the combination of measures was effective, no study as yet has directly estimated the impact of different aspects of these control policies. One factor that makes this difficult for the evaluation of any such national level policy is the impact of variation in space of the disease transmission itself. For COVID-19, one critical issue is the known substantial variation in COVID-19 infection and mortality risk associated with deprivation 1 . The impact of deprivation is likely due to a combination of factors influencing both exposure (e.g. more crowded housing and working conditions) and mortality once exposed (e.g. due to already poorer health) and this should be accounted for in any assessment of the impact of restrictions.
Here, we develop an explicitly spatial agent-based simulation model that accounts for recorded movements-to-work (i.e. "commuter" patterns), modulated by recorded time-varying mobility statistics, and geographically explicit population age structures. Whilst this does not capture all human movement, it is expected to capture a large proportion of the long-range mobility especially affected by lockdown. We also modify the modelled epidemic using deprivation metrics. In these early epidemic stages case data were unevenly recorded, therefore we use this model to estimate transmission characteristics by fitting it to the observed number of COVID-19 related deaths using Approximate Bayesian inference. Our aim is to estimate the impact of travel restrictions and transmission reduction on the spread of COVID-19, and to assess the impact of measures of deprivation in order to estimate its impact on COVID-19 related mortality.
The impact of NPIs is of considerable interest and has been the subject of several analyses, however previous results have focused on statistical comparisons of interventions across multiple countries [2][3][4] . Such approaches are extremely useful for identifying the impact of events, but do not directly allow for disentangling the impact of simultaneous events. Long distance movement and local restrictions are typically applied contemporaneously, and our approach that includes processes at the two scales (local and regional) allows us to estimate the impact of each.

Data sources
We use a combination of publicly available Scottish census data from National Records for Scotland (NRS) 1 and data on COVID-19 held by Public Health Scotland (PHS) and made available via the PHS Electronic Data Research and Innovation Service (eDRIS) 2 . We used data zone (DZ) level resolution where DZs are population census units of approximately 500 to 1,000 residents. The data for assignment of individuals to work locations is drawn from the NRS Census Flows data 3 , Table WU01UK, which provides origin/destination workplace data for the population from the 2011 census. We adjust these with respect to the 2018 population estimates.
Age demographics and movement to work patterns are available at the level of Census Output Areas (OA), each of which contains approximately 20 households or 50 people 4 .
Census data on the Scottish Index of Multiple Deprivation (SIMD) 5 considers multiple relative deprivation measures and combines them into a single value. Deprivation data are publicly available at the DZ level.
We also used publicly available data from Google to estimate mobility levels over time, with respect to commuting patterns 6 .

Model overview
The model we present here breaks down into a number of distinct parts, each of which we describe over the following sections. The code is available from GitHub and is archived with Zenodo. First, to account for the variation in mortality risk due to deprivation factors we develop a statistical model that provides a mortality parameter adjustment to the main model. The core of the simulation model then breaks down into the following parts: • Local transmission-a homogeneous mixing compartmental model for each unit area of the country; • National transmission-a network-based simulation of the movement of individuals between unit areas; • Lockdown simulation-the reduction of both local and national transmission to simulate non-pharmaceutical interventions; • Parameter inference-a Bayesian estimation of the parameters for local transmission.
The simulation framework (Scotland Coronavirus transmission Model, or SCoVMod) considers key aspects of COVID-19  epidemiology including phases for latent infection, infectious  and mildly infected (showing few or no clinical signs) and severely infected (with substantial clinical signs) individuals, hospitalised, recovered and died, similar to other investigations 5,6 . These epidemiological processes are captured as individual disease states (Figure 1). Individuals are also stratified into three age groups: young (0-15), adult  and elderly (65+). Within-OA transmission is assumed to be homogeneously mixed, while between-OA transmission is determined by the empirical age-specific patterns of home and work contact (creating day/night patterns of contact). We do not consider overnight shifts in location or introductions from outside Scotland beyond the impact on the initial seeding. Death rates for all age classes are assumed to be the same, but with differing recovery rates, resulting in age-dependent differences in the infection-related mortality, as is consistent with the data at the time (see Table 1). Deprivation is also known to influence COVID-19 mortality 7 ; we therefore adjust mortality in the model with the average health index in the local council area (see below).
We make a number of simplifying assumptions regarding transmission pathways. First, infections in care homes are not modelled at this stage of the epidemic, as they are assumed to result in few additional infections outside of these locations. Second, hospital-acquired infections are considered 'dead ends', following evidence that they do not result in substantial community outbreaks 8 . Third, population mobility patterns are determined by the patterns of movements to work, recorded in Scottish Census data. We assume that only adults contribute to commuter movement, in the daytime; the remaining proportion of adults and all young and elderly individuals are assumed to move primarily within their local OAs, which also account for non-work activities. Finally, commuting is restricted to healthy and exposed or mildly symptomatic individuals; severely infected and hospitalised individuals do not commute. The day/night patterns also result in two transmission rates.

Health index and mortality
Using the SIMD we constructed two initial statistical models, with a view to guiding the parameters of the core simulation model. These consisted of:    (mildly infected), I S (severely infected), H (hospitalised). Model equations for individuals residing in one OA i and for age class a are therefore: The number of deaths is determined by: The rates of infection are detailed in Figure 2 and other state transition rates are given by: where μ CA is the COVID-19 related mortality rate for a given CA, μ av is the average across CAs, k CA is the CA mean health index value (from the SIMD), and μ mod is a fitted parameter given a prior range of [0, 0.08] in order to preclude negative values for low values of k CA .
The values for all parameters are either established from the literature ( Table 1) or fit (see below).

National transmission (between-OA movements)
Between OAs individuals move daily across a network of locations defined by Scottish Census data adjusted by Google mobility data.
From the current population estimates we draw the number of individuals whose primary residence is mapped onto an OA, with their age group. An individual's workplace is assigned by distributing a proportion of the population of each location to each work location, weighted by the proportion of individuals from each home location in the census flows data who work in another location. For each origin o and destination d we assign a weight w od from the census flow data: where n od is the total number of people who move from o to d to work, and t o is the total number who move from origin o to any location for work. We take the individuals of each home location if they are eligible to work (total n o ); in this case we assume all individuals of adult age 16-65. Each destination is assigned to n o × w od of these individuals. The individuals who remain have no assigned workplace-either they do not work, or they work within their home location.
For each day of the simulation we consider two time steps: a day step where individuals can move to their place of work, and a night step where those individuals move back to their home location. In each day step, we take each destination location d. In order to improve the computational efficiency of the simulation, movements of commuters between OAs were batched into groups of five, with movements between OAs of fewer than five individuals per day retained at a proportionate rate by drawing from a binomial distribution: s mt ~ B(s m , 1 5 ).
If the sampled number of workers s mt is less than or equal to the number of workers who may normally move to destination d, then those who move are sampled randomly from those who may normally move. However, if s mt is greater than the number of workers who may normally move to d, then the additional workers are drawn randomly from workers who have no assigned destination location. While this reduces the overall network link density, the effect on transmission dynamics is expected to be small. We note that this means that interpretation of the combined β D and β N must be made with caution.
For each night of the simulation, the workers who moved in the day step are moved back to their origin location.

Modelling lockdown
To model reductions in activity that are the aim of lockdown, we consider two factors. First, we thin movements in the simulation (mobility reduction) in proportion to the observed changes in mobility 7 (as above), also checking consistency with social contact surveys 25 . This was corroborated for remote areas, where mobility data are few, with an independent dataset (see Figure 3). Second, physical distancing is incorporated via a reduction in contacts applied to both daytime and nighttime transmission rates (transmission reduction). For the post-lockdown period, we contingently fit a reduction in transmission only, assuming that posterior distributions for all other parameters estimated based on the pre-lockdown fit remain relevant, and considering a range of values consistent with independent estimates of the reproduction number R t 8 .

Model inference and computation
Simulated epidemics are compared to the spatio-temporal pattern of COVID-19 spread in Scotland. Non-observable parameters were estimated using the number of deaths where COVID-19 is recorded, considering all weeks beginning 9th March and ending on the 12th April 2020. Recorded deaths are from the weekly PHS records identified by the DZ of residence 9 . These are the most complete and unbiased indicators of infection available during the early epidemic. However, they differ from other official sources as they record the date of registration of death. We assume that for each reported case, death occurs in the week prior to registration, which is allowed up to 8 days postmortem.
Estimation was performed using a Sequential Monte Carlo implementation of Approximate Bayesian Computation (ABC-SMC) 26,27 . We calibrated the model output to the cumulative weekly number of deaths due to COVID-19 aggregated at the level of CAs, using this spatial variation in deaths across Scotland to provide the necessary signature to properly calibrate the role of human mobility. Preliminary attempts to fit the data using weekly incidence had poor results (results not shown). We hypothesise that this is a result of a combination of very few deaths per CA to constrain the fits in the early stages of the epidemic, and also the complications from the multiple changes in the control efforts especially through March. Thus we chose to emphasise the later stages using as our observation the cumulative number of all COVID-19 related deaths per CA.
Simulated and observed summary statistics are compared via a score equal to a sum of squared errors, recorded weekly: where D sim is the cumulative number of deaths per CA simulated and D obs its observed value, Z sim D is the cumulative number of data zones with deaths per CA simulated and Z obs D its observed value. max(x) is the maximum number of x over all CAs, used to balance the relative contributions of the two elements Comparison to Recreation mobility (right). The comparison is relative to the mean value prior to lockdown on March 23rd, 2020. In order to provide confidence that inclusion of mobility reductions across regions is appropriate, we assume that urban areas such as Glasgow and Edinburgh are likely to be well represented, but that rural areas may be less so. To check this, we compare an independent dataset on independent sailings and passenger numbers for ferry services run by Caledonian MacBrayne, who operate all ferry services in the west of Scotland. A comparison of data from 2019 to 2020 and to Google Mobility data, shows a strong fidelity between the two datasets, as well as a substantial reduction in activity at point of lockdown. The similarity prior to lockdown between 2019 and 2020 also suggests that patterns of increased summer activity are unlikely to have had strong influences on our assumptions regarding commuter movements, at least in this area. Data provided by eDRIS is subject to data sharing agreement due to potential for individual disclosure.
of the metrics to the inference. The total number of infected individuals at the start of the simulation (the seeds) are fitted as part of the inference. The seeds are randomly assigned a disease state from E, I M , and I S . Seed locations are stochastically assigned according to the cumulative proportion of deaths registered per Intermediate Zone, as recorded up to the week of March 23rd, 2020. Intermediate zones are aggregates of approximately five DZs; this scale is chosen to account for clustering of infections in areas near to identified deaths.
Uniform prior distributions constrain all parameter values to plausible ranges based on the available literature relevant to the early, pre-lockdown period. Infection dynamics are simulated via a τ-leap algorithm using halfday timesteps 28 . All parameters are listed in Table 1. The inference framework is run on a distributed application framework (Akka) 10 running on a cloud computing infrastructure (Amazon AWS2) 11 . The model code has been written using industry grade software engineering practices including agile development for project task planning, test driven development, pair programming and code reviews to produce unit tested, robust, and reusable software components. The majority of the code has been reviewed by a second software developer.

Health index and mortality
The results of univariable analyses of the SIMD variables for the population level is presented in Table 2. As access has the lowest Akaike information criterion (AIC) and Bayesian information criterion (BIC) and is orthogonal to the other variables we tested it against the remaining variables in turn to check for a multivariable model. The model with access and health had the lowest AIC (276.3), BIC (282.2) and low variance inflation factor (1.291) and was not overdispersed (Table 3). We analysed death records from Public Health Scotland (PHS) data on non-care home deaths due to COVID-19. We examined the impact of different deprivation factors relevant to this period; whilst deprivation overall is significantly associated with increased COVID-19 mortality, this could be further disaggregated: 2. Risk of excess COVID-19 mortality (COVID-19 deaths as a fraction of all deaths) is most closely associated with the access indicator component of SIMD. This indicates that the areas that have good local connectivity and transport will have higher rates of COVID-19 transmission.

Population level risk of COVID-19 mortality is associated with the SIMD indicator that describes (good) accessibility and orthogonally with the SIMD indicator that
3. We hypothesise that access is a proxy for earlier introduction of infection and the mobility patterns in our data. Thus the observed differences in mortality due to access is already accounted for in our model (via the initial seeding plus transmission dynamics). We therefore fit a modifier to COVID-19 mortality considering the average health index (only) in each CA.

Movement and network patterns
The pattern of movement generated at the OA level shows substantial differences in both the average distance travelled, and the connectedness between OAs across the country. Individuals in remote areas move the farthest to work on average, and individuals in the densely populated "Central Belt" are the most connected ( Figure 4).

Parameter posteriors and model fit
In the fitted model, parameter posteriors are strongly unimodal ( Figure 5) with only weak pair-wise correlations. The most likely (mean) estimate for the mortality modifier is μ mod = 0.03, resulting in COVID-19 mortality rates that are 2.45× higher in the CA with the worst health index, compared to the best. Post-lockdown, the best fit value (lowest score) occurs when transmission rates are reduced to 0.12× the pre-lockdown value, starting from 28th March. The number of observed deaths per CA per week mostly lie well within the range of 95% of the fitted simulations (parameters drawn in sets from the posterior parameter distributions; see Figure 6). The number of DZs with deaths is more likely to exceed these limits ( Figure 7) though with low numbers and therefore greater stochastic variability in the data.

Impact of distance reduction and transmission reduction on COVID-19 spread
Lockdown restrictions impact both the local spread of infection via physical distancing measures that change transmission rates and geographical spread (e.g. reductions in travel-to-work).
In only a few cases over 25% of infection occurs outside the OA of residence (Figure 8).
These are typically urban areas, most likely with considerable inward commuting traffic ( Figure 9). However, disease spread is predominantly local, with the most likely outcome being more than 90% of infection occurring within OA.
We consider counterfactuals where some lockdown restrictions are imposed on 9th March, two weeks prior to the actual date, and just before the first death due to COVID-19 in Scotland was reported. The reduction in transmission rates is attributable to physical distancing, or reduction in mobility, or both. With earlier lockdown, we predict a median 581 deaths (95% of simulations within 377 to 1,010 deaths) by 26th April 2020, compared to 2,722 (95% of simulations within 1,294 to 4,050) in the baseline scenario (observed number is 2,795, assuming that all deaths occur in the week prior to the week the death is registered in). Most of this difference is due to the reduction in transmission rates ( Figure 10). The reduction in the total number of deaths also results in a reduction in geographical spread with many fewer DZs affected by COVID-19 mortality in the early lockdown scenario (Figure 11).

Discussion
Exploiting our estimates of the relative contribution of longer distance movements to work and local transmission that include  work and non-work interactions, we show that at this early stage physical distancing had a dominant impact on the death rate, with little evidence that longer range movement restrictions were important for reducing deaths.
Our model is fitted to the spatio-temporal pattern of deaths due to early stage COVID-19 under the assumption that it maps well onto the patterns of infection at this spatial scale and at this stage of the epidemic. While an inference done on case distributions would be a more direct approach to estimating transmission parameters, testing and hence case ascertainment in the early stages of the pandemic was poor. Our modelling approach shows that, even though the spatio-temporal pattern of deaths is subject to substantial stochastic variability, it can be used to generate meaningful opinions in our posterior parameter estimates and form the basis for prediction of counterfactuals.
Movement restrictions are of course critical for preventing the seeding of new areas and therefore a direct estimate of its role in preventing spread at these early stages of the epidemic can provide vital insights into how important rapidly applied measures are and what elements of those measures are most important. Here we show that a movement ban on its own, while it would restrict the important outcome of geographical spread would on its own have had minimal impact at this stage, on the scale of the epidemic. Instead, local physical distancing measures would have had the greatest impact.
Our analysis does not preclude increased transmission due to deprivation at the local level (e.g. due to greater use of public transport, for example) as these factors are not embedded in the model. Of course over the longer term, both elements of long distance travel restrictions and transmission reduction are needed; reductions in both the areas being affected (allowing for potentially more targeted control efforts) and total deaths (both directly saving lives, and reducing impact on the health care system). More localised or clustered infections also will ultimately result in a higher proportion of recovered individuals in hotspot areas, promoting localised herd immunity   effects. Further, as these hotspots were typically in more urban areas, this also has the important consideration of reducing infection rates in remote areas, where access to hospitals and intensive care units (ICUs) is typically poor. Of course, not all long distance movements occur for work purposes and should the patterns of non-work movements be substantially different to movements-to-work patterns, and be of sufficient volume, then it could have an influence on our results. However, the available mobility data suggest that national lockdown restrictions would have had a greater effect on at least the volume of non-work movements.
The best parameter estimates give a high probability that mildly infected individuals have higher transmission rates than the severely infected, though with shorter duration. This does not imply that they are necessarily producing more virus, as those who have severe symptoms may be self-isolating or be more strictly physically distancing. Our model also estimates the number of individuals with infection in early March, with the most likely case being on the order of 1500 infected individuals in Scotland at this time. While our approach is crude and does not take account of continuing importation of individuals over time, it is substantially greater than the estimate of 113 introductions of COVID-19 into Scotland based on viral sequence data 29 , but is plausible if one takes into account the additional infections these introductions would have caused by the start of our simulations.
By explicitly modelling transmission dynamics, this allows us with relatively few data, to infer in our model that health-related deprivation results in a 2.45× difference in the death rate due to COVID-19 across CAs, rather than because these areas carried a greater burden of infections earlier in the pandemic. Deprivation (including health) in Scotland varies substantially at a much finer scale, with zones with the highest deprivation often neighbouring zones with the lowest. This suggests that, were the data available, a deeper interrogation could provide a much more refined assessment of potential health burdens and risks associated with geographical spread.
The effectiveness of lockdown will vary in space and time, due to differences in human behaviour, and also because of non-linear relationships between the numbers of cases, probability of spread, and logistical burdens on care homes, hospitals and ICUs reduced. These potentially counterbalancing factors would of course have to be considered in more detail for a full understanding of lockdown effectiveness. Despite these caveats, our simple approach is useful to strategically examine trade-offs between travel related restrictions, and physical distancing when evaluating future releases from lockdown. Here however, most of the protection from new infections is a result of transmission reduction; long distance restrictions only have a minimal impact because in our model, as overall infection pressure is conserved-i.e. if it does not occur in work locations, it is assumed to occur at home. Thus only extreme travel restrictions are likely to have an impact, at least until build-up of immunity levels is more substantial than observed at this early stage.
Since this analysis, the successful deployment of COVID vaccines in the UK and elsewhere has meant that the response to the current pandemic need not include such extreme measures. However, where vaccine escape mutants become a problem of sufficient concern, and booster vaccines not be readily available, the need for future extreme controls cannot be ruled out, with the recent spread of the Omicron variant one example, although the most extreme measures have thus far proven unnecessary. In situations where such measures are considered, our results indicate that even extreme lockdown measures do not entirely prevent geographical spread between regions; thus any restrictions must include measures to not just reduce mobility but also reduce transmission if, in the future, we are to prevent the spread of COVID-19 to areas that have successfully eradicated local COVID-19 cases.

Data availability
Underlying data Data on underlying population and movements: NRS Census Flows data,

Cameron Zachreson
School of Computing and Information Systems, The University of Melbourne, Melbourne, Vic, Australia This paper presents a spatially explicit model of COVID-19 in Scotland. The population model incorporates medium to long-range mobility patterns as well as local contact rates and the study investigates the effectiveness of lockdown measures with respect to these two model components. The study makes the case that, in Scotland, mobility restrictions were relatively ineffective compared to the influence of local transmission mitigation (i.e., reduction in contact rates within regions). The results indicate that this is because mobility restrictions were implemented too late, and were not stringent enough to prevent the spread of SARS-CoV-2 between areas. The conceptual case is clear and the argument is reasonable, but there are some strong assumptions that could be tested to make the conclusions more robust. There are also some minor clarity issues with the model description, and the presentation of data could be more thorough/professional. Overall I found the paper well-written, the arguments clearly articulated, and the conclusions relatively convincing. Furthermore, the code provided appears to have been developed according to a high professional standard (though it could benefit from additional explanatory documentation and inline comments).
Main issues: strong assumption about initial conditions; posterior correlations not presented.
Suggest analysis of sensitivity to strong assumptions: (i) distribution of initial cases in proportion to observed cumulative deaths (this one is the most important): Models like this are very sensitive to initial conditions (location and number of seed cases).
Since fitting is performed based on the spatial distribution of deaths, the assumption that the initial conditions are themselves strongly dependent on the observed outcomes could be an artificially strong constraint. The true initial conditions were not observed with high levels of accuracy, so some well-reasoned relaxation of the strong constraint used in the design of initial conditions should be explored. For example, a null model could be compared in which the seed cases are distributed with respect to population only (which ○ could be age-stratified etc.).
(ii) Distribution of mobility reduction uniformly in proportion to global mobility reduction: The model defines mobility patterns based on a statistical snapshot. The mobility dynamics under lockdown are treated as uniform adjustments to the snapshot of mobility patterns in the base model. This is a reasonable treatment given the limitations of the data, but represents a strong assumption against which the robustness of the results should be tested. For example, a perturbation approach could be used to systematically randomise the pairwise reductions in mobility, or an alternate assumption could be tested, in which relative mobility reductions depend on distance between origin and destination or some other set of region pair properties. The results in Figure 10 could be presented more clearly (a 1-d histogram for each observable may be preferable given the very high correlation).

○
Equations should be numbered.

○
Model could be more clearly described as an agent-based model, and specified according to agent state variables (Note: this is especially confusing in the discussion of how individuals are assigned work locations and how the mobility flows are computed -it needs to be clear that each working adult is assigned only one location, I had to read the section four times to understand that origins were not stochastically assigned each day).

Is the work clearly and accurately presented and does it cite the current literature? Yes
Is the study design appropriate and is the work technically sound? Yes

Are sufficient details of methods and analysis provided to allow replication by others? Yes
If applicable, is the statistical analysis and its interpretation appropriate?
those used for the posteriors in Figure 5 (e.g. the transmission rate multiplier for mild infectors is given by both "y" and "phi_beta" and differences in capitalisation for day and night).
For the equation for the number of deaths, should it also include the number of susceptible people?

○
The definition of the degree in Figure 4 is not completely clear.

○
In Figure 10, the y axis label should have "data zones" as two words. The caption should be split into two sentences, and the number of simulations done also stated here. The methodology for drawing the rings is not obviously explained. In addition, most of the points look red (all the red points were plotted last?), and so the figure is difficult to read. Could some jittering be applied or the order in which the points are plotted randomised to provide a more representative picture? ○ In Figure 11, the first of the two figure legends should be removed. ○ Also, Figure 3 and Figure 4 are currently presented as low resolution images.

If applicable, is the statistical analysis and its interpretation appropriate? Yes
Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Yes