Flexible modeling based on copulas in nonparametric median regression

Consider the model Y=m(X)+ε, where m(⋅)=med(Y|⋅) is unknown but smooth. It is often assumed that ε and X are independent. However, in practice this assumption is violated in many cases. In this paper we propose modeling the dependence between ε and X by means of a copula model, i.e. (ε,X)∼C_θ(F_ε(⋅),F_X(⋅)), where C_θ is a copula function depending on an unknown parameter θ, and F_ε and F_X are the marginals of ε and X. Since many parametric copula families contain the independent copula as a special case, the so-obtained regression model is more flexible than the ‘classical’ regression model.We estimate the parameter θ via a pseudo-likelihood method and prove the asymptotic normality of the estimator, based on delicate empirical process theory. We also study the estimation of the conditional distribution of Y given X. The procedure is illustrated by means of a simulation study, and the method is applied to data on food expenditures in households.     

On the block thresholding wavelet estimators with censored data

We consider block thresholding wavelet-based density estimators with randomly right-censored data and investigate their asymptotic convergence rates. Unlike for the complete data case, the empirical wavelet coefficients are constructed through the Kaplan-Meier estimators of the distribution functions in the censored data case. On the basis of a result of Stute [W. Stute, The central limit theorem under random censorship, Ann. Statist. 23 (1995) 422-439] that approximates the Kaplan-Meier integrals as averages of i.i.d. random variables with a certain rate in probability, we can show that these wavelet empirical coefficients can be approximated by averages of i.i.d. random variables with a certain error rate in L². Therefore we can show that these estimators, based on block thresholding of empirical wavelet coefficients, achieve optimal convergence rates over a large range of Besov function classes , p≥2, q≥1 and nearly optimal convergence rates when 1≤p<2. We also show that these estimators achieve optimal convergence rates over a large class of functions that involve many irregularities of a wide variety of types, including chirp and Doppler functions, and jump discontinuities. Therefore, in the presence of random censoring, wavelet estimators still provide extensive adaptivity to many irregularities of large function classes. The performance of the estimators is tested via a modest simulation study.     

Nonconcave penalized inverse regression in single-index models with high dimensional predictors

In this paper we aim to estimate the direction in general single-index models and to select important variables simultaneously when a diverging number of predictors are involved in regressions. Towards this end, we propose the nonconcave penalized inverse regression method. Specifically, the resulting estimation with the SCAD penalty enjoys an oracle property in semi-parametric models even when the dimension, p_n, of predictors goes to infinity. Under regularity conditions we also achieve the asymptotic normality when the dimension of predictor vector goes to infinity at the rate of p_n=o(n^1/3) where n is sample size, which enables us to construct confidence interval/region for the estimated index. The asymptotic results are augmented by simulations, and illustrated by analysis of an air pollution dataset.     

Microstructure noise in the continuous case: The pre-averaging approach

This paper presents a generalized pre-averaging approach for estimating the integrated volatility, in the presence of noise. This approach also provides consistent estimators of other powers of volatility — in particular, it gives feasible ways to consistently estimate the asymptotic variance of the estimator of the integrated volatility. We show that our approach, which possesses an intuitive transparency, can generate rate optimal estimators (with convergence rate n^−1/4$n^{-1/4}$).     

Copula-Type Estimators for Flexible Multivariate Density Modeling Using Mixtures

Copulas are popular as models for multivariate dependence because they allow the marginal densities and the joint dependence to be modeled separately. However, they usually require that the transformation from uniform marginals to the marginals of the joint dependence structure is known. This can only be done for a restricted set of copulas, for example, a normal copula. Our article introduces copula-type estimators for flexible multivariate density estimation which also allow the marginal densities to be modeled separately from the joint dependence, as in copula modeling, but overcomes the lack of flexibility of most popular copula estimators. An iterative scheme is proposed for estimating copula-type estimators and its usefulness is demonstrated through simulation and real examples. The joint dependence is modeled by mixture of normals and mixture of normal factor analyzer models, and mixture of t and mixture of t-factor analyzer models. We develop efficient variational Bayes algorithms for fitting these in which model selection is performed automatically. Based on these mixture models, we construct four classes of copula-type densities which are far more flexible than current popular copula densities, and outperform them in a simulated dataset and several real datasets. Supplementary material for this article is available online.     

Semi-parametric estimation of partially linear single-index models

One of the most difficult problems in applications of semi-parametric partially linear single-index models (PLSIM) is the choice of pilot estimators and complexity parameters which may result in radically different estimators. Pilot estimators are often assumed to be root-n consistent, although they are not given in a constructible way. Complexity parameters, such as a smoothing bandwidth are constrained to a certain speed, which is rarely determinable in practical situations.In this paper, efficient, constructible and practicable estimators of PLSIMs are designed with applications to time series. The proposed technique answers two questions from Carroll et al. [Generalized partially linear single-index models, J. Amer. Statist. Assoc. 92 (1997) 477-489]: no root-n pilot estimator for the single-index part of the model is needed and complexity parameters can be selected at the optimal smoothing rate. The asymptotic distribution is derived and the corresponding algorithm is easily implemented. Examples from real data sets (credit-scoring and environmental statistics) illustrate the technique and the proposed methodology of minimum average variance estimation (MAVE).     

Optimal rates of convergence in the Weibull model based on kernel-type estimators

Let F be a distribution function in the maximal domain of attraction of the Gumbel distribution such that −log(1−F(x))=x^1/θL(x) for a positive real number θ, called the Weibull tail index, and a slowly varying function L. It is well known that the estimators of θ have a very slow rate of convergence. We establish here a sharp optimality result in the minimax sense, that is when L is treated as an infinite dimensional nuisance parameter belonging to some functional class. We also establish the rate optimal asymptotic property of a data-driven choice of the sample fraction that is used for estimation.     

Statistical inference in a panel data semiparametric regression model with serially correlated errors

We consider a panel data semiparametric partially linear regression model with an unknown vector β of regression coefficients, an unknown nonparametric function g(·) for nonlinear component, and unobservable serially correlated errors. The correlated errors are modeled by a vector autoregressive process which involves a constant intraclass correlation. Applying the pilot estimators of β and g(·), we construct estimators of the autoregressive coefficients, the intraclass correlation and the error variance, and investigate their asymptotic properties. Fitting the error structure results in a new semiparametric two-step estimator of β, which is shown to be asymptotically more efficient than the usual semiparametric least squares estimator in terms of asymptotic covariance matrix. Asymptotic normality of this new estimator is established, and a consistent estimator of its asymptotic covariance matrix is presented. Furthermore, a corresponding estimator of g(·) is also provided. These results can be used to make asymptotically efficient statistical inference. Some simulation studies are conducted to illustrate the finite sample performances of these proposed estimators.     

Characterization of dependence of multidimensional Lévy processes using Lévy copulas

This paper suggests Lévy copulas in order to characterize the dependence among components of multidimensional Lévy processes. This concept parallels the notion of a copula on the level of Lévy measures. As for random vectors, a version of Sklar's theorem states that the law of a general multivariate Lévy process is obtained by combining arbitrary univariate Lévy processes with an arbitrary Lévy copula. We construct parametric families of Lévy copulas and prove a limit theorem, which indicates how to obtain the Lévy copula of a multivariate Lévy process X from the ordinary copula of the random vector X_t for small t.     

Adaptive estimation of the transition density of a particular hidden Markov chain

We study the following model of hidden Markov chain: with (X_i) a real-valued positive recurrent and stationary Markov chain, and (?_i)_1?i?n+1 a noise independent of the sequence (X_i) having a known distribution. We present an adaptive estimator of the transition density based on the quotient of a deconvolution estimator of the density of X_i and an estimator of the density of (X_i,X_i+1). These estimators are obtained by contrast minimization and model selection. We evaluate the L² risk and its rate of convergence for ordinary smooth and supersmooth noise with regard to ordinary smooth and supersmooth chains. Some examples are also detailed.     

Estimation of location and scale parameters for the Burr XII distribution using generalized order statistics

The minimum variance linear unbiased estimators (MVLUE), the best linear invariant estimators (BLIE) and the maximum likelihood estimators (MLE) based on n-selected generalized order statistics are presented for the parameters of the Burr XII distribution.     

Smooth estimation of survival and density functions for a stationary associated process using Poisson weights

Let {X_n,n≥1} be a sequence of stationary non-negative associated random variables with common marginal density f(x). Here we use the empirical survival function as studied in Bagai and Prakasa Rao (1991) and apply the smoothing technique proposed by Gawronski (1980) (see also Chaubey and Sen, 1996) in proposing a smooth estimator of the density function f and that of the corresponding survival function. Some asymptotic properties of the resulting estimators, similar to those obtained in Chaubey and Sen (1996) for the i.i.d. case, are derived. A simulation study has been carried out to compare the new estimator to the kernel estimator of a density function given in Bagai and Prakasa Rao (1996) and the estimator in Buch-Larsen et al. (2005).     

Berry-Esseen bounds for wavelet estimator in a regression model with linear process errors

In this paper, we derive the Berry-Esseen bounds of the wavelet estimator for a nonparametric regression model with linear process errors generated by φ-mixing sequences. As application, by the suitable choice of some constants, the convergence rate O(n^−1/6) of uniformly asymptotic normality of the wavelet estimator is obtained. Our results generalize some known results in the literature.     

Nonparametric estimation of mixed partial derivatives of a multivariate density

On the basis of a random sample of size n on an m-dimensional random vector X, this note proposes a class of estimators f_n^(p) of f^(p), where f is a density of X w.r.t. a σ-finite measure dominated by the Lebesgue measure on R^m, p = (p₁,…,p_m), p_j ≥ 0, fixed integers, and for x = (x₁,…,x_m) in R^m, f^(p)(x) = ?^p₁+…+p_m f(x)/(?^p₁x₁ … ?^p_mx_m). Asymptotic unbiasedness as well as both almost sure and mean square consistencies of f_n^(p) are examined. Further, a necessary and sufficient condition for uniform asymptotic unbisedness or for uniform mean square consistency of f_n^(p) is given. Finally, applications of estimators of this note to certain statistical problems are pointed out.     

A multidimensional Newton-Raphson method and its applications to the existence of asymptotic Fn-estimators and their stochastic expansions

Given a suitable function F_n we define a class of estimators called asymptotic F_n-estimators (i.e., estimators which approximate the solution of F_n(θ) = 0). It is proved that this class is nonvoid if appropriate regularity conditions are fulfilled and if one has at hand a suitable initial estimator. Furthermore, it is shown that F_n-estimators admit a stochastic expansion (which enables to give results on asymptotic expansions for the distribution of these estimators).     

Local invariance principles and their application to density estimation

Summary Letx ₁,...,x _n be independent random variables with uniform distribution over [0, 1] ^d , andX( ⁿ ) be the centered and normalized empirical process associated tox ₁,...,x _n . Given a Vapnik-Chervonenkis classL of bounded functions from [0, 1] ^d intoR of bounded variation, we apply the one-dimensional dyadic scheme of Komlós, Major and Tusnády to get the best possible rate in Dudley's uniform central limit theorem for the empirical process {E ⁽ⁿ⁾(h):hL}. WhenL fulfills some extra condition, we prove there exists some sequenceB _n of Brownian bridges indexed byL such that whereK (L) denotes the maximal variation of the elements ofL. This result is then applied to maximal deviations distributions for kernel density estimators under minimal assumptions on the sequence of bandwith parameters. We also derive some results concerning strong approximations for empirical processes indexed by classes of sets with uniformly small perimeter. For example, it follows from Beck's paper that the above result is optimal, up to a possible factor , whenL is the class of Euclidean balls with radius less thanr.     

Three general multivariate semi-Pareto distributions and their characterizations

Three general multivariate semi-Pareto distributions are developed in this paper. First one—GMP^(k)(III) has univariate Pareto (III) marginals, it is characterized by the minimum of two independent and identically distributed random vectors. Second one—GMSP has univariate semi-Pareto marginals and it is characterized by finite sample minima. Third one—MSP is characterized through a geometric minimization procedure. All these three characterizations are based on the general and the particular solutions of the Euler's functional equations of k-variates.     

Empirical likelihood confidence region for parameter in the errors-in-variables models

This paper proposes a constrained empirical likelihood confidence region for a parameter β₀ in the linear errors-in-variables model: Y_i=x_i^τβ₀+ε_i,X_i=x_i+u_i,(1?i?n), which is constructed by combining the score function corresponding to the squared orthogonal distance with a constrained region of β₀. It is shown that the coverage error of the confidence region is of order n⁻¹, and Bartlett corrections can reduce the coverage errors to n⁻². An empirical Bartlett correction is given for practical implementation. Simulations show that the proposed confidence region has satisfactory coverage not only for large samples, but also for small to medium samples.     

On the risk of estimates for block decreasing densities

A density f=f(x₁,…,x_d) on [0,∞)^d is block decreasing if for each j∈{1,…,d}, it is a decreasing function of x_j, when all other components are held fixed. Let us consider the class of all block decreasing densities on [0,1]^d bounded by B. We shall study the minimax risk over this class using n i.i.d. observations, the loss being measured by the L₁ distance between the estimate and the true density. We prove that if S=log(1+B), lower bounds for the risk are of the form C(S^d/n)^1/(d+2), where C is a function of d only. We also prove that a suitable histogram with unequal bin widths as well as a variable kernel estimate achieve the optimal multivariate rate. We present a procedure for choosing all parameters in the kernel estimate automatically without loosing the minimax optimality, even if B and the support of f are unknown.