title: Poor fit to the multispecies coalescent is widely detectable in empirical data.
abstract: Model checking is a critical part of Bayesian data analysis, yet it remains largely unused in systematic studies. Phylogeny estimation has recently moved into an era of increasingly complex models that simultaneously account for multiple evolutionary processes, the statistical fit of these models to the data has rarely been tested. Here we develop a posterior predictive simulation-based model check for a commonly used multispecies coalescent model, implemented in *BEAST, and apply it to 25 published data sets. We show that poor model fit is detectable in the majority of data sets; that this poor fit can mislead phylogenetic estimation; and that in some cases it stems from processes of inherent interest to systematists. We suggest that as systematists scale up to phylogenomic data sets, which will be subject to a heterogeneous array of evolutionary processes, critically evaluating the fit of models to data is an analytical step that can no longer be ignored.
authors.affiliation: Department of Biological Sciences, Louisiana State University, Baton Rouge, LA 70803, USA; and Department of Evolution, Ecology & Organismal Biology, Ohio State University, Columbus, OH 43210, USA.
authors.forename: Noah M
authors.initials: NM
authors.lastname: Reid
authors.forename: Sarah M
authors.initials: SM
authors.lastname: Hird
authors.forename: Jeremy M
authors.initials: JM
authors.lastname: Brown
authors.forename: Tara A
authors.initials: TA
authors.lastname: Pelletier
authors.forename: John D
authors.initials: JD
authors.lastname: McVay
authors.forename: Jordan D
authors.initials: JD
authors.lastname: Satler
authors.forename: Bryan C
authors.initials: BC
authors.lastname: Carstens
mesh_terms: D002965:Classification; D003198:Computer Simulation / Q000592:standards*; D056890:Eukaryota / Q000145:classification / Q000235:genetics; D008954:Models, Biological; D010802:Phylogeny