Subsampling,symmetrization, and robust interpolation

The recently developed subsampling methodology has been shown to be valid for the construction of large-sample confidence regions for a general unknown parameter 9 under very minimal conditions. Nevertheless, in some specific cases—e.g. in the case of the sample mean of i.i.d. data—it has been noted that the subsampling distribution estimator underperforms as compared to alternative estimators such as the bootstrap or the asymptotic normal distribution (with estimated variance). In the present report we introduce a (partially) symmetrized.     

Extrapolation of subsampling distribution estimators: The i.i.d. and strong mixing cases

Politis & Romano (1994) proposed a general subsampling methodology for the construction of large‐sample confidence regions for an arbitrary parameter under minimal conditions. Nevertheless, the subsampling distribution estimators may sometimes be inefficient (in the case of the sample mean of i.i.d. data, for instance) as compared to alternative estimators such as the bootstrap and/or the asymptotic normal distribution (with estimated variance). The authors investigate here the extent to which the performance of subsampling distribution estimators can be improved by interpolation and extrapolation techniques, while at the same time retaining the robustness property of consistent distribution estimation even in nonregular cases; both i.i.d. and weakly dependent (mixing) observations are considered.     

Tolerance regions in multivariate linear regression

The construction of tolerance regions is investigated for a multivariate linear regression model under the multivariate normality assumption. In the context of such a model, a tolerance region is a region that will contain, with a certain confidence, at least a specified proportion of the population distribution, for a fixed value of the independent variable in the regression model. The necessary framework is developed for the Monte Carlo estimation of the tolerance factor. Some approximations are developed for the tolerance factor, and the accuracy of the approximations is numerically investigated. The approximations provide tolerance factors that are quite easy to compute, and the numerical results indicate the situations where the approximations are satisfactory. The computations are illustrated using an example.     

Generalized Scan Statistics for Disease Surveillance

In applying scan statistics for disease surveillance, it would be valuable to have an integrated model that simultaneously includes environmental covariates and spatial correlation. In this paper, a generalized scan statistics under quasi‐likelihood functions is proposed to address this issue. We use a two‐step estimation process to obtain estimates of coefficients and adapt a bootstrapping method for the minimal p‐value to address the multiple‐testing problem. Under suitable conditions, the proposed method is consistent and can control the type I error rate. Simulations and applications to real data sets are used to evaluate the method.     

Assessing parameter uncertainty via bootstrap likelihood ratio confidence regions

In this paper, we show that, under certain regularity conditions, constructing likelihood ratio confidence regions using a boostrap estimate of the distribution of the likelihood ratio statistic-instead of the usual chi 2 approximation-leads to regions which have a coverage error of O(n- 2), which is the same as that achieved using a Bartlett-corrected likelihood ratio statistic. We use the boostrap method to assess the uncertainty associated with dose-response parameters that arise in models for the Japanese atomic bomb survivors data.     

Network reliability for multipath TCP networks with a retransmission mechanism under the time constraint

It is essential to reduce data latency and guarantee quality of service for modern computer networks. The emerging networking protocol, Multipath Transmission Control Protocol, can reduce data latency by transmitting data through multiple minimal paths (MPs) and ensure data integrity by the packets retransmission mechanism. The bandwidth of each edge can be considered as multi-state in computer networks because different situations, such as failures, partial failures and maintenance, exist. We evaluate network reliability for a multi-state retransmission flow network through which the data can be successfully transmitted by means of multiple MPs under the time constraint. By generating all minimal bandwidth patterns, the proposed algorithm can satisfy these requirements to calculate network reliability. An example and a practical case of the Pan-European Research and Education Network are applied to demonstrate the proposed algorithm.     

Tolerance regions for random vectors and best linear predictors

In the case that vectors X and Y have a joint multivariate normal distribution, tolerance regions are found for the best linear predictor of Y using X if samples are used to estimate the regression coeffierante. Tolerance regions are also found for Y. In addition, simultaneous tolerance intervals for all linear functions of Y or of the best linear predictor of Y using X are found.     

Regression Tolerance Intervals

In this article, we discuss the utility of tolerance intervals for various regression models. We begin with a discussion of tolerance intervals for linear and nonlinear regression models. We then introduce a novel method for constructing nonparametric regression tolerance intervals by extending the well-established procedure for univariate data. Simulation results and application to real datasets are presented to help visualize regression tolerance intervals and to demonstrate that the methods we discuss have coverage probabilities very close to the specified nominal confidence level.     

Generalized Empirical Likelihood Inference in Generalized Linear Models for Longitudinal Data

In this article, the generalized linear model for longitudinal data is studied. A generalized empirical likelihood method is proposed by combining generalized estimating equations and quadratic inference functions based on the working correlation matrix. It is proved that the proposed generalized empirical likelihood ratios are asymptotically chi-squared under some suitable conditions, and hence it can be used to construct the confidence regions of the parameters. In addition, the maximum empirical likelihood estimates of parameters are obtained, and their asymptotic normalities are proved. Some simulations are undertaken to compare the generalized empirical likelihood and normal approximation-based method in terms of coverage accuracies and average areas/lengths of confidence regions/intervals. An example of a real data is used for illustrating our methods.     

Maximum likelihood estimation for directional conditionally autoregressive models

A spatial process observed over a lattice or a set of irregular regions is usually modeled using a conditionally autoregressive (CAR) model. The neighborhoods within a CAR model are generally formed using only the inter-distances or boundaries between the regions. To accommodate directional spatial variation, a new class of spatial models is proposed using different weights given to neighbors in different directions. The proposed model generalizes the usual CAR model by accounting for spatial anisotropy. Maximum likelihood estimators are derived and shown to be consistent under some regularity conditions. Simulation studies are presented to evaluate the finite sample performance of the new model as compared to the CAR model. Finally, the method is illustrated using a data set on the crime rates of Columbus, OH and on the elevated blood lead levels of children under the age of 72 months observed in Virginia in the year of 2000.     

Prediction of record values

Given a number of record values from independent and identically distributed random variables with a continuous distribution function F, our aim is to predict future record values under suitable assumptions on the tail of F. In this paper, we are primarily concerned with finding reasonable tolerance regions for future record values. Two methods are proposed. The first one deals with the case where we observe only record values. The second one makes use of the information brought by the complete sample.     

Bayesian density regression for discrete outcomes

We develop Bayesian models for density regression with emphasis on discrete outcomes. The problem of density regression is approached by considering methods for multivariate density estimation of mixed scale variables, and obtaining conditional densities from the multivariate ones. The approach to multivariate mixed scale outcome density estimation that we describe represents discrete variables, either responses or covariates, as discretised versions of continuous latent variables. We present and compare several models for obtaining these thresholds in the challenging context of count data analysis where the response may be over‐ and/or under‐dispersed in some of the regions of the covariate space. We utilise a nonparametric mixture of multivariate Gaussians to model the directly observed and the latent continuous variables. The paper presents a Markov chain Monte Carlo algorithm for posterior sampling, sufficient conditions for weak consistency, and illustrations on density, mean and quantile regression utilising simulated and real datasets.     

On multistage group screening designs

In this article we discuss multistage group screening in which group-factors contain differing number of factors. We describe a procedure for grouping the factors in the absence of concrete prior information, so that the relative testing cost is minimal. It Is shown that under quite general conditions, these designs will require fewer runs than the equivalent designs in which the group-factors contain same number of factors.     

Semiparametric Estimation of a Two-component Mixture Model where One Component is known

Abstract.  We consider a two-component mixture model where one component distribution is known while the mixing proportion and the other component distribution are unknown. These kinds of models were first introduced in biology to study the differences in expression between genes. The various estimation methods proposed till now have all assumed that the unknown distribution belongs to a parametric family. In this paper, we show how this assumption can be relaxed. First, we note that generally the above model is not identifiable, but we show that under moment and symmetry conditions some 'almost everywhere' identifiability results can be obtained. Where such identifiability conditions are fulfilled we propose an estimation method for the unknown parameters which is shown to be strongly consistent under mild conditions. We discuss applications of our method to microarray data analysis and to the training data problem. We compare our method to the parametric approach using simulated data and, finally, we apply our method to real data from microarray experiments.     

Asymptotic properties of one-step M-estimators

We study the asymptotic behavior of one-step M-estimators based on not necessarily independent identically distributed observations. In particular, we find conditions for asymptotic normality of these estimators. Asymptotic normality of one-step M-estimators is proven under a wide spectrum of constraints on the exactness of initial estimators. We discuss the question of minimal restrictions on the exactness of initial estimators. We also discuss the asymptotic behavior of the solution to an M-equation closest to the parameter under consideration. As an application, we consider some examples of one-step approximation of quasi-likelihood estimators in nonlinear regression.     

Asymptotic normality of numbers of observations in random regions determined by order statistics

In this paper, we study the joint limiting behaviour of numbers of observations that fall into regions determined by order statistics and Borel sets. We show that suitably centred and normed versions of these numbers are asymptotically multivariate normal under some conditions. We consider two cases: one where the population distribution function is discontinuous and the other where it is continuous and the order statistics are extreme. Finally, we compare results obtained for the two cases with their analogues for absolutely continuous distribution function and central-order statistics.     

A Spatio-Temporal Model for Functional Magnetic Resonance Imaging Data – with a View to Resting State Networks

Abstract.  Functional magnetic resonance imaging (fMRI) is a technique for studying the active human brain. During the fMRI experiment, a sequence of MR images is obtained, where the brain is represented as a set of voxels. The data obtained are a realization of a complex spatio-temporal process with many sources of variation, both biological and technical. We present a spatio-temporal point process model approach for fMRI data where the temporal and spatial activation are modelled simultaneously. It is possible to analyse other characteristics of the data than just the locations of active brain regions, such as the interaction between the active regions. We discuss both classical statistical inference and Bayesian inference in the model. We analyse simulated data without repeated stimuli both for location of the activated regions and for interactions between the activated regions. An example of analysis of fMRI data, using this approach, is presented.     

Empirical likelihood for varying-coefficient semiparametric mixed-effects errors-in-variables models with longitudinal data

In this paper, the empirical likelihood inferences for varying-coefficient semiparametric mixed-effects errors-in-variables models with longitudinal data are investigated. We construct the empirical log-likelihood ratio function for the fixed-effects parameters and the mean parameters of random-effects. The empirical log-likelihood ratio at the true parameters is proven to be asymptotically $\chi ^2_{q+r}$ , where $q$ and $r$ are dimensions of the fixed and random effects respectively, and the corresponding confidence regions for them are then constructed. We also obtain the maximum empirical likelihood estimator of the parameters of interest, and prove it is the asymptotically normal under some suitable conditions. A simulation study and a real data application are undertaken to assess the finite sample performance of the proposed method.     

Accuracy of edge detection methods with local information in speckled imagery

We compare the accuracy of five approaches for contour detection in speckled imagery. Some of these methods take advantage of the statistical properties of speckled data, and all of them employ active contours using B-spline curves. Images obtained with coherent illumination are affected by a noise called speckle, which is inherent to the imaging process. These data have been statistically modeled by a multiplicative model using the G0 distribution, under which regions with different degrees of roughness can be characterized by the value of a parameter. We use this information to find boundaries between regions with different textures. We propose and compare five strategies for boundary detection: three based on the data (maximum discontinuity on raw data, fractal dimension and maximum likelihood) and two based on estimates of the roughness parameter (maximum discontinuity and anisotropic smoothed roughness estimates). In order to compare these strategies, a Monte Carlo experience was performed to assess the accuracy of fitting a curve to a region. The probability of finding the correct edge with less than a specified error is estimated and used to compare the techniques. The two best procedures are then compared in terms of their computational cost and, finally, we show that the maximum likelihood approach on the raw data using the G0 law is the best technique.     

FDR control of detected regions by multiscale matched filtering

Feature extraction from observed noisy samples is a common important problem in statistics and engineering. This paper presents a novel general statistical approach to the region detection problem in long data sequences. The proposed technique is a multiscale kernel regression in conjunction with statistical multiple testing for region detection while controlling the false discovery rate (FDR) and maximizing the signal-to-noise ratio via matched filtering. This is achieved by considering a one-dimensional region detection problem as its equivalent zero-dimensional peak detection problem. The detection method does not require a priori knowledge of the shape of the nonzero regions. However, if the shape of the nonzero regions is known a priori, e.g., rectangular pulse, the signal regions can also be reconstructed from the detected peaks, seen as their topological point representatives. Simulations show that the method can effectively perform signal detection and reconstruction in the simulated data under high noise conditions, while controlling the FDR of detected regions and their reconstructed length.