{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,4,8]],"date-time":"2026-04-08T22:35:59Z","timestamp":1775687759946,"version":"3.50.1"},"reference-count":43,"publisher":"Oxford University Press (OUP)","issue":"11","license":[{"start":{"date-parts":[[2023,11,14]],"date-time":"2023-11-14T00:00:00Z","timestamp":1699920000000},"content-version":"vor","delay-in-days":13,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"name":"Academy of Finland R\u2019Life program funding","award":["329439"],"award-info":[{"award-number":["329439"]}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2023,11,1]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:sec><jats:title>Motivation<\/jats:title><jats:p>Due to advances in measuring technology, many new phenotype, gene expression, and other omics time-course datasets are now commonly available. Cluster analysis may provide useful information about the structure of such data.<\/jats:p><\/jats:sec><jats:sec><jats:title>Results<\/jats:title><jats:p>In this work, we propose BELMM (Bayesian Estimation of Latent Mixture Models): a flexible framework for analysing, clustering, and modelling time-series data in a Bayesian setting. The framework is built on mixture modelling: first, the mean curves of the mixture components are assumed to follow random walk smoothing priors. Second, we choose the most plausible model and the number of mixture components using the Reversible-jump Markov chain Monte Carlo. Last, we assign the individual time series into clusters based on the similarity to the cluster-specific trend curves determined by the latent random walk processes. We demonstrate the use of fast and slow implementations of our approach on both simulated and real time-series data using widely available software R, Stan, and CU-MSDSp.<\/jats:p><\/jats:sec><jats:sec><jats:title>Availability and implementation<\/jats:title><jats:p>The French mortality dataset is available at http:\/\/www.mortality.org, the Drosophila melanogaster embryogenesis gene expression data at https:\/\/www.ncbi.nlm.nih.gov\/geo\/query\/acc.cgi?acc=GSE121160. Details on our simulated datasets are available in the Supplementary Material, and R scripts and a detailed tutorial on GitHub at https:\/\/github.com\/ollisa\/BELMM. The software CU-MSDSp is available on GitHub at https:\/\/github.com\/jtchavisIII\/CU-MSDSp.<\/jats:p><\/jats:sec>","DOI":"10.1093\/bioinformatics\/btad686","type":"journal-article","created":{"date-parts":[[2023,11,14]],"date-time":"2023-11-14T18:45:04Z","timestamp":1699987504000},"source":"Crossref","is-referenced-by-count":3,"title":["BELMM: Bayesian model selection and random walk smoothing in time-series clustering"],"prefix":"10.1093","volume":"39","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-7330-3577","authenticated-orcid":false,"given":"Olli","family":"Sarala","sequence":"first","affiliation":[{"name":"Research Unit of Mathematical Sciences, University of Oulu , FI-90014 Oulu, Finland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-6958-5172","authenticated-orcid":false,"given":"Tanja","family":"Pyh\u00e4j\u00e4rvi","sequence":"additional","affiliation":[{"name":"Department of Forest Sciences, University of Helsinki , FI-00014 Helsinki, Finland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-2808-2768","authenticated-orcid":false,"given":"Mikko J","family":"Sillanp\u00e4\u00e4","sequence":"additional","affiliation":[{"name":"Research Unit of Mathematical Sciences, University of Oulu , FI-90014 Oulu, Finland"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2023,11,14]]},"reference":[{"key":"2023113004051825100_btad686-B1","doi-asserted-by":"crossref","DOI":"10.1007\/978-0-387-92298-0","volume-title":"Bayesian Computation with R. Use R!","author":"Albert","year":"2009"},{"key":"2023113004051825100_btad686-B2","doi-asserted-by":"crossref","first-page":"4970","DOI":"10.1038\/s41467-018-07455-9","article-title":"Quantifying post-transcriptional regulation in the development of Drosophila melanogaster","volume":"9","author":"Becker","year":"2018","journal-title":"Nat Commun"},{"key":"2023113004051825100_btad686-B3","doi-asserted-by":"crossref","first-page":"679","DOI":"10.1111\/1467-9868.00409","article-title":"Reversible jump, birth-and-death and more general continuous time Markov chain Monte Carlo samplers","volume":"65","author":"Capp\u00e9","year":"2003","journal-title":"J R Stat Soc Ser B (Stat Methodol)"},{"key":"2023113004051825100_btad686-B4","doi-asserted-by":"crossref","first-page":"1","DOI":"10.18637\/jss.v061.i06","article-title":"Nbclust: an R package for determining the relevant number of clusters in a data set","volume":"61","author":"Charrad","year":"2014","journal-title":"J Stat Soft"},{"key":"2023113004051825100_btad686-B5","author":"Chavis","year":"2020"},{"key":"2023113004051825100_btad686-B6","doi-asserted-by":"crossref","first-page":"100664","DOI":"10.1016\/j.softx.2021.100664","article-title":"CU-MSDSp: a flexible parallelized Reversible jump Markov chain Monte Carlo method","volume":"14","author":"Chavis","year":"2021","journal-title":"SoftwareX"},{"key":"2023113004051825100_btad686-B7","doi-asserted-by":"crossref","first-page":"e59795","DOI":"10.1371\/journal.pone.0059795","article-title":"Accelerating Bayesian hierarchical clustering of time series data with a randomised algorithm","volume":"8","author":"Darkins","year":"2013","journal-title":"PLoS One"},{"key":"2023113004051825100_btad686-B8","doi-asserted-by":"crossref","first-page":"197","DOI":"10.1023\/A:1008935410038","article-title":"On sequential Monte Carlo sampling methods for Bayesian filtering","volume":"10","author":"Doucet","year":"2000","journal-title":"Stat Comput"},{"key":"2023113004051825100_btad686-B9","doi-asserted-by":"crossref","first-page":"765","DOI":"10.21105\/joss.00765","article-title":"Philentropy: information theory and distance quantification with R","volume":"3","author":"Drost","year":"2018","journal-title":"J Open Source Softw"},{"key":"2023113004051825100_btad686-B10","doi-asserted-by":"crossref","first-page":"162","DOI":"10.1007\/978-3-642-41398-8_15","volume-title":"Advances in Intelligent Data Analysis XII","author":"Eirola","year":"2013"},{"key":"2023113004051825100_btad686-B11","volume-title":"Finite Mixture and Markov Switching Models. Springer Series in Statistics","author":"Fr\u00fchwirth-Schnatter","year":"2006"},{"key":"2023113004051825100_btad686-B12","author":"Golumbeanu","year":"2020"},{"key":"2023113004051825100_btad686-B13","doi-asserted-by":"crossref","first-page":"213","DOI":"10.1038\/s41598-018-36135-3","article-title":"Proteo-transcriptomic dynamics of cellular response to HIV-1 infection","volume":"9","author":"Golumbeanu","year":"2019","journal-title":"Sci Rep"},{"key":"2023113004051825100_btad686-B14","doi-asserted-by":"crossref","first-page":"711","DOI":"10.1093\/biomet\/82.4.711","article-title":"Reversible jump Markov chain Monte Carlo computation and Bayesian model determination","volume":"82","author":"Green","year":"1995","journal-title":"Biometrika"},{"key":"2023113004051825100_btad686-B15","doi-asserted-by":"crossref","first-page":"309","DOI":"10.1111\/j.1467-9574.2012.00516.x","article-title":"Model choice using reversible jump Markov chain Monte Carlo","volume":"66","author":"Hastie","year":"2012","journal-title":"Stat. Neerl"},{"key":"2023113004051825100_btad686-B16","doi-asserted-by":"crossref","first-page":"978","DOI":"10.1007\/11499145_99","volume-title":"Image Analysis, Lecture Notes in Computer Science","author":"Hautam\u00e4ki","year":"2005"},{"key":"2023113004051825100_btad686-B17","author":"Human Mortality Database","year":"2020"},{"key":"2023113004051825100_btad686-B18","doi-asserted-by":"crossref","first-page":"183","DOI":"10.1198\/1061860043010","article-title":"Bayesian P-splines","volume":"13","author":"Lang","year":"2004","journal-title":"J Comput Graph Stat"},{"key":"2023113004051825100_btad686-B19","first-page":"321","article-title":"Geodesic clustering for covariance matrices","volume":"22","author":"Lee","year":"2015","journal-title":"Commun Stat Appl Methods"},{"key":"2023113004051825100_btad686-B20","doi-asserted-by":"crossref","first-page":"589","DOI":"10.1016\/j.cels.2019.10.005","article-title":"Detecting, categorizing, and correcting coverage anomalies of RNA-seq quantification","volume":"9","author":"Ma","year":"2019","journal-title":"Cell Syst"},{"key":"2023113004051825100_btad686-B21","author":"Masarotto"},{"key":"2023113004051825100_btad686-B22","doi-asserted-by":"crossref","first-page":"355","DOI":"10.1146\/annurev-statistics-031017-100325","article-title":"Finite mixture models","volume":"6","author":"McLachlan","year":"2019","journal-title":"Annu Rev Stat Appl"},{"key":"2023113004051825100_btad686-B23","doi-asserted-by":"crossref","first-page":"563","DOI":"10.1007\/s11634-016-0264-8","article-title":"An effective strategy for initializing the EM algorithm in finite mixture models","volume":"10","author":"Michael","year":"2016","journal-title":"Adv Data Anal Classif"},{"key":"2023113004051825100_btad686-B24","doi-asserted-by":"crossref","first-page":"159","DOI":"10.1007\/BF02294245","article-title":"An examination of procedures for determining the number of clusters in a data set","volume":"50","author":"Milligan","year":"1985","journal-title":"Psychometrika"},{"key":"2023113004051825100_btad686-B25","doi-asserted-by":"crossref","first-page":"599","DOI":"10.1016\/j.physa.2006.10.073","article-title":"Tests of the random walk hypothesis for financial data","volume":"377","author":"Nakamura","year":"2007","journal-title":"Physica A"},{"key":"2023113004051825100_btad686-B26","doi-asserted-by":"crossref","first-page":"2159","DOI":"10.1007\/s10994-019-05818-x","article-title":"Covariance-based dissimilarity measures applied to clustering wide-sense stationary ergodic processes","volume":"108","author":"Peng","year":"2019","journal-title":"Mach Learn"},{"key":"2023113004051825100_btad686-B27","doi-asserted-by":"crossref","first-page":"409","DOI":"10.1093\/biomet\/asu008","article-title":"Distances and inference for covariance operators","volume":"101","author":"Pigoli","year":"2014","journal-title":"Biometrika"},{"key":"2023113004051825100_btad686-B28","doi-asserted-by":"crossref","DOI":"10.1002\/9781118573617","volume-title":"High-Dimensional Covariance Estimation. Wiley Series in Probability and Statistics","author":"Pourahmadi","year":"2013"},{"key":"2023113004051825100_btad686-B29","doi-asserted-by":"crossref","first-page":"945","DOI":"10.1093\/genetics\/155.2.945","article-title":"Inference of population structure using multilocus genotype data","volume":"155","author":"Pritchard","year":"2000","journal-title":"Genetics"},{"key":"2023113004051825100_btad686-B30","doi-asserted-by":"crossref","first-page":"e84955","DOI":"10.1371\/journal.pone.0084955","article-title":"Clustering gene expression regulators: new approach to disease subtyping","volume":"9","author":"Pyatnitskiy","year":"2014","journal-title":"PLoS One"},{"key":"2023113004051825100_btad686-B31","doi-asserted-by":"crossref","first-page":"9121","DOI":"10.1073\/pnas.132656399","article-title":"Cluster analysis of gene expression dynamics","volume":"99","author":"Ramoni","year":"2002","journal-title":"Proc Natl Acad Sci USA"},{"key":"2023113004051825100_btad686-B32","doi-asserted-by":"crossref","first-page":"731","DOI":"10.1111\/1467-9868.00095","article-title":"On Bayesian analysis of mixtures with an unknown number of components (with discussion)","volume":"59","author":"Richardson","year":"1997","journal-title":"J R Stat Soc Ser B (Stat Methodol)"},{"key":"2023113004051825100_btad686-B33","doi-asserted-by":"crossref","first-page":"461","DOI":"10.1214\/aos\/1176344136","article-title":"Estimating the dimension of a model","volume":"6","author":"Schwarz","year":"1978","journal-title":"Ann Statist"},{"key":"2023113004051825100_btad686-B34","doi-asserted-by":"crossref","first-page":"289","DOI":"10.32614\/RJ-2016-021","article-title":"mclust 5: clustering, classification and density estimation using Gaussian finite mixture models","volume":"8","author":"Scrucca","year":"2016","journal-title":"R J"},{"key":"2023113004051825100_btad686-B35","doi-asserted-by":"crossref","DOI":"10.1017\/CBO9781139344203","volume-title":"Bayesian Filtering and Smoothing","author":"S\u00e4rkk\u00e4","year":"2013"},{"key":"2023113004051825100_btad686-B36","author":"Stan Development Team","year":"2020"},{"key":"2023113004051825100_btad686-B37","doi-asserted-by":"crossref","first-page":"1020","DOI":"10.1080\/01621459.2016.1147355","article-title":"Detecting and localizing differences in functional time series dynamics: a case study in molecular biophysics","volume":"111","author":"Tavakoli","year":"2016","journal-title":"J Am Stat Assoc"},{"key":"2023113004051825100_btad686-B38","doi-asserted-by":"crossref","first-page":"411","DOI":"10.1111\/1467-9868.00293","article-title":"Estimating the number of clusters in a data set via the gap statistic","volume":"63","author":"Tibshirani","year":"2001","journal-title":"J R Stat Soc Ser B (Stat Methodol)"},{"key":"2023113004051825100_btad686-B39","doi-asserted-by":"crossref","first-page":"100323","DOI":"10.1016\/j.alcr.2019.100323","article-title":"An overview of mixture modelling for latent evolutions in longitudinal data: modelling approaches, fit statistics and software","volume":"43","author":"van der Nest","year":"2020","journal-title":"Adv Life Course Res"},{"key":"2023113004051825100_btad686-B40","doi-asserted-by":"crossref","first-page":"e0131739","DOI":"10.1371\/journal.pone.0131739","article-title":"Overfitting Bayesian mixture models with an unknown number of components","volume":"10","author":"van Havre","year":"2015","journal-title":"PLoS One"},{"key":"2023113004051825100_btad686-B41","doi-asserted-by":"crossref","DOI":"10.7554\/eLife.27451","article-title":"Transcription factor clusters regulate genes in eukaryotic cells","volume":"6","author":"Wollman","year":"2017","journal-title":"eLife"},{"key":"2023113004051825100_btad686-B42","article-title":"Functional data clustering using principal curve methods","volume":"51","author":"Wu","year":"2021","journal-title":"Commun Stat Theory Methods"},{"key":"2023113004051825100_btad686-B43","doi-asserted-by":"crossref","first-page":"106901","DOI":"10.1016\/j.csda.2019.106901","article-title":"Borrowing strength and borrowing index for Bayesian hierarchical models","volume":"144","author":"Xu","year":"2020","journal-title":"Comput Stat Data Anal"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btad686\/53842744\/btad686.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/11\/btad686\/53941376\/btad686.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/11\/btad686\/53941376\/btad686.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2024,11,1]],"date-time":"2024-11-01T22:30:25Z","timestamp":1730500225000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/doi\/10.1093\/bioinformatics\/btad686\/7420213"}},"subtitle":[],"editor":[{"given":"Peter","family":"Robinson","sequence":"additional","affiliation":[],"role":[{"role":"editor","vocabulary":"crossref"}]}],"short-title":[],"issued":{"date-parts":[[2023,11,1]]},"references-count":43,"journal-issue":{"issue":"11","published-print":{"date-parts":[[2023,11,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btad686","relation":{},"ISSN":["1367-4811"],"issn-type":[{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2023,11,1]]},"published":{"date-parts":[[2023,11,1]]},"article-number":"btad686"}}