Galaxy formation : a Bayesian uncertainty analysis.

In many scientific disciplines complex computer models are used to understand the behaviour of large scale physical systems. An uncertainty anal- ysis of such a computer model known as Galform is presented. Galform models the creation and evolution of approximately one million galaxies from the begin- ning of the Universe until the current day, and is regarded as a state-of-the-art model within the cosmology community. It requires the specification of many in- put parameters in order to run the simulation, takes significant time to run, and provides various outputs that can be compared with real world data. A Bayes Linear approach is presented in order to identify the subset of the input space that could give rise to acceptable matches between model output and measured data. This approach takes account of the major sources of uncertainty in a consistent and unified manner, including input parameter uncertainty, function uncertainty, observational error, forcing function uncertainty and structural uncertainty. The approach is known as History Matching, and involves the use of an iterative suc- cession of emulators (stochastic belief specifications detailing beliefs about the Galform function), which are used to cut down the input parameter space. The analysis was successful in producing a large collection of model evaluations that exhibit good fits to the observed data.

[1]  L. M. M.-T. Theory of Probability , 1929, Nature.

[2]  T. J. Mitchell,et al.  Bayesian Prediction of Deterministic Functions, with Applications to the Design and Analysis of Computer Experiments , 1991 .

[3]  Robert Haining,et al.  Statistics for spatial data: by Noel Cressie, 1991, John Wiley & Sons, New York, 900 p., ISBN 0-471-84336-9, US $89.95 , 1993 .

[4]  F. Pukelsheim The Three Sigma Rule , 1994 .

[5]  Adrian E. Raftery,et al.  Inference from a Deterministic Population Dynamics Model for Bowhead Whales , 1995 .

[6]  A. Seheult,et al.  Pressure Matching for Hydrocarbon Reservoirs: A Case Study in the Use of Bayes Linear Strategies for Large Computer Experiments , 1997 .

[7]  R. Ellis,et al.  The 2dF galaxy redshift survey: near-infrared galaxy luminosity functions , 2000, astro-ph/0012429.

[8]  S.Cole,et al.  The 2dF Galaxy Redshift Survey: spectra and redshifts , 2001, astro-ph/0106498.

[9]  A. O'Hagan,et al.  Bayesian calibration of computer models , 2001 .

[10]  D. Madgwick,et al.  The 2dF Galaxy Redshift Survey: The bJ-band galaxy luminosity function and survey selection function , 2001, astro-ph/0111011.

[11]  A. O'Hagan,et al.  Bayesian inference for the uncertainty distribution of computer model outputs , 2002 .

[12]  Edward J. Wollack,et al.  First year Wilkinson Microwave Anisotropy Probe (WMAP) observations: Determination of cosmological parameters , 2003, astro-ph/0302209.

[13]  M. Halpern,et al.  SUBMITTED TO The Astrophysical Journal Preprint typeset using L ATEX style emulateapj v. 11/12/01 FIRST YEAR WILKINSON MICROWAVE ANISOTROPY PROBE (WMAP) OBSERVATIONS: INTERPRETATION OF THE TT AND TE ANGULAR POWER SPECTRUM PEAKS , 2022 .

[14]  Dave Higdon,et al.  Combining Field Data and Computer Simulations for Calibration and Prediction , 2005, SIAM J. Sci. Comput..

[15]  Edinburgh,et al.  Simulating the joint evolution of quasars, galaxies and their large-scale distribution , 2005, astro-ph/0504097.

[16]  J. Peacock,et al.  Simulations of the formation, evolution and clustering of galaxies and quasars , 2005, Nature.

[17]  Michael Goldstein,et al.  Subjective Bayesian Analysis: Principles and Practice , 2006 .

[18]  M. Goldstein Bayes Linear Analysis , 2006 .

[19]  A. OHagan,et al.  Bayesian analysis of computer code outputs: A tutorial , 2006, Reliab. Eng. Syst. Saf..

[20]  C. Baugh,et al.  A primer on hierarchical galaxy formation: the semi-analytical approach , 2006, astro-ph/0610031.

[21]  J. Rougier,et al.  Bayes Linear Calibrated Prediction for Complex Systems , 2006 .

[22]  David Wooff,et al.  Bayes Linear Statistics: Theory and Methods , 2007 .

[23]  D. Higdon,et al.  Computer Model Calibration Using High-Dimensional Output , 2008 .

[24]  J. Rougier Efficient Emulators for Multivariate Deterministic Functions , 2008 .

[25]  David Higdon,et al.  THE COYOTE UNIVERSE. II. COSMOLOGICAL MODELS AND PRECISION EMULATION OF THE NONLINEAR MATTER POWER SPECTRUM , 2009, 0902.0429.

[26]  Michael Goldstein,et al.  Reified Bayesian modelling and inference for physical systems , 2009 .

[27]  Michael Goldstein,et al.  Bayes linear analysis of imprecision in computer models, with application to understanding galaxy formation. , 2009 .

[28]  Michael Goldstein,et al.  Exchangeable computer models. , 2009 .

[29]  Anthony O'Hagan,et al.  Diagnostics for Gaussian Process Emulators , 2009, Technometrics.

[30]  A. O'Hagan,et al.  Gaussian process emulation of dynamic computer codes , 2009 .

[31]  Stanley H. Cohen,et al.  Design and Analysis , 2010 .

[32]  James O. Berger,et al.  External Bayesian analysis for computer simulators , 2010 .

[33]  Sonja Kuhnt,et al.  Design and analysis of computer experiments , 2010 .

[34]  Michael Goldstein,et al.  Bayes linear uncertainty analysis for oil reservoirs based on multiscale computer experiments. , 2010 .

[35]  Mathematics,et al.  The Parameter Space of Galaxy Formation , 2010, 1004.0711.