1
|
Miranda MF. A canonical polyadic tensor basis for fast Bayesian estimation of multi-subject brain activation patterns. Front Neuroinform 2024; 18:1399391. [PMID: 39188665 PMCID: PMC11345152 DOI: 10.3389/fninf.2024.1399391] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2024] [Accepted: 07/29/2024] [Indexed: 08/28/2024] Open
Abstract
Task-evoked functional magnetic resonance imaging studies, such as the Human Connectome Project (HCP), are a powerful tool for exploring how brain activity is influenced by cognitive tasks like memory retention, decision-making, and language processing. A fast Bayesian function-on-scalar model is proposed for estimating population-level activation maps linked to the working memory task. The model is based on the canonical polyadic (CP) tensor decomposition of coefficient maps obtained for each subject. This decomposition effectively yields a tensor basis capable of extracting both common features and subject-specific features from the coefficient maps. These subject-specific features, in turn, are modeled as a function of covariates of interest using a Bayesian model that accounts for the correlation of the CP-extracted features. The dimensionality reduction achieved with the tensor basis allows for a fast MCMC estimation of population-level activation maps. This model is applied to one hundred unrelated subjects from the HCP dataset, yielding significant insights into brain signatures associated with working memory.
Collapse
Affiliation(s)
- Michelle F. Miranda
- Department of Mathematics and Statistics, University of Victoria, Victoria, BC, Canada
| |
Collapse
|
2
|
Zhou F, He K, Wang K, Xu Y, Ni Y. Functional Bayesian networks for discovering causality from multivariate functional data. Biometrics 2023; 79:3279-3293. [PMID: 37635676 PMCID: PMC10840881 DOI: 10.1111/biom.13922] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2022] [Accepted: 08/10/2023] [Indexed: 08/29/2023]
Abstract
Multivariate functional data arise in a wide range of applications. One fundamental task is to understand the causal relationships among these functional objects of interest. In this paper, we develop a novel Bayesian network (BN) model for multivariate functional data where conditional independencies and causal structure are encoded by a directed acyclic graph. Specifically, we allow the functional objects to deviate from Gaussian processes, which is the key to unique causal structure identification even when the functions are measured with noises. A fully Bayesian framework is designed to infer the functional BN model with natural uncertainty quantification through posterior summaries. Simulation studies and real data examples demonstrate the practical utility of the proposed model.
Collapse
Affiliation(s)
- Fangting Zhou
- Department of Statistics, Texas A&M University, College Station, Texas, USA
- Center for Applied Statistics, Institute of Statistics and Big Data, Renmin University of China, Beijing, China
| | - Kejun He
- Center for Applied Statistics, Institute of Statistics and Big Data, Renmin University of China, Beijing, China
| | - Kunbo Wang
- Department of Applied Mathematics and Statistics, Johns Hopkins University, Baltimore, Maryland, USA
| | - Yanxun Xu
- Department of Applied Mathematics and Statistics, Johns Hopkins University, Baltimore, Maryland, USA
| | - Yang Ni
- Department of Statistics, Texas A&M University, College Station, Texas, USA
| |
Collapse
|
3
|
Ren B, Hwang WT. Modeling post-holiday surge in COVID-19 cases in Pennsylvania counties. PLoS One 2022; 17:e0279371. [PMID: 36534663 PMCID: PMC9762594 DOI: 10.1371/journal.pone.0279371] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2022] [Accepted: 12/06/2022] [Indexed: 12/23/2022] Open
Abstract
COVID-19 arrived in the United States in early 2020, with cases quickly being reported in many states including Pennsylvania. Many statistical models have been proposed to understand the trends of the COVID-19 pandemic and factors associated with increasing cases. While Poisson regression is a natural choice to model case counts, this approach fails to account for correlation due to spatial locations. Being a contagious disease and often spreading through community infections, the number of COVID-19 cases are inevitably spatially correlated as locations neighboring counties with a high COVID-19 case count are more likely to have a high case count. In this analysis, we combine generalized estimating equations (GEEs) for Poisson regression, a popular method for analyzing correlated data, with a semivariogram to model daily COVID-19 case counts in 67 Pennsylvania counties between March 20, 2020 to January 23, 2021 in order to study infection dynamics during the beginning of the pandemic. We use a semivariogram that describes the spatial correlation as a function of the distance between two counties as the working correlation. We further incorporate a zero-inflated model in our spatial GEE to accommodate excess zeros in reported cases due to logistical challenges associated with disease monitoring. By modeling time-varying holiday covariates, we estimated the effect of holiday timing on case count. Our analysis showed that the incidence rate ratio was significantly greater than one, 6-8 days after a holiday suggesting a surge in COVID-19 cases approximately one week after a holiday.
Collapse
Affiliation(s)
- Benny Ren
- Department of Biostatistics, Epidemiology, and Informatics, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
- * E-mail:
| | - Wei-Ting Hwang
- Department of Biostatistics, Epidemiology, and Informatics, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
| |
Collapse
|
4
|
Cui E, Li R, Crainiceanu CM, Xiao L. Fast Multilevel Functional Principal Component Analysis. J Comput Graph Stat 2022; 32:366-377. [PMID: 37313008 PMCID: PMC10260118 DOI: 10.1080/10618600.2022.2115500] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2021] [Accepted: 08/06/2022] [Indexed: 10/15/2022]
Abstract
We introduce fast multilevel functional principal component analysis (fast MFPCA), which scales up to high dimensional functional data measured at multiple visits. The new approach is orders of magnitude faster than and achieves comparable estimation accuracy with the original MFPCA (Di et al., 2009). Methods are motivated by the National Health and Nutritional Examination Survey (NHANES), which contains minute-level physical activity information of more than 10000 participants over multiple days and 1440 observations per day. While MFPCA takes more than five days to analyze these data, fast MFPCA takes less than five minutes. A theoretical study of the proposed method is also provided. The associated function mfpca.face() is available in the R package refund.
Collapse
Affiliation(s)
- Erjia Cui
- Department of Biostatistics, Johns Hopkins Bloomberg School of Public Health, 615 N. Wolfe Street, Baltimore, MD 21205
| | - Ruonan Li
- Department of Statistics, North Carolina State University, 2311 Stinson Dr, Raleigh, NC 27607
| | - Ciprian M. Crainiceanu
- Department of Biostatistics, Johns Hopkins Bloomberg School of Public Health, 615 N. Wolfe Street, Baltimore, MD 21205
| | - Luo Xiao
- Department of Statistics, North Carolina State University, 2311 Stinson Dr, Raleigh, NC 27607
| |
Collapse
|
5
|
Huo S, Morris JS, Zhu H. Ultra-Fast Approximate Inference Using Variational Functional Mixed Models. J Comput Graph Stat 2022; 32:353-365. [PMID: 37608921 PMCID: PMC10441618 DOI: 10.1080/10618600.2022.2107532] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2021] [Accepted: 07/23/2022] [Indexed: 10/16/2022]
Abstract
While Bayesian functional mixed models have been shown effective to model functional data with various complex structures, their application to extremely high-dimensional data is limited due to computational challenges involved in posterior sampling. We introduce a new computational framework that enables ultra-fast approximate inference for high-dimensional data in functional form. This framework adopts parsimonious basis to represent functional observations, which facilitates efficient compression and parallel computing in basis space. Instead of performing expensive Markov chain Monte Carlo sampling, we approximate the posterior distribution using variational Bayes and adopt a fast iterative algorithm to estimate parameters of the approximate distribution. Our approach facilitates a fast multiple testing procedure in basis space, which can be used to identify significant local regions that reflect differences across groups of samples. We perform two simulation studies to assess the performance of approximate inference, and demonstrate applications of the proposed approach by using a proteomic mass spectrometry dataset and a brain imaging dataset. Supplementary materials are available online.
Collapse
Affiliation(s)
| | - Jeffrey S Morris
- Department of Biostatistics, Epidemiology and Informatics, Department of Statistics, University of Pennsylvania
| | | |
Collapse
|
6
|
Park Y, Li B, Li Y. Crop Yield Prediction Using Bayesian Spatially Varying Coefficient Models with Functional Predictors. J Am Stat Assoc 2022. [DOI: 10.1080/01621459.2022.2123333] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Affiliation(s)
- Yeonjoo Park
- Management Science and Statistics, University of Texas at San Antonio
| | - Bo Li
- Department of Statistics, University of Illinois at Urbana-Champaign
| | - Yehua Li
- Department of Statistics, University of California at Riverside
| |
Collapse
|
7
|
Li R, Xiao L, Smirnova E, Cui E, Leroux A, Crainiceanu CM. Fixed-effects inference and tests of correlation for longitudinal functional data. Stat Med 2022; 41:3349-3364. [PMID: 35491388 PMCID: PMC9283332 DOI: 10.1002/sim.9421] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2021] [Revised: 01/30/2022] [Accepted: 03/05/2022] [Indexed: 11/19/2022]
Abstract
We propose an inferential framework for fixed effects in longitudinal functional models and introduce tests for the correlation structures induced by the longitudinal sampling procedure. The framework provides a natural extension of standard longitudinal correlation models for scalar observations to functional observations. Using simulation studies, we compare fixed effects estimation under correctly and incorrectly specified correlation structures and also test the longitudinal correlation structure. Finally, we apply the proposed methods to a longitudinal functional dataset on physical activity. The computer code for the proposed method is available at https://github.com/rli20ST758/FILF.
Collapse
Affiliation(s)
- Ruonan Li
- Department of StatisticsNorth Carolina State UniversityRaleighNorth CarolinaUSA
| | - Luo Xiao
- Department of StatisticsNorth Carolina State UniversityRaleighNorth CarolinaUSA
| | - Ekaterina Smirnova
- Department of BiostatisticsVirginia Commonwealth UniversityRichmondVirginiaUSA
| | - Erjia Cui
- Department of BiostatisticsJohns Hopkins Bloomberg School of Public HealthBaltimoreMarylandUSA
| | - Andrew Leroux
- Department of Biostatistics and InformaticsColorado School of Public HealthAuroraColoradoUSA
| | - Ciprian M. Crainiceanu
- Department of BiostatisticsJohns Hopkins Bloomberg School of Public HealthBaltimoreMarylandUSA
| |
Collapse
|
8
|
Li Y, Qiu Y, Xu Y. From multivariate to functional data analysis: fundamentals, recent developments, and emerging areas. J MULTIVARIATE ANAL 2022; 188:104806. [PMID: 39040141 PMCID: PMC11261241 DOI: 10.1016/j.jmva.2021.104806] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Functional data analysis (FDA), which is a branch of statistics on modeling infinite dimensional random vectors resided in functional spaces, has become a major research area for Journal of Multivariate Analysis. We review some fundamental concepts of FDA, their origins and connections from multivariate analysis, and some of its recent developments, including multi-level functional data analysis, high-dimensional functional regression, and dependent functional data analysis. We also discuss the impact of these new methodology developments on genetics, plant science, wearable device data analysis, image data analysis, and business analytics. Two real data examples are provided to motivate our discussions.
Collapse
Affiliation(s)
- Yehua Li
- University of California - Riverside, Riverside, CA 92521, USA
| | - Yumou Qiu
- Iowa State University, Ames, IA 50011, USA
| | - Yuhang Xu
- Bowling Green State University, Bowling Green, OH 43403, USA
| |
Collapse
|
9
|
Yarger D, Stoev S, Hsing T. A functional-data approach to the Argo data. Ann Appl Stat 2022. [DOI: 10.1214/21-aoas1477] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Affiliation(s)
- Drew Yarger
- Department of Statistics, University of Michigan
| | | | - Tailen Hsing
- Department of Statistics, University of Michigan
| |
Collapse
|
10
|
Yu D, He Z. Digital twin-driven intelligence disaster prevention and mitigation for infrastructure: advances, challenges, and opportunities. NATURAL HAZARDS (DORDRECHT, NETHERLANDS) 2022; 112:1-36. [PMID: 35125651 PMCID: PMC8801275 DOI: 10.1007/s11069-021-05190-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/06/2021] [Accepted: 12/19/2021] [Indexed: 06/14/2023]
Abstract
Natural hazards, which have the potential to cause catastrophic damage and loss to infrastructure, have increased significantly in recent decades. Thus, the construction demand for disaster prevention and mitigation for infrastructure (DPMI) systems is increasing. Many studies have applied intelligence technologies to solve key aspects of infrastructure, such as design, construction, disaster prevention and mitigation, and rescue and recovery; however, systematic construction is still lacking. Digital twin (DT) is one of the most promising technologies for multi-stage management which has great potential to solve the above challenges. This paper initially puts forward a scientific concept, in which DT drives the construction of intelligent disaster prevention and mitigation for infrastructure (IDPMI) systematically. To begin with, a scientific review of DT and IDPMI is performed, where the development of DT is summarized and a DT-based life cycle of infrastructures is defined. In addition, the intelligence technologies used in disaster management are key reviewed and their relative merits are illustrated. Furthermore, the development and technical feasibility of DT-driven IDPMI are illustrated by reviewing the relevant practice of DT in infrastructure. In conclusion, a scientific framework of DT-IDPMI is programmed, which not only provides some guidance for the deep integration between DT and IDPMI but also identifies the challenges that inspire the professional community to advance these techniques to address them in future research.
Collapse
Affiliation(s)
- Dianyou Yu
- Department of Civil Engineering, Dalian University of Technology, Dalian, China
| | - Zheng He
- Department of Civil Engineering, Dalian University of Technology, Dalian, China
- State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian, China
| |
Collapse
|
11
|
Cui E, Leroux A, Smirnova E, Crainiceanu CM. Fast Univariate Inference for Longitudinal Functional Models. J Comput Graph Stat 2022; 31:219-230. [PMID: 35712524 PMCID: PMC9197085 DOI: 10.1080/10618600.2021.1950006] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
We propose fast univariate inferential approaches for longitudinal Gaussian and non-Gaussian functional data. The approach consists of three steps: (1) fit massively univariate pointwise mixed effects models; (2) apply any smoother along the functional domain; and (3) obtain joint confidence bands using analytic approaches for Gaussian data or a bootstrap of study participants for non-Gaussian data. Methods are motivated by two applications: (1) Diffusion Tensor Imaging (DTI) measured at multiple visits along the corpus callosum of multiple sclerosis (MS) patients; and (2) physical activity data measured by body-worn accelerometers for multiple days. An extensive simulation study indicates that model fitting and inference are accurate and much faster than existing approaches. Moreover, the proposed approach was the only one that was computationally feasible for the physical activity data application. Methods are accompanied by R software, though the method is "read-and-use", as it can be implemented by any analyst who is familiar with mixed effects model software.
Collapse
Affiliation(s)
- Erjia Cui
- Department of Biostatistics, Johns Hopkins Bloomberg School of Public Health, USA
| | - Andrew Leroux
- Department of Biostatistics and Informatics, University of Colorado, USA
| | | | | |
Collapse
|
12
|
Li Y, Nguyen DV, Kürüm E, Rhee CM, Banerjee S, Şentürk D. Multilevel Varying Coefficient Spatiotemporal Model. Stat (Int Stat Inst) 2021; 11. [DOI: 10.1002/sta4.438] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Yihao Li
- Department of Biostatistics University of California Los Angeles CA USA
| | - Danh V. Nguyen
- Department of Medicine University of California Irvine Orange CA USA
| | - Esra Kürüm
- Department of Statistics University of California Riverside CA USA
| | - Connie M. Rhee
- Department of Medicine University of California Irvine Orange CA USA
- Harold Simmons Center for Chronic Disease Research and Epidemiology University of California Irvine School of Medicine Orange CA USA
| | - Sudipto Banerjee
- Department of Biostatistics University of California Los Angeles CA USA
| | - Damla Şentürk
- Department of Biostatistics University of California Los Angeles CA USA
| |
Collapse
|
13
|
Zhao J, Zhang Q, Fuentes M, Qian Y, Ma L, Moeller G. A Spatio-temporal Model for Detecting the Effect of Cocaine Use Disorder on Functional Connectivity. SPATIAL STATISTICS 2021; 45:100530. [PMID: 34804784 PMCID: PMC8598107 DOI: 10.1016/j.spasta.2021.100530] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Drug addiction can lead to many health-related problems and social concerns. Researchers are interested in the association between long-term drug usage and abnormal functional connectivity. Functional connectivity obtained from functional magnetic resonance imaging data promotes a variety of fundamental understandings in such association. Due to the complex correlation structure and large dimensionality, the modeling and analysis of the functional connectivity from neuroimage are challenging. By proposing a spatio-temporal model for multi-subject neuroimage data, we incorporate voxel-level spatio-temporal dependencies of whole-brain measurements to improve the accuracy of statistical inference. To tackle large-scale spatio-temporal neuroimage data, we develop a computational efficient algorithm to estimate the parameters. Our method is used to first identify functional connectivity, and then detect the effect of cocaine use disorder (CUD) on functional connectivity between different brain regions. The functional connectivity identified by our spatio-temporal model matches existing studies on brain networks, and further indicates that CUD may alter the functional connectivity in the medial orbitofrontal cortex subregions and the supplementary motor areas.
Collapse
Affiliation(s)
- Jifang Zhao
- Department of Statistical Sciences and Operations Research, Virginia Commonwealth University, Richmond, VA
| | - Qiong Zhang
- School of Mathematical and Statistical Sciences, Clemson University, Clemson, SC
| | - Montserrat Fuentes
- Office of the Executive Vice President and Provost, University of Iowa, Iowa City, IA
| | - Yanjun Qian
- Department of Statistical Sciences and Operations Research, Virginia Commonwealth University, Richmond, VA
| | - Liangsuo Ma
- Institute for Drug and Alcohol Studies, Virginia Commonwealth University, Richmond, VA
| | - Gerard Moeller
- Institute for Drug and Alcohol Studies, Virginia Commonwealth University, Richmond, VA
| |
Collapse
|
14
|
Zhang H, Li Y. Unified Principal Component Analysis for Sparse and Dense Functional Data under Spatial Dependency. JOURNAL OF BUSINESS & ECONOMIC STATISTICS : A PUBLICATION OF THE AMERICAN STATISTICAL ASSOCIATION 2021; 40:1523-1537. [PMID: 36582252 PMCID: PMC9793858 DOI: 10.1080/07350015.2021.1938085] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
We consider spatially dependent functional data collected under a geostatistics setting, where locations are sampled from a spatial point process. The functional response is the sum of a spatially dependent functional effect and a spatially independent functional nugget effect. Observations on each function are made on discrete time points and contaminated with measurement errors. Under the assumption of spatial stationarity and isotropy, we propose a tensor product spline estimator for the spatio-temporal covariance function. When a coregionalization covariance structure is further assumed, we propose a new functional principal component analysis method that borrows information from neighboring functions. The proposed method also generates nonparametric estimators for the spatial covariance functions, which can be used for functional kriging. Under a unified framework for sparse and dense functional data, infill and increasing domain asymptotic paradigms, we develop the asymptotic convergence rates for the proposed estimators. Advantages of the proposed approach are demonstrated through simulation studies and two real data applications representing sparse and dense functional data, respectively.
Collapse
|
15
|
Li Y, Nguyen DV, Banerjee S, Rhee CM, Kalantar-Zadeh K, Kürüm E, Şentürk D. Multilevel modeling of spatially nested functional data: Spatiotemporal patterns of hospitalization rates in the US dialysis population. Stat Med 2021; 40:3937-3952. [PMID: 33902165 DOI: 10.1002/sim.9007] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2020] [Revised: 02/15/2021] [Accepted: 04/08/2021] [Indexed: 11/12/2022]
Abstract
End-stage renal disease patients on dialysis experience frequent hospitalizations. In addition to known temporal patterns of hospitalizations over the life span on dialysis, where poor outcomes are typically exacerbated during the first year on dialysis, variations in hospitalizations among dialysis facilities across the US contribute to spatial variation. Utilizing national data from the United States Renal Data System (USRDS), we propose a novel multilevel spatiotemporal functional model to study spatiotemporal patterns of hospitalization rates among dialysis facilities. Hospitalization rates of dialysis facilities are considered as spatially nested functional data (FD) with longitudinal hospitalizations nested in dialysis facilities and dialysis facilities nested in geographic regions. A multilevel Karhunen-Loéve expansion is utilized to model the two-level (facility and region) FD, where spatial correlations are induced among region-specific principal component scores accounting for regional variation. A new efficient algorithm based on functional principal component analysis and Markov Chain Monte Carlo is proposed for estimation and inference. We report a novel application using USRDS data to characterize spatiotemporal patterns of hospitalization rates for over 400 health service areas across the US and over the posttransition time on dialysis. Finite sample performance of the proposed method is studied through simulations.
Collapse
Affiliation(s)
- Yihao Li
- Department of Biostatistics, University of California, Los Angeles, California
| | - Danh V Nguyen
- Department of Medicine, UC Irvine School of Medicine, Orange, California
| | - Sudipto Banerjee
- Department of Biostatistics, University of California, Los Angeles, California
| | - Connie M Rhee
- Department of Medicine, UC Irvine School of Medicine, Orange, California
| | | | - Esra Kürüm
- Department of Statistics, University of California, Riverside, California
| | - Damla Şentürk
- Department of Biostatistics, University of California, Los Angeles, California
| |
Collapse
|
16
|
Shamshoian J, Şentürk D, Jeste S, Telesca D. Bayesian analysis of longitudinal and multidimensional functional data. Biostatistics 2020; 23:558-573. [PMID: 33017019 DOI: 10.1093/biostatistics/kxaa041] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2019] [Revised: 08/31/2020] [Accepted: 09/04/2020] [Indexed: 11/13/2022] Open
Abstract
Multi-dimensional functional data arises in numerous modern scientific experimental and observational studies. In this article, we focus on longitudinal functional data, a structured form of multidimensional functional data. Operating within a longitudinal functional framework we aim to capture low dimensional interpretable features. We propose a computationally efficient nonparametric Bayesian method to simultaneously smooth observed data, estimate conditional functional means and functional covariance surfaces. Statistical inference is based on Monte Carlo samples from the posterior measure through adaptive blocked Gibbs sampling. Several operative characteristics associated with the proposed modeling framework are assessed comparatively in a simulated environment. We illustrate the application of our work in two case studies. The first case study involves age-specific fertility collected over time for various countries. The second case study is an implicit learning experiment in children with autism spectrum disorder.
Collapse
Affiliation(s)
- John Shamshoian
- Department of Biostatistics, University of California, Los Angeles, CA, USA
| | - Damla Şentürk
- Department of Biostatistics, University of California, Los Angeles, CA, USA
| | - Shafali Jeste
- Department of Psychiatry and Biobehavioral Sciences, University of California, Los Angeles, CA, USA
| | - Donatello Telesca
- Department of Biostatistics, University of California, Los Angeles, CA, USA
| |
Collapse
|
17
|
Martínez-Hernández I, Genton MG. Recent developments in complex and spatially correlated functional data. BRAZ J PROBAB STAT 2020. [DOI: 10.1214/20-bjps466] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
18
|
Huang T, Saporta G, Wang H, Wang S. A robust spatial autoregressive scalar-on-function regression with t-distribution. ADV DATA ANAL CLASSI 2020. [DOI: 10.1007/s11634-020-00384-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
19
|
Yang H, Baladandayuthapani V, Rao AUK, Morris JS. Quantile Function on Scalar Regression Analysis for Distributional Data. J Am Stat Assoc 2019; 115:90-106. [PMID: 32981991 PMCID: PMC7517594 DOI: 10.1080/01621459.2019.1609969] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2018] [Revised: 03/08/2019] [Accepted: 04/07/2019] [Indexed: 02/05/2023]
Abstract
Radiomics involves the study of tumor images to identify quantitative markers explaining cancer heterogeneity. The predominant approach is to extract hundreds to thousands of image features, including histogram features comprised of summaries of the marginal distribution of pixel intensities, which leads to multiple testing problems and can miss out on insights not contained in the selected features. In this paper, we present methods to model the entire marginal distribution of pixel intensities via the quantile function as functional data, regressed on a set of demographic, clinical, and genetic predictors to investigate their effects of imaging-based cancer heterogeneity. We call this approach quantile functional regression, regressing subject-specific marginal distributions across repeated measurements on a set of covariates, allowing us to assess which covariates are associated with the distribution in a global sense, as well as to identify distributional features characterizing these differences, including mean, variance, skewness, heavy-tailedness, and various upper and lower quantiles. To account for smoothness in the quantile functions, account for intrafunctional correlation, and gain statistical power, we introduce custom basis functions we call quantlets that are sparse, regularized, near-lossless, and empirically defined, adapting to the features of a given data set and containing a Gaussian subspace so non-Gaussianness can be assessed. We fit this model using a Bayesian framework that uses nonlinear shrinkage of quantlet coefficients to regularize the functional regression coefficients and provides fully Bayesian inference after fitting a Markov chain Monte Carlo. We demonstrate the benefit of the basis space modeling through simulation studies, and apply the method to Magnetic resonance imaging (MRI) based radiomic dataset from Glioblastoma Multiforme to relate imaging-based quantile functions to various demographic, clinical, and genetic predictors, finding specific differences in tumor pixel intensity distribution between males and females and between tumors with and without DDIT3 mutations.
Collapse
Affiliation(s)
- Hojin Yang
- Department of Biostatistics, The University of Texas MD Anderson Cancer Center, Houston, TX 77030
| | | | - Arvind U K Rao
- Department of Bioinformatics and Computational Biology, The University of Texas MD Anderson Cancer Center, Houston, TX 77030
| | - Jeffrey S Morris
- Department of Biostatistics, The University of Texas MD Anderson Cancer Center, Houston, TX 77030
| |
Collapse
|
20
|
Wang Y, Hu J, Do KA, Hobbs BP. An Efficient Nonparametric Estimate for Spatially Correlated Functional Data. STATISTICS IN BIOSCIENCES 2019. [DOI: 10.1007/s12561-019-09233-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
21
|
Fontanella L, Ippoliti L, Valentini P. Predictive functional ANOVA models for longitudinal analysis of mandibular shape changes. Biom J 2019; 61:918-933. [PMID: 30865334 DOI: 10.1002/bimj.201800228] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2018] [Revised: 12/17/2018] [Accepted: 12/19/2018] [Indexed: 11/10/2022]
Abstract
In this paper, we introduce a Bayesian statistical model for the analysis of functional data observed at several time points. Examples of such data include the Michigan growth study where we wish to characterize the shape changes of human mandible profiles. The form of the mandible is often used by clinicians as an aid in predicting the mandibular growth. However, whereas many studies have demonstrated the changes in size that may occur during the period of pubertal growth spurt, shape changes have been less well investigated. Considering a group of subjects presenting normal occlusion, in this paper we thus describe a Bayesian functional ANOVA model that provides information about where and when the shape changes of the mandible occur during different stages of development. The model is developed by defining the notion of predictive process models for Gaussian process (GP) distributions used as priors over the random functional effects. We show that the predictive approach is computationally appealing and that it is useful to analyze multivariate functional data with unequally spaced observations that differ among subjects and times. Graphical posterior summaries show that our model is able to provide a biological interpretation of the morphometric findings and that they comprehensively describe the shape changes of the human mandible profiles. Compared with classical cephalometric analysis, this paper represents a significant methodological advance for the study of mandibular shape changes in two dimensions.
Collapse
Affiliation(s)
- Lara Fontanella
- Department of Legal and Social Sciences, University G. d'Annunzio, Chieti-Pescara, Italy
| | - Luigi Ippoliti
- Department of Economics, University G. d'Annunzio, Chieti-Pescara, Italy
| | - Pasquale Valentini
- Department of Economics, University G. d'Annunzio, Chieti-Pescara, Italy
| |
Collapse
|
22
|
Zhu H, Versace F, Cinciripini PM, Rausch P, Morris JS. Robust and Gaussian spatial functional regression models for analysis of event-related potentials. Neuroimage 2018; 181:501-512. [PMID: 30057352 DOI: 10.1016/j.neuroimage.2018.07.006] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2018] [Revised: 06/01/2018] [Accepted: 07/03/2018] [Indexed: 10/28/2022] Open
Abstract
Event-related potentials (ERPs) summarize electrophysiological brain response to specific stimuli. They can be considered as correlated functions of time with both spatial correlation across electrodes and nested correlations within subjects. Commonly used analytical methods for ERPs often focus on pre-determined extracted components and/or ignore the correlation among electrodes or subjects, which can miss important insights, and tend to be sensitive to outlying subjects, time points or electrodes. Motivated by ERP data in a smoking cessation study, we introduce a Bayesian spatial functional regression framework that models the entire ERPs as spatially correlated functional responses and the stimulus types as covariates. This novel framework relies on mixed models to characterize the effects of stimuli while simultaneously accounting for the multilevel correlation structure. The spatial correlation among the ERP profiles is captured through basis-space Matérn assumptions that allow either separable or nonseparable spatial correlations over time. We induce both adaptive regularization over time and spatial smoothness across electrodes via a correlated normal-exponential-gamma (CNEG) prior on the fixed effect coefficient functions. Our proposed framework includes both Gaussian models as well as robust models using heavier-tailed distributions to make the regression automatically robust to outliers. We introduce predictive methods to select among Gaussian vs. robust models and models with separable vs. non-separable spatiotemporal correlation structures. Our proposed analysis produces global tests for stimuli effects across entire time (or time-frequency) and electrode domains, plus multiplicity-adjusted pointwise inference based on experiment-wise error rate or false discovery rate to flag spatiotemporal (or spatio-temporal-frequency) regions that characterize stimuli differences, and can also produce inference for any prespecified waveform components. Our analysis of the smoking cessation ERP data set reveals numerous effects across different types of visual stimuli.
Collapse
Affiliation(s)
- Hongxiao Zhu
- Department of Statistics, Virginia Tech, Blacksburg, VA, USA.
| | - Francesco Versace
- Department of Behavioral Science, The University of Texas M.D. Anderson Cancer Center, Houston, TX, USA
| | - Paul M Cinciripini
- Department of Behavioral Science, The University of Texas M.D. Anderson Cancer Center, Houston, TX, USA
| | - Philip Rausch
- Department of Psychology, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Jeffrey S Morris
- Department of Biostatistics, The University of Texas M.D. Anderson Cancer Center, Houston, TX, USA
| |
Collapse
|
23
|
|
24
|
Lee W, Miranda MF, Rausch P, Baladandayuthapani V, Fazio M, Downs JC, Morris JS. Bayesian Semiparametric Functional Mixed Models for Serially Correlated Functional Data, with Application to Glaucoma Data. J Am Stat Assoc 2018; 114:495-513. [PMID: 31235987 PMCID: PMC6590079 DOI: 10.1080/01621459.2018.1476242] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2016] [Revised: 12/01/2017] [Indexed: 10/14/2022]
Abstract
Glaucoma, a leading cause of blindness, is characterized by optic nerve damage related to intraocular pressure (IOP), but its full etiology is unknown. Researchers at UAB have devised a custom device to measure scleral strain continuously around the eye under fixed levels of IOP, which here is used to assess how strain varies around the posterior pole, with IOP, and across glaucoma risk factors such as age. The hypothesis is that scleral strain decreases with age, which could alter biomechanics of the optic nerve head and cause damage that could eventually lead to glaucoma. To evaluate this hypothesis, we adapted Bayesian Functional Mixed Models to model these complex data consisting of correlated functions on spherical scleral surface, with nonparametric age effects allowed to vary in magnitude and smoothness across the scleral surface, multi-level random effect functions to capture within-subject correlation, and functional growth curve terms to capture serial correlation across IOPs that can vary around the scleral surface. Our method yields fully Bayesian inference on the scleral surface or any aggregation or transformation thereof, and reveals interesting insights into the biomechanical etiology of glaucoma. The general modeling framework described is very flexible and applicable to many complex, high-dimensional functional data.
Collapse
Affiliation(s)
- Wonyul Lee
- Department of Biostatistics, University of Texas M.D. Anderson Cancer Center, Houston, TX 77230
| | - Michelle F Miranda
- Department of Biostatistics, University of Texas M.D. Anderson Cancer Center, Houston, TX 77230
| | - Philip Rausch
- Department of Psychology, Institut für Psychologie, Humboldt-Universität zu Berlin, Germany
| | | | - Massimo Fazio
- Department of Ophthalmology, University of Alabama at Birmingham, Birmingham, AL 35294
| | - J Crawford Downs
- Department of Ophthalmology, University of Alabama at Birmingham, Birmingham, AL 35294
| | - Jeffrey S Morris
- Department of Biostatistics, University of Texas M.D. Anderson Cancer Center, Houston, TX 77230
| |
Collapse
|
25
|
Park SY, Staicu AM, Xiao L, Crainiceanu CM. Simple fixed-effects inference for complex functional models. Biostatistics 2018; 19:137-152. [PMID: 29036541 PMCID: PMC5862370 DOI: 10.1093/biostatistics/kxx026] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2016] [Revised: 04/09/2017] [Accepted: 05/07/2017] [Indexed: 11/14/2022] Open
Abstract
We propose simple inferential approaches for the fixed effects in complex functional mixed effects models. We estimate the fixed effects under the independence of functional residuals assumption and then bootstrap independent units (e.g. subjects) to conduct inference on the fixed effects parameters. Simulations show excellent coverage probability of the confidence intervals and size of tests for the fixed effects model parameters. Methods are motivated by and applied to the Baltimore Longitudinal Study of Aging, though they are applicable to other studies that collect correlated functional data.
Collapse
Affiliation(s)
- So Young Park
- Department of Statistics, North Carolina State University, Raleigh, NC, USA
| | - Ana-Maria Staicu
- Department of Statistics, North Carolina State University, Raleigh, NC, USA
| | - Luo Xiao
- Department of Statistics, North Carolina State University, Raleigh, NC, USA
| | | |
Collapse
|
26
|
Peterson GCL, Li D, Reich BJ, Brenner D. Spatial prediction of crystalline defects observed in molecular dynamic simulations of plastic damage. J Appl Stat 2017. [DOI: 10.1080/02664763.2016.1221915] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
| | - Dong Li
- Department of Material Sciences, North Carolina State University, Raleigh, NC, USA
| | - Brian J. Reich
- Department of Statistics, North Carolina State University, Raleigh, NC, USA
| | - Donald Brenner
- Department of Material Sciences, North Carolina State University, Raleigh, NC, USA
| |
Collapse
|
27
|
Morris JS, Baladandayuthapani V. Statistical Contributions to Bioinformatics: Design, Modeling, Structure Learning, and Integration. STAT MODEL 2017; 17:245-289. [PMID: 29129969 PMCID: PMC5679480 DOI: 10.1177/1471082x17698255] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The advent of high-throughput multi-platform genomics technologies providing whole-genome molecular summaries of biological samples has revolutionalized biomedical research. These technologiees yield highly structured big data, whose analysis poses significant quantitative challenges. The field of Bioinformatics has emerged to deal with these challenges, and is comprised of many quantitative and biological scientists working together to effectively process these data and extract the treasure trove of information they contain. Statisticians, with their deep understanding of variability and uncertainty quantification, play a key role in these efforts. In this article, we attempt to summarize some of the key contributions of statisticians to bioinformatics, focusing on four areas: (1) experimental design and reproducibility, (2) preprocessing and feature extraction, (3) unified modeling, and (4) structure learning and integration. In each of these areas, we highlight some key contributions and try to elucidate the key statistical principles underlying these methods and approaches. Our goals are to demonstrate major ways in which statisticians have contributed to bioinformatics, encourage statisticians to get involved early in methods development as new technologies emerge, and to stimulate future methodological work based on the statistical principles elucidated in this article and utilizing all availble information to uncover new biological insights.
Collapse
Affiliation(s)
- Jeffrey S Morris
- Department of Biostatistics, The University of Texas M.D. Anderson Cancer Center, Houston, Texas, USA
| | | |
Collapse
|
28
|
Abstract
In this article, Greven and Scheipl describe an impressively general framework for performing functional regression that builds upon the generalized additive modeling framework. Over the past number of years, my collaborators and I have also been developing a general framework for functional regression, functional mixed models, which shares many similarities with this framework, but has many differences as well. In this discussion, I compare and contrast these two frameworks, to hopefully illuminate characteristics of each, highlighting their respecitve strengths and weaknesses, and providing recommendations regarding the settings in which each approach might be preferable.
Collapse
Affiliation(s)
- Jeffrey S Morris
- The University of Texas, MD Anderson Cancer Center, Unit 1411, PO Box 301402, Houston, TX 77230-1402
| |
Collapse
|
29
|
Chen H, Zhao B, Porges EC, Cohen RA, Ebner NC. Edgewise and subgraph-level tests for brain networks. Stat Med 2016; 35:4994-5008. [PMID: 27397632 DOI: 10.1002/sim.7039] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2015] [Revised: 05/27/2016] [Accepted: 06/16/2016] [Indexed: 11/10/2022]
Abstract
Resting-state functional magnetic resonance image is a useful technique for investigating brain functional connectivity at rest. In this work, we develop flexible regression models and methods for determining differences in resting-state functional connectivity as a function of age, gender, drug intervention, or neuropsychiatric disorders. We propose two complementary methods for identifying changes of edges and subgraphs. (i) For detecting changes of edges, we select the optimal model at each edge and then conduct contrast tests to identify the effects of the important variables while controlling the familywise error rate. (ii) We adopt the network-based statistics method to improve power by incorporating the graph topological structure. Both methods have wide applications for low signal-to-noise ratio data. We propose stability criteria for the choice of threshold in the network-based statistics procedure and utilize efficient massive parallel procedure to speed up the estimation and inference procedure. Results from our simulation studies show that the thresholds chosen by the proposed stability criteria outperform the Bonferroni threshold. To demonstrate applicability, we use both methods in the context of the Oxytocin and Aging Study to determine effects of age, gender, and drug treatment on resting-state functional connectivity, as well as in the context of the Autism Brain Imaging Data Exchange Study to determine effects of autism spectrum disorder on functional connectivity at rest. Copyright © 2016 John Wiley & Sons, Ltd.
Collapse
Affiliation(s)
- Huaihou Chen
- Department of Biostatistics, University of Florida, Gainesville, FL, U.S.A.. .,Center for Cognitive Aging and Memory (CAM), Institute on Aging, McKnight Brain Institute, Department of Aging Geriatric Research, University of Florida, Gainesville, FL, U.S.A..
| | - Bingxin Zhao
- Department of Biostatistics, University of Florida, Gainesville, FL, U.S.A
| | - Eric C Porges
- Center for Cognitive Aging and Memory (CAM), Institute on Aging, McKnight Brain Institute, Department of Aging Geriatric Research, University of Florida, Gainesville, FL, U.S.A
| | - Ronald A Cohen
- Center for Cognitive Aging and Memory (CAM), Institute on Aging, McKnight Brain Institute, Department of Aging Geriatric Research, University of Florida, Gainesville, FL, U.S.A
| | - Natalie C Ebner
- Center for Cognitive Aging and Memory (CAM), Institute on Aging, McKnight Brain Institute, Department of Aging Geriatric Research, University of Florida, Gainesville, FL, U.S.A.,Department of Psychology, University of Florida, Gainesville, FL, U.S.A
| |
Collapse
|