Search

Evaluating normalising constants is important across a range of topics in statistical learning, notably Bayesian model selection. However, in many realistic problems this involves the integration of analytically intractable, high-dimensional distributions, and therefore requires the use of stochastic methods such as thermodynamic integration (TI). In this paper we apply a simple but under-appreciated variation of the TI method, here referred to as referenced TI, which computes a single model's normalising constant in an efficient way by using a judiciously chosen reference density. The advantages of the approach and theoretical considerations are set out, along with pedagogical 1 and 2D examples. The approach is shown to be useful in practice when applied to a real problem —to perform model selection for a semi-mechanistic hierarchical Bayesian model of COVID-19 transmission in South Korea involving the integration of a 200D density. [ABSTRACT FROM AUTHOR] more...

The UK and Sweden have among the worst per-capita COVID-19 mortality in Europe. Sweden stands out for its greater reliance on voluntary, rather than mandatory, control measures. We explore how the timing and effectiveness of control measures in the UK, Sweden and Denmark shaped COVID-19 mortality in each country, using a counterfactual assessment: what would the impact have been, had each country adopted the others' policies? Using a Bayesian semi-mechanistic model without prior assumptions on the mechanism or effectiveness of interventions, we estimate the time-varying reproduction number for the UK, Sweden and Denmark from daily mortality data. We use two approaches to evaluate counterfactuals which transpose the transmission profile from one country onto another, in each country's first wave from 13th March (when stringent interventions began) until 1st July 2020. UK mortality would have approximately doubled had Swedish policy been adopted, while Swedish mortality would have more than halved had Sweden adopted UK or Danish strategies. Danish policies were most effective, although differences between the UK and Denmark were significant for one counterfactual approach only. Our analysis shows that small changes in the timing or effectiveness of interventions have disproportionately large effects on total mortality within a rapidly growing epidemic. [ABSTRACT FROM AUTHOR] more...

In the present study, an attempt is made to understand the behavior and pattern of tunnel damage subjected to the static and dynamic loading conditions, through the simulation of natural as well as artificial stress states such as loads due to overburden, building, earthquake and blast. The analysis is carried out by using FEM based numerical tool i.e. ABAQUS for different geo-materials. The Mohr-Coulomb plasticity model is used for modeling the rock mass and Johnson–Cook plasticity model have been used for the simulation of stress–strain response of drop hammer. The results of the analysis are studied for stresses, deformation and safety zones in surrounding rock. The displacement of the crown along the axis of tunnel is quantified. It is observed that the strength properties play a crucial role in quantifying the yielded zone in rock mass under both static and impact loading. The effect of weathering is investigated and the deformations for different weathering stages is measured. The crack dimensions reduce as the brittleness of the material increases. Different zones of influence are investigated for all the materials. It is noticed that the zone is significantly crushed in the case of dynamic loading, due to the high velocity of impact and the length of crack is large in the case of static loading as compared to the dynamic loading. [ABSTRACT FROM AUTHOR] more...

This study is focused on the effect of vibration induced by moving trains in tunnels on the surrounding ground and structures. A three-dimensional finite element model is established for a one-track railway tunnel and an adjacent twelve-storey building frame by using commercial software Midas GTS-NX (2019) and Midas Gen. This study considered the moving load effect of a complete train, which varies with space as well as with time. The effect of factors such as train speed, overburden pressure on the tunnel and variation in soil properties are studied in the time domain. As a result, the variations in horizontal and vertical acceleration for two different sites, i.e., the free ground surface (without structure) and the area containing the structure, are compared. Also, the displacement pattern of the raft foundation is plotted for different train velocities. At lower speeds, the heaving phenomenon is negligible, but as the speed increases, both the heaving and differential settlement increase in the foundation. This study demonstrates that the effect of moving train vibrations should be considered in the design of new nearby structures and proper ground improvement should be considered for existing structures. [ABSTRACT FROM AUTHOR] more...

Background: The COVID-19 pandemic has caused over 7.02 million deaths as of January 2024 and profoundly affected most countries' Gross Domestic Product (GDP). Here, we study the interaction of SARS-CoV-2 transmission, mortality, and economic output between January 2020 and December 2022 across 25 European countries. Methods: We use a Bayesian mixed effects model with auto-regressive terms to estimate the temporal relationships between disease transmission, excess deaths, changes in economic output, transit mobility and non-pharmaceutical interventions (NPIs) across countries. Results: Disease transmission intensity (logRt) decreases GDP and increases excess deaths, where the latter association is longer-lasting. Changes in GDP as well as prior week transmission intensity are both negatively associated with each other (-0.241, 95% CrI: -0.295 - -0.189). We find evidence of risk-averse behaviour, as changes in transit and prior week transmission intensity are negatively associated (-0.055, 95% CrI: -0.074 to -0.036). Our results highlight a complex cost-benefit trade-off from individual NPIs. For example, banning international travel is associated with both increases in GDP (0.014, 0.002—0.025) and decreases in excess deaths (-0.014, 95% CrI: -0.028 - -0.001). Country-specific random effects, such as the poverty rate, are positively associated with excess deaths while the UN government effectiveness index is negatively associated with excess deaths. Interpretation: The interplay between transmission intensity, excess deaths, population mobility and economic output is highly complex, and none of these factors can be considered in isolation. Our results reinforce the intuitive idea that significant economic activity arises from diverse person-to-person interactions. Our analysis quantifies and highlights that the impact of disease on a given country is complex and multifaceted. Long-term economic impairments are not fully captured by our model, as well as long-term disease effects (Long COVID). [ABSTRACT FROM AUTHOR] more...

Background: The ability to accurately predict survival in older adults is crucial as it guides clinical decision making. The added value of using health care usage for predicting mortality remains unexplored. The aim of this study was to investigate if temporal patterns of healthcare expenditures, can improve the predictive performance for mortality, in spousal bereaved older adults, next to other widely used sociodemographic variables. Methods: This is a population-based cohort study of 48,944 Danish citizens 65 years of age and older suffering bereavement within 2013–2016. Individuals were followed from date of spousal loss until death from all causes or 31st of December 2016, whichever came first. Healthcare expenditures were available on weekly basis for each person during the follow-up and used as predictors for mortality risk in Extreme Gradient Boosting models. The extent to which medical spending trajectories improved mortality predictions compared to models with sociodemographics, was assessed with respect to discrimination (AUC), overall prediction error (Brier score), calibration, and clinical benefit (decision curve analysis). Results: The AUC of age and sex for mortality the year after spousal loss was 70.8% [95% CI 68.8, 72.8]. The addition of sociodemographic variables led to an increase of AUC ranging from 0.9% to 3.1% but did not significantly reduce the overall prediction error. The AUC of the model combining the variables above plus medical spending usage was 80.8% [79.3, 82.4] also exhibiting smaller Brier score and better calibration. Overall, patterns of healthcare expenditures improved mortality predictions the most, also exhibiting the highest clinical benefit among the rest of the models. Conclusion: Temporal patterns of medical spending have the potential to significantly improve our assessment on who is at high risk of dying after suffering spousal loss. The proposed methodology can assist in a more efficient risk profiling and prognosis of bereaved individuals. [ABSTRACT FROM AUTHOR] more...

Renewal equations are a popular approach used in modelling the number of new infections, i.e., incidence, in an outbreak. We develop a stochastic model of an outbreak based on a time-varying variant of the Crump–Mode–Jagers branching process. This model accommodates a time-varying reproduction number and a time-varying distribution for the generation interval. We then derive renewal-like integral equations for incidence, cumulative incidence and prevalence under this model. We show that the equations for incidence and prevalence are consistent with the so-called back-calculation relationship. We analyse two particular cases of these integral equations, one that arises from a Bellman–Harris process and one that arises from an inhomogeneous Poisson process model of transmission. We also show that the incidence integral equations that arise from both of these specific models agree with the renewal equation used ubiquitously in infectious disease modelling. We present a numerical discretisation scheme to solve these equations, and use this scheme to estimate rates of transmission from serological prevalence of SARS-CoV-2 in the UK and historical incidence data on Influenza, Measles, SARS and Smallpox. [ABSTRACT FROM AUTHOR] more...

Objective To develop a prognostic model of 1-year mortality for individuals aged 65+ presenting at the emergency department (ED) with a fall based on health care spending patterns to guide clinical decision-making. Design Population-based cohort study (n  = 35,997) included with a fall in 2013 and followed 1 year. Methods Health care spending indicators (dynamical indicators of resilience, DIORs) 2 years before admission were evaluated as potential predictors, along with age, sex and other clinical and sociodemographic covariates. Multivariable logistic regression models were developed and internally validated (10-fold cross-validation). Performance was assessed via discrimination (area under the receiver operating characteristic curve, AUC), Brier scores, calibration and decision curve analysis. Results The AUC of age and sex for mortality was 72.5% [95% confidence interval 71.8 to 73.2]. The best model included age, sex, number of medications and health care spending DIORs. It exhibited high discrimination (AUC: 81.1 [80.5 to 81.6]), good calibration and potential clinical benefit for various threshold probabilities. Overall, health care spending patterns improved predictive accuracy the most while also exhibiting superior performance and clinical benefit. Conclusions Patterns of health care spending have the potential to significantly improve assessments on who is at high risk of dying following admission to the ED with a fall. The proposed methodology can assist in predicting the prognosis of fallers, emphasising the added predictive value of longitudinal health-related information next to clinical and sociodemographic predictors. [ABSTRACT FROM AUTHOR] more...

Uncertainty can be classified as either aleatoric (intrinsic randomness) or epistemic (imperfect knowledge of parameters). The majority of frameworks assessing infectious disease risk consider only epistemic uncertainty. We only ever observe a single epidemic, and therefore cannot empirically determine aleatoric uncertainty. Here, we characterise both epistemic and aleatoric uncertainty using a time-varying general branching process. Our framework explicitly decomposes aleatoric variance into mechanistic components, quantifying the contribution to uncertainty produced by each factor in the epidemic process, and how these contributions vary over time. The aleatoric variance of an outbreak is itself a renewal equation where past variance affects future variance. We find that, superspreading is not necessary for substantial uncertainty, and profound variation in outbreak size can occur even without overdispersion in the offspring distribution (i.e. the distribution of the number of secondary infections an infected person produces). Aleatoric forecasting uncertainty grows dynamically and rapidly, and so forecasting using only epistemic uncertainty is a significant underestimate. Therefore, failure to account for aleatoric uncertainty will ensure that policymakers are misled about the substantially higher true extent of potential risk. We demonstrate our method, and the extent to which potential risk is underestimated, using two historical examples. Intrinsic randomness is a critical source of uncertainty in infectious disease outbreaks. The authors show in a series of analytical results how this source of uncertainty can be better characterised. [ABSTRACT FROM AUTHOR] more...

Background: Spousal bereavement is a life event that affects older people differently. We investigated the impact of spousal bereavement on medical expenditures and mortality in the general population, emphasizing on age and sex. Methods: Data are from a population-based, retrospective cohort study following 924,958 Danish citizens over the age of 65 years, within 2011–2016. Changes in health care expenditures in those who suffer bereavement were compared with time matched changes among those who did not. Mortality hazards were analysed with time to event analysis. Results: A total of 77,722 (~8.4%) individuals experienced bereavement, 65.8% being females. Among males, bereavement was associated with increase of expenditures the year after, that was 42 Euros per week (95% CI, 36 to 48) larger than the non-bereaved group. The corresponding increase for females was 35 Euros per week (95% CI, 30 to 40). The increase of mortality hazards was highest in the first year after bereavement, higher in males than females, in young old and almost absent in the oldest old. Compared with the reference, mortality the year after spousal loss was 70% higher (HR 1.70 [95% CI 1.40 to 2.08]) for males aged 65–69 years and remained elevated for a period of six years. Mortality for females aged 65–69 years was 27% higher in the first year (HR 1.27, [1.07 to 1.52]), normalizing thereafter. Conclusion: Bereavement affects older people differently with younger males being most frail with limited recovery potential. [ABSTRACT FROM AUTHOR] more...

Delhi, the national capital of India, experienced multiple severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreaks in 2020 and reached population seropositivity of >50% by 2021. During April 2021, the city became overwhelmed by COVID-19 cases and fatalities, as a new variant, B.1.617.2 (Delta), replaced B.1.1.7 (Alpha). A Bayesian model explains the growth advantage of Delta through a combination of increased transmissibility and reduced sensitivity to immune responses generated against earlier variants (median estimates: 1.5-fold greater transmissibility and 20% reduction in sensitivity). Seropositivity of an employee and family cohort increased from 42% to 87.5% between March and July 2021, with 27% reinfections, as judged by increased antibody concentration after a previous decline. The likely high transmissibility and partial evasion of immunity by the Delta variant contributed to an overwhelming surge in Delhi. [ABSTRACT FROM AUTHOR] more...

Background: Many popular disease transmission models have helped nations respond to the COVID-19 pandemic by informing decisions about pandemic planning, resource allocation, implementation of social distancing measures, lockdowns, and other non-pharmaceutical interventions. We study how five epidemiological models forecast and assess the course of the pandemic in India: a baseline curve-fitting model, an extended SIR (eSIR) model, two extended SEIR (SAPHIRE and SEIR-fansy) models, and a semi-mechanistic Bayesian hierarchical model (ICM).Methods: Using COVID-19 case-recovery-death count data reported in India from March 15 to October 15 to train the models, we generate predictions from each of the five models from October 16 to December 31. To compare prediction accuracy with respect to reported cumulative and active case counts and reported cumulative death counts, we compute the symmetric mean absolute prediction error (SMAPE) for each of the five models. For reported cumulative cases and deaths, we compute Pearson's and Lin's correlation coefficients to investigate how well the projected and observed reported counts agree. We also present underreporting factors when available, and comment on uncertainty of projections from each model.Results: For active case counts, SMAPE values are 35.14% (SEIR-fansy) and 37.96% (eSIR). For cumulative case counts, SMAPE values are 6.89% (baseline), 6.59% (eSIR), 2.25% (SAPHIRE) and 2.29% (SEIR-fansy). For cumulative death counts, the SMAPE values are 4.74% (SEIR-fansy), 8.94% (eSIR) and 0.77% (ICM). Three models (SAPHIRE, SEIR-fansy and ICM) return total (sum of reported and unreported) cumulative case counts as well. We compute underreporting factors as of October 31 and note that for cumulative cases, the SEIR-fansy model yields an underreporting factor of 7.25 and ICM model yields 4.54 for the same quantity. For total (sum of reported and unreported) cumulative deaths the SEIR-fansy model reports an underreporting factor of 2.97. On October 31, we observe 8.18 million cumulative reported cases, while the projections (in millions) from the baseline model are 8.71 (95% credible interval: 8.63-8.80), while eSIR yields 8.35 (7.19-9.60), SAPHIRE returns 8.17 (7.90-8.52) and SEIR-fansy projects 8.51 (8.18-8.85) million cases. Cumulative case projections from the eSIR model have the highest uncertainty in terms of width of 95% credible intervals, followed by those from SAPHIRE, the baseline model and finally SEIR-fansy.Conclusions: In this comparative paper, we describe five different models used to study the transmission dynamics of the SARS-Cov-2 virus in India. While simulation studies are the only gold standard way to compare the accuracy of the models, here we were uniquely poised to compare the projected case-counts against observed data on a test period. The largest variability across models is observed in predicting the "total" number of infections including reported and unreported cases (on which we have no validation data). The degree of under-reporting has been a major concern in India and is characterized in this report. Overall, the SEIR-fansy model appeared to be a good choice with publicly available R-package and desired flexibility plus accuracy. [ABSTRACT FROM AUTHOR] more...

Cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in Manaus, Brazil, resurged in late 2020 despite previously high levels of infection. Genome sequencing of viruses sampled in Manaus between November 2020 and January 2021 revealed the emergence and circulation of a novel SARS-CoV-2 variant of concern. Lineage P.1 acquired 17 mutations, including a trio in the spike protein (K417T, E484K, and N501Y) associated with increased binding to the human ACE2 (angiotensin-converting enzyme 2) receptor. Molecular clock analysis shows that P.1 emergence occurred around mid-November 2020 and was preceded by a period of faster molecular evolution. Using a two-category dynamical model that integrates genomic and mortality data, we estimate that P.1 may be 1.7- to 2.4-fold more transmissible and that previous (non-P.1) infection provides 54 to 79% of the protection against infection with P.1 that it provides against non-P.1 lineages. Enhanced global genomic surveillance of variants of concern, which may exhibit increased transmissibility and/or immune evasion, is critical to accelerate pandemic responsiveness. [ABSTRACT FROM AUTHOR] more...

The ammonium salt of diisopropyl dithiophosphate (A1), O,O′- diisopropanediyl S-(N-phthalimido methyl) dithiophosphate (L1) and zinc diisopropyl dithiophosphate (L2) were isolated in non-aqueous medium and well-characterized by UV–visible spectroscopy, infrared (IR) spectroscopy, NMR (1H, 13C and 31P) spectra and single-crystal X-ray diffraction analysis. The isopropyl group attached to the metal in compound L2 functions as the inter- and intra-chelating group which deviate slightly from planarity with the metal atoms. The central eight-membered ring possesses the "cradle" configuration. The fluorescence spectra of the three compounds A1, L1 and L2 were compared and it was observed that derivatives show excitation due to electron-donor and -acceptor groups. The intramolecular charge transfer (ICT) occurs from P = S to one of the C = O groups in the phthalimide ring and thiophosphate when these molecules are excited electronically in a polar solvent and, therefore, ICT state displays a fluorescence property. The molecular docking study was also carried out for the three compounds. Based on docking studies, a zinc-based metal complex of isopropyl dithiophosphate showed good antifungal activity whereas phthalimide-based isopropyl dithiophosphate derivative showed better antibacterial activity. [ABSTRACT FROM AUTHOR] more...

Developing new methods for modelling infectious diseases outbreaks is important for monitoring transmission and developing policy. In this paper we propose using semi-mechanistic Hawkes Processes for modelling malaria transmission in near-elimination settings. Hawkes Processes are well founded mathematical methods that enable us to combine the benefits of both statistical and mechanistic models to recreate and forecast disease transmission beyond just malaria outbreak scenarios. These methods have been successfully used in numerous applications such as social media and earthquake modelling, but are not yet widespread in epidemiology. By using domain-specific knowledge, we can both recreate transmission curves for malaria in China and Eswatini and disentangle the proportion of cases which are imported from those that are community based. Author summary: This paper introduces a mathematically well-founded method for infectious disease outbreaks known as Hawkes Processes. These semi-mechanistic models are relatively new to the infectious diseases toolkit and enable us to combine disease specific information such as the infectious profile with statistical rigour to recreate temporal disease transmission. We show that these methods are very suited to modelling malaria in communities close to eliminating malaria—in particular China and Eswatini—where we are able to disentangle the contribution of exogenous (external) transmission and endogenous (person-to-person) transmission. This is particularly important for developing policies when counties are approaching elimination. [ABSTRACT FROM AUTHOR] more...

Background: There were 25.6 million attendances at Emergency Departments (EDs) in England in 2019 corresponding to an increase of 12 million attendances over the past ten years. The steadily rising demand at EDs creates a constant challenge to provide adequate quality of care while maintaining standards and productivity. Managing hospital demand effectively requires an adequate knowledge of the future rate of admission. We develop a novel predictive framework to understand the temporal dynamics of hospital demand.Methods: We compare and combine state-of-the-art forecasting methods to predict hospital demand 1, 3 or 7 days into the future. In particular, our analysis compares machine learning algorithms to more traditional linear models as measured in a mean absolute error (MAE) and we consider two different hyperparameter tuning methods, enabling a faster deployment of our models without compromising performance. We believe our framework can readily be used to forecast a wide range of policy relevant indicators.Results: We find that linear models often outperform machine learning methods and that the quality of our predictions for any of the forecasting horizons of 1, 3 or 7 days are comparable as measured in MAE. Our approach is able to predict attendances at these emergency departments one day in advance up to a mean absolute error of ±14 and ±10 patients corresponding to a mean absolute percentage error of 6.8% and 8.6% respectively.Conclusions: Simple linear methods like generalized linear models are often better or at least as good as ensemble learning methods like the gradient boosting or random forest algorithm. However, though sophisticated machine learning methods are not necessarily better than linear models, they improve the diversity of model predictions so that stacked predictions can be more robust than any single model including the best performing one. [ABSTRACT FROM AUTHOR] more...

In order to monitor progress towards malaria elimination, it is crucial to be able to measure changes in spatio-temporal transmission. However, common metrics of malaria transmission such as parasite prevalence are under powered in elimination contexts. China has achieved major reductions in malaria incidence and is on track to eliminate, having reporting zero locally-acquired malaria cases in 2017 and 2018. Understanding the spatio-temporal pattern underlying this decline, especially the relationship between locally-acquired and imported cases, can inform efforts to maintain elimination and prevent re-emergence. This is particularly pertinent in Yunnan province, where the potential for local transmission is highest. Using a geo-located individual-level dataset of cases recorded in Yunnan province between 2011 and 2016, we introduce a novel Bayesian framework to model a latent diffusion process and estimate the joint likelihood of transmission between cases and the number of cases with unobserved sources of infection. This is used to estimate the case reproduction number, Rc. We use these estimates within spatio-temporal geostatistical models to map how transmission varied over time and space, estimate the timeline to elimination and the risk of resurgence. We estimate the mean Rc between 2011 and 2016 to be 0.171 (95% CI = 0.165, 0.178) for P. vivax cases and 0.089 (95% CI = 0.076, 0.103) for P. falciparum cases. From 2014 onwards, no cases were estimated to have a Rc value above one. An unobserved source of infection was estimated to be moderately likely (p>0.5) for 19/ 611 cases and high (p>0.8) for 2 cases, suggesting very high levels of case ascertainment. Our estimates suggest that, maintaining current intervention efforts, Yunnan is unlikely to experience sustained local transmission up to 2020. However, even with a mean of 0.005 projected up to 2020, locally-acquired cases are possible due to high levels of importation. Author summary: Although malaria is still responsible for a great deal of death and illness in many parts of the world, many national control programmes have made great strides in controlling malaria and now are in a position to aim for elimination. However, in order to monitor progress towards elimination and plan interventions, it is crucial to measure malaria transmission and how it varies over space and time. However, traditional metrics used to measure malaria transmission are not suitable in elimination settings. China is one example of a country approaching elimination, with aims to eliminate the disease by 2020. Using a detailed individual level dataset of the times and locations of people showing symptoms of malaria, we use approaches adapted from the study of how information spreads through social networks to estimate the likelihood of transmission occurring between cases. This information is used to estimate how many people we expect each case to go on to infect. In elimination settings, this number is an indication of how quickly elimination will be reached and how likely we are to see a resurgence in cases once elimination is achieved. Our results show a decline in this metric over time, as well as seasonal changes in transmission which are different to the patterns in when the most cases were observed. [ABSTRACT FROM AUTHOR] more...

If herd immunity had already been reached, we would expect no correlation, or even a negative correlation, as lockdown would not alter the herd immunity threshold in the population or the ultimate death rate per capita. Taking Spain as an example, for the country to have achieved herd immunity, one would have to assume that the herd immunity threshold differs by a factor of ten between regions. [Extracted from the article] more...