Partial covariate adjusted regression

Covariate adjusted regression (CAR) is a recently proposed adjustment method for regression analysis where both the response and predictors are not directly observed [?entürk, D., Müller, H.G., 2005. Covariate adjusted regression. Biometrika 92, 75–89]. The available data have been distorted by unknown functions of an observable confounding covariate. CAR provides consistent estimators for the coefficients of the regression between the variables of interest, adjusted for the confounder. We develop a broader class of partial covariate adjusted regression (PCAR) models to accommodate both distorted and undistorted (adjusted/unadjusted) predictors. The PCAR model allows for unadjusted predictors, such as age, gender and demographic variables, which are common in the analysis of biomedical and epidemiological data. The available estimation and inference procedures for CAR are shown to be invalid for the proposed PCAR model. We propose new estimators and develop new inference tools for the more general PCAR setting. In particular, we establish the asymptotic normality of the proposed estimators and propose consistent estimators of their asymptotic variances. Finite sample properties of the proposed estimators are investigated using simulation studies and the method is also illustrated with a Pima Indians diabetes data set.     

Bayesian inverse problems with heterogeneous variance

We consider inverse problems in Hilbert spaces under correlated Gaussian noise, and use a Bayesian approach to find their regularized solution. We focus on mildly ill-posed inverse problems with fractional noise, using a novel wavelet-based vaguelette–vaguelette approach. It allows us to apply sequence space methods without assuming that all operators are simultaneously diagonalizable. The results are proved for more general bases and covariance operators. Our primary aim is to study posterior contraction rate in such inverse problems over Sobolev classes and compare it to the derived minimax rate. Secondly, we study effect of plugging in a consistent estimator of variances in sequence space on the posterior contraction rate. This result is applied to the problem with error in forward operator. Thirdly, we show that empirical Bayes posterior distribution with a plugged-in maximum marginal likelihood estimator of the prior scale contracts at the optimal rate, adaptively, in the minimax sense.     

Rates of convergence of autocorrelation estimates for periodically correlated autoregressive Hilbertian processes

Autoregressive Hilbertian (ARH) processes are of great importance in the analysis of functional time series data and estimation of the autocorrelation operators attracts the attention of various researchers. In this paper, we study estimators of the autocorrelation operators of periodically correlated autoregressive Hilbertian processes of order one (PCARH(1)), which is an extension of ARH(1) processes. The estimation method is based on the spectral decomposition of the covariance operator and considers two main cases: known and unknown eigenvectors. We show the consistency in the mean integrated quadratic sense of the estimators of the autocorrelation operators and present upper bounds for the corresponding rates.     

Analysis of multivariate arma processes with non-stationary innovations

The paper considers vector ARMA processes with nonstationary innovations. It is suggested that this class of models provide a very efficient framework for nonstationary problems. A generalization of the Yule-Walker equations relating the underlying process is obtained. Identification procedures are discussed. The associated prediction problem is solved using the Hilbert space approach.     

General linear processes in Hilbert spaces and prediction

We discuss a general definition of linear processes in Hilbert spaces that takes into account the outstanding role played by this model in prediction theory.     

Chebyshev's inequality for Hilbert-space valued random elements with estimated mean and covariance

We obtain a generalization of the Chebyshev's inequality for random elements taking values in a separable Hilbert space with estimated mean and covariance.     

Estimation of the Sobol indices in a linear functional multidimensional model

We consider a functional linear model where the explicative variables are known stochastic processes taking values in a Hilbert space, the main example is given by Gaussian processes in L²([0,1])${\mathbb{L}}^2([0,1])$. We propose estimators of the Sobol indices in this functional linear model. Our estimators are based on U-statistics. We prove the asymptotic normality and the efficiency of our estimators and we compare them from a theoretical and practical point of view with classical estimators of Sobol indices.     

Forecasting contemporal time series aggregates

Forecast of a contemporal aggregate of several time series can be obtained from ‘1’ an aggregate series, ‘2’ individual component processes, or ‘3’ a joint multiple forecasting model. Through general Hilbert space theory and some illustrative examples, this paper establishes the relative efficiencies among the three methods     

A New Type of Discrete Self-Decomposability and Its Application to Continuous-Time Markov Processes for Modeling Count Data Time Series

We propose a family of extended thinning operators, indexed by a parameter γ in [0, 1), with the boundary case of γ=0 corresponding to the well-known binomial thinning operator. The extended thinning operators can be used to construct a class of continuous-time Markov processes for modeling count time series data. The class of stationary distributions of these processes is called generalized discrete self-decomposable, denoted by DSD (γ). We obtain characterization results for the DSD (γ) class and investigate relationships among the classes for different γ's.     

Bayesian Inference and Forecasting in Dynamic Neural Networks with Fully Markov Switching ARCH Noises

We deal with one-layer feed-forward neural network for the Bayesian analysis of nonlinear time series. Noises are modeled nonlinearly and nonnormally, by means of ARCH models whose parameters are all dependent on a hidden Markov chain. Parameter estimation is performed by sampling from the posterior distribution via Evolutionary Monte Carlo algorithm, in which two new crossover operators have been introduced. Unknown parameters of the model also include the missing values which can occur within the observed series, so, considering future values as missing, it is also possible to compute point and interval multi-step-ahead predictions.     

Bayes Hilbert Spaces

A Bayes linear space is a linear space of equivalence classes of proportional σ‐finite measures, including probability measures. Measures are identified with their density functions. Addition is given by Bayes' rule and substraction by Radon–Nikodym derivatives. The present contribution shows the subspace of square‐log‐integrable densities to be a Hilbert space, which can include probability and infinite measures, measures on the whole real line or discrete measures. It extends the ideas from the Hilbert space of densities on a finite support towards Hilbert spaces on general measure spaces. It is also a generalisation of the Euclidean structure of the simplex, the sample space of random compositions. In this framework, basic notions of mathematical statistics get a simple algebraic interpretation. A key tool is the centred‐log‐ratio transformation, a generalization of that used in compositional data analysis, which maps the Hilbert space of measures into a subspace of square‐integrable functions. As a consequence of this structure, distances between densities, orthonormal bases, and Fourier series representing measures become available. As an application, Fourier series of normal distributions and distances between them are derived, and an example related to grain size distributions is presented. The geometry of the sample space of random compositions, known as Aitchison geometry of the simplex, is obtained as a particular case of the Hilbert space when the measures have discrete and finite support.     

Moments and generating functions for the absortion distribution and its negative binomial analogue

The moments of the absorption are difficult to obtain. The generating functions are basic hypergeometric functions. This paper shows how to define two shift operators that allow elementary arguments to be used to develop recursions for the expected values of general functions. The exact moments of the distribution follow. The generating function for the negative binomial analogue gives the moments directly.     

A Note on Estimation in Hilbertian Linear Models

We study estimation and prediction in linear models where the response and the regressor variable both take values in some Hilbert space. Our main objective is to obtain consistency of a principal component‐based estimator for the regression operator under minimal assumptions. In particular, we avoid some inconvenient technical restrictions that have been used throughout the literature. We develop our theory in a time‐dependent setup that comprises as important special case the autoregressive Hilbertian model.     

Fractional approaches for the distribution of innovation sequence of INAR(1) processes

Abstract

In this paper, we present a fractional decomposition of the probability generating function of the innovation process of the first-order non-negative integer-valued autoregressive [INAR(1)] process to obtain the corresponding probability mass function. We also provide a comprehensive review of integer-valued time series models, based on the concept of thinning operators with geometric-type marginals. In particular, we develop two fractional approaches to obtain the distribution of innovation processes of the INAR(1) model and show that the distribution of the innovations sequence has geometric-type distribution. These approaches are discussed in detail and illustrated through a few examples.     

Sparse Principal Component Analysis in Hilbert Space

Technical advances in many areas have produced more complicated high‐dimensional data sets than the usual high‐dimensional data matrix, such as the fMRI data collected in a period for independent trials, or expression levels of genes measured in different tissues. Multiple measurements exist for each variable in each sample unit of these data. Regarding the multiple measurements as an element in a Hilbert space, we propose Principal Component Analysis (PCA) in Hilbert space. The principal components (PCs) thus defined carry information about not only the patterns of variations in individual variables but also the relationships between variables. To extract the features with greatest contributions to the explained variations in PCs for high‐dimensional data, we also propose sparse PCA in Hilbert space by imposing a generalized elastic‐net constraint. Efficient algorithms to solve the optimization problems in our methods are provided. We also propose a criterion for selecting the tuning parameter.     

Strong convergence properties for partial sums of asymptotically negatively associated random vectors in Hilbert spaces

Abstract

In this paper, we investigate the almost sure convergence for partial sums of asymptotically negatively associated (ANA, for short) random vectors in Hilbert spaces. The Khintchine-Kolmogorov type convergence theorem, three series theorem and the Kolmogorov type strong law of large numbers for partial sums of ANA random vectors in Hilbert spaces are obtained. The results obtained in the paper generalize some corresponding ones for independent random vectors and negatively associated random vectors in Hilbert spaces.     

Young bitableaux,lattice paths and Hilbert functions

A recent result on the enumeration of p-tuples of nonintersecting lattice paths in an integral rectangle is used to deduce a formula of Abhyankar for standard Young bitableaux of certain type, which gives the Hilbert function of a class of determinantal ideals. The lattice path formula is also shown to yield the numerator of the Hilbert series of these determinantal ideals and the h-vectors of the associated simplicial complexes. As a consequence, the a-invariant of these determinantal ideals is obtained in some cases, extending an earlier result of Gräbe. Some problems concerning generalizations of these results to ‘higher dimensions’ are also discussed. In an appendix, the equivalence of Abhyankar's formula for unitableaux of a given shape and a formula of Hodge, obtained in connection with his determination of Hilbert functions of Schubert varieties in Grassmannians, is outlined.     

Signed binomial approximation of binomial mixtures via differential calculus for linear operators

We obtain sharp estimates in signed binomial approximation of binomial mixtures with respect to the total variation distance. We provide closed form expressions for the leading terms, and show that the corresponding leading coefficients depend on the zeros of appropriate Krawtchouk polynomials. The special case of Pólya–Eggenberger distributions is discussed in detail. Our approach is based on a differential calculus for linear operators represented by stochastic processes, which allows us to give unified proofs.     

再生核Hilbert空间展开的函数型回归模型变量选择

针对协变量是函数型、响应变量是标量的多元函数型回归模型,文章提出了函数系数基于再生核Hilbert空间展开的变量选择方法。首先,利用带积分余项的泰勒展开式和再生核Hilbert空间内积性质将模型转化为结构化形式,其次,通过自适应弹性网惩罚对结构化模型中的组间和组内系数同时进行压缩。结果证明了这种压缩估计具有Oracle性质,蒙特卡罗模拟结果也显示新方法在不同样本量、不同噪声和变量相关性干扰下均优于基于普通基函数展开的变量选择方法,且尤其适用于原始协变量高度相关的情形。最后,通过分析一个商品房平均销售价格影响因素数据演示了新方法的应用。     

Pricing of American options in discrete time using least squares estimates with complexity penalties

Pricing of American options in discrete time is considered, where the option is allowed to be based on several underlying stocks. It is assumed that the price processes of the underlying stocks are given by Markov processes. We use the Monte Carlo approach to generate artificial sample paths of these price processes, and then we use nonparametric regression estimates to estimate from this data so-called continuation values, which are defined as mean values of the American option for given values of the underlying stocks at time t subject to the constraint that the option is not exercised at time t. As nonparametric regression estimates we use least squares estimates with complexity penalties, which include as special cases least squares spline estimates, least squares neural networks, smoothing splines and orthogonal series estimates. General results concerning rate of convergence are presented and applied to derive results for the special cases mentioned above. Furthermore the pricing of American options is illustrated by simulated data.