skip to content
 

Optimal dimension-reduced calibration and the terminal case for spatial models

Presented by: 
Daniel Williamson University of Exeter
Date: 
Wednesday 7th March 2018 - 11:45 to 12:30
Venue: 
INI Seminar Room 1
Abstract: 
Since the seminal paper by Kennedy and O’Hagan in 2001, the calibration of computer models using Gaussian process emulators has represented a gold standard for scientists and statisticians working to quantify uncertainty using complex computer codes. When the output of such codes is high dimensional, such as with the spatial fields routinely produced by climate models, the standard approach (attributed to Higdon in 2008) is to take principal components across the model output, and use the loadings on these as a lower dimensional representation of the model output that can be used within the Kennedy and O’Hagan framework. In this talk I will argue that, in general, we should not expect this to work. I will introduce what we term a “terminal case analysis” for general computer model calibration, show the implications for inference of a terminal case analysis and argue that though a high dimensional computer model may not lead to a terminal case analysis, the standard statistical treatment outlined above invariably leads to one artificially. I will then present our solution to this which uses rotation ideas to fix the search directions of our lower dimensional representations so that general calibration of spatio-temporal models is possible. We apply our method to idealised examples and to the output of the state of the art Canadian atmosphere model CanAGCM4. We will see that the problem of calibrating climate models requires a great deal of novel statistical thinking before we, as a community, can claim to have a solution ready for this important application area. This is work done with (and by) Dr James Salter.
The video for this talk should appear here if JavaScript is enabled.
If it doesn't, something may have gone wrong with our embedded player.
We'll get it fixed as soon as possible.
University of Cambridge Research Councils UK
    Clay Mathematics Institute London Mathematical Society NM Rothschild and Sons