build_species_tree
Build a species or guide tree from a join formula and Imap, then write PREFIX.stree and PREFIX.nwk for use as BPP stree_file.
Instructions
Build the species (or guide) tree from a join formula (bpp-tree).
Writes PREFIX.stree (the species&tree block, with individual counts taken
from the Imap, plus any migration block) and PREFIX.nwk. Pass PREFIX.stree
to make_control_file as stree_file.
joins: comma-separated joins, each 'X+Y' or 'X+Y=name', e.g. 'chimp+bonobo=pan, pan+human, gorilla+pan_human'. A label like A_B refers to the clade of A and B. Tip names are the Imap's species names. Ask the user for the tree (or guide tree for delimitation); do not invent one.migration(MSC-M): bands 'SRC->DST, ...'. Mutually exclusive withintrogression. The control file then also needs a wprior.introgression(MSC-I): events 'DONOR->RECIP phi=0.1, ...'. The control file then also needs a phiprior. Look up those priors with lookup_docs; lint_control_file reports what is missing.
Read in the report: newick, taxa, species_counts,
species_and_tree_block, and warnings (e.g. ROOT_AUTO_JOINED: two
clades were joined at the root automatically) and errors.
server.diagram is an ASCII drawing of the tree. ALWAYS show the user the
diagram and ask them to confirm the topology before going on.
Input Schema
| Name | Required | Description | Default |
|---|---|---|---|
| imap | Yes | ||
| joins | Yes | ||
| migration | No | ||
| out_prefix | Yes | ||
| introgression | No |