Piecewise Approximate Bayesian Computation: fast inference for discretely observed Markov models using a factorised posterior distribution

Many modern statistical applications involve inference for complicated stochastic models for which the likelihood function is difficult or even impossible to calculate, and hence conventional likelihood-based inferential techniques cannot be used. In such settings, Bayesian inference can be performed using Approximate Bayesian Computation (ABC). However, in spite of many recent developments to ABC methodology, in many applications the computational cost of ABC necessitates the choice of summary statistics and tolerances that can potentially severely bias the estimate of the posterior. We propose a new “piecewise” ABC approach suitable for discretely observed Markov models that involves writing the posterior density of the parameters as a product of factors, each a function of only a subset of the data, and then using ABC within each factor. The approach has the advantage of side-stepping the need to choose a summary statistic and it enables a stringent tolerance to be set, making the posterior “less approximate”. We investigate two methods for estimating the posterior density based on ABC samples for each of the factors: the first is to use a Gaussian approximation for each factor, and the second is to use a kernel density estimate. Both methods have their merits. The Gaussian approximation is simple, fast, and probably adequate for many applications. On the other hand, using instead a kernel density estimate has the benefit of consistently estimating the true piecewise ABC posterior as the number of ABC samples tends to infinity. We illustrate the piecewise ABC approach with four examples; in each case, the approach offers fast and accurate inference

Statistical Mechanics of Horizontal Gene Transfer in Evolutionary Ecology

The biological world, especially its majority microbial component, is strongly interacting and may be dominated by collective effects. In this review, we provide a brief introduction for statistical physicists of the way in which living cells communicate genetically through transferred genes, as well as the ways in which they can reorganize their genomes in response to environmental pressure. We discuss how genome evolution can be thought of as related to the physical phenomenon of annealing, and describe the sense in which genomes can be said to exhibit an analogue of information entropy. As a direct application of these ideas, we analyze the variation with ocean depth of transposons in marine microbial genomes, predicting trends that are consistent with recent observations using metagenomic surveys.Comment: Accepted by Journal of Statistical Physic

A multi-ancestry genome-wide study incorporating gene-smoking interactions identifies multiple new loci for pulse pressure and mean arterial pressure

Elevated blood pressure (BP), a leading cause of global morbidity and mortality, is influenced by both genetic and lifestyle factors. Cigarette smoking is one such lifestyle factor. Across five ancestries, we performed a genome-wide gene–smoking interaction study of mean arterial pressure (MAP) and pulse pressure (PP) in 129 913 individuals in stage 1 and follow-up analysis in 480 178 additional individuals in stage 2. We report here 136 loci significantly associated with MAP and/or PP. Of these, 61 were previously published through main-effect analysis of BP traits, 37 were recently reported by us for systolic BP and/or diastolic BP through gene–smoking interaction analysis and 38 were newly identified (P < 5 × 10−8, false discovery rate < 0.05). We also identified nine new signals near known loci. Of the 136 loci, 8 showed significant interaction with smoking status. They include CSMD1 previously reported for insulin resistance and BP in the spontaneously hypertensive rats. Many of the 38 new loci show biologic plausibility for a role in BP regulation. SLC26A7 encodes a chloride/bicarbonate exchanger expressed in the renal outer medullary collecting duct. AVPR1A is widely expressed, including in vascular smooth muscle cells, kidney, myocardium and brain. FHAD1 is a long non-coding RNA overexpressed in heart failure. TMEM51 was associated with contractile function in cardiomyocytes. CASP9 plays a central role in cardiomyocyte apoptosis. Identified only in African ancestry were 30 novel loci. Our findings highlight the value of multi-ancestry investigations, particularly in studies of interaction with lifestyle factors, where genomic and lifestyle differences may contribute to novel findings

Inclusive-photon production and its dependence on photon isolation in pp collisions at s√ = 13 TeV using 139 fb−1 of ATLAS data

Measurements of differential cross sections are presented for inclusive isolated-photon production in pp collisions at a centre-of-mass energy of 13 TeV provided by the LHC and using 139 fb−1 of data recorded by the ATLAS experiment. The cross sections are measured as functions of the photon transverse energy in different regions of photon pseudorapidity. The photons are required to be isolated by means of a fixed-cone method with two different cone radii. The dependence of the inclusive-photon production on the photon isolation is investigated by measuring the fiducial cross sections as functions of the isolation-cone radius and the ratios of the differential cross sections with different radii in different regions of photon pseudorapidity. The results presented in this paper constitute an improvement with respect to those published by ATLAS earlier: the measurements are provided for different isolation radii and with a more granular segmentation in photon pseudorapidity that can be exploited in improving the determination of the proton parton distribution functions. These improvements provide a more in-depth test of the theoretical predictions. Next-to-leading-order QCD predictions from JETPHOX and SHERPA and next-to-next-to-leading-order QCD predictions from NNLOJET are compared to the measurements, using several parameterisations of the proton parton distribution functions. The measured cross sections are well described by the fixed-order QCD predictions within the experimental and theoretical uncertainties in most of the investigated phase-space region

Measurements of Zγ+jets differential cross sections in pp collisions at s√ = 13 TeV with the ATLAS detector

Differential cross-section measurements of Zγ production in association with hadronic jets are presented, using the full 139 fb−1 dataset of s√ = 13 TeV proton–proton collisions collected by the ATLAS detector during Run 2 of the LHC. Distributions are measured using events in which the Z boson decays leptonically and the photon is usually radiated from an initial-state quark. Measurements are made in both one and two observables, including those sensitive to the hard scattering in the event and others which probe additional soft and collinear radiation. Different Standard Model predictions, from both parton-shower Monte Carlo simulation and fixed-order QCD calculations, are compared with the measurements. In general, good agreement is observed between data and predictions from MATRIX and MiNNLOPS, as well as next-to-leading-order predictions from MADGRAPH5_AMC@NLO and SHERPA

Evidence for the charge asymmetry in pp → tt¯ production at s√ = 13 TeV with the ATLAS detector

Inclusive and differential measurements of the top–antitop (tt¯) charge asymmetry Att¯C and the leptonic asymmetry Aℓℓ¯C are presented in proton–proton collisions at s√ = 13 TeV recorded by the ATLAS experiment at the CERN Large Hadron Collider. The measurement uses the complete Run 2 dataset, corresponding to an integrated luminosity of 139 fb−1, combines data in the single-lepton and dilepton channels, and employs reconstruction techniques adapted to both the resolved and boosted topologies. A Bayesian unfolding procedure is performed to correct for detector resolution and acceptance effects. The combined inclusive tt¯ charge asymmetry is measured to be Att¯C = 0.0068 ± 0.0015, which differs from zero by 4.7 standard deviations. Differential measurements are performed as a function of the invariant mass, transverse momentum and longitudinal boost of the tt¯ system. Both the inclusive and differential measurements are found to be compatible with the Standard Model predictions, at next-to-next-to-leading order in quantum chromodynamics perturbation theory with next-to-leading-order electroweak corrections. The measurements are interpreted in the framework of the Standard Model effective field theory, placing competitive bounds on several Wilson coefficients

Search for heavy resonances decaying into a Z or W boson and a Higgs boson in final states with leptons and b-jets in 139 fb−1 of pp collisions at s√ = 13 TeV with the ATLAS detector

This article presents a search for new resonances decaying into a Z or W boson and a 125 GeV Higgs boson h, and it targets the νν¯¯¯bb¯¯, ℓ+ℓ−bb¯¯, or ℓ±νbb¯¯ final states, where ℓ = e or μ, in proton-proton collisions at s√ = 13 TeV. The data used correspond to a total integrated luminosity of 139 fb−1 collected by the ATLAS detector during Run 2 of the LHC at CERN. The search is conducted by examining the reconstructed invariant or transverse mass distributions of Zh or Wh candidates for evidence of a localised excess in the mass range from 220 GeV to 5 TeV. No significant excess is observed and 95% confidence-level upper limits between 1.3 pb and 0.3 fb are placed on the production cross section times branching fraction of neutral and charged spin-1 resonances and CP-odd scalar bosons. These limits are converted into constraints on the parameter space of the Heavy Vector Triplet model and the two-Higgs-doublet model

The ATLAS trigger system for LHC Run 3 and trigger performance in 2022

The ATLAS trigger system is a crucial component of the ATLAS experiment at the LHC. It is responsible for selecting events in line with the ATLAS physics programme. This paper presents an overview of the changes to the trigger and data acquisition system during the second long shutdown of the LHC, and shows the performance of the trigger system and its components in the proton-proton collisions during the 2022 commissioning period as well as its expected performance in proton-proton and heavy-ion collisions for the remainder of the third LHC data-taking period (2022–2025)