Nonparametric Bayes factors based on empirical likelihood ratios

Bayes methodology provides posterior distribution functions based on parametric likelihoods adjusted for prior distributions. A distribution-free alternative to the parametric likelihood is use of empirical likelihood (EL) techniques, well known in the context of nonparametric testing of statistical hypotheses. Empirical likelihoods have been shown to exhibit many of the properties of conventional parametric likelihoods. In this paper, we propose and examine Bayes factors (BF) methods that are derived via the EL ratio approach. Following Kass and Wasserman (1995), we consider Bayes factors type decision rules in the context of standard statistical testing techniques. We show that the asymptotic properties of the proposed procedure are similar to the classical BF's asymptotic operating characteristics. Although we focus on hypothesis testing, the proposed approach also yields confidence interval estimators of unknown parameters. Monte Carlo simulations were conducted to evaluate the theoretical results as well as to demonstrate the power of the proposed test.     

Iterative bias reduction: a comparative study

Multivariate nonparametric smoothers, such as kernel based smoothers and thin plate splines smoothers, are adversely impacted by the sparseness of data in high dimension, also known as the curse of dimensionality. Adaptive smoothers, that can exploit the underlying smoothness of the regression function, may partially mitigate this effect. This paper presents a comparative simulation study of a novel adaptive smoother (IBR) with competing multivariate smoothers available as package or function within the R language and environment for statistical computing. Comparison between the methods are made on simulated datasets of moderate size, from 50 to 200 observations, with two, five or 10 potential explanatory variables, and on a real dataset. The results show that the good asymptotic properties of IBR are complemented by a very good behavior on moderate sized datasets, results which are similar to those obtained with Duchon low rank splines.     

Progressive type II censored order statistics from exponential distributions

In the model of progressive type II censoring, point and interval estimation as well as relations for single and product moments are considered. Based on two-parameter exponential distributions, maximum likelihood estimators (MLEs), uniformly minimum variance unbiased estimators (UMVUEs) and best linear unbiased estimators (BLUEs) are derived for both location and scale parameters. Some properties of these estimators are shown. Moreover, results for single and product moments of progressive type II censored order statistics are presented to obtain recurrence relations from exponential and truncated exponential distributions. These relations may then be used to compute all the means, variances and covariances of progressive type II censored order statistics based on exponential distributions for arbitrary censoring schemes. The presented recurrence relations simplify those given by Aggarwala and Balakrishnan (1996)     

Random coefficients regression models. a review 1

A review is given of random regression coefficients models. The emphasis is put on the problem of estimating the mean regression coefficients and the covariance matrix of the coefficients. Prediction of the individual random coefficients is not discussed. The main purpose of the review is to point to the practical aspects of the models and the problem of statistical inference in finite samples. Some problems for future research are indicated.     

Charts for finding upper percentage points of Student's t in the range .01 to .00001

The Bonferroni t-statistic is a versatile tool in multiple comparisons problems. The need for "oddball percentage points" may lead to extensive tables or heavy computation. Charts of t_p as a function of log p enable near two-decimal accuracy for any percentage point between .01 and .00001     

A Unified Way to Represent and Compute the Cumulative Distribution Function for Noncentral t,F, and χ2

This short article shows an unified approach to representing and computing the cumulative distribution function for noncentral t, F, and χ². Unlike the existing algorithms, which involve different expansion and/or recurrence, the new approach consistently represents all the three noncentral cumulative distribution functions as the integral of the normal cumulative distribution function and χ² density function.     

Computation of the central and noncentral f distributions

Recursion relations suitable for rapid computation are derived for the cumulative distribution of F′ = (X/m)/(Y/n) where X is χ²(λ, m) and Y is independently χ²(n). When n is even no complicated function evaluations are needed. For n odd, a special doubly noncentral t distribution is needed to start the computation. Series representations for this t distribution are given with rigorous bounds on truncation errors. Proper recursion techniques for numerical evaluation of the special functions are given.     

Aspects of Statistical Computing: What Packages for the 1980's Ought to Do

A wish list of desirable statistical computing capabilities is presented. This may help one question which of these capabilities can be satisfied by existing packages, which might be met through reasonable extensions to these packages, which might require substantial new development, and which ought to be supplied by the computing environment rather than the packages. These questions are explored, taking into account the nature of the statistical work and the choices presented by technology. Attention is given to the barriers to be overcome if future statistical packages are to take full advantage of new technology.     

A Theory for Coloring Bivariate Statistical Maps

Consideration of some practical uses of statistical bivariate maps—for example, display of association between variables—leads to principles for making effective use of color to represent data values. Effective color schemes for bivariate maps are viewed as continuous transformations from color models to the unit square with appropriate restrictions involving hue, saturation, and brightness. Several schemes, including those used by the U.S. Census Bureau, are criticized on the basis of this theory.