Remarks on between estimator in the intraclass correlation model with missing data

In this paper, we consider the between estimator under the intraclass correlation model with missing data. We give a necessary and sufficient condition for existing exact simultaneous confidence intervals for all contrasts in the means under the between transformed model, which indicates the F-test statistic and simultaneous confidence intervals, constructed by Seo et al. [T. Seo, J. Kikuchi, K. Koizumi, On simultaneous confidence intervals for all contracts in the means of the intraclass correlation model with missing data, J. Multivariate Anal. 97 (2006) 1976–1983] based on the between estimator, is invalid. Furthermore, using the distribution of the between estimator, we present the exact test statistics and confidence intervals for partial contrasts.     

Exact inference on contrasts in means of intraclass correlation models with missing responses

Intraclass correlation models with missing data at random are considered. With a properly reduced model, a general method, which allows repeated observations with missing data in a non-monotone pattern, is proposed to construct exact test statistics and simultaneous confidence intervals for linear contrasts in the means. Simulation results are given to compare exact and asymptotic simultaneous confidence intervals. A real example is provided for the illustration of the proposed method.     

Fiducial inference on the largest mean of a multivariate normal distribution

Inference on the largest mean of a multivariate normal distribution is a surprisingly difficult and unexplored topic. Difficulties arise when two or more of the means are simultaneously the largest mean. Our proposed solution is based on an extension of R.A. Fisher’s fiducial inference methods termed generalized fiducial inference. We use a model selection technique along with the generalized fiducial distribution to allow for equal largest means and alleviate the overestimation that commonly occurs. Our proposed confidence intervals for the largest mean have asymptotically correct frequentist coverage and simulation results suggest that they possess promising small sample empirical properties. In addition to the theoretical calculations and simulations we also applied this approach to the air quality index of the four largest cities in the northeastern United States (Baltimore, Boston, New York, and Philadelphia).     

Empirical likelihood for semiparametric regression model with missing response data

A bias-corrected technique for constructing the empirical likelihood ratio is used to study a semiparametric regression model with missing response data. We are interested in inference for the regression coefficients, the baseline function and the response mean. A class of empirical likelihood ratio functions for the parameters of interest is defined so that undersmoothing for estimating the baseline function is avoided. The existing data-driven algorithm is also valid for selecting an optimal bandwidth. Our approach is to directly calibrate the empirical log-likelihood ratio so that the resulting ratio is asymptotically chi-squared. Also, a class of estimators for the parameters of interest is constructed, their asymptotic distributions are obtained, and consistent estimators of asymptotic bias and variance are provided. Our results can be used to construct confidence intervals and bands for the parameters of interest. A simulation study is undertaken to compare the empirical likelihood with the normal approximation-based method in terms of coverage accuracies and average lengths of confidence intervals. An example for an AIDS clinical trial data set is used for illustrating our methods.     

Exact confidence coefficients of simultaneous confidence intervals for multinomial proportions

Simultaneous confidence intervals for multinomial proportions are useful in many areas of science. Since 1964, approximate simultaneous 1-α confidence intervals have been proposed for multinomial proportions. Although at each point in the parameter space, these confidence sets have asymptotic 1-α coverage probability, the exact confidence coefficients of these simultaneous confidence intervals for a fixed sample size are unknown before.In this paper, we propose a procedure for calculating exact confidence coefficients for simultaneous confidence intervals of multinomial proportions for any fixed sample size. With this methodology, exact confidence coefficients can be clearly derived, and the point at which the infimum of the coverage probability occurs can be clearly identified.     

Empirical likelihood inference in partially linear single-index models for longitudinal data

The empirical likelihood method is especially useful for constructing confidence intervals or regions of parameters of interest. Yet, the technique cannot be directly applied to partially linear single-index models for longitudinal data due to the within-subject correlation. In this paper, a bias-corrected block empirical likelihood (BCBEL) method is suggested to study the models by accounting for the within-subject correlation. BCBEL shares some desired features: unlike any normal approximation based method for confidence region, the estimation of parameters with the iterative algorithm is avoided and a consistent estimator of the asymptotic covariance matrix is not needed. Because of bias correction, the BCBEL ratio is asymptotically chi-squared, and hence it can be directly used to construct confidence regions of the parameters without any extra Monte Carlo approximation that is needed when bias correction is not applied. The proposed method can naturally be applied to deal with pure single-index models and partially linear models for longitudinal data. Some simulation studies are carried out and an example in epidemiology is given for illustration.     

More on the inadmissibility of step-up

Cohen and Sackrowitz [Characterization of Bayes procedures for multiple endpoint problems and inadmissibility of the step-up procedure, Ann. Statist. 33 (2005) 145-158] proved that the step-up multiple testing procedure is inadmissible for a multivariate normal model with unknown mean vector and known intraclass covariance matrix. The hypotheses tested are each mean is zero vs. each mean is positive. The risk function is a 2×1 vector where one component is average size and the other component is one minus average power. In this paper, we extend the inadmissibility result to several different models, to two-sided alternatives, and to other risk functions. The models include one-parameter exponential families, independent t-variables, independent χ²-variables, t-tests arising from the analysis of variance, and t-tests arising from testing treatments against a control. The additional risk functions are linear combinations where one component is the false discovery rate (FDR).     

Nonparametric inference on the difference of location parameters of correlated variables from fragmentary samples

Summary In this paper, two types of robust estimators and approximate confidence intervals for the difference of location parameters of correlated random variables are proposed and investigated when some observations are missing. It is shown that the suggested estimators are consistent and asymptotically normally distributed. In addition, the proposed approximate confidence intervals are also shown to enjoy some nice asymptotic properties.     

The Factor Analytic Approach to Simultaneous Inference in the General Linear Model

Abstract

Consider the general linear model (GLM) Y = Xβ + ε. Suppose Θ₁,…, Θ_k, a subset of the β's, are of interest; Θ₁,…, Θ_k may be treatment contrasts in an ANOVA setting or regression coefficients in a response surface setting. Existing simultaneous confidence intervals for Θ₁,…, Θ_k are relatively conservative or, in the case of the MEANS option in PROC GLM of SAS, possibly misleading. The difficulty is with the multidimensionality of the integration required to compute exact coverage probability when X does not follow a nice textbook design. Noting that such exact coverage probabilities can be computed if the correlation matrix R of the estimators of Θ₁, …, Θ_k has a one-factor structure in the factor analytic sense, it is proposed that approximate simultaneous confidence intervals be computed by approximating R with the closest one-factor structure correlation matrix. Computer simulations of hundreds of randomly generated designs in the settings of regression, analysis of covariance, and unbalanced block designs show that the coverage probabilities are practically exact, more so than can be anticipated by even second-order probability bounds.     

平衡单向分类模型中几种同时置信区间的比较

在平衡单向分类模型中,通过对对照的三种同时置信区间长度表达式的分析,可以看到,在不同情况下,三种区间的优劣.     

On the simultaneous anova and manova tests

Summary In the present paper, the simultaneous confidence intervals associated with the Simultaneous Multivariate Analysis of Variance (SMANOVA) tests are derived when the hypotheses can be tested in a “quasi-independentrd manner. With respect to the Simultaneous Analysis of Variance (SANOVA) test, the quasi-independence is weakened and the distribution problems are investigated when the joint distribution of the sums of squares associated with various hypotheses is a “multivariate chi-square ” distribution. It is shown that the lengths of the confidence intervals associated with the SANOVA test (dependent case) are shorter than the lengths of the corresponding confidence intervals obtained by using Scheffé’s method when the hypotheses on several contrasts are tested simultaneously. It is also shown that the SMANOVA test yields narrower confidence intervals than the MANOVA test when several quasi-independent hypotheses are tested simultaneously. Part of the work was done under U.S. Air Force Contract No. AF 49 (638) 213 when the author was at the University of North Carolina. Part of the author’s Ph. D. thesis submitted to the University of Minnesota.     

Analysis of variance with general errors and grouped and non-grouped data: Some iterative algorithms

In this paper we consider some iterative estimation algorithms, which are valid to analyse the variance of data, which may be either non-grouped or grouped with different classification intervals. This situation appears, for instance, when data is collected from different sources and the grouping intervals differ from one source to another. The analysis of variance is carried out by means of general linear models, whose error terms may be general. An initial procedure in the line of the EM, although it does not necessarily agree with it, opens the paper and gives rise to a simplified version where we avoid the double iteration, which implicitly appears in the EM and, also, in the initial procedure mentioned above. The asymptotic stochastic properties of the resulting estimates have been investigated in depth and used to test ANOVA hypothesis.     

Minimax least squares and quasiminimax estimation in linear models with missing values

We consider the problem of estimating the parameter vector in the linear model when observations on the independent variables are partially missing or incorrect. New estimators are developed, which systematically combine prior information with the incomplete data. We compare these methods with the alternative strategy of deleting incomplete observations.Support by Deutsche Forschungsgemeinschaft, Grant No. 284/1-2 is gratefully acknowledged.     

Empirical likelihood for single-index models

The empirical likelihood method is especially useful for constructing confidence intervals or regions of the parameter of interest. This method has been extensively applied to linear regression and generalized linear regression models. In this paper, the empirical likelihood method for single-index regression models is studied. An estimated empirical log-likelihood approach to construct the confidence region of the regression parameter is developed. An adjusted empirical log-likelihood ratio is proved to be asymptotically standard chi-square. A simulation study indicates that compared with a normal approximation-based approach, the proposed method described herein works better in terms of coverage probabilities and areas (lengths) of confidence regions (intervals).     

Smoothed jackknife empirical likelihood method for ROC curve

In this paper we propose a smoothed jackknife empirical likelihood method to construct confidence intervals for the receiver operating characteristic (ROC) curve. By applying the standard empirical likelihood method for a mean to the jackknife sample, the empirical likelihood ratio statistic can be calculated by simply solving a single equation. Therefore, this procedure is easy to implement. Wilks’ theorem for the empirical likelihood ratio statistic is proved and a simulation study is conducted to compare the performance of the proposed method with other methods.     

Second-order accurate inference on eigenvalues of covariance and correlation matrices

Edgeworth expansions and saddlepoint approximations for the distributions of estimators of certain eigenfunctions of covariance and correlation matrices are developed. These expansions depend on second-, third-, and fourth-order moments of the sample covariance matrix. Expressions for and estimators of these moments are obtained. The expansions and moment expressions are used to construct second-order accurate confidence intervals for the eigenfunctions. The expansions are illustrated and the results of a small simulation study that evaluates the finite-sample performance of the confidence intervals are reported.     

Oracally efficient spline smoothing of nonlinear additive autoregression models with simultaneous confidence band

Under weak conditions of smoothness and mixing, we propose spline-backfitted spline (SBS) estimators of the component functions for a nonlinear additive autoregression model that is both computationally expedient for analyzing high dimensional large time series data, and theoretically reliable as the estimator is oracally efficient and comes with asymptotically simultaneous confidence band. Simulation evidence strongly corroborates with the asymptotic theory.     

Bartlett correctable two-sample adjusted empirical likelihood

We propose a two-sample adjusted empirical likelihood (AEL) to construct confidence regions for the difference of two d-dimensional population means. This method eliminates the non-definition of the usual two-sample empirical likelihood (EL) and is shown to be Bartlett correctable. We further show that when the adjustment level is half the Bartlett correction factor for the usual two-sample EL, the two-sample AEL has the same high-order precision as the EL with Bartlett correction. To enhance the performance of the two-sample AEL with adjustment level being half the Bartlett correction factor, we propose a less biased estimate of the Bartlett correction factor. The efficiency of the proposed method is illustrated by simulations and a real data example.     

Nonparametric variance function estimation with missing data

In this paper, a fixed design regression model where the errors follow a strictly stationary process is considered. In this model the conditional mean function and the conditional variance function are unknown curves. Correlated errors when observations are missing in the response variable are assumed. Four nonparametric estimators of the conditional variance function based on local polynomial fitting are proposed. Expressions of the asymptotic bias and variance of these estimators are obtained. A simulation study illustrates the behavior of the proposed estimators.     

Mean–Mean Multiple Comparison Displays for Families of Linear Contrasts

Traditional tabular and graphical displays of results of simultaneous confidence intervals or hypothesis tests are deficient in several respects. Expanding on earlier work, we present new mean–mean multiple comparison graphs that succinctly and compactly display the results of traditional procedures for multiple comparisons of population means or linear contrasts involving means. The MMC plot can be used with unbalanced, multifactor designs with covariates. After reviewing the construction of these displays in the S language (S-Plus and R), we demonstrate their application to four multiple comparison scenarios.