Skip to main content
Dryad logo

Generalized hidden Markov models for phylogenetic comparative datasets

Citation

Boyko, James (2020), Generalized hidden Markov models for phylogenetic comparative datasets, Dryad, Dataset, https://doi.org/10.5061/dryad.vx0k6djpg

Abstract

  1. Hidden Markov models (HMM) have emerged as an important tool for understanding the evolution of characters that take on discrete states. Their flexibility and biological sensibility make them appealing for many phylogenetic comparative applications.
  2. Previously available packages placed unnecessary limits on the number of observed and hidden states that can be considered when estimating transition rates and inferring ancestral states on a phylogeny.
  3. To address these issues, we expanded the capabilities of the R package corHMM to handle n-state and n-character problems and provide users with a streamlined set of functions to create custom HMMs for any biological question of arbitrary complexity.
  4. We show that increasing the number of observed states increases the accuracy of ancestral state reconstruction. We also explore the conditions for when an HMM is most effective, finding that an HMM is an appropriate model when the degree of rate heterogeneity is moderate to high.
  5. Finally, we demonstrate the importance of these generalizations by reconstructing the phyllotaxy of the ancestral angiosperm flower. Partially contradicting previous results, we find the most likely state to be a whorled perianth, whorled androecium, whorled gynoecium. The difference between our analysis and previous studies was that our modeling explicitly allowed for the correlated evolution of several flower characters.

Usage Notes

CorHMMSimsData: includes std_rates (which compares 2,3,4 state models), state-process-diff (results of the state-dependent model fits), param-process-bias (results of the parameter-process bias model fits)

CorHMMCaseStudyDataAndResults: includes the data necessary and model fits from the case study

CorHMMScripts: the scripts needed to recreate our simulations, figures, and analysis