Abstract
Symmetric positive definite matrix data commonly appear in computer vision and medical imaging, such as diffusion tensor imaging. The aim of this paper is to develop a nonparametric estimation method for a symmetric positive definite matrix regression function given covariates. By obtaining a suitable parametrization based on the Cholesky decomposition, we make it possible to apply univariate smoothing methods to the matrix regression problem. The parametrization also guarantees that the proposed estimator is symmetric positive definite over the entire domain. We adopt the Wishart log-likelihood and a smoothing technique using the basis methodology to define our estimator. The rate of convergence of the proposed estimator is obtained under some regularity conditions. Simulations are performed to investigate the finite sample properties of our proposed method using natural splines. Moreover, we present the results of an analysis of real diffusion tensor imaging data where the estimated fractional anisotropy is provided using 3 × 3 symmetric positive definite matrices measured at consecutive positions along a fiber tract in the brain of subjects.
Original language | English |
---|---|
Pages (from-to) | 795-817 |
Number of pages | 23 |
Journal | Journal of the Korean Statistical Society |
Volume | 50 |
Issue number | 3 |
DOIs | |
State | Published - Sep 2021 |
Keywords
- Cholesky factorization
- Diffusion tensor imaging
- Matrix regression function
- Minimaxity
- Stein loss
- Wishart distribution