Mixture of Species Sampling Models

We introduce mixtures of species sampling sequences (mSSS) and discuss how these sequences are related to various types of Bayesian models. As a particular case, we recover species sampling sequences with general (not necessarily diffuse) base measures. These models include some “spike-and-slab” non-parametric priors recently introduced to provide sparsity. Furthermore, we show how mSSS arise while considering hierarchical species sampling random probabilities (e.g., the hierarchical Dirichlet process). Extending previous results, we prove that mSSS are obtained by assigning the values of an exchangeable sequence to the classes of a latent exchangeable random partition. Using this representation, we give an explicit expression of the Exchangeable Partition Probability Function of the partition generated by an mSSS. Some special cases are discussed in detail—in particular, species sampling sequences with general base measures and a mixture of species sampling sequences with Gibbs-type latent partition. Finally, we give explicit expressions of the predictive distributions of an mSSS.

Propagation of epistemic uncertainty in magnitude–frequency relations through PSHA using predictive distributions

The current practice of Probabilistic Seismic Hazard Analysis (PSHA) considers the inclusion of epistemic uncertainties involved in different parts of the analysis via the logic-tree approach. Given the complexity of modern PSHA models, numerous branches are needed, which in some cases leads to concerns regarding performance issues. We introduce the use of a magnitude exceedance rate which, following Bayesian conventions, we call predictive exceedance rate. This rate is the original Gutenberg–Richter relation after having included the effect of the epistemic uncertainty in parameter β. The predictive exceedance rate was first proposed by Campbell but to our best knowledge is seldom used in current PSHA. We show that the predictive exceedance rate is as accurate as the typical logic-tree approach but allows for much faster computations, a very useful property given the complexity of some modern PSHA models.

Predictive Constructions Based on Measure-Valued Pólya Urn Processes

Measure-valued Pólya urn processes (MVPP) are Markov chains with an additive structure that serve as an extension of the generalized k-color Pólya urn model towards a continuum of possible colors. We prove that, for any MVPP (μn)n≥0 on a Polish space X, the normalized sequence (μn/μn(X))n≥0 agrees with the marginal predictive distributions of some random process (Xn)n≥1. Moreover, μn=μn−1+RXn, n≥1, where x↦Rx is a random transition kernel on X; thus, if μn−1 represents the contents of an urn, then Xn denotes the color of the ball drawn with distribution μn−1/μn−1(X) and RXn—the subsequent reinforcement. In the case RXn=WnδXn, for some non-negative random weights W1,W2,…, the process (Xn)n≥1 is better understood as a randomly reinforced extension of Blackwell and MacQueen’s Pólya sequence. We study the asymptotic properties of the predictive distributions and the empirical frequencies of (Xn)n≥1 under different assumptions on the weights. We also investigate a generalization of the above models via a randomization of the law of the reinforcement.

Evaluating epidemic forecasts in an interval format

For practical reasons, many forecasts of case, hospitalization, and death counts in the context of the current Coronavirus Disease 2019 (COVID-19) pandemic are issued in the form of central predictive intervals at various levels. This is also the case for the forecasts collected in the COVID-19 Forecast Hub (https://covid19forecasthub.org/). Forecast evaluation metrics like the logarithmic score, which has been applied in several infectious disease forecasting challenges, are then not available as they require full predictive distributions. This article provides an overview of how established methods for the evaluation of quantile and interval forecasts can be applied to epidemic forecasts in this format. Specifically, we discuss the computation and interpretation of the weighted interval score, which is a proper score that approximates the continuous ranked probability score. It can be interpreted as a generalization of the absolute error to probabilistic forecasts and allows for a decomposition into a measure of sharpness and penalties for over- and underprediction.