Confidence sets for asset correlations in portfolio credit risk

Asset correlations are of critical importance in quantifying portfolio credit risk and economic capitalin financial institutions. Estimation of asset correlation with rating transition data has focusedon the point estimation of the correlation without giving any consideration to the uncertaintyaround these point estimates. In this article we use Bayesian methods to estimate a dynamicfactor model for default risk using rating data (McNeil et al., 2005; McNeil and Wendin, 2007).Bayesian methods allow us to formally incorporate human judgement in the estimation of assetcorrelation, through the prior distribution and fully characterize a confidence set for the correlations.Results indicate: i) a two factor model rather than the one factor model, as proposed bythe Basel II framework, better represents the historical default data. ii) importance of unobservedfactors in this type of models is reinforced and point out that the levels of the implied asset correlationscritically depend on the latent state variable used to capture the dynamics of default,as well as other assumptions on the statistical model. iii) the posterior distributions of the assetcorrelations show that the Basel recommended bounds, for this parameter, undermine the levelof systemic risk.

Stratospheric Memory and Skill of Extended-Range Weather Forecasts

We use an empirical statistical model to demonstrate significant skill in making extended-range forecasts of the monthly-mean Arctic Oscillation (AO). Forecast skill derives from persistent circulation anomalies in the lowermost stratosphere and is greatest during boreal winter. A comparison to the Southern Hemisphere provides evidence that both the time scale and predictability of the AO depend on the presence of persistent circulation anomalies just above the tropopause. These circulation anomalies most likely affect the troposphere through changes to waves in the upper troposphere, which induce surface pressure changes that correspond to the AO.

Teleconnections between Ethiopian summer rainfall and sea surface temperature: part II. Seasonal forecasting

A seasonal forecasting system that is capable of skilfully predicting rainfall totals on a regional scale would be of great value to Ethiopia. Here, we describe how a statistical model can exploit the teleconnections described in part 1 of this pair of papers to develop such a system. We show that, in most cases, the predictors selected objectively by the statistical model can be interpreted in the light of physical teleconnections with Ethiopian rainfall, and discuss why, in some cases, unexpected regions are chosen as predictors. We show that the forecast has skill in all parts of Ethiopia, and argue that this method could provide the basis of an operational seasonal forecasting system for Ethiopia.

Technical note: A proposed method to determine the extent of degradation of a feed in the rumen from the degradation profile obtained with the in vitro gas production technique using feces as the inoculum

A method is proposed to determine the extent of degradation in the rumen involving a two-stage mathematical modeling process. In the first stage, a statistical model shifts (or maps) the gas accumulation profile obtained using a fecal inoculum to a ruminal gas profile. Then, a kinetic model determines the extent of degradation in the rumen from the shifted profile. The kinetic model is presented as a generalized mathematical function, allowing any one of a number of alternative equation forms to be selected. This method might allow the gas production technique to become an approach for determining extent of degradation in the rumen, decreasing the need for surgically modified animals while still maintaining the link with the animal. Further research is needed before the proposed methodology can be used as a standard method across a range of feeds.

Bayesian estimation of ancestral character states on phylogenies

Biologists frequently attempt to infer the character states at ancestral nodes of a phylogeny from the distribution of traits observed in contemporary organisms. Because phylogenies are normally inferences from data, it is desirable to account for the uncertainty in estimates of the tree and its branch lengths when making inferences about ancestral states or other comparative parameters. Here we present a general Bayesian approach for testing comparative hypotheses across statistically justified samples of phylogenies, focusing on the specific issue of reconstructing ancestral states. The method uses Markov chain Monte Carlo techniques for sampling phylogenetic trees and for investigating the parameters of a statistical model of trait evolution. We describe how to combine information about the uncertainty of the phylogeny with uncertainty in the estimate of the ancestral state. Our approach does not constrain the sample of trees only to those that contain the ancestral node or nodes of interest, and we show how to reconstruct ancestral states of uncertain nodes using a most-recent-common-ancestor approach. We illustrate the methods with data on ribonuclease evolution in the Artiodactyla. Software implementing the methods ( BayesMultiState) is available from the authors.

Declining survival rates in a reintroduced population of the Mauritius kestrel: evidence for non-linear density dependence and environmental stochasticity

1. We studied a reintroduced population of the formerly critically endangered Mauritius kestrel Falco punctatus Temmink from its inception in 1987 until 2002, by which time the population had attained carrying capacity for the study area. Post-1994 the population received minimal management other than the provision of nestboxes. 2. We analysed data collected on survival (1987-2002) using program MARK to explore the influence of density-dependent and independent processes on survival over the course of the population's development. 3.We found evidence for non-linear, threshold density dependence in juvenile survival rates. Juvenile survival was also strongly influenced by climate, with the temporal distribution of rainfall during the cyclone season being the most influential climatic variable. Adult survival remained constant throughout. 4. Our most parsimonious capture-mark-recapture statistical model, which was constrained by density and climate, explained 75.4% of the temporal variation exhibited in juvenile survival rates over the course of the population's development. 5. This study is an example of how data collected as part of a threatened species recovery programme can be used to explore the role and functional form of natural population regulatory processes. With the improvements in conservation management techniques and the resulting success stories, formerly threatened species offer unique opportunities to further our understanding of the fundamental principles of population ecology.

Quantifying the residual creep life of polyester mooring ropes

This paper assesses the ARELIS (Assured Residual Life Span) method for estimating residual creep life of polyester rope used in deepwater mooring lines. A statistical model has been developed to quantify the uncertainties in the method, such as the scatter in creep rupture test data and load sharing between sub-ropes. This model can be used to determine the required test load, duration and number of ARELIS tests, in order to guarantee a minimum creep life for a mooring line at its service load. Creep rupture tests have been performed to provide input for the statistical model.

Influence of drought-induced acidification on the mobility of dissolved organic carbon in peat soils

A strong relationship between dissolved organic carbon (DOC) and sulphate (SO42−) dynamics under drought conditions has been revealed from analysis of a 10-year time series (1993–2002). Soil solution from a blanket peat at 10 cm depth and stream water were collected at biweekly and weekly intervals, respectively, by the Environmental Change Network at Moor House-Upper Teesdale National Nature Reserve in the North Pennine uplands of Britain. DOC concentrations in soil solution and stream water were closely coupled, displaying a strong seasonal cycle with lowest concentrations in early spring and highest in late summer/early autumn. Soil solution DOC correlated strongly with seasonal variations in soil temperature at the same depth 4-weeks prior to sampling. Deviation from this relationship was seen, however, in years with significant water table drawdown (>−25 cm), such that DOC concentrations were up to 60% lower than expected. Periods of drought also resulted in the release of SO42−, because of the oxidation of inorganic/organic sulphur stored in the peat, which was accompanied by a decrease in pH and increase in ionic strength. As both pH and ionic strength are known to control the solubility of DOC, inclusion of a function to account for DOC suppression because of drought-induced acidification accounted for more of the variability of DOC in soil solution (R2=0.81) than temperature alone (R2=0.58). This statistical model of peat soil solution DOC at 10 cm depth was extended to reproduce 74% of the variation in stream DOC over this period. Analysis of annual budgets showed that the soil was the main source of SO42− during droughts, while atmospheric deposition was the main source in other years. Mass balance calculations also showed that most of the DOC originated from the peat. The DOC flux was also lower in the drought years of 1994 and 1995, reflecting low DOC concentrations in soil and stream water. The analysis presented in this paper suggests that lower concentrations of DOC in both soil and stream waters during drought years can be explained in terms of drought-induced acidification. As future climate change scenarios suggest an increase in the magnitude and frequency of drought events, these results imply potential for a related increase in DOC suppression by episodic acidification.

Regularized logistic models for probabilistic forecasting and diagnostics

Logistic models are studied as a tool to convert dynamical forecast information (deterministic and ensemble) into probability forecasts. A logistic model is obtained by setting the logarithmic odds ratio equal to a linear combination of the inputs. As with any statistical model, logistic models will suffer from overfitting if the number of inputs is comparable to the number of forecast instances. Computational approaches to avoid overfitting by regularization are discussed, and efficient techniques for model assessment and selection are presented. A logit version of the lasso (originally a linear regression technique), is discussed. In lasso models, less important inputs are identified and the corresponding coefficient is set to zero, providing an efficient and automatic model reduction procedure. For the same reason, lasso models are particularly appealing for diagnostic purposes.

Decline and recovery of total column ozone using a multimodel time series analysis

Simulations of 15 coupled chemistry climate models, for the period 1960–2100, are presented. The models include a detailed stratosphere, as well as including a realistic representation of the tropospheric climate. The simulations assume a consistent set of changing greenhouse gas concentrations, as well as temporally varying chlorofluorocarbon concentrations in accordance with observations for the past and expectations for the future. The ozone results are analyzed using a nonparametric additive statistical model. Comparisons are made with observations for the recent past, and the recovery of ozone, indicated by a return to 1960 and 1980 values, is investigated as a function of latitude. Although chlorine amounts are simulated to return to 1980 values by about 2050, with only weak latitudinal variations, column ozone amounts recover at different rates due to the influence of greenhouse gas changes. In the tropics, simulated peak ozone amounts occur by about 2050 and thereafter total ozone column declines. Consequently, simulated ozone does not recover to values which existed prior to the early 1980s. The results also show a distinct hemispheric asymmetry, with recovery to 1980 values in the Northern Hemisphere extratropics ahead of the chlorine return by about 20 years. In the Southern Hemisphere midlatitudes, ozone is simulated to return to 1980 levels only 10 years ahead of chlorine. In the Antarctic, annually averaged ozone recovers at about the same rate as chlorine in high latitudes and hence does not return to 1960s values until the last decade of the simulations.

Ultraconserved words point to deep language ancestry across Eurasia

The search for ever deeper relationships among the World’s languages is bedeviled by the fact that most words evolve too rapidly to preserve evidence of their ancestry beyond 5,000 to 9,000 y. On the other hand, quantitative modeling indicates that some “ultraconserved” words exist that might be used to find evidence for deep linguistic relationships beyond that time barrier. Here we use a statistical model, which takes into account the frequency with which words are used in common everyday speech, to predict the existence of a set of such highly conserved words among seven language families of Eurasia postulated to form a linguistic superfamily that evolved from a common ancestor around 15,000 y ago. We derive a dated phylogenetic tree of this proposed superfamily with a time-depth of ∼14,450 y, implying that some frequently used words have been retained in related forms since the end of the last ice age. Words used more than once per 1,000 in everyday speech were 7- to 10-times more likely to show deep ancestry on this tree. Our results suggest a remarkable fidelity in the transmission of some words and give theoretical justification to the search for features of language that might be preserved across wide spans of time and geography.

Oblique h.f. radiowave propagation in the main trough region of the ionosphere

The propagation of 7.335 MHz, c.w. signals over a 5212 km sub-auroral, west-east path is studied. Measurements and semi-empirical predictions are made of the amplitude distributions and Doppler shifts of the received signals. The observed amplitude distribution is fitted with one produced by a numerical fading model, yielding the power losses suffered by the signals during propagation via the predominating modes. The signals are found to suffer exceptionally low losses at certain local times under geomagnetically quiet conditions. The mid-latitude trough in the F2 peak ionization density is predicted by a statistical model to be at the latitudes of this path at these times and at low Kp values. A sharp cut-off in low-power losses at a mean Kp of 2.75 strongly implicates the trough in the propagation of these signals. The Doppler shifts observed at these times cannot be explained by a simple ray-tracing model. It is shown however, that a simple extension of this model to allow for the trough can reproduce the form of the observed diurnal variation.

Detecting regular sound changes in linguistics as events of concerted evolution

Background: Concerted evolution is normally used to describe parallel changes at different sites in a genome, but it is also observed in languages where a specific phoneme changes to the same other phoneme in many words in the lexicon—a phenomenon known as regular sound change. We develop a general statistical model that can detect concerted changes in aligned sequence data and apply it to study regular sound changes in the Turkic language family. Results: Linguistic evolution, unlike the genetic substitutional process, is dominated by events of concerted evolutionary change. Our model identified more than 70 historical events of regular sound change that occurred throughout the evolution of the Turkic language family, while simultaneously inferring a dated phylogenetic tree. Including regular sound changes yielded an approximately 4-fold improvement in the characterization of linguistic change over a simpler model of sporadic change, improved phylogenetic inference, and returned more reliable and plausible dates for events on the phylogenies. The historical timings of the concerted changes closely follow a Poisson process model, and the sound transition networks derived from our model mirror linguistic expectations. Conclusions: We demonstrate that a model with no prior knowledge of complex concerted or regular changes can nevertheless infer the historical timings and genealogical placements of events of concerted change from the signals left in contemporary data. Our model can be applied wherever discrete elements—such as genes, words, cultural trends, technologies, or morphological traits—can change in parallel within an organism or other evolving group.

Forecasting elections

In this paper we assess opinion polls, prediction markets, expert opinion and statistical modelling over a large number of US elections in order to determine which perform better in terms of forecasting outcomes. In line with existing literature, we bias-correct opinion polls. We consider accuracy, bias and precision over different time horizons before an election, and we conclude that prediction markets appear to provide the most precise forecasts and are similar in terms of bias to opinion polls. We find that our statistical model struggles to provide competitive forecasts, while expert opinion appears to be of value. Finally we note that the forecast horizon matters; whereas prediction market forecasts tend to improve the nearer an election is, opinion polls appear to perform worse, while expert opinion performs consistently throughout. We thus contribute to the growing literature comparing election forecasts of polls and prediction markets.

Mixture of changing uniaxial micellar forms in lyotropic biaxial nematics

Lyotropic nematics consisting of surfactant-cosurfactant water solutions may present a biaxial phase or direct U(+) <-> U(-) transitions, in different regions of the temperature-relative concentration phase diagram, for different systems and compositions. We propose that these may be related to changes of uniaxial micellar form, which may occur either smoothly or abruptly. Smooth change of cylinder-like into disc-like shapes requires a distribution of Maier-Saupe interaction constants and we consider two limiting cases for the distribution of forms: a single Gaussian and a double Gaussian. Alternatively, an abrupt change of form is described by a discontinuous distribution of interaction constants. Our results show that the dispersive distributions yield a biaxial phase, while an abrupt change of shape leads to coexistence of uniaxial phases. Fitting the theory to the experiment for the ternary system KL/decanol/D2O leads to transition lines in very good agreement with experimental results. In order to rationalise the results of the comparison, we analyse temperature and concentration form dependence, which connects micellar and experimental macroscopic parameters. Physically consistent variations of micellar asymmetry, amphiphile partitioning and volume are obtained. To the best of the authors` knowledge, this is the first truly statistical microscopic approach that is able to model experimentally observed lyotropic biaxial nematic phases.