Estimating spatial quantile regression with functional coefficients: A robust semiparametric framework

This paper considers an estimation of semiparametric functional (varying)-coefficient quantile regression with spatial data. A general robust framework is developed that treats quantile regression for spatial data in a natural semiparametric way. The local M-estimators of the unknown functional-coefficient functions are proposed by using local linear approximation, and their asymptotic distributions are then established under weak spatial mixing conditions allowing the data processes to be either stationary or nonstationary with spatial trends. Application to a soil data set is demonstrated with interesting findings that go beyond traditional analysis.Comment: Published in at http://dx.doi.org/10.3150/12-BEJ480 the Bernoulli (http://isi.cbs.nl/bernoulli/) by the International Statistical Institute/Bernoulli Society (http://isi.cbs.nl/BS/bshome.htm

Local Adaptive Grouped Regularization and its Oracle Properties for Varying Coefficient Regression

Varying coefficient regression is a flexible technique for modeling data where the coefficients are functions of some effect-modifying parameter, often time or location in a certain domain. While there are a number of methods for variable selection in a varying coefficient regression model, the existing methods are mostly for global selection, which includes or excludes each covariate over the entire domain. Presented here is a new local adaptive grouped regularization (LAGR) method for local variable selection in spatially varying coefficient linear and generalized linear regression. LAGR selects the covariates that are associated with the response at any point in space, and simultaneously estimates the coefficients of those covariates by tailoring the adaptive group Lasso toward a local regression model with locally linear coefficient estimates. Oracle properties of the proposed method are established under local linear regression and local generalized linear regression. The finite sample properties of LAGR are assessed in a simulation study and for illustration, the Boston housing price data set is analyzed.Comment: 30 pages, one technical appendix, two figure

Estimation in semiparametric spatial regression

Nonparametric methods have been very popular in the last couple of decades in time series and regression, but no such development has taken place for spatial models. A rather obvious reason for this is the curse of dimensionality. For spatial data on a grid evaluating the conditional mean given its closest neighbors requires a four-dimensional nonparametric regression. In this paper a semiparametric spatial regression approach is proposed to avoid this problem. An estimation procedure based on combining the so-called marginal integration technique with local linear kernel estimation is developed in the semiparametric spatial regression setting. Asymptotic distributions are established under some mild conditions. The same convergence rates as in the one-dimensional regression case are established. An application of the methodology to the classical Mercer and Hall wheat data set is given and indicates that one directional component appears to be nonlinear, which has gone unnoticed in earlier analyses.Additive approximation; asymptotic theory; conditional autoregression; local linear kernel estimate; marginal integration; semiparametric regression; spatial mixing process

Local linear spatial regression

A local linear kernel estimator of the regression function x\mapsto g(x):=E[Y_i|X_i=x], x\in R^d, of a stationary (d+1)-dimensional spatial process {(Y_i,X_i),i\in Z^N} observed over a rectangular domain of the form I_n:={i=(i_1,...,i_N)\in Z^N| 1\leq i_k\leq n_k,k=1,...,N}, n=(n_1,...,n_N)\in Z^N, is proposed and investigated. Under mild regularity assumptions, asymptotic normality of the estimators of g(x) and its derivatives is established. Appropriate choices of the bandwidths are proposed. The spatial process is assumed to satisfy some very general mixing conditions, generalizing classical time-series strong mixing concepts. The size of the rectangular domain I_n is allowed to tend to infinity at different rates depending on the direction in Z^N.Comment: Published at http://dx.doi.org/10.1214/009053604000000850 in the Annals of Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical Statistics (http://www.imstat.org

Local Linear Fitting Under Near Epoch Dependence: Uniform consistency with Convergence Rates

Local linear fitting is a popular nonparametric method in statistical and econometric modelling. Lu and Linton (2007) established the pointwise asymptotic distribution for the local linear estimator of a nonparametric regression function under the condition of near epoch dependence. In this paper, we further investigate the uniform consistency of this estimator. The uniform strong and weak consistencies with convergence rates for the local linear fitting are established under mild conditions. Furthermore, general results regarding uniform convergence rates for nonparametric kernel-based estimators are provided. The results of this paper will be of wide potential interest in time series semiparametric modelling.α-mixing, local linear fitting, near epoch dependence, convergence rates, uniform consistency

Specification testing in nonlinear and nonstationary time series autoregression

This paper considers a class of nonparametric autoregressive models with nonstationarity. We propose a nonparametric kernel test for the conditional mean and then establish an asymptotic distribution of the proposed test. Both the setting and the results differ from earlier work on nonparametric autoregression with stationarity. In addition, we develop a new bootstrap simulation scheme for the selection of a suitable bandwidth parameter involved in the kernel test as well as the choice of a simulated critical value. The finite-sample performance of the proposed test is assessed using one simulated example and one real data example.Comment: Published in at http://dx.doi.org/10.1214/09-AOS698 the Annals of Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical Statistics (http://www.imstat.org

Local bilinear multiple-output quantile/depth regression

A new quantile regression concept, based on a directional version of Koenker and Bassett's traditional single-output one, has been introduced in [Ann. Statist. (2010) 38 635-669] for multiple-output location/linear regression problems. The polyhedral contours provided by the empirical counterpart of that concept, however, cannot adapt to unknown nonlinear and/or heteroskedastic dependencies. This paper therefore introduces local constant and local linear (actually, bilinear) versions of those contours, which both allow to asymptotically recover the conditional halfspace depth contours that completely characterize the response's conditional distributions. Bahadur representation and asymptotic normality results are established. Illustrations are provided both on simulated and real data.Comment: Published at http://dx.doi.org/10.3150/14-BEJ610 in the Bernoulli (http://isi.cbs.nl/bernoulli/) by the International Statistical Institute/Bernoulli Society (http://isi.cbs.nl/BS/bshome.htm

Nonparametric Specification Testing for Nonlinear Time Series with Nonstationarity

This paper considers a nonparametric time series regression model with a nonstationary regressor. We construct a nonparametric test for testing whether the regression is of a known parametric form indexed by a vector of unknown parameters. We establish the asymptotic distribution of the proposed test statistic. Both the setting and the results differ from earlier work on nonparametric time series regression with stationarity. In addition, we develop a bootstrap simulation scheme for the selection of suitable bandwidth parameters involved in the kernel test as well as the choice of simulated critical values. An example of implementation is given to show that the proposed test works in practice.integrated regressor, kernel test, nonparametric regression, nonstationary time series, random walk

Estimating Value At Risk

Significantly driven by JP Morgan's RiskMetrics system with EWMA (exponentially weighted moving average) forecasting technique, value-at-risk (VaR) has turned to be a popular measure of the degree of various risks in financial risk management. In this paper we propose a new approach termed skewed-EWMA to forecast the changing volatility and formulate an adaptively efficient procedure to estimate the VaR. Differently from the JP Morgan's standard-EWMA, which is derived from a Gaussian distribution, and the Guermat and Harris (2001)'s robust-EWMA, from a Laplace distribution, we motivate and derive our skewed-EWMA procedure from an asymmetric Laplace distribution, where both skewness and heavy tails in return distribution and the time-varying nature of them in practice are taken into account. An EWMA-based procedure that adaptively adjusts the shape parameter controlling the skewness and kurtosis in the distribution is suggested. Backtesting results show that our proposed skewed-EWMA method offers a viable improvement in forecasting VaR