{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,1,28]],"date-time":"2026-01-28T21:20:14Z","timestamp":1769635214629,"version":"3.49.0"},"reference-count":25,"publisher":"Oxford University Press (OUP)","issue":"5","license":[{"start":{"date-parts":[[2018,8,8]],"date-time":"2018-08-08T00:00:00Z","timestamp":1533686400000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/journals\/pages\/open_access\/funder_policies\/chorus\/standard_publication_model"}],"funder":[{"DOI":"10.13039\/100000069","name":"National Institute of Arthritis and Musculoskeletal and Skin Diseases","doi-asserted-by":"publisher","award":["RO1AR05974"],"award-info":[{"award-number":["RO1AR05974"]}],"id":[{"id":"10.13039\/100000069","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100006108","name":"National Center for Advancing Translational Sciences","doi-asserted-by":"publisher","award":["(BU-CTSI) 1UL1TR001430"],"award-info":[{"award-number":["(BU-CTSI) 1UL1TR001430"]}],"id":[{"id":"10.13039\/100006108","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2019,3,1]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:sec><jats:title>Motivation<\/jats:title><jats:p>Clustering algorithms like K-Means and standard Gaussian mixture models (GMM) fail to account for the structure of variability of replicated data or repeated measures over time. Additionally, a priori cluster number assumptions add an additional complexity to the process. Current methods to optimize cluster labels and number can be inaccurate or computationally intensive for temporal gene expression data with this additional variability.<\/jats:p><\/jats:sec><jats:sec><jats:title>Results<\/jats:title><jats:p>An extension to a model-based clustering algorithm is proposed using mixtures of mixed effects polynomial regression models and the EM algorithm with an entropy penalized log-likelihood function (EPEM). The EPEM is used to cluster temporal gene expression data with this additional variability. The addition of random effects in our model decreased the misclassification error when compared to mixtures of fixed effects models or other methods such as K-Means and GMM. Applying our method to microarray data from a fracture healing study revealed distinct temporal patterns of gene expression.<\/jats:p><\/jats:sec><jats:sec><jats:title>Availability and Implementation<\/jats:title><jats:p>https:\/\/github.com\/darlenelu72\/EPEM-GMM.<\/jats:p><\/jats:sec><jats:sec><jats:title>Supplementary information<\/jats:title><jats:p>Supplementary data are available at Bioinformatics online.<\/jats:p><\/jats:sec>","DOI":"10.1093\/bioinformatics\/bty696","type":"journal-article","created":{"date-parts":[[2018,8,7]],"date-time":"2018-08-07T19:12:59Z","timestamp":1533669179000},"page":"778-786","source":"Crossref","is-referenced-by-count":3,"title":["Clustering of temporal gene expression data with mixtures of mixed effects models with a penalized likelihood"],"prefix":"10.1093","volume":"35","author":[{"given":"Darlene","family":"Lu","sequence":"first","affiliation":[{"name":"Biostatistics, Boston University, Boston, MA, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Yorghos","family":"Tripodis","sequence":"additional","affiliation":[{"name":"Biostatistics, Boston University, Boston, MA, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Louis C","family":"Gerstenfeld","sequence":"additional","affiliation":[{"name":"Orthopaedic Surgery, Boston University, Boston, MA, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Serkalem","family":"Demissie","sequence":"additional","affiliation":[{"name":"Biostatistics, Boston University, Boston, MA, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2018,8,8]]},"reference":[{"key":"2023013107250632700_bty696-B1","doi-asserted-by":"crossref","first-page":"716","DOI":"10.1109\/TAC.1974.1100705","article-title":"A new look at the statistical model identification","volume":"19","author":"Aike","year":"1974","journal-title":"IEEE Trans. Automatic Control"},{"key":"2023013107250632700_bty696-B2","doi-asserted-by":"crossref","first-page":"e5393.","DOI":"10.1371\/journal.pone.0005393","article-title":"Transcriptional analysis of fracture healing and the induction of embryonic stem cell-related genes","volume":"4","author":"Bais","year":"2009","journal-title":"PloS One"},{"key":"2023013107250632700_bty696-B4","doi-asserted-by":"crossref","first-page":"719","DOI":"10.1109\/34.865189","article-title":"Assessing a mixture model for clustering with the integrated classification likelihood","volume":"22","author":"Biernacki","year":"2000","journal-title":"IEEE Trans. Pattern Anal. Mach. Intell"},{"key":"2023013107250632700_bty696-B3","doi-asserted-by":"crossref","first-page":"561","DOI":"10.1016\/S0167-9473(02)00163-9","article-title":"Choosing starting values for the EM algorithm for getting the highest likelihood in multivariate Gaussian mixture models","volume":"41","author":"Biernacki","year":"2003","journal-title":"Comput. Stat. Data Anal"},{"key":"2023013107250632700_bty696-B5","doi-asserted-by":"crossref","first-page":"243","DOI":"10.1191\/1471082X05st096oa","article-title":"Mixture of linear mixed models for clustering gene expression profiles from repeated microarray experiments","volume":"5","author":"Celeux","year":"2005","journal-title":"Stat. Model"},{"key":"2023013107250632700_bty696-B6","doi-asserted-by":"crossref","first-page":"2308","DOI":"10.1080\/00949655.2015.1109096","article-title":"Unsupervised learning of regression mixture models with unknown number of components","volume":"86","author":"Chamroukhi","year":"2016","journal-title":"J. Stat. Comput. Simul"},{"key":"2023013107250632700_bty696-B7","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1111\/j.2517-6161.1977.tb01600.x","article-title":"Maximum likelihood from incomplete data via the EM Algorithm","volume":"39","author":"Dempster","year":"1977","journal-title":"J. R. Stat. Soc"},{"key":"2023013107250632700_bty696-B8","doi-asserted-by":"crossref","first-page":"381","DOI":"10.1109\/34.990138","article-title":"Unsupervised learning of finite mixture models","volume":"24","author":"Figueiredo","year":"2002","journal-title":"IEEE Trans. Pattern Anal. Mach. Intell"},{"key":"2023013107250632700_bty696-B9","doi-asserted-by":"crossref","first-page":"155","DOI":"10.3233\/IDA-2007-11204","article-title":"Mixture-model cluster analysis using information theoretical criteria","volume":"11","author":"Fonesca","year":"2007","journal-title":"Intell. Data Anal"},{"key":"2023013107250632700_bty696-B10","article-title":"mclust Version 4 for R: normal mixture modeling for model-based clustering, classification, and density estimation","author":"Fraley","year":"2012","journal-title":"Technical Report"},{"key":"2023013107250632700_bty696-B11","author":"Gaffney","year":"2004"},{"key":"2023013107250632700_bty696-B12","doi-asserted-by":"crossref","first-page":"63","DOI":"10.1145\/312129.312198","volume-title":"Proceedings of the Fifth ACM SIGKDD International Conference on Knowledge Discovery and Data Mining","author":"Gaffney","year":"1999"},{"key":"2023013107250632700_bty696-B13","doi-asserted-by":"crossref","first-page":"2597","DOI":"10.1002\/jbmr.486","article-title":"The transcriptome of fracture healing defines mechanisms of coordination of skeletal and vascular development during endochondral bone formation","volume":"26","author":"Grimes","year":"2011","journal-title":"J. Bone Min. Res"},{"key":"2023013107250632700_bty696-B14","first-page":"100","article-title":"Algorithm AS 136 A K-means clustering algorithm","volume":"28","author":"Hartigan","year":"1979","journal-title":"J. R. Stat. Soc. Ser. C (Appl. Stat.)"},{"key":"2023013107250632700_bty696-B15","doi-asserted-by":"crossref","first-page":"193","DOI":"10.1007\/BF01908075","article-title":"Comparing partitions","volume":"2","author":"Hubert","year":"1985","journal-title":"J. Classif"},{"key":"2023013107250632700_bty696-B16","doi-asserted-by":"crossref","first-page":"1204","DOI":"10.1359\/jbmr.080317","article-title":"Genetic variation in the patterns of skeletal progenitor cell differentiation and progression during endochondral bone formation affects the rate of fracture healing","volume":"23","author":"Jepsen","year":"2008","journal-title":"J. Bone Min. Res"},{"key":"2023013107250632700_bty696-B17","doi-asserted-by":"crossref","first-page":"25","DOI":"10.2991\/jsta.2016.15.1.3","article-title":"Parameter estimation in gamma mixture model using normal-based approximation","volume":"15","author":"Lakshmi","year":"2016","journal-title":"J. Stat. Theory Appl"},{"key":"2023013107250632700_bty696-B18","doi-asserted-by":"crossref","DOI":"10.1002\/9780470191613","volume-title":"The EM Algorithm and Extensions","author":"McLachlan","year":"2008"},{"key":"2023013107250632700_bty696-B19","doi-asserted-by":"crossref","first-page":"175","DOI":"10.1002\/0471721182","volume-title":"Finite Mixture Models","author":"McLachlan","year":"2000"},{"key":"2023013107250632700_bty696-B20","doi-asserted-by":"crossref","first-page":"1745","DOI":"10.1093\/bioinformatics\/btl165","article-title":"A mixture model with random-effects components for clustering correlated gene-expression profiles","volume":"22","author":"Ng","year":"2006","journal-title":"Bioinformatics"},{"key":"2023013107250632700_bty696-B21","doi-asserted-by":"crossref","first-page":"894","DOI":"10.1080\/01621459.1997.10474044","article-title":"Practical Bayesian density estimation using mixtures of normals","volume":"92","author":"Roeder","year":"1997","journal-title":"J. Am. Stat. Assoc"},{"key":"2023013107250632700_bty696-B22","volume-title":"Numerical Taxonomy. The Principles and Practice of Numerical Classification","author":"Sneath","year":"1973","edition":"1st edn."},{"key":"2023013107250632700_bty696-B23","doi-asserted-by":"crossref","first-page":"12837","DOI":"10.1073\/pnas.0504609102","article-title":"Significance analysis of time course microarray experiments","volume":"102","author":"Storey","year":"2005","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"2023013107250632700_bty696-B24","doi-asserted-by":"crossref","first-page":"724","DOI":"10.1359\/jbmr.091021","article-title":"Acute phosphate restriction leads to impaired fracture healing and resistance to bmp-2","volume":"25","author":"Wigner","year":"2010","journal-title":"J. Bone Min. Res"},{"key":"2023013107250632700_bty696-B25","doi-asserted-by":"crossref","first-page":"3950","DOI":"10.1016\/j.patcog.2012.04.031","article-title":"A robust EM clustering algorithm for Gaussian mixture models","volume":"45","author":"Yang","year":"2012","journal-title":"Pattern Recogn"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/35\/5\/778\/48965861\/bioinformatics_35_5_778.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/35\/5\/778\/48965861\/bioinformatics_35_5_778.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2024,7,9]],"date-time":"2024-07-09T00:51:12Z","timestamp":1720486272000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/35\/5\/778\/5068161"}},"subtitle":[],"editor":[{"given":"Jonathan","family":"Wren","sequence":"additional","affiliation":[],"role":[{"role":"editor","vocabulary":"crossref"}]}],"short-title":[],"issued":{"date-parts":[[2018,8,8]]},"references-count":25,"journal-issue":{"issue":"5","published-print":{"date-parts":[[2019,3,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/bty696","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2019,3,1]]},"published":{"date-parts":[[2018,8,8]]}}}