Statistical Practice . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1, 124, 254, 297 History Corner . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .20, 179 Teacher’s Corner . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26, 173, 263, 335 General . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40, 147, 211, 366 Statistical Computing and Graphics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56 Statistical Computing and Software Reviews . . . . . . . . . . . . . . . . . . . . . . . . 75, 187 Reviews of Books and Teaching Materials . . . . . . . . . . . . . . . . . 92, 189, 281, 401 Brief Reviews . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100, 195, 292, 404 Letters to the Editor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 102, 197, 294, 406 Special Section: Statistical Training and Curricular Revision . . . . . . . . . . . . . 105 Errata . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 200 Special Section: Opportunities and Challenges for the Discipline . . . . . . . . . 201 Software Review . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40

The American statistician

The HI 21cm absorption line is masked by bright foregrounds and systematic distortions that arise due to the chromaticity of the antenna used to make the observation coupling to the spectral inhomogeneity of these foregrounds. We demonstrate that these distortions are sufficient to conceal the 21cm signal when the antenna is not perfectly achromatic and that simple corrections assuming a constant spatial distribution of foreground power are insufficient to overcome them. We then propose a new physics-motivated method of modelling the foregrounds of 21cm experiments in order to fit the chromatic distortions as part of the foregrounds. This is done by generating a simulated sky model across the observing band by dividing the sky into $N$ regions and scaling a base map assuming a distinct uniform spectral index in each region. The resulting sky map can then be convolved with a model of the antenna beam to give a model of foregrounds and chromaticity parameterised by the spectral indices of the $N$ regions. We demonstrate that fitting this model for varying $N$ using a Bayesian nested sampling algorithm and comparing the results using the evidence allows the 21cm signal to be reliably detected in data of a relatively smooth conical log spiral antenna. We also test a much more chromatic conical sinuous antenna and find this model will not produce a reliable signal detection, but in a manner that is easily distinguishable from a true detection.

A general Bayesian framework for foreground modelling and chromaticity correction for global 21 cm experiments

We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the application of NS in three scientific fields in which the algorithm has proved to be useful: cosmology, gravitational-wave astronomy, and materials science. We close by making recommendations for best practice when using NS and by summarizing potential limitations and optimizations of NS. 

/pdf/nested-sampling-for-physical-scientists-2v0bgpsa.pdf

Nested sampling for physical scientists

Maximally Smooth Functions (MSFs) are a form of constrained functions in which there are no inflection points or zero crossings in high order derivatives. Consequently, they have applications to signal recovery in experiments where signals of interest are expected to be non-smooth features masked by larger smooth signals or foregrounds. They can also act as a powerful tool for diagnosing the presence of systematics. The constrained nature of MSFs makes fitting these functions a non-trivial task. We introduce maxsmooth, an open source package that uses quadratic programming to rapidly fit MSFs. We demonstrate the efficiency and reliability of maxsmooth by comparison to commonly used fitting routines. We show that by using quadratic programming we can reduce the fitting time by approximately two orders of magnitude. maxsmooth features a built-in library of MSF models and allows the user to define their own. We also introduce and implement with maxsmooth Partially Smooth Functions, which are useful for describing elements of non-smooth structure in foregrounds. This work has been motivated by the problem of foreground modelling in 21-cm cosmology for which MSFs have been shown to be a viable alternative to polynomials. We discuss applications of maxsmooth to 21-cm cosmology and highlight this with examples using data from the Experiment to Detect the Global Epoch of Reionization Signature (EDGES) and the Large-aperture Experiment to Detect the Dark Ages (LEDA) experiments. MSFs are applied to data from LEDA for the first time in this paper. maxsmooth is pip installable and available for download at:this https URL

MAXSMOOTH: rapid maximally smooth function fitting with applications in Global 21-cm cosmology

Observations of the 21-cm line from primordial hydrogen promise to be one of the best tools to study the early epochs of the Universe: the dark ages, the cosmic dawn and the subsequent epoch of reionization. In 2018, the Experiment to Detect the Global Epoch of Reionization Signature (EDGES) caught the attention of the cosmology community with a potential detection of an absorption feature in the sky-averaged radio spectrum centred at 78 MHz. The feature is deeper than expected, and, if confirmed, would call for new physics. However, different groups have re-analysed the EDGES data and questioned the reliability of the signal. The Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) is a sky-averaged 21-cm experiment aiming at improving the current observations by tackling the issues faced by current instruments related to residual systematic signals in the data. The novel experimental approach focuses on detecting and jointly explaining these systematics together with the foregrounds and the cosmological signal using Bayesian statistics. To achieve this, REACH features simultaneous observations with two different antennas, an ultra-wideband system (redshift range about 7.5 to 28) and a receiver calibrator based on in-field measurements. Simulated observations forecast percent-level constraints on astrophysical parameters, potentially opening up a new window to the infant Universe. The Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) aims to detect the sky-averaged 21-cm neutral hydrogen line from the early Universe by jointly analysing the cosmological signal, foreground emission and systematic effects. 

/pdf/the-reach-radiometer-for-detecting-the-21-cm-hydrogen-signal-1qrrv5lq.pdf

The REACH radiometer for detecting the 21-cm hydrogen signal from redshift z ≈ 7.5–28

We modelled the two major layer of Earth's ionosphere, the F-layer and the D-layer, by a simplified spatial model with temporal variance to study the chromatic ionospheric effects on global 21-cm observations. From the analyses, we found that the magnitude of the ionospheric disruptions due to ionospheric refraction and absorption can be greater than the expected global 21-cm signal, and the variation of its magnitude can differ, depending on the ionospheric conditions. Within the parameter space adopted in the model, the shape of the global 21-cm signal is distorted after propagating through the ionosphere, while its amplitude is weakened. It is observed that the ionospheric effects do not cancel out over time, and thus should be accounted for in the foreground calibration at each timestep to account for the chromaticity introduced by the ionosphere.

Quantifying ionospheric effects on global 21-cm observations

Global 21cm experiments aim to measure the sky averaged HI absorption signal from cosmic dawn and the epoch of reionisation. However, antenna chromaticity coupling to bright foregrounds can introduce distortions into the observational data of such experiments. We demonstrate a method for guiding the antenna design of a global experiment through data analysis simulations. This is done by performing simulated observations for a range of inserted 21cm signals, then attempting to identify the signals with a data analysis pipeline. We demonstrate this method on five antennae that were considered as potential designs for the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH); a conical log spiral antenna, an inverted conical sinuous antenna and polygonal-, square- and elliptical-bladed dipoles. We find that the log spiral performs significantly better than the other antennae tested, able to correctly and confidently identify every inserted 21cm signal. In second place is the polygonal dipole antenna, which was only unable to detect signals with both very low amplitudes of 0.05K and low centre frequency of 80MHz. The conical sinuous antenna was found to perform least accurately, only able to detect the highest amplitude 21cm signals, and even then with biases. We also demonstrate that, due to the non-trivial nature of chromatic distortion and the processes of correcting for it, these are not the results that could have been expected superficially from the extent of chromatic variation in each antenna.

Informing antenna design for sky-averaged 21-cm experiments using a simulated Bayesian data analysis pipeline

Following the reported detection of an absorption profile associated with the 21~cm sky-averaged signal from the Cosmic Dawn by the EDGES experiment in 2018, a number of experiments have been set up to verify this result. This paper discusses the design process used for global 21~cm experiments, focusing specifically on the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH). This experiment will seek to understand and compensate for systematic errors present using detailed modelling and characterization of the instrumentation. There is detailed the quantitative figures of merit and numerical modelling used to assist the design process of the REACH dipole antenna (one of the 2 antenna designs for REACH Phase I). This design process produced a 2.5:1 frequency bandwidth dipole. The aim of this design was to balance spectral smoothness and low impedance reflections with the ability to describe and understand the antenna response to the sky signal to inform the critically important calibration during observation and data analysis.

Dominic Anstey

Papers

A general Bayesian framework for foreground modelling and chromaticity correction for global 21 cm experiments

Quantifying ionospheric effects on global 21-cm observations

Informing antenna design for sky-averaged 21-cm experiments using a simulated Bayesian data analysis pipeline

Radio antenna design for sky-averaged 21 cm cosmology experiments: the REACH case