package dr.evomodelxml.tree; import dr.evolution.tree.Tree; import dr.evolution.tree.TreeTrait; import dr.evolution.tree.TreeTraitProvider; import dr.evolution.util.Taxa; import dr.evolution.util.TaxonList; import dr.evomodel.tree.AncestralTrait; import dr.evomodel.tree.TreeModel; import dr.xml.*; /** */ public class AncestralTraitParser extends AbstractXMLObjectParser { public static final String ANCESTRAL_TRAIT = "ancestralTrait"; public static final String ANCESTRAL_STATE = "ancestralState"; public static final String NAME = "name"; public static final String MRCA = "mrca"; public static final String TRAIT_NAME = "traitName"; public static final String STATES = "states"; public String getParserName() { return ANCESTRAL_TRAIT; } public String[] getParserNames() { // provide a synonym to maintain backwards compatibility return new String[]{ getParserName(), ANCESTRAL_STATE }; } public Object parseXMLObject(XMLObject xo) throws XMLParseException { String traitName = xo.getAttribute(TRAIT_NAME, STATES); String name = xo.getAttribute(NAME, traitName); Tree tree = (Tree) xo.getChild(Tree.class); TreeTraitProvider treeTraitProvider = (TreeTraitProvider) xo.getChild(TreeTraitProvider.class); TaxonList taxa = null; if (xo.hasChildNamed(MRCA)) { taxa = (TaxonList) xo.getElementFirstChild(MRCA); } TreeTrait trait = treeTraitProvider.getTreeTrait(traitName); if (trait == null) { throw new XMLParseException("A trait called, " + traitName + ", was not available from the TreeTraitProvider supplied to " + getParserName() + (xo.hasId() ? ", with ID " + xo.getId() : "")); } try { return new AncestralTrait(name, trait, tree, taxa); } catch (Tree.MissingTaxonException mte) { throw new XMLParseException("Taxon, " + mte + ", in " + getParserName() + "was not found in the tree."); } } //************************************************************************ // AbstractXMLObjectParser implementation //************************************************************************ public String getParserDescription() { return "A statistic that has as its value the height of the most recent common ancestor of a set of taxa in a given tree"; } public Class getReturnType() { return AncestralTrait.class; } public XMLSyntaxRule[] getSyntaxRules() { return rules; } private XMLSyntaxRule[] rules = new XMLSyntaxRule[]{ new StringAttributeRule(NAME, "A name for this statistic primarily for the purposes of logging", true), new StringAttributeRule(TRAIT_NAME, "The name of the trait to log", true), new ElementRule(TreeModel.class), new ElementRule(TreeTraitProvider.class), new ElementRule(MRCA, new XMLSyntaxRule[]{new ElementRule(Taxa.class)}, "The MRCA to reconstruct the trait at (default root node)", true) }; }