Longitudinal functional principal component modelling via Stochastic Approximation Monte Carlo
Type
ArticleKAUST Grant Number
KUS-CI-016-04Date
2010-06Permanent link to this record
http://hdl.handle.net/10754/598736
Metadata
Show full item recordAbstract
The authors consider the analysis of hierarchical longitudinal functional data based upon a functional principal components approach. In contrast to standard frequentist approaches to selecting the number of principal components, the authors do model averaging using a Bayesian formulation. A relatively straightforward reversible jump Markov Chain Monte Carlo formulation has poor mixing properties and in simulated data often becomes trapped at the wrong number of principal components. In order to overcome this, the authors show how to apply Stochastic Approximation Monte Carlo (SAMC) to this problem, a method that has the potential to explore the entire space and does not become trapped in local extrema. The combination of reversible jump methods and SAMC in hierarchical longitudinal functional data is simplified by a polar coordinate representation of the principal components. The approach is easy to implement and does well in simulated data in determining the distribution of the number of principal components, and in terms of its frequentist estimation properties. Empirical applications are also presented.Citation
Martinez JG, Liang F, Zhou L, Carroll RJ (2010) Longitudinal functional principal component modelling via Stochastic Approximation Monte Carlo. Can J Statistics 38: 256–270. Available: http://dx.doi.org/10.1002/cjs.10062.Sponsors
We wish to thank the editor and a referee for their detailed reading of the original manuscript and their many suggestions that led to a significant improvement in the paper. We also thank the participants of the International Symposium in Statistics on Inferences in Generalized Linear Longitudinal Mixed Models, whose comments led us to rethink important parts of this work. Martinez was supported by a post-doctoral training grant from the National Cancer Institute (R25T-CA90301), Carroll's research was supported by a grant from the National Cancer Institute (CA57030), Liang's research was supported by a grant from the National Science Foundation (DMS-0706755) and Zhou was supported by a grant from the National Science Foundation (DMS-0907170). Both Carroll and Liang were also supported by Award Number KUS-CI-016-04, made by King Abdullah University of Science and Technology (KAUST).Publisher
WileyJournal
Canadian Journal of StatisticsPubMed ID
20689648PubMed Central ID
PMC2915799ae974a485f413a2113503eed53cd6c53
10.1002/cjs.10062
Scopus Count
Collections
Publications Acknowledging KAUST SupportRelated articles
- Bayesian phylogeny analysis via stochastic approximation Monte Carlo.
- Authors: Cheon S, Liang F
- Issue date: 2009 Nov
- Phylogenetic tree construction using sequential stochastic approximation Monte Carlo.
- Authors: Cheon S, Liang F
- Issue date: 2008 Jan
- An efficient interpolation technique for jump proposals in reversible-jump Markov chain Monte Carlo calculations.
- Authors: Farr WM, Mandel I, Stevens D
- Issue date: 2015 Jun
- Estimation of under-reporting in epidemics using approximations.
- Authors: Gamado K, Streftaris G, Zachary S
- Issue date: 2017 Jun
- Abrupt motion tracking via intensively adaptive Markov-chain Monte Carlo sampling.
- Authors: Zhou X, Lu Y, Lu J, Zhou J
- Issue date: 2012 Feb