A Functional Estimate of Covariation
本文用有限元方法构造函数型Choleski因子来估计协方差函数，得到的逆协方差矩阵估计具有强带状结构，在模拟和Fels女性生长数据平滑中表现优于经典方法。
The analysis of functional data calls for a bivariate functional covariance function σ(s, t) that may be evaluated at any discrete set of points to define a variance-covariance matrix Σ. This article uses finite element methodology to construct a representation of a functional Choleski factor λ(w, s) to define σ(s, t) = ∫λ(w, s)λ(w, t) dw. An estimate of Σ-1 is especially important for applications and, where the eigenstructure of the covariance permits, this is readily available since the resulting Σ is almost always positive definite. A simulation study compares the performance of estimates of Σ and Σ-1 to those from the classic covariance matrix estimate and an estimate using glasso package in R. The method’s capability of constraining estimates of Σ-1 to be strongly band-structured resulted in superior estimates. The real data application is to the smoothing of the Fels female growth data where σ(s, t) estimates the residual covariance structure in the presence of sampling points varying from one case to another. Supplementary materials are available online.