-
Intruder band mixing in an ab initio description of 12Be
Authors:
Anna E. McCoy,
Mark A. Caprio,
Pieter Maris,
Patrick J. Fasano
Abstract:
The spectrum of 12Be exhibits exotic features, e.g., an intruder ground state and shape coexistence, normally associated with the breakdown of a shell closure. While previous phenomenological treatments indicated the ground state has substantial contributions from intruder configurations, it is only with advances in computational abilities and improved interactions that this intruder mixing is obs…
▽ More
The spectrum of 12Be exhibits exotic features, e.g., an intruder ground state and shape coexistence, normally associated with the breakdown of a shell closure. While previous phenomenological treatments indicated the ground state has substantial contributions from intruder configurations, it is only with advances in computational abilities and improved interactions that this intruder mixing is observed in ab initio no-core shell model (NCSM) predictions. In this work, we extract electromagnetic observables and symmetry decompositions from the NCSM wave functions to demonstrate that the low-lying positive parity spectrum can be explained in terms of mixing of rotational bands with very different intrinsic structure coexisting within the low-lying spectrum. These observed bands exhibit an approximate SU(3) symmetry and are qualitatively consistent with Elliott model predictions.
△ Less
Submitted 19 February, 2024;
originally announced February 2024.
-
Explainable Machine Learning for Categorical and Mixed Data with Lossless Visualization
Authors:
Boris Kovalerchuk,
Elijah McCoy
Abstract:
Building accurate and interpretable Machine Learning (ML) models for heterogeneous/mixed data is a long-standing challenge for algorithms designed for numeric data. This work focuses on developing numeric coding schemes for non-numeric attributes for ML algorithms to support accurate and explainable ML models, methods for lossless visualization of n-D non-numeric categorical data with visual rule…
▽ More
Building accurate and interpretable Machine Learning (ML) models for heterogeneous/mixed data is a long-standing challenge for algorithms designed for numeric data. This work focuses on developing numeric coding schemes for non-numeric attributes for ML algorithms to support accurate and explainable ML models, methods for lossless visualization of n-D non-numeric categorical data with visual rule discovery in these visualizations, and accurate and explainable ML models for categorical data. This study proposes a classification of mixed data types and analyzes their important role in Machine Learning. It presents a toolkit for enforcing interpretability of all internal operations of ML algorithms on mixed data with a visual data exploration on mixed data. A new Sequential Rule Generation (SRG) algorithm for explainable rule generation with categorical data is proposed and successfully evaluated in multiple computational experiments. This work is one of the steps to the full scope ML algorithms for mixed data supported by lossless visualization of n-D data in General Line Coordinates beyond Parallel Coordinates.
△ Less
Submitted 22 November, 2023; v1 submitted 28 May, 2023;
originally announced May 2023.
-
Nuclear $β$ decay as a probe for physics beyond the Standard Model
Authors:
M. Brodeur,
N. Buzinsky,
M. A. Caprio,
V. Cirigliano,
J. A. Clark,
P. J. Fasano,
J. A. Formaggio,
A. T. Gallant,
A. Garcia,
S. Gandolfi,
S. Gardner,
A. Glick-Magid,
L. Hayen,
H. Hergert,
J. D. Holt,
M. Horoi,
M. Y. Huang,
K. D. Launey,
K. G. Leach,
B. Longfellow,
A. Lovato,
A. E. McCoy,
D. Melconian,
P. Mohanmurthy,
D. C. Moore
, et al. (21 additional authors not shown)
Abstract:
This white paper was submitted to the 2022 Fundamental Symmetries, Neutrons, and Neutrinos (FSNN) Town Hall Meeting in preparation for the next NSAC Long Range Plan. We advocate to support current and future theoretical and experimental searches for physics beyond the Standard Model using nuclear $β$ decay.
This white paper was submitted to the 2022 Fundamental Symmetries, Neutrons, and Neutrinos (FSNN) Town Hall Meeting in preparation for the next NSAC Long Range Plan. We advocate to support current and future theoretical and experimental searches for physics beyond the Standard Model using nuclear $β$ decay.
△ Less
Submitted 10 January, 2023;
originally announced January 2023.
-
Measuring e-Commerce Metric Changes in Online Experiments
Authors:
C. H. Bryan Liu,
Emma J. McCoy
Abstract:
Digital technology organizations routinely use online experiments (e.g. A/B tests) to guide their product and business decisions. In e-commerce, we often measure changes to transaction- or item-based business metrics such as Average Basket Value (ABV), Average Basket Size (ABS), and Average Selling Price (ASP); yet it remains a common pitfall to ignore the dependency between the value/size of tran…
▽ More
Digital technology organizations routinely use online experiments (e.g. A/B tests) to guide their product and business decisions. In e-commerce, we often measure changes to transaction- or item-based business metrics such as Average Basket Value (ABV), Average Basket Size (ABS), and Average Selling Price (ASP); yet it remains a common pitfall to ignore the dependency between the value/size of transactions/items during experiment design and analysis. We present empirical evidence on such dependency, its impact on measurement uncertainty, and practical implications on A/B test outcomes if left unmitigated. By making the evidence available, we hope to drive awareness of the pitfall among experimenters in e-commerce and hence encourage the adoption of established mitigation approaches. We also share lessons learned when incorporating selected mitigation approaches into our experimentation analysis platform currently in production.
△ Less
Submitted 17 April, 2023; v1 submitted 31 October, 2022;
originally announced October 2022.
-
Can a latent Hawkes process be used for epidemiological modelling?
Authors:
Stamatina Lamprinakou,
Axel Gandy,
Emma McCoy
Abstract:
Understanding the spread of COVID-19 has been the subject of numerous studies, highlighting the significance of reliable epidemic models. Here, we introduce a novel epidemic model using a latent Hawkes process with temporal covariates for modelling the infections. Unlike other models, we model the reported cases via a probability distribution driven by the underlying Hawkes process. Modelling the…
▽ More
Understanding the spread of COVID-19 has been the subject of numerous studies, highlighting the significance of reliable epidemic models. Here, we introduce a novel epidemic model using a latent Hawkes process with temporal covariates for modelling the infections. Unlike other models, we model the reported cases via a probability distribution driven by the underlying Hawkes process. Modelling the infections via a Hawkes process allows us to estimate by whom an infected individual was infected. We propose a Kernel Density Particle Filter (KDPF) for inference of both latent cases and reproduction number and for predicting the new cases in the near future. The computational effort is proportional to the number of infections making it possible to use particle filter type algorithms, such as the KDPF. We demonstrate the performance of the proposed algorithm on synthetic data sets and COVID-19 reported cases in various local authorities in the UK, and benchmark our model to alternative approaches.
△ Less
Submitted 15 August, 2022;
originally announced August 2022.
-
Explainable Mixed Data Representation and Lossless Visualization Toolkit for Knowledge Discovery
Authors:
Boris Kovalerchuk,
Elijah McCoy
Abstract:
Developing Machine Learning (ML) algorithms for heterogeneous/mixed data is a longstanding problem. Many ML algorithms are not applicable to mixed data, which include numeric and non-numeric data, text, graphs and so on to generate interpretable models. Another longstanding problem is developing algorithms for lossless visualization of multidimensional mixed data. The further progress in ML heavil…
▽ More
Developing Machine Learning (ML) algorithms for heterogeneous/mixed data is a longstanding problem. Many ML algorithms are not applicable to mixed data, which include numeric and non-numeric data, text, graphs and so on to generate interpretable models. Another longstanding problem is developing algorithms for lossless visualization of multidimensional mixed data. The further progress in ML heavily depends on success interpretable ML algorithms for mixed data and lossless interpretable visualization of multidimensional data. The later allows developing interpretable ML models using visual knowledge discovery by end-users, who can bring valuable domain knowledge which is absent in the training data. The challenges for mixed data include: (1) generating numeric coding schemes for non-numeric attributes for numeric ML algorithms to provide accurate and interpretable ML models, (2) generating methods for lossless visualization of n-D non-numeric data and visual rule discovery in these visualizations. This paper presents a classification of mixed data types, analyzes their importance for ML and present the developed experimental toolkit to deal with mixed data. It combines the Data Types Editor, VisCanvas data visualization and rule discovery system which is available on GitHub.
△ Less
Submitted 13 June, 2022;
originally announced June 2022.
-
Symmetry and shape coexistence in 10Be
Authors:
M. A. Caprio,
A. E. McCoy,
P. J. Fasano,
T. Dytrych
Abstract:
Within the low-lying spectrum of 10Be, multiple rotational bands are found, with strikingly different moments of inertia. A proposed interpretation has been that these bands variously represent triaxial rotation and prolate axially-deformed rotation. The bands are well-reproduced in ab initio no-core configuration interaction (NCCI) calculations. We use the calculated wave functions to elucidate t…
▽ More
Within the low-lying spectrum of 10Be, multiple rotational bands are found, with strikingly different moments of inertia. A proposed interpretation has been that these bands variously represent triaxial rotation and prolate axially-deformed rotation. The bands are well-reproduced in ab initio no-core configuration interaction (NCCI) calculations. We use the calculated wave functions to elucidate the nuclear shapes underlying these bands, by examining the Elliott SU(3) symmetry content of these wave functions. The ab initio results support an interpretation in which the ground-state band, along with an accompanying K=2 side band, represent a triaxial rotor, arising from an SU(3) irreducible representation in the 0 hbar omega space. Then, the lowest excited K=0 band represents a prolate rotor, arising from an SU(3) irreducible representation in the 2 hbar omega space.
△ Less
Submitted 7 December, 2021;
originally announced December 2021.
-
Datasets for Online Controlled Experiments
Authors:
C. H. Bryan Liu,
Ângelo Cardoso,
Paul Couturier,
Emma J. McCoy
Abstract:
Online Controlled Experiments (OCE) are the gold standard to measure impact and guide decisions for digital products and services. Despite many methodological advances in this area, the scarcity of public datasets and the lack of a systematic review and categorization hinder its development. We present the first survey and taxonomy for OCE datasets, which highlight the lack of a public dataset to…
▽ More
Online Controlled Experiments (OCE) are the gold standard to measure impact and guide decisions for digital products and services. Despite many methodological advances in this area, the scarcity of public datasets and the lack of a systematic review and categorization hinder its development. We present the first survey and taxonomy for OCE datasets, which highlight the lack of a public dataset to support the design and running of experiments with adaptive stopping, an increasingly popular approach to enable quickly deploying improvements or rolling back degrading changes. We release the first such dataset, containing daily checkpoints of decision metrics from multiple, real experiments run on a global e-commerce platform. The dataset design is guided by a broader discussion on data requirements for common statistical tests used in digital experimentation. We demonstrate how to use the dataset in the adaptive stopping scenario using sequential and Bayesian hypothesis tests and learn the relevant parameters for each approach.
△ Less
Submitted 14 January, 2022; v1 submitted 19 November, 2021;
originally announced November 2021.
-
The First Interstellar Astronauts Will Not Be Human
Authors:
Stephen Lantin,
Sophie Mendell,
Ghassan Akkad,
Alexander N. Cohen,
Xander Apicella,
Emma McCoy,
Eliana Beltran-Pardo,
Michael Waltemathe,
Prasanna Srinivasan,
Pradeep M. Joshi,
Joel H. Rothman,
Philip Lubin
Abstract:
Our ability to explore the cosmos by direct contact has been limited to a small number of lunar and interplanetary missions. However, the NASA Starlight program points a path forward to send small, relativistic spacecraft far outside our solar system via standoff directed-energy propulsion. These miniaturized spacecraft are capable of robotic exploration but can also transport seeds and organisms,…
▽ More
Our ability to explore the cosmos by direct contact has been limited to a small number of lunar and interplanetary missions. However, the NASA Starlight program points a path forward to send small, relativistic spacecraft far outside our solar system via standoff directed-energy propulsion. These miniaturized spacecraft are capable of robotic exploration but can also transport seeds and organisms, marking a profound change in our ability to both characterize and expand the reach of known life. Here we explore the biological and technological challenges of interstellar space biology, focusing on radiation-tolerant microorganisms capable of cryptobiosis. Additionally, we discuss planetary protection concerns and other ethical considerations of sending life to the stars.
△ Less
Submitted 25 October, 2021;
originally announced October 2021.
-
First Measurement of the $B(E2; 3/2^- \rightarrow 1/2^-)$ Transition Strength in $^7$Be: Testing Ab Initio Predictions for $A=7$ Nuclei
Authors:
S. L. Henderson,
T. Ahn,
M. A. Caprio,
P. J. Fasano,
A. Simon,
W. Tan,
P. O'Malley,
J. Allen,
D. W. Bardayan,
D. Blankstein,
B. Frentz,
M. R. Hall,
J. J. Kolata,
A. E. McCoy,
S. Moylan,
C. S. Reingold,
S. Y. Strauss,
R. O. Torres-Isea
Abstract:
Electromagnetic observables are able to give insight into collective and emergent features in nuclei, including nuclear clustering. These observables also provide strong constraints for ab initio theory, but comparison of these observables between theory and experiment can be difficult due to the lack of convergence for relevant calculated values, such as $E2$ transition strengths. By comparing th…
▽ More
Electromagnetic observables are able to give insight into collective and emergent features in nuclei, including nuclear clustering. These observables also provide strong constraints for ab initio theory, but comparison of these observables between theory and experiment can be difficult due to the lack of convergence for relevant calculated values, such as $E2$ transition strengths. By comparing the ratios of $E2$ transition strengths for mirror transitions, we find that a wide range of ab initio calculations give robust and consistent predictions for this ratio. To experimentally test the validity of these ab initio predictions, we performed a Coulomb excitation experiment to measure the $B(E2; 3/2^- \rightarrow 1/2^-)$ transition strength in $^7$Be for the first time. A $B(E2; 3/2^- \rightarrow 1/2^-)$ value of $26(6)(3) \, e^2 \mathrm{fm}^4$ was deduced from the measured Coulomb excitation cross section. This result is used with the experimentally known $^7$Li $B(E2; 3/2^- \rightarrow 1/2^-)$ value to provide an experimental ratio to compare with the ab initio predictions. Our experimental value is consistent with the theoretical ratios within $1 σ$ uncertainty, giving experimental support for the value of these ratios. Further work in both theory and experiment can give insight into the robustness of these ratios and their physical meaning.
△ Less
Submitted 15 September, 2021;
originally announced September 2021.
-
Remeasuring the anomalously enhanced $B(E2; 2^+ \rightarrow 1^+)$ in $^8\mathrm{Li}$
Authors:
S. L. Henderson,
T. Ahn,
P. J. Fasano,
A. E. McCoy,
S. Aguilar,
D. T. Blankstein,
L. Caves,
A. C. Dombos,
R. K. Grzywacz,
K. L. Jones,
S. Jin,
R. Kelmar,
J. J. Kolata,
P. D. O'Malley,
C. S. Reingold,
A. Simon,
K. Smith
Abstract:
The large reported $E2$ strength between the $2^+$ ground state and $1^+$ first excited state of $\isotope[8]{Li}$, $B(E2; 2^+ \rightarrow 1^+)= 55(15)\,e^2\fm^4$, presents a puzzle. Unlike in neighboring $A=7\text{--}9$ isotopes, where enhanced $E2$ strengths may be understood to arise from deformation as rotational in-band transitions, the $2^+\rightarrow1^+$ transition in $^8$Li cannot be under…
▽ More
The large reported $E2$ strength between the $2^+$ ground state and $1^+$ first excited state of $\isotope[8]{Li}$, $B(E2; 2^+ \rightarrow 1^+)= 55(15)\,e^2\fm^4$, presents a puzzle. Unlike in neighboring $A=7\text{--}9$ isotopes, where enhanced $E2$ strengths may be understood to arise from deformation as rotational in-band transitions, the $2^+\rightarrow1^+$ transition in $^8$Li cannot be understood in any simple way as a rotational in-band transition. Moreover, the reported strength exceeds \textit{ab initio} predictions by an order of magnitude. In light of this discrepancy, we revisited the Coulomb excitation measurement of this strength, now using particle-$γ$ coincidences, yielding a revised $B(E2; 2^+ \rightarrow 1^+)$ of $19(^{+7}_{-6})(2)$~e$^2$fm$^4$. We explore how this value compares to what might be expected in the limits of rotational models. While the present value is about a factor of three smaller than previously reported, it remains anomalously enhanced.
△ Less
Submitted 5 May, 2023; v1 submitted 13 September, 2021;
originally announced September 2021.
-
Causal Analysis at Extreme Quantiles with Application to London Traffic Flow Data
Authors:
Prajamitra Bhuyan,
Kaushik Jana,
Emma J. McCoy
Abstract:
Transport engineers employ various interventions to enhance traffic-network performance. Quantifying the impacts of Cycle Superhighways is complicated due to the non-random assignment of such an intervention over the transport network. Treatment effects on asymmetric and heavy-tailed distributions are better reflected at extreme tails rather than at the median. We propose a novel method to estimat…
▽ More
Transport engineers employ various interventions to enhance traffic-network performance. Quantifying the impacts of Cycle Superhighways is complicated due to the non-random assignment of such an intervention over the transport network. Treatment effects on asymmetric and heavy-tailed distributions are better reflected at extreme tails rather than at the median. We propose a novel method to estimate the treatment effect at extreme tails incorporating heavy-tailed features in the outcome distribution. The analysis of London transport data using the proposed method indicates that the extreme traffic flow increased substantially after Cycle Superhighways came into operation.
△ Less
Submitted 7 August, 2023; v1 submitted 23 August, 2021;
originally announced August 2021.
-
Quadrupole moments and proton-neutron structure in p-shell mirror nuclei
Authors:
M. A. Caprio,
P. J. Fasano,
P. Maris,
A. E. McCoy
Abstract:
Electric quadrupole (E2) matrix elements provide a measure of nuclear deformation and related collective structure. Ground-state quadrupole moments in particular are known to high precision in many p-shell nuclei. While the experimental electric quadrupole moment only measures the proton distribution, both proton and neutron quadrupole moments are needed to probe proton-neutron asymmetry in the nu…
▽ More
Electric quadrupole (E2) matrix elements provide a measure of nuclear deformation and related collective structure. Ground-state quadrupole moments in particular are known to high precision in many p-shell nuclei. While the experimental electric quadrupole moment only measures the proton distribution, both proton and neutron quadrupole moments are needed to probe proton-neutron asymmetry in the nuclear deformation. We seek insight into the relation between these moments through the ab initio no-core configuration interaction (NCCI), or no-core shell model (NCSM), approach. Converged ab initio calculations for quadrupole moments are particularly challenging, due to sensitivity to long-range behavior of the wave functions. We therefore study more robustly-converged ratios of quadrupole moments: across mirror nuclides, or of proton and neutron quadrupole moments within the same nuclide. In calculations for mirror pairs in the p shell, we explore how well the predictions for mirror quadrupole moments agree with experiment and how well isospin (mirror) symmetry holds for quadrupole moments across a mirror pair.
△ Less
Submitted 18 September, 2021; v1 submitted 22 June, 2021;
originally announced June 2021.
-
Assessing the validity of Bayesian inference using loss functions
Authors:
Yu Luo,
David A. Stephens,
Daniel J. Graham,
Emma J. McCoy
Abstract:
In the usual Bayesian setting, a full probabilistic model is required to link the data and parameters, and the form of this model and the inference and prediction mechanisms are specified via de Finetti's representation. In general, such a formulation is not robust to model mis-specification of its component parts. An alternative approach is to draw inference based on loss functions, where the qua…
▽ More
In the usual Bayesian setting, a full probabilistic model is required to link the data and parameters, and the form of this model and the inference and prediction mechanisms are specified via de Finetti's representation. In general, such a formulation is not robust to model mis-specification of its component parts. An alternative approach is to draw inference based on loss functions, where the quantity of interest is defined as a minimizer of some expected loss, and to construct posterior distributions based on the loss-based formulation; this strategy underpins the construction of the Gibbs posterior. We develop a Bayesian non-parametric approach; specifically, we generalize the Bayesian bootstrap, and specify a Dirichlet process model for the distribution of the observables. We implement this using direct prior-to-posterior calculations, but also using predictive sampling. We also study the assessment of posterior validity for non-standard Bayesian calculations, and provide an efficient way to calibrate the scaling parameter in the Gibbs posterior so that it can achieve the desired coverage rate. We show that the developed non-standard Bayesian updating procedures yield valid posterior distributions in terms of consistency and asymptotic normality under model mis-specification. Simulation studies show that the proposed methods can recover the true value of the parameter efficiently and achieve frequentist coverage even when the sample size is small. Finally, we apply our methods to evaluate the causal impact of speed cameras on traffic collisions in England.
△ Less
Submitted 9 February, 2023; v1 submitted 6 March, 2021;
originally announced March 2021.
-
Rotational bands beyond the Elliott model
Authors:
Ryan Zbikowski,
Calvin W. Johnson,
Anna E. McCoy,
Mark A. Caprio,
Patrick J. Fasano
Abstract:
Rotational bands are commonplace in the spectra of atomic nuclei. Inspired by early descriptions of these bands by quadrupole deformations of a liquid drop, Elliott constructed a discrete nucleon representations of $\mathrm{SU}(3)$ from fermionic creation and annihilation operators. Ever since, Elliott's model has been foundational to descriptions of rotation in nuclei. Later work, however, sugges…
▽ More
Rotational bands are commonplace in the spectra of atomic nuclei. Inspired by early descriptions of these bands by quadrupole deformations of a liquid drop, Elliott constructed a discrete nucleon representations of $\mathrm{SU}(3)$ from fermionic creation and annihilation operators. Ever since, Elliott's model has been foundational to descriptions of rotation in nuclei. Later work, however, suggested the symplectic extension $\mathrm{Sp}(3,R)$ provides a more unified picture. We decompose no-core shell-model nuclear wave functions into symmetry-defined subspaces for several beryllium isotopes, as well as $^{20}$Ne, using the quadratic Casimirs of both Elliott's $\mathrm{SU}(3)$ and $\mathrm{Sp}(3,R)$. The band structure, delineated by strong $B(E2)$ values, has a more consistent description in $\mathrm{Sp}(3,R)$ rather than $\mathrm{SU}(3)$. {In particular, we confirm previous work finding in some nuclides strongly connected upper and lower bands with the same underlying symplectic structure.
△ Less
Submitted 16 November, 2020;
originally announced November 2020.
-
Emergent Sp(3,R) dynamical symmetry in the nuclear many-body system from an ab initio description
Authors:
A. E. McCoy,
M. A. Caprio,
T. Dytrych,
P. J. Fasano
Abstract:
Ab initio nuclear theory provides not only a microscopic framework for quantitative description of the nuclear many-body system, but also a foundation for deeper understanding of emergent collective correlations. A symplectic Sp(3,R)$\supset$U(3) dynamical symmetry is identified in ab initio predictions, from a no-core configuration interaction approach, and found to provide a qualitative understa…
▽ More
Ab initio nuclear theory provides not only a microscopic framework for quantitative description of the nuclear many-body system, but also a foundation for deeper understanding of emergent collective correlations. A symplectic Sp(3,R)$\supset$U(3) dynamical symmetry is identified in ab initio predictions, from a no-core configuration interaction approach, and found to provide a qualitative understanding of the spectrum of 7Be. Low-lying states form an Elliott SU(3) spectrum, while an Sp(3,R) excitation gives rise to an excited rotational band with strong quadrupole connections to the ground state band.
△ Less
Submitted 12 August, 2020;
originally announced August 2020.
-
An Evaluation Framework for Personalization Strategy Experiment Designs
Authors:
C. H. Bryan Liu,
Emma J. McCoy
Abstract:
Online Controlled Experiments (OCEs) are the gold standard in evaluating the effectiveness of changes to websites. An important type of OCE evaluates different personalization strategies, which present challenges in low test power and lack of full control in group assignment. We argue that getting the right experiment setup -- the allocation of users to treatment/analysis groups -- should take pre…
▽ More
Online Controlled Experiments (OCEs) are the gold standard in evaluating the effectiveness of changes to websites. An important type of OCE evaluates different personalization strategies, which present challenges in low test power and lack of full control in group assignment. We argue that getting the right experiment setup -- the allocation of users to treatment/analysis groups -- should take precedence of post-hoc variance reduction techniques in order to enable the scaling of the number of experiments. We present an evaluation framework that, along with a few simple rule of thumbs, allow experimenters to quickly compare which experiment setup will lead to the highest probability of detecting a treatment effect under their particular circumstance.
△ Less
Submitted 9 May, 2023; v1 submitted 22 July, 2020;
originally announced July 2020.
-
BART-based inference for Poisson processes
Authors:
Stamatina Lamprinakou,
Mauricio Barahona,
Seth Flaxman,
Sarah Filippi,
Axel Gandy,
Emma McCoy
Abstract:
The effectiveness of Bayesian Additive Regression Trees (BART) has been demonstrated in a variety of contexts including non-parametric regression and classification. A BART scheme for estimating the intensity of inhomogeneous Poisson processes is introduced. Poisson intensity estimation is a vital task in various applications including medical imaging, astrophysics and network traffic analysis. Th…
▽ More
The effectiveness of Bayesian Additive Regression Trees (BART) has been demonstrated in a variety of contexts including non-parametric regression and classification. A BART scheme for estimating the intensity of inhomogeneous Poisson processes is introduced. Poisson intensity estimation is a vital task in various applications including medical imaging, astrophysics and network traffic analysis. The new approach enables full posterior inference of the intensity in a non-parametric regression setting. The performance of the novel scheme is demonstrated through simulation studies on synthetic and real datasets up to five dimensions, and the new scheme is compared with alternative approaches.
△ Less
Submitted 12 November, 2022; v1 submitted 16 May, 2020;
originally announced May 2020.
-
Intrinsic operators for the translationally-invariant many-body problem
Authors:
M. A. Caprio,
A. E. McCoy,
P. J. Fasano
Abstract:
The need to enforce fermionic antisymmetry in the nuclear many-body problem commonly requires use of single-particle coordinates, defined relative to some fixed origin. To obtain physical operators which nonetheless act on the nuclear many-body system in a Galilean-invariant fashion, thereby avoiding spurious center-of-mass contributions to observables, it is necessary to express these operators w…
▽ More
The need to enforce fermionic antisymmetry in the nuclear many-body problem commonly requires use of single-particle coordinates, defined relative to some fixed origin. To obtain physical operators which nonetheless act on the nuclear many-body system in a Galilean-invariant fashion, thereby avoiding spurious center-of-mass contributions to observables, it is necessary to express these operators with respect to the translational intrinsic frame. Several commonly-encountered operators in nuclear many-body calculations, including the magnetic dipole and electric quadrupole operators (in the impulse approximation) and generators of SU(3) and Sp(3,R) symmetry groups, are bilinear in the coordinates and momenta of the nucleons and, when expressed in intrinsic form, become two-body operators. To work with such operators in a second-quantized many-body calculation, it is necessary to relate three distinct forms: the defining intrinsic-frame expression, an explicitly two-body expression in terms of two-particle relative coordinates, and a decomposition into one-body and separable two-body parts. We establish the relations between these forms, for general (non-scalar and non-isoscalar) operators bilinear in coordinates and momenta.
△ Less
Submitted 18 August, 2020; v1 submitted 24 April, 2020;
originally announced April 2020.
-
Analysing the causal effect of London cycle superhighways on traffic congestion
Authors:
Prajamitra Bhuyan,
Emma J. McCoy,
Haojie Li,
Daniel J. Graham
Abstract:
Transport operators have a range of intervention options available to improve or enhance their networks. Such interventions are often made in the absence of sound evidence on resulting outcomes. Cycling superhighways were promoted as a sustainable and healthy travel mode, one of the aims of which was to reduce traffic congestion. Estimating the impacts that cycle superhighways have on congestion i…
▽ More
Transport operators have a range of intervention options available to improve or enhance their networks. Such interventions are often made in the absence of sound evidence on resulting outcomes. Cycling superhighways were promoted as a sustainable and healthy travel mode, one of the aims of which was to reduce traffic congestion. Estimating the impacts that cycle superhighways have on congestion is complicated due to the non-random assignment of such intervention over the transport network. In this paper, we analyse the causal effect of cycle superhighways utilising pre-intervention and post-intervention information on traffic and road characteristics along with socio-economic factors. We propose a modeling framework based on the propensity score and outcome regression model. The method is also extended to the doubly robust set-up. Simulation results show the superiority of the performance of the proposed method over existing competitors. The method is applied to analyse a real dataset on the London transport network. The methodology proposed can assist in effective decision making to improve network performance.
△ Less
Submitted 23 February, 2021; v1 submitted 19 March, 2020;
originally announced March 2020.
-
Ab initio rotation in 10Be
Authors:
M. A. Caprio,
P. J. Fasano,
A. E. McCoy,
P. Maris,
J. P. Vary
Abstract:
Ab initio theory describes nuclei from a fully microscopic formulation, with no presupposition of collective degrees of freedom, yet signatures of clustering and rotation nonetheless arise. We can therefore look to ab initio theory for an understanding of the nature of these emergent phenomena. To probe the nature of rotation in 10Be, we examine the predicted rotational spectroscopy from no-core c…
▽ More
Ab initio theory describes nuclei from a fully microscopic formulation, with no presupposition of collective degrees of freedom, yet signatures of clustering and rotation nonetheless arise. We can therefore look to ab initio theory for an understanding of the nature of these emergent phenomena. To probe the nature of rotation in 10Be, we examine the predicted rotational spectroscopy from no-core configuration interaction (NCCI) calculations with the Daejeon16 internucleon interaction, and find spectra suggestive of coexisting rotational structures having qualitatively different intrinsic deformations: one triaxial and the other with large axial deformation arising primarily from the neutrons.
△ Less
Submitted 12 December, 2019;
originally announced December 2019.
-
Probing ab initio emergence of nuclear rotation
Authors:
M. A. Caprio,
P. J. Fasano,
P. Maris,
A. E. McCoy,
J. P. Vary
Abstract:
Structural phenomena in nuclei, from shell structure and clustering to superfluidity and collective rotations and vibrations, reflect emergent degrees of freedom. Ab initio theory describes nuclei directly from a fully microscopic formulation. We can therefore look to ab initio theory as a means of exploring the emergence of effective degrees of freedom in nuclei. For the illustrative case of emer…
▽ More
Structural phenomena in nuclei, from shell structure and clustering to superfluidity and collective rotations and vibrations, reflect emergent degrees of freedom. Ab initio theory describes nuclei directly from a fully microscopic formulation. We can therefore look to ab initio theory as a means of exploring the emergence of effective degrees of freedom in nuclei. For the illustrative case of emergent rotational bands in the Be isotopes, we establish an understanding of the underlying oscillator space and angular momentum (orbital and spin) structure. We consider no-core configuration interaction (NCCI) calculations for 7,9,11Be with the Daejeon16 internucleon interaction. Although shell model or rotational degrees of freedom are not assumed in the ab initio theory, the NCCI results are suggestive of the emergence of effective shell model degrees of freedom (0 hbar-omega and 2 hbar-omega excitations) and LS-scheme rotational degrees of freedom, consistent with an Elliott-Wilsdon SU(3) description. These results provide some basic insight into the connection between emergent effective collective rotational and shell model degrees of freedom in these light nuclei and the underlying ab initio microscopic description.
△ Less
Submitted 29 November, 2019;
originally announced December 2019.
-
Symplectic no-core configuration interaction framework for ab initio nuclear structure
Authors:
Anna E. McCoy,
Mark A. Caprio,
Tomas Dytrych
Abstract:
We introduce a symplectic no-core configuration interaction (SpNCCI) framework for ab initio nuclear structure calculations, in a correlated many-body basis which encodes an approximate Sp(3,R) symmetry of the nucleus. Such a scheme potentially provides a means of restricting the many-body space to include only those highly-excited configurations which dominantly contribute to the nuclear wave fun…
▽ More
We introduce a symplectic no-core configuration interaction (SpNCCI) framework for ab initio nuclear structure calculations, in a correlated many-body basis which encodes an approximate Sp(3,R) symmetry of the nucleus. Such a scheme potentially provides a means of restricting the many-body space to include only those highly-excited configurations which dominantly contribute to the nuclear wave function. We examine the symplectic symmetry structure arising in an illustrative ab initio SpNCCI calculation for 6Li. We observe both the dominance of symplectic symmetry in individual wave functions and the emergence of families of states related by symplectic symmetry.
△ Less
Submitted 5 February, 2018;
originally announced February 2018.
-
Quantifying the causal effect of speed cameras on road traffic accidents via an approximate Bayesian doubly robust estimator
Authors:
Daniel J Graham,
Cian Naik,
Emma J McCoy,
Haojie Li
Abstract:
This paper quantifies the effect of speed cameras on road traffic collisions using an approximate Bayesian doubly-robust (DR) causal inference estimation method. Previous empirical work on this topic, which shows a diverse range of estimated effects, is based largely on outcome regression (OR) models using the Empirical Bayes approach or on simple before and after comparisons. Issues of causality…
▽ More
This paper quantifies the effect of speed cameras on road traffic collisions using an approximate Bayesian doubly-robust (DR) causal inference estimation method. Previous empirical work on this topic, which shows a diverse range of estimated effects, is based largely on outcome regression (OR) models using the Empirical Bayes approach or on simple before and after comparisons. Issues of causality and confounding have received little formal attention. A causal DR approach combines propensity score (PS) and OR models to give an average treatment effect (ATE) estimator that is consistent and asymptotically normal under correct specification of either of the two component models. We develop this approach within a novel approximate Bayesian framework to derive posterior predictive distributions for the ATE of speed cameras on road traffic collisions. Our results for England indicate significant reductions in the number of collisions at speed cameras sites (mean ATE = -15%). Our proposed method offers a promising approach for evaluation of transport safety interventions.
△ Less
Submitted 16 August, 2019; v1 submitted 17 March, 2017;
originally announced March 2017.
-
The many-nucleon theory of nuclear collective structure and its macroscopic limits: an algebraic perspective
Authors:
D. J. Rowe,
A. E. McCoy,
M. A. Caprio
Abstract:
The nuclear collective models introduced by Bohr, Mottelson and Rainwater, together with the Mayer-Jensen shell model, have provided the central framework for the development of nuclear physics. This paper reviews the microscopic evolution of the collective models and their underlying foundations. In particular, it is shown that the Bohr-Mottelson models have expressions as macroscopic limits of m…
▽ More
The nuclear collective models introduced by Bohr, Mottelson and Rainwater, together with the Mayer-Jensen shell model, have provided the central framework for the development of nuclear physics. This paper reviews the microscopic evolution of the collective models and their underlying foundations. In particular, it is shown that the Bohr-Mottelson models have expressions as macroscopic limits of microscopic models that have precisely-defined expressions in many-nucleon quantum mechanics. Understanding collective models in this way is especially useful because it enables the analysis of nuclear properties in terms of them to be revisited and reassessed in the light of their microscopic foundations.
△ Less
Submitted 14 March, 2017;
originally announced March 2017.
-
Multiply robust dose-response estimation for multivalued causal inference problems
Authors:
Cian Naik,
Emma J. McCoy,
Daniel J. Graham
Abstract:
This paper develops a multiply robust (MR) dose-response estimator for causal inference problems involving multivalued treatments. We combine a family of generalised propensity score (GPS) models and a family of outcome regression (OR) models to achieve an average potential outcomes estimator that is consistent if just one of the GPS or OR models in each family is correctly specified. We provide p…
▽ More
This paper develops a multiply robust (MR) dose-response estimator for causal inference problems involving multivalued treatments. We combine a family of generalised propensity score (GPS) models and a family of outcome regression (OR) models to achieve an average potential outcomes estimator that is consistent if just one of the GPS or OR models in each family is correctly specified. We provide proofs and simulations that demonstrate multiple robustness in the context of multivalued causal inference problems.
△ Less
Submitted 17 May, 2017; v1 submitted 8 November, 2016;
originally announced November 2016.
-
Algebraic evaluation of matrix elements in the Laguerre function basis
Authors:
A. E. McCoy,
M. A. Caprio
Abstract:
The Laguerre functions constitute one of the fundamental basis sets for calculations in atomic and molecular electron-structure theory, with applications in hadronic and nuclear theory as well. While similar in form to the Coulomb bound-state eigenfunctions (from the Schroedinger eigenproblem) or the Coulomb-Sturmian functions (from a related Sturm-Liouville problem), the Laguerre functions, unlik…
▽ More
The Laguerre functions constitute one of the fundamental basis sets for calculations in atomic and molecular electron-structure theory, with applications in hadronic and nuclear theory as well. While similar in form to the Coulomb bound-state eigenfunctions (from the Schroedinger eigenproblem) or the Coulomb-Sturmian functions (from a related Sturm-Liouville problem), the Laguerre functions, unlike these former functions, constitute a complete, discrete, orthonormal set for square-integrable functions in three dimensions. We construct the SU(1,1)xSO(3) dynamical algebra for the Laguerre functions and apply the ideas of factorization (or supersymmetric quantum mechanics) to derive shift operators for these functions. We use the resulting algebraic framework to derive analytic expressions for matrix elements of several basic radial operators (involving powers of the radial coordinate and radial derivative) in the Laguerre function basis. We illustrate how matrix elements for more general spherical tensor operators in three dimensional space, such as the gradient, may then be constructed from these radial matrix elements.
△ Less
Submitted 16 May, 2016;
originally announced May 2016.
-
Doubly robust dose-response estimation for continuous treatments via generalized propensity score augmented outcome regression
Authors:
Daniel J. Graham,
Emma J. McCoy,
David A. Stephens
Abstract:
This paper constructs a doubly robust estimator for continuous dose-response estimation. An outcome regression model is augmented with a set of inverse generalized propensity score covariates to correct for potential misspecification bias. From the augmented model we can obtain consistent estimates of mean average potential outcomes for distinct strata of the treatment. A polynomial regression is…
▽ More
This paper constructs a doubly robust estimator for continuous dose-response estimation. An outcome regression model is augmented with a set of inverse generalized propensity score covariates to correct for potential misspecification bias. From the augmented model we can obtain consistent estimates of mean average potential outcomes for distinct strata of the treatment. A polynomial regression is then fitted to these point estimates to derive a Taylor approximation to the continuous dose-response function. The bootstrap is used for variance estimation. Analytical results and simulations show that our approach can provide a good approximation to linear or nonlinear dose-response functions under various sources of misspecification of the outcome regression or propensity score models. Efficiency in finite samples is good relative to minimum variance consistent estimators.
△ Less
Submitted 16 June, 2015;
originally announced June 2015.
-
Approximate Bayesian Computation for Smoothing
Authors:
James S. Martin,
Ajay Jasra,
Sumeetpal S. Singh,
Nick Whiteley,
Emma McCoy
Abstract:
We consider a method for approximate inference in hidden Markov models (HMMs). The method circumvents the need to evaluate conditional densities of observations given the hidden states. It may be considered an instance of Approximate Bayesian Computation (ABC) and it involves the introduction of auxiliary variables valued in the same space as the observations. The quality of the approximation may…
▽ More
We consider a method for approximate inference in hidden Markov models (HMMs). The method circumvents the need to evaluate conditional densities of observations given the hidden states. It may be considered an instance of Approximate Bayesian Computation (ABC) and it involves the introduction of auxiliary variables valued in the same space as the observations. The quality of the approximation may be controlled to arbitrary precision through a parameter ε>0 . We provide theoretical results which quantify, in terms of ε, the ABC error in approximation of expectations of additive functionals with respect to the smoothing distributions. Under regularity assumptions, this error is O(nε), where n is the number of time steps over which smoothing is performed. For numerical implementation we adopt the forward-only sequential Monte Carlo (SMC) scheme of [16] and quantify the combined error from the ABC and SMC approximations. This forms some of the first quantitative results for ABC methods which jointly treat the ABC and simulation errors, with a finite number of data and simulated samples. When the HMM has unknown static parameters, we consider particle Markov chain Monte Carlo [2] (PMCMC) methods for batch statistical inference.
△ Less
Submitted 22 June, 2012;
originally announced June 2012.
-
Inference for a Class of Partially Observed Point Process Models
Authors:
James S. Martin,
Ajay Jasra,
Emma McCoy
Abstract:
This paper presents a simulation-based framework for sequential inference from partially and discretely observed point process (PP's) models with static parameters. Taking on a Bayesian perspective for the static parameters, we build upon sequential Monte Carlo (SMC) methods, investigating the problems of performing sequential filtering and smoothing in complex examples, where current methods ofte…
▽ More
This paper presents a simulation-based framework for sequential inference from partially and discretely observed point process (PP's) models with static parameters. Taking on a Bayesian perspective for the static parameters, we build upon sequential Monte Carlo (SMC) methods, investigating the problems of performing sequential filtering and smoothing in complex examples, where current methods often fail. We consider various approaches for approximating posterior distributions using SMC. Our approaches, with some theoretical discussion are illustrated on a doubly stochastic point process applied in the context of finance.
△ Less
Submitted 22 January, 2012;
originally announced January 2012.
-
Gamma Ray Bursts, The Principle of Relative Locality and Connection Normal Coordinates
Authors:
A. E. McCoy
Abstract:
The launch of the Fermi telescope in 2008 opened up the possibility of measuring the energy dependence of the speed of light by considering the time delay in the arrival of gamma ray bursts emitted simultaneously from very distant sources.The expected time delay between the arrival of gamma rays of significantly different energies as predicted by the framework of relative locality has already been…
▽ More
The launch of the Fermi telescope in 2008 opened up the possibility of measuring the energy dependence of the speed of light by considering the time delay in the arrival of gamma ray bursts emitted simultaneously from very distant sources.The expected time delay between the arrival of gamma rays of significantly different energies as predicted by the framework of relative locality has already been calculated in Riemann normal coordinates. In the following, we calculate the time delay in more generality and then specialize to the connection normal coordinate system as a check that the results are coordinate independent. We also show that this result does not depend on the presence of torsion.
△ Less
Submitted 3 January, 2012;
originally announced January 2012.
-
Unitary representations of the Galilean line group: Quantum mechanical principle of equivalence
Authors:
B. R. MacGregor,
A. E. McCoy,
S. Wickramasekara
Abstract:
We present a formalism of Galilean quantum mechanics in non-inertial reference frames and discuss its implications for the equivalence principle. This extension of quantum mechanics rests on the Galilean line group, the semidirect product of the real line and the group of analytic functions from the real line to the Euclidean group in three dimensions. This group provides transformations between a…
▽ More
We present a formalism of Galilean quantum mechanics in non-inertial reference frames and discuss its implications for the equivalence principle. This extension of quantum mechanics rests on the Galilean line group, the semidirect product of the real line and the group of analytic functions from the real line to the Euclidean group in three dimensions. This group provides transformations between all inertial and non-inertial reference frames and contains the Galilei group as a subgroup. We construct a certain class of unitary representations of the Galilean line group and show that these representations determine the structure of quantum mechanics in non-inertial reference frames. Our representations of the Galilean line group contain the usual unitary projective representations of the Galilei group, but have a more intricate cocycle structure. The transformation formula for the Hamiltonian under the Galilean line group shows that in a non-inertial reference frame it acquires a fictitious potential energy term that is proportional to the inertial mass, suggesting the equivalence of inertial mass and gravitational mass in quantum mechanics.
△ Less
Submitted 12 July, 2011;
originally announced July 2011.
-
Racah's method for general subalgebra chains: Coupling coefficients of SO(5) in canonical and physical bases
Authors:
M. A. Caprio,
K. D. Sviratcheva,
A. E. McCoy
Abstract:
It is shown that the method of infinitesimal generators ("Racah's method") can be broadly and systematically formulated as a method applicable to the calculation of reduced coupling coefficients for a generic subalgebra chain G>H, provided the reduced matrix elements of the generators of G and the recoupling coefficients of H are known. The calculation of SO(5)>SO(4) reduced coupling coefficients…
▽ More
It is shown that the method of infinitesimal generators ("Racah's method") can be broadly and systematically formulated as a method applicable to the calculation of reduced coupling coefficients for a generic subalgebra chain G>H, provided the reduced matrix elements of the generators of G and the recoupling coefficients of H are known. The calculation of SO(5)>SO(4) reduced coupling coefficients is considered as an example, and a procedure for transformation of reduced coupling coefficients between canonical and physical subalegebra chains is presented. The problem of calculating coupling coefficients for generic irreps of SO(5), reduced with respect to any of its subalgebra chains, is completely resolved by this approach.
△ Less
Submitted 14 June, 2010;
originally announced June 2010.
-
Non-Regular Likelihood Inference for Seasonally Persistent Processes
Authors:
Emma J. McCoy,
Sofia C. Olhede,
David A. Stephens
Abstract:
The estimation of parameters in the frequency spectrum of a seasonally persistent stationary stochastic process is addressed. For seasonal persistence associated with a pole in the spectrum located away from frequency zero, a new Whittle-type likelihood is developed that explicitly acknowledges the location of the pole. This Whittle likelihood is a large sample approximation to the distribution…
▽ More
The estimation of parameters in the frequency spectrum of a seasonally persistent stationary stochastic process is addressed. For seasonal persistence associated with a pole in the spectrum located away from frequency zero, a new Whittle-type likelihood is developed that explicitly acknowledges the location of the pole. This Whittle likelihood is a large sample approximation to the distribution of the periodogram over a chosen grid of frequencies, and constitutes an approximation to the time-domain likelihood of the data, via the linear transformation of an inverse discrete Fourier transform combined with a demodulation. The new likelihood is straightforward to compute, and as will be demonstrated has good, yet non-standard, properties. The asymptotic behaviour of the proposed likelihood estimators is studied; in particular, $N$-consistency of the estimator of the spectral pole location is established. Large finite sample and asymptotic distributions of the score and observed Fisher information are given, and the corresponding distributions of the maximum likelihood estimators are deduced. A study of the small sample properties of the likelihood approximation is provided, and its superior performance to previously suggested methods is shown, as well as agreement with the developed distributional approximations.
△ Less
Submitted 2 September, 2007;
originally announced September 2007.