Geostatistical inference in the presence of geomasking:A composite-likelihood approach

In almost any geostatistical analysis, one of the underlying, often implicit, modelling assumptions is that the spatial locations, where measurements are taken, are recorded without error. In this study we develop geostatistical inference when this assumption is not valid. This is often the case when, for example, individual address information is randomly altered to provide privacy protection or imprecisions are induced by geocoding processes and measurement devices. Our objective is to develop a method of inference based on the composite likelihood that overcomes the inherent computational limits of the full likelihood method as set out in Fanshawe and Diggle (2011). Through a simulation study, we then compare the performance of our proposed approach with an N-weighted least squares estimation procedure, based on a corrected version of the empirical variogram. Our results indicate that the composite-likelihood approach outperforms the latter, leading to smaller root-mean-square-errors in the parameter estimates. Finally, we illustrate an application of our method to analyse data on malnutrition from a Demographic and Health Survey conducted in Senegal in 2011, where locations were randomly perturbed to protect the privacy of respondents. We conclude that our approach based on the composite likelihood is a feasible and computationally more efficient alternative option to existing likelihood-based methods that deal with positional error in a geostatistical context

Measuring the spatial heterogeneity on the reduction of vaginal fistula burden in Ethiopia between 2005 and 2016

Vaginal fistula is a shattering maternal complication characterized by an anomalous opening between the bladder and/or rectum and vagina resulting in continuous leakage of urine or stool. Although prevalent in Ethiopia, its magnitude and distribution is not well studied. We used statistical mapping models using 2005 and 2016 Ethiopia Demographic Health Surveys data combined with a suite of potential risk factors to estimate the burden of vaginal fistula among women of childbearing age. The estimated number of women of childbearing age with lifetime and untreated vaginal fistula in 2016 were 72,533 (95% CI 38,235–124,103) and 31,961 (95% CI 11,596–70,309) respectively. These figures show reduction from the 2005 estimates: 98,098 (95% CI 49,819–170,737) lifetime and 59,114 (95% CI 26,580–118,158) untreated cases of vaginal fistula. The number of districts having more than 200 untreated cases declined drastically from 54 in 2005 to 6 in 2016. Our results show a significant subnational variation in the burden of vaginal fistula. Overall, between 2005 and 2016 there was substantial reduction in the prevalence of vaginal fistula in Ethiopia. Our results help guide local level tracking, planning, spatial targeting of resources and implementation of interventions against vaginal fistula

Estimating the burden of mycetoma in Sudan for the period 1991-2018 using a model-based geostatistical approach

Mycetoma is widespread in tropical and subtropical regions favouring arid areas with low humidity and a short rainy season. Sudan is one of the highly endemic countries for mycetoma. Estimating the population at risk and the number of cases is critical for delivering targeted and equitable prevention and treatment services. In this study, we have combined a large dataset of mycetoma cases recorded by the Mycetoma Research Centre (MRC) in Sudan over 28 years (1991-2018) with a collection of environmental and water and hygiene-related datasets in a geostatistical framework to produce estimates of the disease burden across the country. We developed geostatistical models to predict the number of cases of actinomycetoma and eumycetoma in areas considered environmentally suitable for the two mycetoma forms. Then used the raster dataset (gridded map) with the population estimates for 2020 to compute the potentially affected population since 1991. The geostatistical models confirmed this heterogeneous and distinct distribution of the estimated cases of eumycetoma and actinomycetoma across Sudan. For eumycetoma, these higher-risk areas were smaller and scattered across Al Jazirah, Khartoum, White Nile and Sennar states, while for actinomycetoma a higher risk for infection is shown across the rural districts of North and West Kurdufan. Nationally, we estimated 63,825 people (95%CI: 13,693 to 197,369) to have been suffering from mycetoma since 1991 in Sudan,51,541 people (95%CI: 9,893 – 166,073) with eumycetoma and 12,284 people (95%CI: 3,800 – 31,296) with actinomycetoma. In conclusion, the risk of mycetoma in Sudan is particularly high in certain restricted areas, but cases are ubiquitous across all states. Both prevention and treatment services are required to address the burden. Such work provides a guide for future control and prevention programs for mycetoma, highly endemic areas are clearly targeted, and resources are directed to areas with high demand

Spatiotemporal Heterogeneity in the Distribution of Chikungunya and Zika Virus Case Incidences and Risk Factors During Their Epidemics in Barranquilla, Colombia, between 2014 and 2016: An Ecological Study

Chikungunya virus (CHIKV) and Zika virus (ZIKV) have recently emerged as global infections with consequential disability adjusted life years (DALYs) and economic burden. This study aimed to explore the spatiotemporal heterogeneity in the occurrence of CHIKV and ZIKV outbreaks throughout Barranquilla, Colombia during 2014 and 2016 and explored the potential for clustering. Incidence data were fitted using multiple Bayesian Poisson models based on a suite of explanatory variables as potential risk factors and multiple options for random effects. A best fit model was used to analyse the case incidence risk for both epidemics to identify any risk factors during their epidemics. Neighbourhoods in the northern region of Barranquilla were hotspots for the outbreaks of both CHIKV and ZIKV. Additional hotspots occurred in the south-western and central regions of the CHIKV and ZIKV outbreaks, respectively. Multivariate conditional autoregressive models strongly identified higher socioeconomic strata (SES) and residing in detached houses as risk factors for ZIKV case incidences. These novel findings challenge the belief that these infections are driven by social vulnerability and merits further study both in Barranquilla and throughout the tropical and subtropical regions of the world.&amp;nbsp;</jats:p

Exponential growth, high prevalence of SARS-CoV-2, and vaccine effectiveness associated with the Delta variant.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections were rising during early summer 2021 in many countries as a result of the Delta variant. We assessed reverse transcription polymerase chain reaction swab positivity in the Real-time Assessment of Community Transmission–1 (REACT-1) study in England. During June and July 2021, we observed sustained exponential growth with an average doubling time of 25 days, driven by complete replacement of the Alpha variant by Delta and by high prevalence at younger, less-vaccinated ages. Prevalence among unvaccinated people [1.21% (95% credible interval 1.03%, 1.41%)] was three times that among double-vaccinated people [0.40% (95% credible interval 0.34%, 0.48%)]. However, after adjusting for age and other variables, vaccine effectiveness for double-vaccinated people was estimated at between ~50% and ~60% during this period in England. Increased social mixing in the presence of Delta had the potential to generate sustained growth in infections, even at high levels of vaccination.The study was funded by the Department of Health and Social Care in England. Sequencing was provided through funding from the COVID-19 Genomics UK (COG-UK) Consortium. P.E. is Director of the Medical Research Council (MRC) Centre for Environment and Health (MR/L01341X/1, MR/S019669/1). P.E. acknowledges support from Health Data Research UK (HDR UK); the National Institute for Health Research (NIHR) Imperial Biomedical Research Centre; NIHR Health Protection Research Units (HPRUs) in Chemical and Radiation Threats and Hazards, and Environmental Exposures and Health; the British Heart Foundation Centre for Research Excellence at Imperial College London (RE/18/4/34215); and the UK Dementia Research Institute at Imperial (MC_PC_17114). S.R., C.A.D. acknowledge support: MRC Centre for Global Infectious Disease Analysis, NIHR HPRU in Modelling and Health Economics, Wellcome Trust (200861/Z/16/Z, 200187/Z/15/Z), and Centres for Disease Control and Prevention (US, U01CK0005-01-02). G.C. is supported by an NIHR Professorship. H.War. acknowledges support from an NIHR Senior Investigator Award and the Wellcome Trust (205456/Z/16/Z). We thank The Huo Family Foundation for their support of our work on COVID-19. Quadram authors gratefully acknowledge the support of the Biotechnology and Biological Sciences Research Council (BBSRC); their research was funded by the BBSRC Institute Strategic Programme Microbes in the Food Chain BB/R012504/1 and its constituent project BBS/E/F/000PR10352. We thank members of the COVID-19 Genomics Consortium UK (COG-UK) for their contributions to generating the genomic data used in this study. COG-UK is supported by funding from the MRC, part of UK Research & Innovation (UKRI), NIHR and Genome Research Limited, operating as the Wellcome Sanger Institute

Using a model-based geostatistical approach to design and analyse the prevalence of schistosomiasis in Kenya

BackgroundInfections caused by both Schistosoma mansoni and Schistosoma haematobium are endemic in Kenya, with over six million children at risk. A national school-based deworming programme was launched in 2012 with the goal of eliminating parasitic worms as a public health problem. This study used a model-based geostatistical (MBG) approach to design and analyse the impact of the programme and inform treatment strategy changes for schistosomiasis (SCH).MethodsA cross-sectional survey of 200 schools across 27 counties of Kenya was utilised. The study design, selection of the schools, and analysis followed the MBG approach, which incorporated historical data on treatment, morbidity, and environmental covariates.ResultsThe overall SCH prevalence was 5.0% (95% CI 4.9%–5.2%) and was estimated, with a high predictive probability of 0.999, to be between 1% and&lt; 10%. The predictive probabilities at county level revealed county heterogeneity, with that of four counties estimated to be between 0% and&lt; 1%, that of 20 counties estimated to be between 1% and&lt; 10%, that of two counties estimated to be between 10% and&lt; 20%, and that of one county estimated to be between 20% and&lt; 50%.ConclusionSCH treatment requirements can now be confidently refined based on the World Health Organization’s guidelines. The four counties with prevalences of between 0% and&lt; 1% may consider suspending treatment only in areas (i.e., sub-counties and wards) where the prevalence is&lt; 1%

An Ensemble Framework for Projecting the Impact of Lymphatic Filariasis Interventions Across Sub-Saharan Africa at a Fine Spatial Scale

BackgroundLymphatic filariasis (LF) is a neglected tropical disease targeted for elimination as a public health problem by 2030. Although mass treatments have led to huge reductions in LF prevalence, some countries or regions may find it difficult to achieve elimination by 2030 owing to various factors, including local differences in transmission. Subnational projections of intervention impact are a useful tool in understanding these dynamics, but correctly characterizing their uncertainty is challenging.MethodsWe developed a computationally feasible framework for providing subnational projections for LF across 44 sub-Saharan African countries using ensemble models, guided by historical control data, to allow assessment of the role of subnational heterogeneities in global goal achievement. Projected scenarios include ongoing annual treatment from 2018 to 2030, enhanced coverage, and biannual treatment.ResultsOur projections suggest that progress is likely to continue well. However, highly endemic locations currently deploying strategies with the lower World Health Organization recommended coverage (65%) and frequency (annual) are expected to have slow decreases in prevalence. Increasing intervention frequency or coverage can accelerate progress by up to 5 or 6 years, respectively.ConclusionsWhile projections based on baseline data have limitations, our methodological advancements provide assessments of potential bottlenecks for the global goals for LF arising from subnational heterogeneities. In particular, areas with high baseline prevalence may face challenges in achieving the 2030 goals, extending the "tail" of interventions. Enhancing intervention frequency and/or coverage will accelerate progress. Our approach facilitates preimplementation assessments of the impact of local interventions and is applicable to other regions and neglected tropical diseases

An Ensemble Framework for Projecting the Impact of Lymphatic Filariasis Interventions Across Sub-Saharan Africa at a Fine Spatial Scale.

Background Lymphatic filariasis (LF) is a neglected tropical disease targeted for elimination as a public health problem by 2030. Although mass treatments have led to huge reductions in LF prevalence, some countries or regions may find it difficult to achieve elimination by 2030 owing to various factors, including local differences in transmission. Subnational projections of intervention impact are a useful tool in understanding these dynamics, but correctly characterizing their uncertainty is challenging. Methods We developed a computationally feasible framework for providing subnational projections for LF across 44 sub-Saharan African countries using ensemble models, guided by historical control data, to allow assessment of the role of subnational heterogeneities in global goal achievement. Projected scenarios include ongoing annual treatment from 2018 to 2030, enhanced coverage, and biannual treatment. Results Our projections suggest that progress is likely to continue well. However, highly endemic locations currently deploying strategies with the lower World Health Organization recommended coverage (65%) and frequency (annual) are expected to have slow decreases in prevalence. Increasing intervention frequency or coverage can accelerate progress by up to 5 or 6 years, respectively. Conclusions While projections based on baseline data have limitations, our methodological advancements provide assessments of potential bottlenecks for the global goals for LF arising from subnational heterogeneities. In particular, areas with high baseline prevalence may face challenges in achieving the 2030 goals, extending the "tail" of interventions. Enhancing intervention frequency and/or coverage will accelerate progress. Our approach facilitates preimplementation assessments of the impact of local interventions and is applicable to other regions and neglected tropical diseases

Predicted Impact of COVID-19 on Neglected Tropical Disease Programs and the Opportunity for Innovation

Due to the COVID-19 pandemic, many key neglected tropical disease (NTD) activities have been postponed. This hindrance comes at a time when the NTDs are progressing towards their ambitious goals for 2030. Mathematical modelling on several NTDs, namely gambiense sleeping sickness, lymphatic filariasis, onchocerciasis, schistosomiasis, soil-transmitted helminthiases (STH), trachoma, and visceral leishmaniasis, shows that the impact of this disruption will vary across the diseases. Programs face a risk of resurgence, which will be fastest in high-transmission areas. Furthermore, of the mass drug administration diseases, schistosomiasis, STH, and trachoma are likely to encounter faster resurgence. The case-finding diseases (gambiense sleeping sickness and visceral leishmaniasis) are likely to have fewer cases being detected but may face an increasing underlying rate of new infections. However, once programs are able to resume, there are ways to mitigate the impact and accelerate progress towards the 2030 goals.</p

Accelerating Progress Towards the 2030 Neglected Tropical Diseases Targets: How Can Quantitative Modeling Support Programmatic Decisions?

Over the past decade, considerable progress has been made in the control, elimination, and eradication of neglected tropical diseases (NTDs). Despite these advances, most NTD programs have recently experienced important setbacks; for example, NTD interventions were some of the most frequently and severely impacted by service disruptions due to the coronavirus disease 2019 (COVID-19) pandemic. Mathematical modeling can help inform selection of interventions to meet the targets set out in the NTD road map 2021-2030, and such studies should prioritize questions that are relevant for decision-makers, especially those designing, implementing, and evaluating national and subnational programs. In September 2022, the World Health Organization hosted a stakeholder meeting to identify such priority modeling questions across a range of NTDs and to consider how modeling could inform local decision making. Here, we summarize the outputs of the meeting, highlight common themes in the questions being asked, and discuss how quantitative modeling can support programmatic decisions that may accelerate progress towards the 2030 targets