Statistics for Spatial Data.

We consider prediction and uncertainty analysis for systems which are approximated using complex mathematical models. Such models, implemented as computer codes, are often generic in the sense that by a suitable choice of some of the model's input parameters the code can be used to predict the behaviour of the system in a variety of specific applications. However, in any specific application the values of necessary parameters may be unknown. In this case, physical observations of the system in the specific context are used to learn about the unknown parameters. The process of fitting the model to the observed data by adjusting the parameters is known as calibration. Calibration is typically effected by ad hoc fitting, and after calibration the model is used, with the fitted input values, to predict the future behaviour of the system. We present a Bayesian calibration technique which improves on this traditional approach in two respects. First, the predictions allow for all sources of uncertainty, including the remaining uncertainty over the fitted parameters. Second, they attempt to correct for any inadequacy of the model which is revealed by a discrepancy between the observed data and the model predictions from even the best-fitting parameter values. The method is illustrated by using data from a nuclear radiation release at Tomsk, and from a more complex simulated nuclear accident exercise.

https://www.jstor.org/stable/pdfplus/2680584.pdf

Bayesian Calibration of computer models

The author's preface gives an outline: "This book is about weakconvergence methods in metric spaces, with applications sufficient to show their power and utility. The Introduction motivates the definitions and indicates how the theory will yield solutions to problems arising outside it. Chapter 1 sets out the basic general theorems, which are then specialized in Chapter 2 to the space C[0, l ] of continuous functions on the unit interval and in Chapter 3 to the space D [0, 1 ] of functions with discontinuities of the first kind. The results of the first three chapters are used in Chapter 4 to derive a variety of limit theorems for dependent sequences of random variables. " The book develops and expands on Donsker's 1951 and 1952 papers on the invariance principle and empirical distributions. The basic random variables remain real-valued although, of course, measures on C[0, l ] and D[0, l ] are vitally used. Within this framework, there are various possibilities for a different and apparently better treatment of the material. More of the general theory of weak convergence of probabilities on separable metric spaces would be useful. Metrizability of the convergence is not brought up until late in the Appendix. The close relation of the Prokhorov metric and a metric for convergence in probability is (hence) not mentioned (see V. Strassen, Ann. Math. Statist. 36 (1965), 423-439; the reviewer, ibid. 39 (1968), 1563-1572). This relation would illuminate and organize such results as Theorems 4.1, 4.2 and 4.4 which give isolated, ad hoc connections between weak convergence of measures and nearness in probability. In the middle of p. 16, it should be noted that C*(S) consists of signed measures which need only be finitely additive if 5 is not compact. On p. 239, where the author twice speaks of separable subsets having nonmeasurable cardinal, he means "discrete" rather than "separable." Theorem 1.4 is Ulam's theorem that a Borel probability on a complete separable metric space is tight. Theorem 1 of Appendix 3 weakens completeness to topological completeness. After mentioning that probabilities on the rationals are tight, the author says it is an

Convergence of probability measures

(2002). Modelling Extremal Events for Insurance and Finance. Journal of the American Statistical Association: Vol. 97, No. 457, pp. 360-360.

Modelling Extremal Events for Insurance and Finance

▶ Addresses a wide range of timely environment, economic and energy topics ▶ A forum to review, analyze and stimulate the development, testing and implementation of mitigation and adaptation strategies at regional, national and global scales ▶ Contributes to real-time policy analysis and development as national and international policies and agreements are discussed and promulgated ▶ 94% of authors who answered a survey reported that they would definitely publish or probably publish in the journal again

Mitigation and Adaptation Strategies for Global Change

Statistics for the Environment.

Recently there has been a lot of effort to model extremes of spatially dependent data. These efforts seem to be divided into two distinct groups: the study of max-stable processes, together with the development of statistical models within this framework; the use of more pragmatic, flexible models using Bayesian hierarchical models (BHM) and simulation based inference techniques. Each modeling strategy has its strong and weak points. While max-stable models capture the local behavior of spatial extremes correctly, hierarchical models based on the conditional independence assumption, lack the asymptotic arguments the max-stable models enjoy. On the other hand, they are very flexible in allowing the introduction of physical plausibility into the model. When the objective of the data analysis is to estimate return levels or kriging of extreme values in space, capturing the correct dependence structure between the extremes is crucial and max-stable processes are better suited for these purposes. However when the primary interest is to explain the sources of variation in extreme events Bayesian hierarchical modeling is a very flexible tool due to the ease with which random effects are incorporated in the model. In this paper we model a data set on Portuguese wildfires to show the flexibility of BHM in incorporating spatial dependencies acting at different resolutions.

/pdf/asymptotic-models-and-inference-for-extremes-of-spatio-2tmzpsbg8m.pdf

Asymptotic models and inference for extremes of spatio-temporal data

Spatial and temporal patterns of large fire (>100 ha) incidence in Portugal over the period 1984–2004 were modeled using extreme value statistics, namely the Peaks Over Threshold approach, which uses the Generalized Pareto Distribution (GPD) as a model. The original dataset includes all fires larger than 5 ha (30 616 fires) that were observed in Portugal during the study period, mapped from Landsat satellite imagery. The country was divided into eight regions, considered internally homogeneous from the perspective of their fire regimes and respective environmental correlates. The temporal analysis showed that there does not appear to be any trend in the incidence of very large fires, but revealed a cyclical behavior in the values of the GPD shape parameter, with a period in the range of 3 to 5 years. Spatial analysis highlighted strong regional differences in the incidence of large fires, and allowed the calculation of return levels for a range of fire sizes. This analysis was affected by the presence of a few outlying observations, which may correspond to clusters of contiguous fire scars, resulting in artificially large burned areas. We discuss some of the implications of our findings in terms of consequences for fire management aimed at preventing the occurrence of extremely large fires, and present ideas for extending the present study.

Kamil Feridun Turkman

Papers

Statistics for the Environment.

Asymptotic models and inference for extremes of spatio-temporal data

Spatial and temporal extremes of wildfire sizes in Portugal (1984–2004)

Statistics for the Environment 2: Water Related Issues.

Non-Linear Time Series