Spatial autoregressive and moving average Hilbertian processes

This paper addresses the introduction and study of structural properties of Hilbert-valued spatial autoregressive processes (SARH(1) processes), and Hilbert-valued spatial moving average processes (SMAH(1) processes), with innovations given by two-parameter (spatial) matingale differences. For inference purposes, the conditions under which the tensorial product of standard autoregressive Hilbertian (ARH(1)) processes (respectively, of standard moving average Hilbertian (MAH(1)) processes) is a standard SARH(1) process (respectively, it is a standard SMAH(1) process) are studied. Examples related to the spatial functional observation of two-parameter Markov and diffusion processes are provided. Some open research lines are described in relation to the formulation of SARMAH processes, as well as General Spatial Linear Processes in Functional Spaces.     

Asymptotic distribution of the OLS estimator for a purely autoregressive spatial model

We derive the asymptotics of the OLS estimator for a purely autoregressive spatial model. Only low-level conditions are used. As the sample size increases, the spatial matrix is assumed to approach a square-integrable function on the square (0,1)². The asymptotic distribution is a ratio of two infinite linear combinations of χ² variables. The formula involves eigenvalues of an integral operator associated with the function approached by the spatial matrices. Under the conditions imposed identification conditions for the maximum likelihood method and method of moments fail. A corrective two-step procedure using the OLS estimator is proposed.     

On unit roots for spatial autoregressive models

In this paper we consider the unit root problem for one rather simple autoregressive model Y_t,s=aY_t-1,s+bY_t,s-1+?_t,s on a two-dimensional lattice. We show that the growth of variance of Y_t,s is essentially different from corresponding growth in the unit root case for AR(1) or AR(2) time series models. We also show that the dimension of the lattice plays an important role: the growth of variance of autoregressive field on a d-dimensional lattice is different for d=2,3 and d≥4.     

Bayesian multivariate spatial models for roadway traffic crash mapping

We consider several Bayesian multivariate spatial models for estimating the crash rates from different kinds of crashes. Multivariate conditional autoregressive (CAR) models are considered to account for the spatial effect. The models considered are fully Bayesian. A general theorem for each case is proved to ensure posterior propriety under noninformative priors. The different models are compared according to some Bayesian criterion. Markov chain Monte Carlo (MCMC) is used for computation. We illustrate these methods with Texas Crash Data.     

Evaluation of reproducibility for paired functional data

Evaluation of reproducibility is important in assessing whether a new method or instrument can reproduce the results from a traditional gold standard approach. In this paper, we propose a measure to assess measurement agreement for functional data which are frequently encountered in medical research and many other research fields. Formulae to compute the standard error of the proposed estimator and confidence intervals for the proposed measure are derived. The estimators and the coverage probabilities of the confidence intervals are empirically tested for small-to-moderate sample sizes via Monte Carlo simulations. A real data example in physiology study is used to illustrate the proposed statistical inference procedures.     

Stochastic optimization algorithms for barrier dividend strategies

This work focuses on finding optimal barrier policy for an insurance risk model when the dividends are paid to the share holders according to a barrier strategy. A new approach based on stochastic optimization methods is developed. Compared with the existing results in the literature, more general surplus processes are considered. Precise models of the surplus need not be known; only noise-corrupted observations of the dividends are used. Using barrier-type strategies, a class of stochastic optimization algorithms are developed. Convergence of the algorithm is analyzed; rate of convergence is also provided. Numerical results are reported to demonstrate the performance of the algorithm.     

A comparative study of Gaussian geostatistical models and Gaussian Markov random field models

Gaussian geostatistical models (GGMs) and Gaussian Markov random fields (GMRFs) are two distinct approaches commonly used in spatial models for modeling point-referenced and areal data, respectively. In this paper, the relations between GGMs and GMRFs are explored based on approximations of GMRFs by GGMs, and approximations of GGMs by GMRFs. Two new metrics of approximation are proposed : (i) the Kullback-Leibler discrepancy of spectral densities and (ii) the chi-squared distance between spectral densities. The distances between the spectral density functions of GGMs and GMRFs measured by these metrics are minimized to obtain the approximations of GGMs and GMRFs. The proposed methodologies are validated through several empirical studies. We compare the performance of our approach to other methods based on covariance functions, in terms of the average mean squared prediction error and also the computational time. A spatial analysis of a dataset on PM_2.5 collected in California is presented to illustrate the proposed method.     

Robust utility maximization with limited downside risk in incomplete markets

In this article we consider the portfolio selection problem of an agent with robust preferences in the sense of Gilboa and Schmeidler [Itzhak Gilboa, David Schmeidler, Maxmin expected utility with non-unique prior, Journal of Mathematical Economics 18 (1989) 141–153] in an incomplete market. Downside risk is constrained by a robust version of utility-based shortfall risk. We derive an explicit representation of the optimal terminal wealth in terms of certain worst case measures which can be characterized as minimizers of a dual problem. This dual problem involves a three-dimensional analogue of f$f$-divergences which generalize the notion of relative entropy.     

Measuring the coupled risks: A copula-based CVaR model

Integrated risk management for financial institutions requires an approach for aggregating risk types (such as market and credit) whose distributional shapes vary considerably. The financial institutions often ignore risks’ coupling influence so as to underestimate the financial risks. We constructed a copula-based Conditional Value-at-Risk (CVaR) model for market and credit risks. This technique allows us to incorporate realistic marginal distributions that capture essential empirical features of these risks, such as skewness and fat-tails while allowing for a rich dependence structure. Finally, the numerical simulation method is used to implement the model. Our results indicate that the coupled risks for the listed company’s stock maybe are undervalued if credit risk is ignored, especially for the listed company with bad credit quality.     

Relational models for contingency tables

The paper considers general multiplicative models for complete and incomplete contingency tables that generalize log-linear and several other models and are entirely coordinate free. Sufficient conditions for the existence of maximum likelihood estimates under these models are given, and it is shown that the usual equivalence between multinomial and Poisson likelihoods holds if and only if an overall effect is present in the model. If such an effect is not assumed, the model becomes a curved exponential family and a related mixed parameterization is given that relies on non-homogeneous odds ratios. Several examples are presented to illustrate the properties and use of such models.     

On Hyers-Ulam stability for a class of functional equations

Summary In this paper we prove some stability theorems for functional equations of the formg[F(x, y)]=H[g(x), g(y), x, y]. As special cases we obtain well known results for Cauchy and Jensen equations and for functional equations in a single variable. Work supported by M.U.R.S.T. Research funds (60%).     

Further results on stability of normal regions for linear homogeneous functional equations

In this paper we consider the stability of normal regions on the plane, determined by continuous solutions of the linear homogeneous functional inequality in the case where continuous solutions of the corresponding linear functional equation do not depend continuously on initial conditions.     

The index of an isolated critical point for a class of non-differentiable elliptic operators in reflexive Banach spaces

This paper is motivated by the problem of computing the index of an isolated critical point for a certain class of nonlinear elliptic operators. Such a class includes operators A which are defined by

     

Estimation in generalized linear models for functional data via penalized likelihood

We analyze in a regression setting the link between a scalar response and a functional predictor by means of a Functional Generalized Linear Model. We first give a theoretical framework and then discuss identifiability of the model. The functional coefficient of the model is estimated via penalized likelihood with spline approximation. The L² rate of convergence of this estimator is given under smoothness assumption on the functional coefficient. Heuristic arguments show how these rates may be improved for some particular frameworks.     

Maps leaving functional values of operator products invariant

Let A be a standard operator algebra on a complex Hilbert space H of dimension greater than 2. By invariants of certain functional values of operator products, we characterize some surjective maps on A. Furthermore, several kinds of general preserver problems on standard operator algebras are solved when we take respectively the functional as, for example, k-numerical radius (k?1), operator norm, Ky Fan k-norm, Schatten p-norm (1?p<∞), and so on.     

Sample covariance shrinkage for high dimensional dependent data

For high dimensional data sets the sample covariance matrix is usually unbiased but noisy if the sample is not large enough. Shrinking the sample covariance towards a constrained, low dimensional estimator can be used to mitigate the sample variability. By doing so, we introduce bias, but reduce variance. In this paper, we give details on feasible optimal shrinkage allowing for time series dependent observations.     

A generalized Mahalanobis distance for mixed data

A distance for mixed nominal, ordinal and continuous data is developed by applying the Kullback-Leibler divergence to the general mixed-data model, an extension of the general location model that allows for ordinal variables to be incorporated in the model. The distance obtained can be considered as a generalization of the Mahalanobis distance to data with a mixture of nominal, ordinal and continuous variables. Moreover, it includes as special cases previous Mahalanobis-type distances developed by Bedrick et al. (Biometrics 56 (2000) 394) and Bar-Hen and Daudin (J. Multivariate Anal. 53 (1995) 332). Asymptotic results regarding the maximum likelihood estimator of the distance are discussed. The results of a simulation study on the level and power of the tests are reported and a real-data example illustrates the method.     

REML estimation for binary data in GLMMs

The restricted maximum likelihood (REML) procedure is useful for inferences about variance components in mixed linear models. However, its extension to hierarchical generalized linear models (HGLMs) is often hampered by analytically intractable integrals. Numerical integration such as Gauss-Hermite quadrature (GHQ) is generally not recommended when the dimensionality of the integral is high. With binary data various extensions of the REML method have been suggested, but they have had unsatisfactory biases in estimation. In this paper we propose a statistically and computationally efficient REML procedure for the analysis of binary data, which is applicable over a wide class of models and design structures. We propose a bias-correction method for models such as binary matched pairs and discuss how the REML estimating equations for mixed linear models can be modified to implement more general models.     

On the equivariant formality of Kähler manifolds with torus group actions

We consider simply connected compact Kähler manifolds which have a holomorphic action of a torus group. We use the existing equivariant models for rational homotopy to show that these spaces satisfy an equivariant formality condition over the complex numbers.