Geostatistical modeling of positive-definite matrices: An application to diffusion tensor imaging
- 1 June 2022
- journal article
- research article
- Published by Oxford University Press (OUP) in Biometrics
- Vol. 78 (2), 548-559
- https://doi.org/10.1111/biom.13445
Abstract
Geostatistical modeling for continuous point-referenced data has extensively been applied to neuroimaging because it produces efficient and valid statistical inference. However, diffusion tensor imaging (DTI), a neuroimaging technique characterizing the brain's anatomical structure, produces a positive-definite (p.d.) matrix for each voxel. Currently, only a few geostatistical models for p.d. matrices have been proposed because introducing spatial dependence among p.d. matrices properly is challenging. In this paper, we use the spatial Wishart process, a spatial stochastic process (random field), where each p.d. matrix-variate random variable marginally follows a Wishart distribution, and spatial dependence between random matrices is induced by latent Gaussian processes. This process is valid on an uncountable collection of spatial locations and is almost-surely continuous, leading to a reasonable way of modeling spatial dependence. Motivated by a DTI data set of cocaine users, we propose a spatial matrix-variate regression model based on the spatial Wishart process. A problematic issue is that the spatial Wishart process has no closed-form density function. Hence, we propose an approximation method to obtain a feasible Cholesky decomposition model, which we show to be asymptotically equivalent to the spatial Wishart process model. A local likelihood approximation method is also applied to achieve fast computation. The simulation studies and real data application demonstrate that the Cholesky decomposition process model produces reliable inference and improved performance, compared to other methods.Funding Information
- Foundation for the National Institutes of Health (P30CA016059, R01DE024984, U54DA038999)
This publication has 38 references indexed in Scilit:
- Mapping the Voxel-Wise Effective Connectome in Resting State fMRIPLOS ONE, 2013
- A hitchhiker's guide to diffusion tensor imagingFrontiers in Neuroscience, 2013
- Diffusion Tensor Imaging and Decision Making in Cocaine DependencePLOS ONE, 2010
- Bayesian Treed Gaussian Process Models With an Application to Computer ModelingJournal of the American Statistical Association, 2008
- Spatial statistical models that use flow and stream distanceEnvironmental and Ecological Statistics, 2006
- Orthogonal tensor invariants and the analysis of diffusion tensor magnetic resonance imagesMagnetic Resonance in Medicine, 2005
- Drugs, the Brain, and Behavior: the Pharmacology of Abuse and Dependence. Edited by John Brick and Carlton Erickson. 186 pp., illustrated. New York: Haworth Medical Press.International Journal of Neuropsychopharmacology, 2000
- Semiparametric Regression for Repeated Outcomes with Nonignorable NonresponseJournal of the American Statistical Association, 1998
- Estimating a Population of Parameter Values Using Bayes and Empirical Bayes MethodsJournal of the American Statistical Association, 1984
- The Folded Normal DistributionTechnometrics, 1961