BEAST 2

Constraints Of Monophyly

Using BEAUti

Constraining a clade to be monphyletic is easiest done in via BEAUti from the “Priors” tab and selecting the “+” button at the bottom of the screen. A window is shown where you can specify the taxa of a clade and the name of the clade.

TaxonSetEditor

 

When closing, a new entry appears in the priors panel. Click the monophyletic checkbox to ensure the clade is forced to be monophyletic.

Monophyletic

 

By editing XML

Firstly create a distribution that contains the group you want to constrain to be monophyletic. By adding the distribution element directly to the XML inside the element with id “prior” like so:

<distribution id="ph.prior" spec="beast.math.distributions.MRCAPrior" tree="@Tree.t:xyz">
      <taxonset id="ph" spec="TaxonSet">
             <taxon id="Homo_sapiens" spec="Taxon"/>
             <taxon id="Pan" spec="Taxon"/>
      </taxonset>
</distribution>

Note 1: the tree constrained here is Tree.t:xyz — you need to change xyz so that it matches the name of the tree you want to constrain.
Note 2: if the taxon is already specified before (because it is already in another constraint, or when you start with a SNAPP analysis), you should refer to that taxon instead of creating a new one, like so: <taxon idref="Homo_sapiens"/>, so leave out the spec attribute and replace id with idref.

To constrain the taxonset to be monophyletic, add a monophyletic="true" attribute to the distribution:

<distribution id="ph.prior" monophyletic="true" spec="beast.math.distributions.MRCAPrior" tree="@Tree.t:xyz">
      <taxonset id="ph" spec="TaxonSet">
             <taxon id="Homo_sapiens" spec="Taxon"/>
             <taxon id="Pan" spec="Taxon"/>
      </taxonset>
</distribution>

If you add this to the log file then it displays the height of the ph clade in the log file.

  <log idref='ph.prior'/>

But you will have to provide a starting tree that is valid to begin with so the default random starting tree won’t work. The easiest way to do this is by using a randomtree as starting tree (the default in BEAUti) or to supply your own initial tree. See this tutorial for details of how to do this.

Rooting trees using an outgroup

You can use this technique to root the sampled trees using an outgroup. First create a distibution, as above, containing your ingroup taxa. Then follow the above instructions to keep your ingroup monophyletic with respect to the outgroup. Note that this is essentially a prior and you may wish to estimate the posterior probability that you outgroup is indeed an outgroup.

Bayesian evolutionary analysis by sampling trees

Served through Jekyll, customised theme based on the twentyfourteen wordpress theme.