{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,10,23]],"date-time":"2026-10-23T22:16:10Z","timestamp":1792793770230,"version":"4.5.0"},"reference-count":25,"publisher":"Oxford University Press (OUP)","issue":"1","license":[{"start":{"date-parts":[[2026,1,1]],"date-time":"2026-01-01T00:00:00Z","timestamp":1767225600000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/pages\/standard-publication-reuse-rights"}],"funder":[{"DOI":"10.13039\/100000065","name":"National Institute of Neurological Disorders and Stroke","doi-asserted-by":"publisher","award":["R01NS097423-01"],"award-info":[{"award-number":["R01NS097423-01"]}],"id":[{"id":"10.13039\/100000065","id-type":"DOI","asserted-by":"publisher"}]},{"name":"European Union and the Next-Generation EU","id":[{"id":"https:\/\/ror.org\/019w4f821","id-type":"ROR","asserted-by":"crossref"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2026,1,20]]},"abstract":"<jats:title>Summary<\/jats:title>\n                  <jats:p>Accelerometer studies typically include repeated 24-h observations over several days and more recently have collected information on participants for weeks, months, or years. Meanwhile, there is growing evidence to suggest that components of physical activity beyond average or expected values are important for quantifying and understanding activity. Motivated by these emerging aspects, we introduce Multilevel Functional Quantile Principal Component Analysis (MFQPCA), a novel dimension-reduction technique that extends Functional Quantile Principal Component Analysis to hierarchical functional data. We examine accelerometry data collected for 4 to 9\u2009months from congestive heart failure patients, where daily activity profiles are nested within individuals and exhibit substantial variability: participants often differ in the timing and intensity of physical activity behaviors across days and from each other, and capturing information beyond expected-value curves is necessary for a robust quantification of diurnal patterns. To this end, MFQPCA decomposes quantile-specific variability into between-participant and within-participant components, estimating shared patterns and corresponding scores at each level and quantile of interest. Our method robustly captures complex distributional features in hierarchical data, disentangles between- and within-participant sources of variability at different quantile levels, and facilitates the study of longitudinal changes in these distributions. We produce participant- and participant-day-level 10%, 30%, 50%, 70% and 90% quantile curves over 24\u2009h, revealing that day-to-day variability dominates sedentary periods, whereas between-participant differences grow in vigorous activity. An efficient alternating minimization algorithm is implemented in the open-source R package FunQ, making MFQPCA readily applicable to a wide range of hierarchical functional data settings.<\/jats:p>","DOI":"10.1093\/biostatistics\/kxag017","type":"journal-article","created":{"date-parts":[[2026,5,22]],"date-time":"2026-05-22T11:42:53Z","timestamp":1779450173000},"source":"Crossref","is-referenced-by-count":0,"title":["Multilevel functional quantile principal component analysis"],"prefix":"10.1093","volume":"27","author":[{"ORCID":"https:\/\/orcid.org\/0000-0003-2044-4170","authenticated-orcid":false,"given":"\u00c1lvaro","family":"M\u00e9ndez-Civieta","sequence":"first","affiliation":[{"name":"Department of Biostatistics, Columbia University , 722 W 178 St , New York,NY 10032,","place":["United States"]},{"name":"uc3m-Santander Big Data Institute, University Carlos III of Madrid , C. Madrid, 126 ,Madrid 28903,","place":["Spain"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Ying","family":"Wei","sequence":"additional","affiliation":[{"name":"Department of Biostatistics, Columbia University , 722 W 178 St , New York,NY 10032,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Jeff","family":"Goldsmith","sequence":"additional","affiliation":[{"name":"Department of Biostatistics, Columbia University , 722 W 178 St , New York,NY 10032,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2026,6,4]]},"reference":[{"key":"2026102110230509600_kxag017-B1","doi-asserted-by":"crossref","first-page":"611","DOI":"10.1111\/1467-9868.00196","article-title":"Probabilistic principal component analysis","volume":"61","author":"Bishop","year":"1999","journal-title":"J R Stat Soc Ser B Method"},{"key":"2026102110230509600_kxag017-B2","doi-asserted-by":"crossref","first-page":"126","DOI":"10.1002\/sta4.50","article-title":"Multilevel sparse functional principal component analysis","volume":"3","author":"Di","year":"2014","journal-title":"Stat"},{"key":"2026102110230509600_kxag017-B3","doi-asserted-by":"crossref","first-page":"458","DOI":"10.1214\/08-AOAS206","article-title":"Multilevel functional principal component analysis","volume":"3","author":"Di","year":"2009","journal-title":"Ann Appl Stat"},{"key":"2026102110230509600_kxag017-B4","doi-asserted-by":"crossref","first-page":"46","DOI":"10.1016\/j.csda.2016.07.010","article-title":"A note on modeling sparse exponential-family functional response curves","volume":"105","author":"Gertheiss","year":"2017","journal-title":"Comput Stat Data Anal"},{"key":"2026102110230509600_kxag017-B5","doi-asserted-by":"crossref","first-page":"41","DOI":"10.1111\/j.1541-0420.2012.01808.x","article-title":"Corrected confidence bands for functional data using principal components","volume":"69","author":"Goldsmith","year":"2013","journal-title":"Biometrics"},{"key":"2026102110230509600_kxag017-B6","doi-asserted-by":"crossref","first-page":"1723","DOI":"10.1249\/MSS.0000000000000968","article-title":"New insights into activity patterns in children, found using functional data analyses","volume":"48","author":"Goldsmith","year":"2016","journal-title":"Med Sci Sports Exerc"},{"key":"2026102110230509600_kxag017-B7","doi-asserted-by":"crossref","first-page":"344","DOI":"10.1111\/biom.12278","article-title":"Generalized multilevel function-on-scalar regression and principal component analysis","volume":"71","author":"Goldsmith","year":"2015","journal-title":"Biometrics"},{"key":"2026102110230509600_kxag017-B8","doi-asserted-by":"crossref","first-page":"1022","DOI":"10.1214\/10-EJS575","article-title":"Longitudinal functional principal component analysis","volume":"4","author":"Greven","year":"2010","journal-title":"Electron J Stat"},{"key":"2026102110230509600_kxag017-B9","doi-asserted-by":"publisher","author":"He","year":"2023","DOI":"10.32614\/CRAN.pacakge.conquer"},{"key":"2026102110230509600_kxag017-B10","doi-asserted-by":"crossref","first-page":"367","DOI":"10.1016\/j.jeconom.2021.07.010","article-title":"Smoothed quantile regression with large-scale inference","volume":"232","author":"He","year":"2023","journal-title":"J Econom"},{"key":"2026102110230509600_kxag017-B11","doi-asserted-by":"publisher","first-page":"553","DOI":"10.1080\/01621459.2018.1482750","article-title":"Multilevel matrix-variate analysis and its application to accelerometry-measured physical activity in clinical populations","volume":"114","author":"Huang","year":"2019","journal-title":"J Am Stat Assoc"},{"key":"2026102110230509600_kxag017-B12","doi-asserted-by":"crossref","first-page":"587","DOI":"10.1093\/biomet\/87.3.587","article-title":"Principal component models for sparse functional data","volume":"87","author":"James","year":"2000","journal-title":"Biometrika"},{"key":"2026102110230509600_kxag017-B13","doi-asserted-by":"publisher","author":"Koenker","year":"2023","DOI":"10.32614\/CRAN.package.quantreg,"},{"key":"2026102110230509600_kxag017-B14","doi-asserted-by":"crossref","first-page":"1299","DOI":"10.1080\/09291016.2021.1929673","article-title":"Registration of 24-hour accelerometric rest-activity profiles and its application to human chronotypes","volume":"53","author":"McDonnell","year":"2022","journal-title":"Biol Rhythm Res"},{"key":"2026102110230509600_kxag017-B15","doi-asserted-by":"crossref","first-page":"kxae040","DOI":"10.1093\/biostatistics\/kxae040","article-title":"Functional quantile principal component analysis","volume":"26","author":"M\u00e9ndez-Civieta","year":"2025","journal-title":"Biostatistics"},{"key":"2026102110230509600_kxag017-B16","doi-asserted-by":"crossref","first-page":"1352","DOI":"10.1198\/016214506000000465","article-title":"Using wavelet-based functional mixed models to characterize population heterogeneity in accelerometer profiles: a case study","volume":"101","author":"Morris","year":"2006","journal-title":"J Am Stat Assoc"},{"key":"2026102110230509600_kxag017-B17","volume-title":"R: a language and environment for statistical computing","author":"R Core Team","year":"2024"},{"key":"2026102110230509600_kxag017-B18","author":"Scheipl","year":"2023"},{"issue":":","key":"2026102110230509600_kxag017-B19","doi-asserted-by":"publisher","first-page":"5222","DOI":"10.1002\/sim.5989","article-title":"An introduction with medical applications to functional data analysis","volume":"32","author":"S\u00f8rensen","year":"2013","journal-title":"Stat Med"},{"key":"2026102110230509600_kxag017-B20","doi-asserted-by":"publisher","first-page":"102","DOI":"10.1016\/j.ypmed.2017.05.030","article-title":"Re-evaluating the effect of age on physical activity over the lifespan","volume":"101","author":"Varma","year":"2017","journal-title":"Prev Med"},{"issue":":","key":"2026102110230509600_kxag017-B21","doi-asserted-by":"publisher","first-page":"48","DOI":"10.1111\/biom.12963","article-title":"Registration for exponential family functional data","volume":"75","author":"Wrobel","year":"2019","journal-title":"Biometrics"},{"issue":":","key":"2026102110230509600_kxag017-B22","doi-asserted-by":"publisher","first-page":"352","DOI":"10.1093\/biostatistics\/kxu045","article-title":"Quantifying the lifetime circadian rhythm of physical activity: a covariate-dependent functional approach","volume":"16","author":"Xiao","year":"2015","journal-title":"Biostatistics"},{"issue":":","key":"2026102110230509600_kxag017-B23","doi-asserted-by":"publisher","first-page":"409","DOI":"10.1007\/s11222-014-9485-x","article-title":"Fast covariance estimation for high-dimensional functional data","volume":"26","author":"Xiao","year":"2016","journal-title":"Stat Comput"},{"issue":":","key":"2026102110230509600_kxag017-B25","doi-asserted-by":"crossref","first-page":"577","DOI":"10.1198\/016214504000001745","article-title":"Functional data analysis for sparse longitudinal data","volume":"100","author":"Yao","year":"2005","journal-title":"J Am Stat Assoc"},{"key":"2026102110230509600_kxag017-B26","doi-asserted-by":"publisher","first-page":"852","DOI":"10.1198\/jcgs.2011.10122","article-title":"Multilevel functional principal component analysis for high-dimensional data","volume":"20","author":"Zipunnikov","year":"2011","journal-title":"J Comput Graph Stat"}],"container-title":["Biostatistics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/biostatistics\/article-pdf\/27\/1\/kxag017\/68530088\/kxag017.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/biostatistics\/article-pdf\/27\/1\/kxag017\/68530088\/kxag017.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,10,21]],"date-time":"2026-10-21T14:23:19Z","timestamp":1792592599000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/biostatistics\/article\/doi\/10.1093\/biostatistics\/kxag017\/8708315"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,1]]},"references-count":25,"journal-issue":{"issue":"1","published-print":{"date-parts":[[2026,1,20]]}},"URL":"https:\/\/doi.org\/10.1093\/biostatistics\/kxag017","relation":{},"ISSN":["1465-4644","1468-4357"],"issn-type":[{"value":"1465-4644","type":"print"},{"value":"1468-4357","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2026,1]]},"published":{"date-parts":[[2026,1]]},"article-number":"kxag017"}}