A frequency domain comparison of estimation methods for gaussian fields

This paper demonstrates that well-known parameter estimation methods for Gaussian fields place different emphasis on the high and low frequency components of the data. As a consequence, the relative importance of the frequencies under the objective of the analysis should be taken into account when selecting an estimation method, in addition to other considerations such as statistical and computational efficiency. The paper also shows that when noise is added to the Gaussian field, maximum pseudolikelihood automatically sets the smoothing parameter of the model equal to one. A simulation study then indicates that generalised cross-validation is more robust than maximum likelihood un-

der model misspecification in smoothing and image restoration problems. This has implications for Bayesian procedures since these use the same weightings of the frequencies as the likelihood.     

Questions and Answers

In addition to his contributions to biostatistics and clinical trials, Paul Meier had a long-term interest in the legal applications of statistics. As part of this, he had extensive experience as a statistical consultant. Legal consulting can be a minefield, but as a result of his background, Paul had excellent advice to give to those starting out on how to function successfully in this environment.     

Automatic locally adaptive smoothing for tree-based set estimation

Tree-based methods similar to CART have recently been utilized for problems in which the main goal is to estimate some set of interest. It is often the case that the boundary of the true set is smooth in some sense, however tree-based estimates will not be smooth, as they will be a union of ‘boxes’. We propose a general methodology for smoothing such sets that allows for varying levels of smoothness on the boundary automatically. The method is similar to the idea underlying support vector machines, which is applying a computationally simple technique to data after a non-linear mapping to produce smooth estimates in the original space. In particular, we consider the problem of level-set estimation for regression functions and the dyadic tree-based method of Willett and Nowak [Minimax optimal level-set estimation, IEEE Trans. Image Process. 16 (2007), pp. 2965–2979].     

A survey of tests for exponentiality

A wide selection of tests for exponentiality is discussed and compared. Power computations, using simulations, were done for each procedure. Certain tests (e.g. Gnedenko (1969), Lin and Mudholkar (1980), Harris (1376), Cox and Oakes (1384), and Deshpande (1983)) performed well for alternative distributions with non-monotonic hazard rates, while others (e.g. Deshpande (1983), Gail and Gastwirth (1978), Kolmogorov-Smirnov (LillViefors (1969)), Hahn and Shapiro (1967), Hollander and Proschan (1972), and Cox and Oakes (1984)) fared well for monotonic hazard rates. Of all the procedures compared, the score test presented in Cox and Oakes (1984) appears to be the best if one does not have a particular alternative in mind.     

Linking aggregate production planning with marketing pricing decisions

In this paper, we link production planning decisions to marketing decisions that involve the price of product groups. The focus of this paper is the development of a closed-loop procedure for aggregate production planning and pricing. We seek to satisfy uncertain demand while minimizing total costs that include material, labour, and inventory holding costs. The procedure is useful for variable demand to update short-term aggregate plans.     

Linearity of regression for non-adjacent record values

Let X₁,X₂,… be a sequence of iid random variables having a continuous distribution; by R₁,R₂,… denote the corresponding record values. All the distributions allowing linearity of regressions either E(R_m+k|R_m) or E(R_m|R_m+k) are identified.     

On the equivalence of some test criteria based on BAN estimators for the multivariate exponential family

For the general multivariate exponential family of distributions it is shown that Rao's test criterion based on efficient scores is algebraically identical to the general chi-squared criterion based on maximum likelihood estimates and, similarly, that the Wald statistic is algebraically identical to the general minimum modified chi-squared statistic using linearization; these results are valid also for the multisample versions. Thus, these are extensions to the general exponential family of the findings due to Silvey (1970) and Bhapkar (1966), respectively, for the special case of the multinomial family.It is also shown that the general forms of the chi-squared and modified chi-squared criteria reduce to their respective well-known forms for the multivariate symmetric power series distribution. This finding is, thus, an extension of results noted by Ferguson (1958) and Clickner (1976) for the special case of the multinomial distribution.     

The estimation of derivatives of a nonparametric regression function when the data are correlated

A kernel estimator of a derivative of arbitrary order of a nonparametric average population curve is considered for a correlated-errors model with balanced replicate measurements at each design point. Asymptotic expansions of the mean squared error are derived for two classes of correlation functions in the model. Consistency, choice of smoothing parameter, and rates of convergence are examined for the important special cases of estimating the first and second derivatives.     

Flexible Copula Density Estimation with Penalized Hierarchical B‐splines

The paper introduces a new method for flexible spline fitting for copula density estimation. Spline coefficients are penalized to achieve a smooth fit. To weaken the curse of dimensionality, instead of a full tensor spline basis, a reduced tensor product based on so called sparse grids (Notes Numer. Fluid Mech. Multidiscip. Des., 31, 1991, 241‐251) is used. To achieve uniform margins of the copula density, linear constraints are placed on the spline coefficients, and quadratic programming is used to fit the model. Simulations and practical examples accompany the presentation.     

Point estimation of lower quantiles based on two-sampling scheme

ABSTRACT

Consider a two-sampling scheme in which an initial sample is first taken from the underlying population and then by assuming a suitable restriction on this sample, some more data points are observed as a new restricted sample. This sampling scheme is used to do inference about the lower quantiles of the underlying distribution. The results are compared with those of simple random sampling in view of mean squared error and Pitman’s measure of closeness criteria for exponential and uniform distributions. It will be shown that the proposed sampling scheme would improve the performance of the point estimators of the lower quantiles of the population.