Multilevel maximum likelihood estimation with application to covariance matrices

The asymptotic variance of the maximum likelihood estimate is proved to decrease when the maximization is restricted to a subspace that contains the true parameter value. Maximum likelihood estimation allows a systematic fitting of covariance models to the sample, which is important in data assimilation. The hierarchical maximum likelihood approach is applied to the spectral diagonal covariance model with different parameterizations of eigenvalue decay, and to the sparse inverse covariance model with specified parameter values on different sets of nonzero entries. It is shown computationally that using smaller sets of parameters can decrease the sampling noise in high dimension substantially.     

Model Checks in Inverse Regression Models with Convolution‐Type Operators

Abstract. We consider the problem of testing parametric assumptions in an inverse regression model with a convolution‐type operator. An L ₂ ‐type goodness‐of‐fit test is proposed which compares the distance between a parametric and a non‐parametric estimate of the regression function. Asymptotic normality of the corresponding test statistic is shown under the null hypothesis and under a general non‐parametric alternative with different rates of convergence in both cases. The feasibility of the proposed test is demonstrated by means of a small simulation study. In particular, the power of the test against certain types of alternative is investigated. Finally, an empirical example is provided, in which the proposed methods are applied to the determination of the shape of the luminosity profile of the elliptical galaxy NGC 5017.     

A bootstrap control chart for inverse Gaussian percentiles

The conventional Shewhart-type control chart is developed essentially on the central limit theorem. Thus, the Shewhart-type control chart performs particularly well when the observed process data come from a near-normal distribution. On the other hand, when the underlying distribution is unknown or non-normal, the sampling distribution of a parameter estimator may not be available theoretically. In this case, the Shewhart-type charts are not available. Thus, in this paper, we propose a parametric bootstrap control chart for monitoring percentiles when process measurements have an inverse Gaussian distribution. Through extensive Monte Carlo simulations, we investigate the behaviour and performance of the proposed bootstrap percentile charts. The average run lengths of the proposed percentage charts are investigated.     

Optimal quantization applied to sliced inverse regression

In this paper we consider a semiparametric regression model involving a d-dimensional quantitative explanatory variable X and including a dimension reduction of X via an index β′X. In this model, the main goal is to estimate the Euclidean parameter β and to predict the real response variable Y conditionally to X. Our approach is based on sliced inverse regression (SIR) method and optimal quantization in L^p-norm. We obtain the convergence of the proposed estimators of β and of the conditional distribution. Simulation studies show the good numerical behavior of the proposed estimators for finite sample size.     

A comparative study of confidence intervals for negative binomial proportion

In this paper, we investigate four existing and three new confidence interval estimators for the negative binomial proportion (i.e., proportion under inverse/negative binomial sampling). An extensive and systematic comparative study among these confidence interval estimators through Monte Carlo simulations is presented. The performance of these confidence intervals are evaluated in terms of their coverage probabilities and expected interval widths. Our simulation studies suggest that the confidence interval estimator based on saddlepoint approximation is more appealing for large coverage levels (e.g., nominal level≤1% ) whereas the score confidence interval estimator is more desirable for those commonly used coverage levels (e.g., nominal level>1% ). We illustrate these confidence interval construction methods with a real data set from a maternal congenital heart disease study.     

A stochastic procedure to solve linear ill-posed problems

In this work, we propose a stochastic procedure of Robbins–Monro type to resolve linear inverse problems in Hilbert space. We study the probability of large deviation between the exact solution and the approximated one and build a confidence domain for the approximated solution while precising the rate of convergence. To check the validity of our work, we give a simulation application into a deconvolution problem.     

Asymptotic expansions for sequential confidence levels in inverse linear regression

Non-linear renewal theory is used to derive second order asymptotic expansions for the coverage probability of a fixed-width sequential confidence interval for an unknown parameter xin the inverse linear regression model. These expansions are obtained for a two-stage sequential procedure, proposed by Perng and Tong (1974) for the construction of a confidence interval for x.     

Effect of a new selection procedure on adaptive estimation

Adaptive estimation of parameters of some failure time distributionsis considered. A new procedure named the F-procedure has beendeveloped for selecting an appropriate model out of two possible models by Pandey et.al. (1991). Applying this F-procedure adaptive estimatorsof parameters of exponential, Wei bull, inverse Gaussian (IG) and Wald failure time distributions have been proposed in this paper. Comparison of these estimators has been undertaken with MLE's of the respective parameters and with some previous adaptiveestimators by simulation of samples using the Monte Carlo method.Adaptive estimation of parameters of some failure time distributions is considered. A new procedure named the F-procedure has been developedfor selecting an appropriate model out of two possible models by Pandey et.al. (1991). Applying this F-procedure adaptive estimators of parameters of exponential, Wei bull, inverse Gaussian (IG) and Wald failure time distributions have been proposed in this paper. Comparison of these estimators has been undertaken with MLE's of the respective parameters and with some previous adaptive estimators by simulation of samples using the Monte Carlo method.     

Application of the Bootstrap Approach to the Choice of Dimension and the α Parameter in the SIRα Method

To reduce the dimensionality of regression problems, sliced inverse regression approaches make it possible to determine linear combinations of a set of explanatory variables X related to the response variable Y in general semiparametric regression context. From a practical point of view, the determination of a suitable dimension (number of the linear combination of X) is important. In the literature, statistical tests based on the nullity of some eigenvalues have been proposed. Another approach is to consider the quality of the estimation of the effective dimension reduction (EDR) space. The square trace correlation between the true EDR space and its estimate can be used as goodness of estimation. In this article, we focus on the SIR_α method and propose a naïve bootstrap estimation of the square trace correlation criterion. Moreover, this criterion could also select the α parameter in the SIR_α method. We indicate how it can be used in practice. A simulation study is performed to illustrate the behavior of this approach.