A BAYESIAN HIERARCHICAL SMALL AREA POPULATION MODEL ACCOUNTING FOR DATA SOURCE SPECIFIC METHODOLOGIES FROM AMERICAN COMMUNITY SURVEY, POPULATION ESTIMATES PROGRAM, AND DECENNIAL CENSUS DATA.

Small area population counts are necessary for many epidemiological studies, yet their quality and accuracy are often not assessed. In the United States, small area population counts are published by the United States Census Bureau (USCB) in the form of the decennial census counts, intercensal population projections (PEP), and American Community Survey (ACS) estimates. Although there are significant relationships between these three data sources, there are important contrasts in data collection, data availability, and processing methodologies such that each set of reported population counts may be subject to different sources and magnitudes of error. Additionally, these data sources do not report identical small area population counts due to post-survey adjustments specific to each data source. Consequently, in public health studies, small area disease/mortality rates may differ depending on which data source is used for denominator data. To accurately estimate annual small area population counts and their associated uncertainties, we present a Bayesian population (BPop) model, which fuses information from all three USCB sources, accounting for data source specific methodologies and associated errors. We produce comprehensive small area race-stratified estimates of the true population, and associated uncertainties, given the observed trends in all three USCB population estimates. The main features of our framework are: (1) a single model integrating multiple data sources, (2) accounting for data source specific data generating mechanisms and specifically accounting for data source specific errors, and (3) prediction of population counts for years without USCB reported data. We focus our study on the Black and White only populations for 159 counties of Georgia and produce estimates for years 2006-2023. We compare BPop population estimates to decennial census counts, PEP annual counts, and ACS multi-year estimates. Additionally, we illustrate and explain the different types of data source specific errors. Lastly, we compare model performance using simulations and validation exercises. Our Bayesian population model can be extended to other applications at smaller spatial granularity and for demographic subpopulations defined further by race, age, and sex, and/or for other geographical regions.

Joint quantile disease mapping with application to malaria and G6PD deficiency.

Statistical analysis based on quantile methods is more comprehensive, flexible and less sensitive to outliers when compared to mean methods. Joint disease mapping is useful for inferring correlation between different diseases. Most studies investigate this link through multiple correlated mean regressions. We propose a joint quantile regression framework for multiple diseases where different quantile levels can be considered. We are motivated by the theorized link between the presence of malaria and the gene deficiency G6PD, where medical scientists have anecdotally discovered a possible link between high levels of G6PD and lower than expected levels of malaria initially pointing towards the occurrence of G6PD inhibiting the occurrence of malaria. Thus, the need for flexible joint quantile regression in a disease mapping framework arises. Our model can be used for linear and nonlinear effects of covariates by stochastic splines since we define it as a latent Gaussian model. We perform Bayesian inference using the R integrated nested Laplace approximation, suitable even for large datasets. Finally, we illustrate the model's applicability by considering data from 21 countries, although better data are needed to prove a significant relationship. The proposed methodology offers a framework for future studies of interrelated disease phenomena.

Interoperability of statistical models in pandemic preparedness: principles and reality.

We present interoperability as a guiding framework for statistical modelling to assist policy makers asking multiple questions using diverse datasets in the face of an evolving pandemic response. Interoperability provides an important set of principles for future pandemic preparedness, through the joint design and deployment of adaptable systems of statistical models for disease surveillance using probabilistic reasoning. We illustrate this through case studies for inferring and characterising spatial-temporal prevalence and reproduction numbers of SARS-CoV-2 infections in England.

Time varying association between deprivation, ethnicity and SARS-CoV-2 infections in England: A population-based ecological study.

BACKGROUND: Ethnically diverse and socio-economically deprived communities have been differentially affected by the COVID-19 pandemic in the UK. METHOD: Using a multilevel regression model we assessed the time-varying association between SARS-CoV-2 infections and areal level deprivation and ethnicity from 1st of June 2020 to the 19th of September 2021. We separately considered weekly test positivity rate and estimated debiased prevalence at the Lower Tier Local Authority (LTLA) level, adjusting for confounders and spatio-temporal correlation structure. FINDINGS: Comparing the least deprived and predominantly White areas with most deprived and predominantly non-White areas over the whole study period, the weekly positivity rate increases from 2·977% (95% CrI 2.913%-3.029%) to 3·347% (95% CrI 3.300%-3.402%). Similarly, prevalence increases from 0·369% (95% CrI 0.361%-0.375%) to 0·405% (95% CrI 0.399%-0.412%). Deprivation has a stronger effect until October 2020, while the effect of ethnicity becomes more pronounced at the peak of the second wave and then again in May-June 2021. In the second wave of the pandemic, LTLAs with large South Asian populations were the most affected, whereas areas with large Black populations did not show increased values for either outcome during the entire period under analysis. INTERPRETATION: Deprivation and proportion of non-White populations are both associated with an increased COVID-19 burden in terms of disease spread and monitoring, but the strength of association varies over the course of the pandemic and for different ethnic subgroups. The consistency of results across the two outcomes suggests that deprivation and ethnicity have a differential impact on disease exposure or susceptibility rather than testing access and habits. FUNDINGS: EPSRC, MRC, The Alan Turing Institute, NIH, UKHSA, DHSC.

Improving local prevalence estimates of SARS-CoV-2 infections using a causal debiasing framework.

Global and national surveillance of SARS-CoV-2 epidemiology is mostly based on targeted schemes focused on testing individuals with symptoms. These tested groups are often unrepresentative of the wider population and exhibit test positivity rates that are biased upwards compared with the true population prevalence. Such data are routinely used to infer infection prevalence and the effective reproduction number, Rt, which affects public health policy. Here, we describe a causal framework that provides debiased fine-scale spatiotemporal estimates by combining targeted test counts with data from a randomized surveillance study in the United Kingdom called REACT. Our probabilistic model includes a bias parameter that captures the increased probability of an infected individual being tested, relative to a non-infected individual, and transforms observed test counts to debiased estimates of the true underlying local prevalence and Rt. We validated our approach on held-out REACT data over a 7-month period. Furthermore, our local estimates of Rt are indicative of 1-week- and 2-week-ahead changes in SARS-CoV-2-positive case numbers. We also observed increases in estimated local prevalence and Rt that reflect the spread of the Alpha and Delta variants. Our results illustrate how randomized surveys can augment targeted testing to improve statistical accuracy in monitoring the spread of emerging and ongoing infectious disease.

Time varying association between deprivation, ethnicity and SARS-CoV-2 infections in England: a space-time study.

BACKGROUND: Ethnically diverse and socio-economically deprived communities have been differentially affected by the COVID-19 pandemic in the UK. METHOD: Using a multilevel regression model we assess the time-varying association between SARS-CoV-2 infections and areal level deprivation and ethnicity. We separately consider weekly test positivity rate (number of positive tests over the total number of tests) and estimated unbiased prevalence (proportion of individuals in the population who would test positive) at the Lower Tier Local Authority (LTLA) level. The model also adjusts for age, urbanicity, vaccine uptake and spatio-temporal correlation structure. FINDINGS: Comparing the least deprived and predominantly White areas with most deprived and predominantly non-White areas over the whole study period, the weekly positivity rate increases by 13% from 297% to 335%. Similarly, prevalence increases by 10% from 037% to 041%. Deprivation has a stronger effect until October 2020, while the effect of ethnicity becomes slightly more pronounced at the peak of the second wave and then again in May-June 2021. Not all BAME groups were equally affected: in the second wave of the pandemic, LTLAs with large South Asian populations were the most affected, whereas areas with large Black populations did not show increased values for either outcome during the entire period under analysis. INTERPRETATION: At the area level, IMD and BAME% are both associated with an increased COVID-19 burden in terms of prevalence (disease spread) and test positivity (disease monitoring), and the strength of association varies over the course of the pandemic. The consistency of results across the two outcome measures suggests that community level characteristics such as deprivation and ethnicity have a differential impact on disease exposure or susceptibility rather than testing access and habits. FUNDINGS: EPSRC, MRC, The Alan Turing Institute, NIH, UKHSA, DHSC, NIHR.

Long-term exposure to air-pollution and COVID-19 mortality in England: A hierarchical spatial analysis.

Recent studies suggested a link between long-term exposure to air-pollution and COVID-19 mortality. However, due to their ecological design based on large spatial units, they neglect the strong localised air-pollution patterns, and potentially lead to inadequate confounding adjustment. We investigated the effect of long-term exposure to NO2 and PM2.5 on COVID-19 mortality in England using high geographical resolution. In this nationwide cross-sectional study in England, we included 38,573 COVID-19 deaths up to June 30, 2020 at the Lower Layer Super Output Area level (n = 32,844 small areas). We retrieved averaged NO2 and PM2.5 concentration during 2014-2018 from the Pollution Climate Mapping. We used Bayesian hierarchical models to quantify the effect of air-pollution while adjusting for a series of confounding and spatial autocorrelation. We find a 0.5% (95% credible interval: -0.2%, 1.2%) and 1.4% (95% CrI: -2.1%, 5.1%) increase in COVID-19 mortality risk for every 1 µg/m3 increase in NO2 and PM2.5 respectively, after adjusting for confounding and spatial autocorrelation. This corresponds to a posterior probability of a positive effect equal to 0.93 and 0.78 respectively. The spatial relative risk at LSOA level revealed strong patterns, similar for the different pollutants. This potentially captures the spread of the disease during the first wave of the epidemic. Our study provides some evidence of an effect of long-term NO2 exposure on COVID-19 mortality, while the effect of PM2.5 remains more uncertain.

Long-term exposure to air-pollution and COVID-19 mortality in England: a hierarchical spatial analysis.

BACKGROUND: Recent studies suggested a link between long-term exposure to air-pollution and COVID-19 mortality. However, due to their ecological design, based on large spatial units, they neglect the strong localised air-pollution patterns, and potentially lead to inadequate confounding adjustment. We investigated the effect of long-term exposure to NO2 and PM2.5 on COVID-19 deaths up to June 30, 2020 in England using high geographical resolution. METHODS: We included 38 573 COVID-19 deaths up to June 30, 2020 at the Lower Layer Super Output Area level in England (n=32 844 small areas). We retrieved averaged NO2 and PM2.5 concentration during 2014-2018 from the Pollution Climate Mapping. We used Bayesian hierarchical models to quantify the effect of air-pollution while adjusting for a series of confounding and spatial autocorrelation. FINDINGS: We find a 0.5% (95% credible interval: -0.2%-1.2%) and 1.4% (-2.1%-5.1%) increase in COVID-19 mortality rate for every 1µg/m3 increase in NO2 and PM2.5 respectively, after adjusting for confounding and spatial autocorrelation. This corresponds to a posterior probability of a positive effect of 0.93 and 0.78 respectively. The spatial relative risk at LSOA level revealed strong patterns, similar for the different pollutants. This potentially captures the spread of the disease during the first wave of the epidemic. INTERPRETATION: Our study provides some evidence of an effect of long-term NO2 exposure on COVID-19 mortality, while the effect of PM2.5 remains more uncertain. FUNDING: Medical Research Council, Wellcome Trust, Environmental Protection Agency and National Institutes of Health.