An empirical Bayes' procedure for ranking players in Ryder Cup golf

We describe a model to obtain strengths and rankings of players appearing in golf's Ryder Cup. Obtaining rankings is complicated because of two reasons. First, competitors do not compete on an equal number of occasions, with some competitors appearing too infrequently for their ranking to be estimated with any degree of certainty, and second, different competitors experience different levels of volatility in results. Our approach is to assume the competitor strengths are drawn from some common distribution. For small numbers of competitors, as is the case here, we fit the model using Monte-Carlo integration. Results suggest there is very little difference between the top performing players, though Scotland's Colin Montgomerie is estimated as the strongest Ryder Cup player.     

Euro area inflation as a predictor of national inflation rates

The stability of inflation differentials is an important condition for the smooth working of a currency area, such as the European Economic and Monetary Union. In the presence of stability, changes in national inflation rates, while holding Euro area inflation fixed contemporaneously, should be only transitory. If this is the case, the rate of inflation of the whole area can also be interpreted as a predictor, at least in the long-run, of the different national inflation rates. However, in this paper we show that this condition is satisfied only for a small number of countries, including France and Italy. Better convergence results for inflation differentials are, instead, found for the USA. Some policy implications are drawn for the Eurozone.     

Reliability of stress–strength model for exponentiated Pareto distributions

In this paper, the reliability of a system is discussed when the strength of the system and the stress imposed on it are independent, non-identical exponentiated Pareto distributed random variables. Different point estimations and interval estimations are proposed. The point estimators obtained are maximum likelihood, uniformly minimum variance unbiased and Bayesian estimators. The interval estimations obtained are approximate, exact, bootstrap-p and bootstrap-t confidence intervals and Bayesian credible interval. Different methods and the corresponding confidence intervals are compared using Monte-carlo simulations.     

Mixed integer (0–1) fractional programming for decision support in paper production industry

This paper presents an effective and efficient method for solving a special class of mixed integer fractional programming (FP) problems. We take a classical reformulation approach for continuous FP as a starting point and extend it for solving a more general class of mixed integer (0–1) fractional programming problems.To stress the practical relevance of the research we focus on a real-life application in paper production industry. The constantly advancing physical knowledge of large scale pulp and paper production did have a substantial impact on an existing DSS in which mixed integer (0–1) fractional programming is introduced. We show that the motivation to solve a real-life fractional programming problem can provide the basis for a new approach in a new context that has an added value of its own, even outside the given application area. We describe the main characteristics of the DSS, the necessity to develop a non-iterative solution procedure and demonstrate both the effectiveness and efficiency of the proposed approach from practical data sets.     

Symptoms of Eating Disorders Among Females in Drug Addiction Treatment

The aim of this study is to measure and describe symptoms of eating disorders among females in treatment for drug addiction in Norway. Previous clinical and epidemiological studies have revealed coprevalence between eating disorders and substance use or abuse. However, few studies have measured eating disorders in drug-using samples and even fewer within the context of drug treatment. In this study, 29 females with drug use disorder in residential treatment were tested with the Eating Disorder Inventory–2. A subgroup of 9 females (31%) with significant symptoms of eating disorders was identified. The characteristics of this group and possible clinical consequences are discussed.     

LAN property for an ergodic Ornstein–Uhlenbeck process with Poisson jumps

In this article, we consider an ergodic Ornstein–Uhlenbeck process with jumps driven by a Brownian motion and a compensated Poisson process, whose drift and diffusion coefficients as well as its jump intensity depend on unknown parameters. Considering the process discretely observed at high frequency, we derive the local asymptotic normality property. To obtain this result, Malliavin calculus and Girsanov’s theorem are applied to write the log-likelihood ratio in terms of sums of conditional expectations, for which a central limit theorem for triangular arrays can be applied.     

Marginal likelihood estimation from the Metropolis output: tips and tricks for efficient implementation in generalized linear latent variable models

The marginal likelihood can be notoriously difficult to compute, and particularly so in high-dimensional problems. Chib and Jeliazkov employed the local reversibility of the Metropolis–Hastings algorithm to construct an estimator in models where full conditional densities are not available analytically. The estimator is free of distributional assumptions and is directly linked to the simulation algorithm. However, it generally requires a sequence of reduced Markov chain Monte Carlo runs which makes the method computationally demanding especially in cases when the parameter space is large. In this article, we study the implementation of this estimator on latent variable models which embed independence of the responses to the observables given the latent variables (conditional or local independence). This property is employed in the construction of a multi-block Metropolis-within-Gibbs algorithm that allows to compute the estimator in a single run, regardless of the dimensionality of the parameter space. The counterpart one-block algorithm is also considered here, by pointing out the difference between the two approaches. The paper closes with the illustration of the estimator in simulated and real-life data sets.     

Nonparametric Bayesian inference for multivariate density functions using Feller priors

Multivariate density estimation plays an important role in investigating the mechanism of high-dimensional data. This article describes a nonparametric Bayesian approach to the estimation of multivariate densities. A general procedure is proposed for constructing Feller priors for multivariate densities and their theoretical properties as nonparametric priors are established. A blocked Gibbs sampling algorithm is devised to sample from the posterior of the multivariate density. A simulation study is conducted to evaluate the performance of the procedure.     

Bias correction for parameter estimates of spatial point process models

When a spatial point process model is fitted to spatial point pattern data using standard software, the parameter estimates are typically biased. Contrary to folklore, the bias does not reflect weaknesses of the underlying mathematical methods, but is mainly due to the effects of discretization of the spatial domain. We investigate two approaches to correcting the bias: a Newton–Raphson-type correction and Richardson extrapolation. In simulation experiments, Richardson extrapolation performs best.     

A pool-adjacent-violators-algorithm approach to detect infinite parameter estimates in one-regressor dose–response models with asymptotes

Binary response models are often applied in dose–response settings where the number of dose levels is limited. Commonly, one can find cases where the maximum likelihood estimation process for these models produces infinite values for at least one of the parameters, often corresponding to the ‘separated data’ issue. Algorithms for detecting such data have been proposed, but are usually incorporated directly into in the parameter estimation. Additionally, they do not consider the use of asymptotes in the model formulation. In order to study this phenomenon in greater detail, we define the class of specifiably degenerate functions where this can occur (including the popular logistic and Weibull models) that allows for asymptotes in the dose–response specification. We demonstrate for this class that the well-known pool-adjacent-violators algorithm can efficiently pre-screen for non-estimable data. A simulation study demonstrates the frequency with which this problem can occur for various response models and conditions.