Semiparametric latent variable regression models for spatiotemporal modelling of mobile source particles in the greater Boston area

Summary.  Traffic particle concentrations show considerable spatial variability within a metropolitan area. We consider latent variable semiparametric regression models for modelling the spatial and temporal variability of black carbon and elemental carbon concentrations in the greater Boston area. Measurements of these pollutants, which are markers of traffic particles, were obtained from several individual exposure studies that were conducted at specific household locations as well as 15 ambient monitoring sites in the area. The models allow for both flexible non-linear effects of covariates and for unexplained spatial and temporal variability in exposure. In addition, the different individual exposure studies recorded different surrogates of traffic particles, with some recording only outdoor concentrations of black or elemental carbon, some recording indoor concentrations of black carbon and others recording both indoor and outdoor concentrations of black carbon. A joint model for outdoor and indoor exposure that specifies a spatially varying latent variable provides greater spatial coverage in the area of interest. We propose a penalized spline formulation of the model that relates to generalized kriging of the latent traffic pollution variable and leads to a natural Bayesian Markov chain Monte Carlo algorithm for model fitting. We propose methods that allow us to control the degrees of freedom of the smoother in a Bayesian framework. Finally, we present results from an analysis that applies the model to data from summer and winter separately.     

Likelihood ratio tests in linear mixed models with one variance component

Summary.  We consider the problem of testing null hypotheses that include restrictions on the variance component in a linear mixed model with one variance component and we derive the finite sample and asymptotic distribution of the likelihood ratio test and the restricted likelihood ratio test. The spectral representations of the likelihood ratio test and the restricted likelihood ratio test statistics are used as the basis of efficient simulation algorithms of their null distributions. The large sample χ ² mixture approximations using the usual asymptotic theory for a null hypothesis on the boundary of the parameter space have been shown to be poor in simulation studies. Our asymptotic calculations explain these empirical results. The theory of Self and Liang applies only to linear mixed models for which the data vector can be partitioned into a large number of independent and identically distributed subvectors. One-way analysis of variance and penalized splines models illustrate the results.     

Small area estimates of labour force participation under a multinomial logit mixed model

Summary.  A new methodology is developed for estimating unemployment or employment characteristics in small areas, based on the assumption that the sample totals of unemployed and employed individuals follow a multinomial logit model with random area effects. The method is illustrated with UK labour force data aggregated by sex–age groups. For these data, the accuracy of direct estimates is poor in comparison with estimates that are derived from the multinomial logit model. Furthermore, two different estimators of the mean-squared errors are given: an analytical approximation obtained by Taylor linearization and an estimator based on bootstrapping. A simulation study for comparison of the two estimators shows the good performance of the bootstrap estimator.     

On the mvue for the distribution in the multivariate normal regression

In the multivariate normal regression setting, the estimability of a distribution is studied generalizing earlier results for the univariate case. The MVUE of an estimable distribution is obtained.     

Estimation of the mean of an exponential distribution in the presence of an outlier

Kale and Sinha (1971) have found an estimator of the mean of an exponential distribution in the présence of an outlying observation with higher expected value. Here an alternative estimator of the mean is proposed and it is compared with the estimator of Kale and Sinha (1971) and the maximum likelihood estimator given by Kale (1975). The proposed estimator is found to be more efficient than the latter two estimators in some cases.     

Autocorrelation and tolerance limits

Tolerance limits are limits that include a specified proportion of the population at a given confidence level. They are used to make sure that the production will not be outside specifications. Tolerance limits are either designed based on the normality assumption, or nonparametric tolerance limits are established. In either case, no provision for autocorrelated processes is made in the available design tables of tolerance limits. It is shown how to construct tolerance limits to cover a specified proportion of the population when autocorrelation is present in the process. A comparison of four different tolerance limits is provided, and recommendations are given for choosing the "best" estimator of the process variability for the construction of tolerance limits.     

Parametric scalp mapping and inference via spatially smooth linear models for mismatch negativity studies

Mismatch negativity (MMN) is a neurophysiological tool that can be used to investigate various facets of comprehension. Subjects are presented with different stimuli to elicit the MMN response, which is derived from electroencephalography (EEG) signals recorded at electrodes across the brain. We propose a methodology to extend single electrode analyses of MMN data by generating smooth scalp maps of estimated experimental effects. It is shown that penalized least squares estimates of effect maps can be produced using a two step procedure involving (a) ANOVA at each electrode and (b) spatial smoothing across electrodes. A Fisher von-Mises kernel is used for smoothing scalp maps with cross-validated bandwidth selection. The methodology is applied to a case control study involving aphasics (language disordered individuals). Analysis of residuals shows possible heteroscedasticity and non-Gaussian tail behavior. For robust inference, a semiparametric multivariate approach is proposed to determine the significance of parametric maps. A variety of global and regional test statistics are developed to investigate the significance of spatial patterns in treatment effects. The methodology is seen to confirm previous findings from single electrode analysis and identifies some new significant spatial patterns of difference between controls and aphasics.     

An iterative sparse algorithm for the penalized maximum likelihood estimator in mixed effects model

In this paper, we propose a new iterative sparse algorithm (ISA) to compute the maximum likelihood estimator (MLE) or penalized MLE of the mixed effects model. The sparse approximation based on the arrow-head (A-H) matrix is one solution which is popularly used in practice. The A-H method provides an easy computation of the inverse of the Hessian matrix and is computationally efficient. However, it often has non-negligible error in approximating the inverse of the Hessian matrix and in the estimation. Unlike the A-H method, in the ISA, the sparse approximation is applied “iteratively” to reduce the approximation error at each Newton Raphson step. The advantages of the ISA over the exact and A-H method are illustrated using several synthetic and real examples.     

Bounds on the empirical Bayes and compound risks of truncated versions of Robbins's estimator of a binomial parameter

Robbins (1956) in his original paper on empirical Bayes methods suggested a method of estimating a binomial success probability. We give explicit bounds for the empirical Bayes risk of natural variants of the Robbins estimator that show convergence to an optimal risk at $\text{O(n}{}^{\mathrm{<mtext>?</mtext><mtext>1</mtext><mtext>2</mtext>}}\text{)}$ rate. Bounds that yield the same convergence rate are also obtained in the related compound estimation problem.