On the Existence of Monotone Pure‐Strategy Equilibria in Bayesian Games

We generalize Athey's (2001) and McAdams' (2003) results on the existence of monotone pure‐strategy equilibria in Bayesian games. We allow action spaces to be compact locally complete metric semilattices and type spaces to be partially ordered probability spaces. Our proof is based on contractibility rather than convexity of best‐reply sets. Several examples illustrate the scope of the result, including new applications to multi‐unit auctions with risk‐averse bidders.     

Spatially varying dynamic coefficient models

In this work we present a flexible class of linear models to treat observations made in discrete time and continuous space, where the regression coefficients vary smoothly in time and space. This kind of model is particularly appealing in situations where the effect of one or more explanatory processes on the response present substantial heterogeneity in both dimensions. We describe how to perform inference for this class of models and also how to perform forecasting in time and interpolation in space, using simulation techniques. The performance of the algorithm to estimate the parameters of the model and to perform prediction in time is investigated with simulated data sets. The proposed methodology is used to model pollution levels in the Northeast of the United States.     

Bayesian locally optimal design of knockout tournaments

The elimination or knockout format is one of the most common designs for pairing competitors in tournaments and leagues. In each round of a knockout tournament, the losers are eliminated while the winners advance to the next round. Typically, the goal of such a design is to identify the overall best player. Using a common probability model for expressing relative player strengths, we develop an adaptive approach to pairing players each round in which the probability that the best player advances to the next round is maximized. We evaluate our method using simulated game outcomes under several data-generating mechanisms, and compare it to random pairings, to the standard knockout format, and to two variants of the standard format.     

Bayesian monitoring of local residual autocorrelations taking into account the run-length

This paper proposes new two-sided monitoring algorithms for detecting the presence of first order residual autocorrelations in Dynamic Normal Models. The methodology uses a Bayesian decision approach with loss function which takes into account the run-length of the process. The power and mean run-length of the proposed algorithms are analysed by Monte Carlo methods. The results obtained improve those corresponding to the monitoring algorithm for residual autocorrelations proposed in Gargallo and Salvador [2003. Monitoring residual autocorrelations in dynamic linear models. Comm. Statist. Simulation Comput. 32(4), 1079–1104.] with respect to the run-length, and also exhibit more homogeneous behaviour.     

Statistical inference for start-up demonstration tests with rejection of units upon observing d failures

In this paper, we consider the statistical inference for the success probability in the case of start-up demonstration tests in which rejection of units is possible when a pre-fixed number of failures is observed before the required number of consecutive successes are achieved for acceptance of the unit. Since the expected value of the stopping time is not a monotone function of the unknown parameter, the method of moments is not useful in this situation. Therefore, we discuss two estimation methods for the success probability: (1) the maximum likelihood estimation (MLE) via the expectation-maximization (EM) algorithm and (2) Bayesian estimation with a beta prior. We examine the small-sample properties of the MLE and Bayesian estimator. Finally, we present an example to illustrate the method of inference discussed here.     

Bayesian analysis of elapsed times in continuous‐time Markov chains

The authors consider Bayesian analysis for continuous‐time Markov chain models based on a conditional reference prior. For such models, inference of the elapsed time between chain observations depends heavily on the rate of decay of the prior as the elapsed time increases. Moreover, improper priors on the elapsed time may lead to improper posterior distributions. In addition, an infinitesimal rate matrix also characterizes this class of models. Experts often have good prior knowledge about the parameters of this matrix. The authors show that the use of a proper prior for the rate matrix parameters together with the conditional reference prior for the elapsed time yields a proper posterior distribution. The authors also demonstrate that, when compared to analyses based on priors previously proposed in the literature, a Bayesian analysis on the elapsed time based on the conditional reference prior possesses better frequentist properties. The type of prior thus represents a better default prior choice for estimation software.     

Properties of Bayes Factors Based on Test Statistics

Abstract.  This article examines the consistency, interpretation and application of Bayes factors constructed from standard test statistics. Primary conclusions are that Bayes factors based on multinomial and normal test statistics are consistent for suitable choices of the hyperparameters used to specify alternative hypotheses, and that such constructions can be extended to obtain consistent Bayes factors based on likelihood ratio statistics. A connection between Bayes factors based on likelihood ratio statistics and the Bayesian information criterion is exposed, as is a connection between Bayes factors based on F statistics and parametric Bayes factors based on normal-inverse gamma models. Similarly, Bayes factors based on chi-squared statistics for multinomial data are shown to provide accurate approximations to Bayes factors based on multinomial/Dirichlet models. An illustration of how the simple form of these Bayes factors can be exploited to generate easily interpretable summaries of the experimental 'weight of evidence' is provided.     

New methods for fitting multiple sinusoids from irregularly sampled data

A novel framework is proposed for the estimation of multiple sinusoids from irregularly sampled time series. This spectral analysis problem is addressed as an under-determined inverse problem, where the spectrum is discretized on an arbitrarily thin frequency grid. As we focus on line spectra estimation, the solution must be sparse, i.e. the amplitude of the spectrum must be zero almost everywhere. Such prior information is taken into account within the Bayesian framework. Two models are used to account for the prior sparseness of the solution, namely a Laplace prior and a Bernoulli–Gaussian prior, associated to optimization and stochastic sampling algorithms, respectively. Such approaches are efficient alternatives to usual sequential prewhitening methods, especially in case of strong sampling aliases perturbating the Fourier spectrum. Both methods should be intensively tested on real data sets by physicists.     

Efficient estimation for the proportional hazards model with bivariate current status data

We consider efficient estimation of regression and association parameters jointly for bivariate current status data with the marginal proportional hazards model. Current status data occur in many fields including demographical studies and tumorigenicity experiments and several approaches have been proposed for regression analysis of univariate current status data. We discuss bivariate current status data and propose an efficient score estimation approach for the problem. In the approach, the copula model is used for joint survival function with the survival times assumed to follow the proportional hazards model marginally. Simulation studies are performed to evaluate the proposed estimates and suggest that the approach works well in practical situations. A real life data application is provided for illustration.     

Modelling heterogeneous space–time occurrences of earthquakes and its residual analysis

Summary. Earthquake intensities are modelled as a function of previous activity whose specific form is based on established empirical laws in seismology, but whose parameter values can vary from place to place. This model is used for characterizing regional features of seismic activities in and around Japan, and also for exploring regions where the actual seismicity rate systematically deviates from that of the modelled rate.