A robust and efficient adaptive reweighted estimator of multivariate location and scatter

This article proposes a reweighted estimator of multivariate location and scatter, with weights adaptively computed from the data. Its breakdown point and asymptotic behavior under elliptical distributions are established. This adaptive estimator is able to attain simultaneously the maximum possible breakdown point for affine equivariant estimators and full asymptotic efficiency at the multivariate normal distribution. For the special case of hard-rejection weights and the MCD as initial estimator, it is shown to be more efficient than its non-adaptive counterpart for a broad range of heavy-tailed elliptical distributions. A Monte Carlo study shows that the adaptive estimator is as robust as its non-adaptive relative for several types of bias-inducing contaminations, while it is remarkably more efficient under normality for sample sizes as small as 200.     

Asymptotic distributions of nonparametric regression estimators for longitudinal or functional data

The estimation of a regression function by kernel method for longitudinal or functional data is considered. In the context of longitudinal data analysis, a random function typically represents a subject that is often observed at a small number of time points, while in the studies of functional data the random realization is usually measured on a dense grid. However, essentially the same methods can be applied to both sampling plans, as well as in a number of settings lying between them. In this paper general results are derived for the asymptotic distributions of real-valued functions with arguments which are functionals formed by weighted averages of longitudinal or functional data. Asymptotic distributions for the estimators of the mean and covariance functions obtained from noisy observations with the presence of within-subject correlation are studied. These asymptotic normality results are comparable to those standard rates obtained from independent data, which is illustrated in a simulation study. Besides, this paper discusses the conditions associated with sampling plans, which are required for the validity of local properties of kernel-based estimators for longitudinal or functional data.     

U-statistics in danach spaces

U-statistics in Banach spaces are considered and thoroughly investigated. The martingale structure, estimates of moments, the law of large numbers, the central limit theorem, the invariance principle, estimates of the rate of convergence, and large deviations are established     

The σ<Emphasis Type="Italic">g</Emphasis>-Drazin Inverse and the Generalized Mbekhta Decomposition

In this paper we define and study an extension of the g-Drazin for elements of a Banach algebra and for bounded linear operators based on an isolated spectral set rather than on an isolated spectral point. We investigate salient properties of the new inverse and its continuity, and illustrate its usefulness with an application to differential equations. Generalized Mbekhta subspaces are introduced and the corresponding extended Mbekhta decomposition gives a characterization of circularly isolated spectral sets.     

Entropy estimate for high-dimensional monotonic functions

We establish upper and lower bounds for the metric entropy and bracketing entropy of the class of d-dimensional bounded monotonic functions under L^p norms. It is interesting to see that both the metric entropy and bracketing entropy have different behaviors for p<d/(d-1) and p>d/(d-1). We apply the new bounds for bracketing entropy to establish a global rate of convergence of the MLE of a d-dimensional monotone density.     

Estimation of intrinsic processes affected by additive fractal noise

Fractal Gaussian models have been widely used to represent the singular behavior of phenomena arising in different applied fields; for example, fractional Brownian motion and fractional Gaussian noise are considered as monofractal models in subsurface hydrology and geophysical studies Mandelbrot [The Fractal Geometry of Nature, Freeman Press, San Francisco, 1982 [13]]. In this paper, we address the problem of least-squares linear estimation of an intrinsic fractal input random field from the observation of an output random field affected by fractal noise (see Angulo et al. [Estimation and filtering of fractional generalised random fields, J. Austral. Math. Soc. A 69 (2000) 1-26 [2]], Ruiz-Medina et al. [Fractional generalized random fields on bounded domains, Stochastic Anal. Appl. 21 (2003a) 465-492], Ruiz-Medina et al. [Fractional-order regularization and wavelet approximation to the inverse estimation problem for random fields, J. Multivariate Anal. 85 (2003b) 192-216]. Conditions on the fractality order of the additive noise are studied to obtain a bounded inversion of the associated Wiener-Hopf equation. A stable solution is then obtained in terms of orthogonal bases of the reproducing kernel Hilbert spaces associated with the random fields involved. Such bases are constructed from orthonormal wavelet bases (see Angulo and Ruiz-Medina [Multiresolution approximation to the stochastic inverse problem, Adv. in Appl. Probab. 31 (1999) 1039-1057], Angulo et al. [Wavelet-based orthogonal expansions of fractional generalized random fields on bounded domains, Theoret. Probab. Math. Stat. (2004), in press]). A simulation study is carried out to illustrate the influence of the fractality orders of the output random field and the fractal additive noise on the stability of the solution derived.     

Regular variation of order 1 nonlinear AR-ARCH models

We prove both geometric ergodicity and regular variation of the stationary distribution for a class of nonlinear stochastic recursions that includes nonlinear AR-ARCH models of order 1. The Lyapounov exponent for the model, the index of regular variation and the spectral measure for the regular variation all are characterized by a simple two-state Markov chain.     

Asymptotic expansions in sequential estimation for the first-order random coefficient autoregressive model: Regenerative approach

Edgeworth expansions for the distribution of a sequential least squares estimator in the random coefficient autoregressive (RCA) model are derived. The regenerative approach to second-order asymptotic analysis of Markov-type statistical models is developed.     

A new joint model for longitudinal and survival data with a cure fraction

We develop a new joint cure rate model for longitudinal and survival data. The model allows for multiple longitudinal markers as well as a cure structure for the survival component based on the promotion time cure rate model, as described in Ibrahim et al. (Bayesian Survival Analysis, Springer, New York, 2001). Several characteristics and properties of the new model are discussed and examined. A real dataset from a melanoma clinical trial is given to demonstrate the methodology.     

Exploring interactions in high-dimensional genomic data: an overview of Logic Regression, with applications

Logic Regression is an adaptive regression methodology mainly developed to explore high-order interactions in genomic data. Logic Regression is intended for situations where most of the covariates in the data to be analyzed are binary. The goal of Logic Regression is to find predictors that are Boolean (logical) combinations of the original predictors. In this article, we give an overview of the methodology and discuss some applications. We also describe the software for Logic Regression, which is available as an R and S-Plus package.