BAli---Phy

Release Notes

See NEWS file for more details, and for older changes.

• Speed & Memory
  • Speed: 0-300% faster.    Memory: 0-300% less memory.
  • [3.6] Speed up fixed-alignment analyses.
  • Site-compression for fixed-alignment partitions.
• Alignments
  • Allow analyzing 1-sequence and 2-sequence alignments.
  • [3.5] Automatically reconstruct ancestors for summary alignment.
• Large Trees
  • Compute likelihoods for very large trees (likelihood rescaling).
  • New treelength prior: don't prefer huge treelengths for trees with lots of leaves.
• General
  • Add unit tests and likelihood testsuite.
  • Allow selecting character sets from a file: sequences.fasta:100-240,300-900
  • Allow recovering from initial -Infinity log-likelihood in MCMC.
  • Install package manager bali-phy-pkg.
• Models
  • Priors:
    • Explicit priors, including default priors
    • Sample initial branch-lengths from prior.
    • substitution model: -S 'hky85[kappa=2]'
    • insertion/deletion model: -I 'rs07[log_rate~laplace[-4,0.707]]'
    • branch lengths: -B '~iid[num_branches[T],gamma[0.5,div[2,num_branches[T]]]]'
    • partition-specific scale factors: -R '1,3:~gamma[0.5,2]'
  • New syntax:
    • let[m=hky85,mixture[List[m,m,m+Rates.gamma]]] (user-defined variables)
    • function[w,gy94[omega=w]] (functions)
  • New models:
    • mixture[List[hky85,jc69]] (rates-across sites mixtures)
    • hky85+multi_rate[beta[2,3]] (rates-across sites with any distribution)
    • +fe for equal frequencies.
    • [3.5] +mut_sel or +mut_sel_aa to add selection on a model.
  • RNA stem (16-state) models
    • RNA models for fixed alignments. (Preliminary, because you have to pair the letters manually.)
    • RNA stem models: RNA.16a, gtr_sym+x2_sym+f, gtr+x2+mut_sel
  • Codon and triplet models
    • Improved codon models: gy94, mg94, f1x4, f3x4, gy94_ext, mg94_ext, x3, dNdS, etc.
    • dNdS mixtures with mg94 and fMutSel: m3[function[w,fMutSel[omega=w]]]
    • gtr+x3+dNdS+mut_sel (syntax for building up codon models piecewise)
  • Heterotachy/Covarion models
    • [3.5] Tuffley-Steel (1998) ON/OFF model.
    • [3.5] Huelsenbeck (2002) ON/OFF + Gamma rates-across-sites model.
    • [3.5] Galtier (2001) Gamma rate-switching model.
    • [3.5] Wang et al (2007) ON/OFF + Gamma rate-switching model.
• Summarization tools
  • Show 2D & 3D topology convergencs figures (MDS).
  • Reorganize logged statistics.
  • Compute alignment summaries over all runs, not just the first one.
  • [3.5] Rewrite bp-analyze in python.
• Help
  • Help for models, functions, distributions, and commands via bali-phy help topic.
  • Rewrite documentation for substitution models and insertion/deletion models.
• Tools
  • cut-range: allow reading alignments from multiple files.
  • alignment-distances: new tool, add accuracy and recal metrics
  • tree-tool: scales trees, prunes trees, computes diameter, etc.
  • alignment-thin: clean up options and man pages
  • bali-subsample: rename from subsample to avoid conflict with phylobayes.
  • [3.5] extract-ancestors: Extract and name ancestral sequences.
  • [3.5] summarize-ancestors: Construct summary alignment with ancestral sequences.
  • [3.5] tree-tool: add option to remove internal node names.
  • [3.5] Speed up tree drawing in bp-analyze.
  • [3.6] bp-analyze: show up to 4 chains in tree mixing plots.
  • [3.6] bp-analyze: allow analyzing bare Newick tree samples.
• Install
  • Transition homebrew package to brewsci/bio.
• Fixes
  • [3.5] Stabilize matrix exponential.
  • [3.5] Fix dirichlet process prior.
  • [3.5] Don't print out Newick trees rooted at a leaf.
  • [3.5] Compute posterior probabilities for MAP and greedy consensus.
  • [3.6] bali-phy: fix crash in MCMC when resampling alignments.
  • [3.6] bali-phy: allow sampling ancestral sequences with compressed, fixed alignments.
  • [3.6] bp-analyze: fix excessive cpu/memory usage when analyzing several chains.
  • [3.6] alignment-distances: fix NxN analysis.
  • [3.6.1] BES: fix overflow in likelihood.
  • [3.6.1] m8,m8a,m8a_test: fix help documentation.
  • [3.6.1] fix mixture distribution.