{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,9,26]],"date-time":"2025-09-26T00:12:02Z","timestamp":1758845522636},"reference-count":54,"publisher":"Springer Science and Business Media LLC","issue":"1","license":[{"start":{"date-parts":[[2020,8,28]],"date-time":"2020-08-28T00:00:00Z","timestamp":1598572800000},"content-version":"tdm","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"},{"start":{"date-parts":[[2020,8,28]],"date-time":"2020-08-28T00:00:00Z","timestamp":1598572800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"}],"content-domain":{"domain":["link.springer.com"],"crossmark-restriction":false},"short-container-title":["BMC Bioinformatics"],"published-print":{"date-parts":[[2020,12]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:sec><jats:title>Background<\/jats:title><jats:p>As the barriers to incorporating RNA sequencing (RNA-Seq) into biomedical studies continue to decrease, the complexity and size of RNA-Seq experiments are rapidly growing. Paired, longitudinal, and other correlated designs are becoming commonplace, and these studies offer immense potential for understanding how transcriptional changes within an individual over time differ depending on treatment or environmental conditions. While several methods have been proposed for dealing with repeated measures within RNA-Seq analyses, they are either restricted to handling only paired measurements, can only test for differences between two groups, and\/or have issues with maintaining nominal false positive and false discovery rates. In this work, we propose a Bayesian hierarchical negative binomial generalized linear mixed model framework that can flexibly model RNA-Seq counts from studies with arbitrarily many repeated observations, can include covariates, and also maintains nominal false positive and false discovery rates in its posterior inference.<\/jats:p><\/jats:sec><jats:sec><jats:title>Results<\/jats:title><jats:p>In simulation studies, we showed that our proposed method (MCMSeq) best combines high statistical power (i.e. sensitivity or recall) with maintenance of nominal false positive and false discovery rates compared the other available strategies, especially at the smaller sample sizes investigated. This behavior was then replicated in an application to real RNA-Seq data where MCMSeq was able to find previously reported genes associated with tuberculosis infection in a cohort with longitudinal measurements.<\/jats:p><\/jats:sec><jats:sec><jats:title>Conclusions<\/jats:title><jats:p>Failing to account for repeated measurements when analyzing RNA-Seq experiments can result in significantly inflated false positive and false discovery rates. Of the methods we investigated, whether they model RNA-Seq counts directly or worked on transformed values, the Bayesian hierarchical model implemented in the mcmseq R package (available at<jats:ext-link xmlns:xlink=\"http:\/\/www.w3.org\/1999\/xlink\" ext-link-type=\"uri\" xlink:href=\"https:\/\/github.com\/stop-pre16\/mcmseq\">https:\/\/github.com\/stop-pre16\/mcmseq<\/jats:ext-link>) best combined sensitivity and nominal error rate control.<\/jats:p><\/jats:sec>","DOI":"10.1186\/s12859-020-03715-y","type":"journal-article","created":{"date-parts":[[2020,8,28]],"date-time":"2020-08-28T14:02:59Z","timestamp":1598623379000},"update-policy":"http:\/\/dx.doi.org\/10.1007\/springer_crossmark_policy","source":"Crossref","is-referenced-by-count":11,"title":["MCMSeq: Bayesian hierarchical modeling of clustered and repeated measures RNA sequencing experiments"],"prefix":"10.1186","volume":"21","author":[{"given":"Brian E.","family":"Vestal","sequence":"first","affiliation":[]},{"given":"Camille M.","family":"Moore","sequence":"additional","affiliation":[]},{"given":"Elizabeth","family":"Wynn","sequence":"additional","affiliation":[]},{"given":"Laura","family":"Saba","sequence":"additional","affiliation":[]},{"given":"Tasha","family":"Fingerlin","sequence":"additional","affiliation":[]},{"given":"Katerina","family":"Kechris","sequence":"additional","affiliation":[]}],"member":"297","published-online":{"date-parts":[[2020,8,28]]},"reference":[{"issue":"1","key":"3715_CR1","doi-asserted-by":"crossref","first-page":"2308","DOI":"10.1038\/s41467-018-04579-w","volume":"9","author":"A Singhania","year":"2018","unstructured":"Singhania A, Verma R, Graham CM, Lee J, Tran T, Richardson M, Lecine P, Leissner P, Berry MP, Wilkinson RJ, et al. A modular transcriptional signature identifies phenotypic heterogeneity of human tuberculosis infection. Nat Commun. 2018; 9(1):2308.","journal-title":"Nat Commun"},{"issue":"9","key":"3715_CR2","doi-asserted-by":"crossref","first-page":"1007290","DOI":"10.1371\/journal.ppat.1007290","volume":"14","author":"BR Rosenberg","year":"2018","unstructured":"Rosenberg BR, Depla M, Freije CA, Gaucher D, Mazouz S, Boisvert M, B\u00e9dard N, Bruneau J, Rice CM, Shoukry NH. Longitudinal transcriptomic characterization of the immune response to acute hepatitis c virus infection in patients with spontaneous viral clearance. PLoS Pathog. 2018; 14(9):1007290.","journal-title":"PLoS Pathog"},{"issue":"2","key":"3715_CR3","first-page":"87","volume":"15","author":"S Cui","year":"2016","unstructured":"Cui S, Ji T, Li J, Cheng J, Qiu J. What if we ignore the random effects when analyzing rna-seq data in a multifactor experiment. SStat Appl Genet Mol Biol. 2016; 15(2):87\u2013105.","journal-title":"SStat Appl Genet Mol Biol"},{"key":"3715_CR4","volume-title":"Applied longitudinal analysis, vol. 998","author":"GM Fitzmaurice","year":"2012","unstructured":"Fitzmaurice GM, Laird NM, Ware JH. Applied longitudinal analysis, vol. 998. Hoboken: John Wiley & Sons; 2012."},{"issue":"10","key":"3715_CR5","doi-asserted-by":"crossref","first-page":"4288","DOI":"10.1093\/nar\/gks042","volume":"40","author":"DJ McCarthy","year":"2012","unstructured":"McCarthy DJ, Chen Y, Smyth GK. Differential expression analysis of multifactor rna-seq experiments with respect to biological variation. Nucleic Acids Res. 2012; 40(10):4288\u201397.","journal-title":"Nucleic Acids Res"},{"key":"3715_CR6","doi-asserted-by":"crossref","unstructured":"Harrell FE. Regression modeling strategies: with applications to linear models, logistic and ordinal regression, and survival analysis: Springer; 2015.","DOI":"10.1007\/978-3-319-19425-7"},{"issue":"21","key":"3715_CR7","doi-asserted-by":"crossref","first-page":"2881","DOI":"10.1093\/bioinformatics\/btm453","volume":"23","author":"MD Robinson","year":"2007","unstructured":"Robinson MD, Smyth GK. Moderated statistical tests for assessing differences in tag abundance. Bioinformatics. 2007; 23(21):2881\u20137.","journal-title":"Bioinformatics"},{"issue":"2","key":"3715_CR8","doi-asserted-by":"crossref","first-page":"321","DOI":"10.1093\/biostatistics\/kxm030","volume":"9","author":"MD Robinson","year":"2007","unstructured":"Robinson MD, Smyth GK. Small-sample estimation of negative binomial dispersion, with applications to sage data. Biostatistics. 2007; 9(2):321\u201332.","journal-title":"Biostatistics"},{"issue":"3","key":"3715_CR9","doi-asserted-by":"crossref","first-page":"863","DOI":"10.2307\/2532104","volume":"46","author":"WW Piegorsch","year":"1990","unstructured":"Piegorsch WW. Maximum likelihood estimation for the negative binomial dispersion parameter. Biometrics. 1990; 46(3):863\u20137.","journal-title":"Biometrics"},{"issue":"1","key":"3715_CR10","doi-asserted-by":"crossref","first-page":"309","DOI":"10.2307\/2532055","volume":"45","author":"SJ Clark","year":"1989","unstructured":"Clark SJ, Perry JN. Estimation of the negative binomial parameter \u03ba by maximum quasi-likelihood. Biometrics. 1989; 45(1):309\u201316.","journal-title":"Biometrics"},{"issue":"3\/4","key":"3715_CR11","doi-asserted-by":"crossref","first-page":"358","DOI":"10.2307\/2332388","volume":"37","author":"FJ Anscombe","year":"1950","unstructured":"Anscombe FJ. Sampling theory of the negative binomial and logarithmic series distributions. Biometrika. 1950; 37(3\/4):358\u201382.","journal-title":"Biometrika"},{"issue":"1","key":"3715_CR12","doi-asserted-by":"crossref","first-page":"109","DOI":"10.2307\/2530749","volume":"40","author":"LJ Willson","year":"1984","unstructured":"Willson LJ, Folks JL, Young J. Multistage estimation compared with fixed-sample-size estimation of the negative binomial parameter k. Biometrics. 1984; 40(1):109\u201317.","journal-title":"Biometrics"},{"issue":"1\/2","key":"3715_CR13","doi-asserted-by":"crossref","first-page":"37","DOI":"10.2307\/2333044","volume":"45","author":"C Bliss","year":"1958","unstructured":"Bliss C, Owen A. Negative binomial distributions with a common k. Biometrika. 1958; 45(1\/2):37\u201358.","journal-title":"Biometrika"},{"issue":"1","key":"3715_CR14","doi-asserted-by":"crossref","first-page":"139","DOI":"10.1093\/bioinformatics\/btp616","volume":"26","author":"MD Robinson","year":"2010","unstructured":"Robinson MD, McCarthy DJ, Smyth GK. Edger: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010; 26(1):139\u201340.","journal-title":"Bioinformatics"},{"issue":"12","key":"3715_CR15","doi-asserted-by":"crossref","first-page":"550","DOI":"10.1186\/s13059-014-0550-8","volume":"15","author":"MI Love","year":"2014","unstructured":"Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for rna-seq data with deseq2. Genome Biol. 2014; 15(12):550.","journal-title":"Genome Biol"},{"issue":"1","key":"3715_CR16","doi-asserted-by":"crossref","first-page":"91","DOI":"10.1186\/1471-2105-14-91","volume":"14","author":"C Soneson","year":"2013","unstructured":"Soneson C, Delorenzi M. A comparison of methods for differential expression analysis of rna-seq data. BMC Bioinformatics. 2013; 14(1):91.","journal-title":"BMC Bioinformatics"},{"issue":"421","key":"3715_CR17","doi-asserted-by":"crossref","first-page":"9","DOI":"10.1080\/01621459.1993.10594284","volume":"88","author":"NE Breslow","year":"1993","unstructured":"Breslow NE, Clayton DG. Approximate inference in generalized linear mixed models. J Am Stat Assoc. 1993; 88(421):9\u201325.","journal-title":"J Am Stat Assoc"},{"key":"3715_CR18","unstructured":"Rocke D. Excess false positives in negative-binomial based analysis of data from rna-seq experiments. http:\/\/dmrocke.ucdavis.edu\/papers\/RNA-Seq%20UVA%20Seminar2016-03.pdf. Accessed 15 Aug 2019."},{"issue":"11","key":"3715_CR19","doi-asserted-by":"crossref","first-page":"106","DOI":"10.1093\/nar\/gkx204","volume":"45","author":"S Sun","year":"2017","unstructured":"Sun S, Hood M, Scott L, Peng Q, Mukherjee S, Tung J, Zhou X. Differential expression analysis for rnaseq using poisson mixed models. Nucleic Acids Res. 2017; 45(11):106.","journal-title":"Nucleic Acids Res"},{"issue":"9","key":"3715_CR20","doi-asserted-by":"crossref","first-page":"1577","DOI":"10.1002\/sim.6449","volume":"34","author":"H Zhang","year":"2015","unstructured":"Zhang H, Xu J, Jiang N, Hu X, Luo Z. Plnseq: a multivariate poisson lognormal distribution for high-throughput matched rna-sequencing read count data. Stat Med. 2015; 34(9):1577\u201389.","journal-title":"Stat Med"},{"issue":"1","key":"3715_CR21","doi-asserted-by":"crossref","first-page":"248","DOI":"10.1186\/s12859-016-1111-9","volume":"17","author":"G Kang","year":"2016","unstructured":"Kang G, Du L, Zhang H. multide: a dimension reduced model based statistical method for differential expression analysis using rna-sequencing data with multiple treatment conditions. BMC Bioinformatics. 2016; 17(1):248.","journal-title":"BMC Bioinformatics"},{"issue":"5","key":"3715_CR22","doi-asserted-by":"crossref","first-page":"787","DOI":"10.1093\/bioinformatics\/bty731","volume":"35","author":"Y Bian","year":"2019","unstructured":"Bian Y, He C, Hou J, Cheng J, Qiu J. PairedFB: a full hierarchical bayesian model for paired rna-seq data with heterogeneous treatment effects. Bioinformatics. 2019; 35(5):787\u201397.","journal-title":"Bioinformatics"},{"issue":"1","key":"3715_CR23","doi-asserted-by":"crossref","first-page":"116","DOI":"10.1186\/1471-2105-15-116","volume":"15","author":"MA van de Wiel","year":"2014","unstructured":"van de Wiel MA, Neerincx M, Buffart TE, Sie D, Verheul HM. Shrinkbayes: a versatile r-package for analysis of count-based sequencing data in complex study designs. BMC Bioinformatics. 2014; 15(1):116.","journal-title":"BMC Bioinformatics"},{"key":"3715_CR24","doi-asserted-by":"crossref","DOI":"10.1007\/978-3-319-07212-8","volume-title":"Statistical analysis of next generation sequencing data","author":"S Datta","year":"2014","unstructured":"Datta S, Nettleton D. Statistical analysis of next generation sequencing data. New York: Springer; 2014."},{"issue":"4","key":"3715_CR25","doi-asserted-by":"crossref","first-page":"589","DOI":"10.1093\/biostatistics\/kxx005","volume":"18","author":"D Agniel","year":"2017","unstructured":"Agniel D, Hejblum BP. Variance component score test for time-course gene set analysis of longitudinal rna-seq data. Biostatistics. 2017; 18(4):589\u2013604.","journal-title":"Biostatistics"},{"issue":"7","key":"3715_CR26","doi-asserted-by":"crossref","first-page":"47","DOI":"10.1093\/nar\/gkv007","volume":"43","author":"ME Ritchie","year":"2015","unstructured":"Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. Limma powers differential expression analyses for rna-sequencing and microarray studies. Nucleic Acids Res. 2015; 43(7):47.","journal-title":"Nucleic Acids Res"},{"issue":"2","key":"3715_CR27","doi-asserted-by":"crossref","first-page":"29","DOI":"10.1186\/gb-2014-15-2-r29","volume":"15","author":"CW Law","year":"2014","unstructured":"Law CW, Chen Y, Shi W, Smyth GK. Voom: Precision weights unlock linear model analysis tools for rna-seq read counts. Genome Biol. 2014; 15(2):29.","journal-title":"Genome Biol"},{"key":"3715_CR28","volume-title":"Bioinformatics and Computational Biology Solutions Using R and Bioconductor","author":"GK Smyth","year":"2005","unstructured":"Smyth GK. Limma: linear models for microarray data. In: Bioinformatics and Computational Biology Solutions Using R and Bioconductor. New York: Springer: 2005. p. 397\u2013420."},{"issue":"1","key":"3715_CR29","doi-asserted-by":"crossref","first-page":"57","DOI":"10.1023\/A:1018509429360","volume":"7","author":"D Gamerman","year":"1997","unstructured":"Gamerman D. Sampling from the posterior distribution in generalized linear mixed models. Stat Comput. 1997; 7(1):57\u201368.","journal-title":"Stat Comput"},{"issue":"1","key":"3715_CR30","doi-asserted-by":"crossref","first-page":"143","DOI":"10.1186\/s12859-017-1539-6","volume":"18","author":"P Rudra","year":"2017","unstructured":"Rudra P, Shi WJ, Vestal B, Russell PH, Odell A, Dowell RD, Radcliffe RA, Saba LM, Kechris K. Model based heritability scores for high-throughput sequencing data. BMC Bioinformatics. 2017; 18(1):143.","journal-title":"BMC Bioinformatics"},{"key":"3715_CR31","volume-title":"Evaluating the accuracy of sampling-based approaches to the calculation of posterior moments, vol. 196","author":"J Geweke","year":"1991","unstructured":"Geweke J, et al. Evaluating the accuracy of sampling-based approaches to the calculation of posterior moments, vol. 196. Minneapolis: Federal Reserve Bank of Minneapolis, Research Department Minneapolis, MN; 1991."},{"issue":"1","key":"3715_CR32","doi-asserted-by":"crossref","first-page":"289","DOI":"10.1111\/j.2517-6161.1995.tb02031.x","volume":"57","author":"Y Benjamini","year":"1995","unstructured":"Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995; 57(1):289\u2013300.","journal-title":"J R Stat Soc Ser B Methodol"},{"key":"3715_CR33","volume-title":"Bayesian inference in statistical analysis, vol. 40","author":"GE Box","year":"2011","unstructured":"Box GE, Tiao GC. Bayesian inference in statistical analysis, vol. 40. Hoboken: John Wiley & Sons; 2011."},{"key":"3715_CR34","doi-asserted-by":"crossref","DOI":"10.1002\/9781119942412","volume-title":"Bayesian biostatistics","author":"E Lesaffre","year":"2012","unstructured":"Lesaffre E, Lawson AB. Bayesian biostatistics. New York: Wiley; 2012."},{"key":"3715_CR35","volume-title":"CORE and Institut de Statistique","author":"BA Turlach","year":"1993","unstructured":"Turlach BA. Bandwidth selection in kernel density estimation: A review. In: CORE and Institut de Statistique. Berlin: Springer: 1993."},{"issue":"2","key":"3715_CR36","doi-asserted-by":"crossref","first-page":"479","DOI":"10.1111\/j.1541-0420.2007.00895.x","volume":"64","author":"JS Morris","year":"2008","unstructured":"Morris JS, Brown PJ, Herrick RC, Baggerly KA, Coombes KR. Bayesian analysis of mass spectrometry proteomic data using wavelet-based functional mixed models. Biometrics. 2008; 64(2):479\u201389.","journal-title":"Biometrics"},{"issue":"397","key":"3715_CR37","doi-asserted-by":"crossref","first-page":"106","DOI":"10.1080\/01621459.1987.10478396","volume":"82","author":"G Casella","year":"1987","unstructured":"Casella G, Berger RL. Reconciling bayesian and frequentist evidence in the one-sided testing problem. J Am Stat Assoc. 1987; 82(397):106\u201311.","journal-title":"J Am Stat Assoc"},{"key":"3715_CR38","doi-asserted-by":"publisher","first-page":"3188","DOI":"10.7287\/peerj.preprints.3188v1","volume":"5","author":"D Eddelbuettel","year":"2017","unstructured":"Eddelbuettel D, Balamuta JJ. Extending extitR with extitC++: a brief introduction to extitRcpp. PeerJ Prepr. 2017; 5:3188\u20131. https:\/\/doi.org\/10.7287\/peerj.preprints.3188v1.","journal-title":"PeerJ Prepr"},{"key":"3715_CR39","unstructured":"Allaire J, Francois R, Ushey K, Vandenbrouck G, Geelnard M, Intel. RcppParallel: Parallel Programming Tools for \u2019Rcpp\u2019. R package version 4.4.1. 2018."},{"key":"3715_CR40","doi-asserted-by":"crossref","first-page":"1054","DOI":"10.1016\/j.csda.2013.02.005","volume":"71","author":"D Eddelbuettel","year":"2014","unstructured":"Eddelbuettel D, Sanderson C. Rcpparmadillo: Accelerating r with high-performance c++ linear algebra. Comput Stat Data Anal. 2014; 71:1054\u201363.","journal-title":"Comput Stat Data Anal"},{"key":"3715_CR41","doi-asserted-by":"crossref","unstructured":"Carpenter B, Gelman A, Hoffman MD, Lee D, Goodrich B, Betancourt M, Brubaker M, Guo J, Li P, Riddell A. Stan: A probabilistic programming language. J Stat Softw. 2017;76(1).","DOI":"10.18637\/jss.v076.i01"},{"issue":"25","key":"3715_CR42","doi-asserted-by":"crossref","first-page":"3049","DOI":"10.1002\/sim.3680","volume":"28","author":"D Lunn","year":"2009","unstructured":"Lunn D, Spiegelhalter D, Thomas A, Best N. The bugs project: Evolution, critique and future directions. Stat Med. 2009; 28(25):3049\u201367.","journal-title":"Stat Med"},{"key":"3715_CR43","unstructured":"Team S. Rstan: the r interface to stan. R package version. 2016;2(1)."},{"key":"3715_CR44","unstructured":"Zhang W. cplm: Monte carlo em algorithms and bayesian methods for fitting tweedie compound poisson linear models. R package version 0.2-1. 2011."},{"issue":"13","key":"3715_CR45","doi-asserted-by":"publisher","first-page":"1","DOI":"10.18637\/jss.v082.i13","volume":"82","author":"A Kuznetsova","year":"2017","unstructured":"Kuznetsova A, Brockhoff PB, Christensen RHB. lmerTest package: Tests in linear mixed effects models. J Stat Softw. 2017; 82(13):1\u201326. https:\/\/doi.org\/10.18637\/jss.v082.i13.","journal-title":"J Stat Softw"},{"key":"3715_CR46","doi-asserted-by":"crossref","first-page":"233","DOI":"10.1080\/10556788.2011.597854","volume":"27","author":"DA Fournier","year":"2012","unstructured":"Fournier DA, Skaug HJ, Ancheta J, Ianelli J, Magnusson A, Maunder MN, Nielsen A, Sibert J. AD Model Builder: using automatic differentiation for statistical inference of highly parameterized complex nonlinear models. Optim Methods Softw. 2012; 27:233\u201349.","journal-title":"Optim Methods Softw"},{"issue":"6","key":"3715_CR47","doi-asserted-by":"crossref","first-page":"1873","DOI":"10.1183\/13993003.02121-2015","volume":"47","author":"S Blankley","year":"2016","unstructured":"Blankley S, Graham CM, Levin J, Turner J, Berry MP, Bloom CI, Xu Z, Pascual V, Banchereau J, Chaussabel D, et al. A 380-gene meta-signature of active tuberculosis compared with healthy controls. Eur Respir J. 2016; 47(6):1873\u20136.","journal-title":"Eur Respir J"},{"issue":"10035","key":"3715_CR48","doi-asserted-by":"crossref","first-page":"2312","DOI":"10.1016\/S0140-6736(15)01316-1","volume":"387","author":"DE Zak","year":"2016","unstructured":"Zak DE, Penn-Nicholson A, Scriba TJ, Thompson E, Suliman S, Amon LM, Mahomed H, Erasmus M, Whatney W, Hussey GD, et al. A blood rna signature for tuberculosis disease risk: a prospective cohort study. Lancet. 2016; 387(10035):2312\u201322.","journal-title":"Lancet"},{"issue":"10","key":"3715_CR49","doi-asserted-by":"crossref","first-page":"1001538","DOI":"10.1371\/journal.pmed.1001538","volume":"10","author":"M Kaforou","year":"2013","unstructured":"Kaforou M, Wright VJ, Oni T, French N, Anderson ST, Bangani N, Banwell CM, Brent AJ, Crampin AC, Dockrell HM, et al. Detection of tuberculosis in hiv-infected and-uninfected african adults using whole blood rna expression signatures: a case-control study. PLoS Med. 2013; 10(10):1001538.","journal-title":"PLoS Med"},{"issue":"9","key":"3715_CR50","doi-asserted-by":"crossref","first-page":"73230","DOI":"10.1371\/journal.pone.0073230","volume":"8","author":"SA Joosten","year":"2013","unstructured":"Joosten SA, Fletcher HA, Ottenhoff TH. A helicopter perspective on tb biomarkers: pathway and process based analysis of gene expression data provides new insight into tb pathogenesis. PloS ONE. 2013; 8(9):73230.","journal-title":"PloS ONE"},{"issue":"1","key":"3715_CR51","doi-asserted-by":"crossref","first-page":"15","DOI":"10.1038\/gene.2010.51","volume":"12","author":"J Maertzdorf","year":"2011","unstructured":"Maertzdorf J, Repsilber D, Parida SK, Stanley K, Roberts T, Black G, Walzl G, Kaufmann SH. Human gene expression profiles of susceptibility and resistance in tuberculosis. Genes Immun. 2011; 12(1):15.","journal-title":"Genes Immun"},{"key":"3715_CR52","doi-asserted-by":"crossref","unstructured":"Nguyen Y, Nettleton D. rmrnaseq Differential expression analysis for repeated-measures rna-seq data. Bioinformatics. 2020.","DOI":"10.1093\/bioinformatics\/btaa525"},{"key":"3715_CR53","volume-title":"Markov Chain Montecarlo in Practice","author":"D Clayton","year":"1996","unstructured":"Clayton D. Generalized linear mixed models In: Gilks W, Richardson S, Spiegelhalter D, editors. Markov Chain Montecarlo in Practice. London: Chapman and Hall: 1996. p. 275\u2013302."},{"issue":"2","key":"3715_CR54","doi-asserted-by":"crossref","first-page":"630","DOI":"10.1177\/0962280212462859","volume":"25","author":"M Chavance","year":"2016","unstructured":"Chavance M, Escolano S. Misspecification of the covariance structure in generalized linear mixed models. Stat Methods Med Res. 2016; 25(2):630\u201343.","journal-title":"Stat Methods Med Res"}],"container-title":["BMC Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1186\/s12859-020-03715-y.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/article\/10.1186\/s12859-020-03715-y\/fulltext.html","content-type":"text\/html","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1186\/s12859-020-03715-y.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2024,8,12]],"date-time":"2024-08-12T19:34:22Z","timestamp":1723491262000},"score":1,"resource":{"primary":{"URL":"https:\/\/bmcbioinformatics.biomedcentral.com\/articles\/10.1186\/s12859-020-03715-y"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2020,8,28]]},"references-count":54,"journal-issue":{"issue":"1","published-print":{"date-parts":[[2020,12]]}},"alternative-id":["3715"],"URL":"https:\/\/doi.org\/10.1186\/s12859-020-03715-y","relation":{},"ISSN":["1471-2105"],"issn-type":[{"type":"electronic","value":"1471-2105"}],"subject":[],"published":{"date-parts":[[2020,8,28]]},"assertion":[{"value":"8 May 2020","order":1,"name":"received","label":"Received","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"18 August 2020","order":2,"name":"accepted","label":"Accepted","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"28 August 2020","order":3,"name":"first_online","label":"First Online","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"No ethics approval was required for this study. All data analyzed in this manuscript was either simulated or downloaded from publicly available sources.","order":1,"name":"Ethics","group":{"name":"EthicsHeading","label":"Ethics approval and consent to participate"}},{"value":"Not Applicable.","order":2,"name":"Ethics","group":{"name":"EthicsHeading","label":"Consent for publication"}},{"value":"Katerina Kechris was a member of the editorial board of <i>BMC Bioinformatics<\/i> from 2015-2019. All other authors declare that they have no competing interests.","order":3,"name":"Ethics","group":{"name":"EthicsHeading","label":"Competing interests"}}],"article-number":"375"}}