package dr.app.beauti.components.tipdatesampling; import dr.app.beauti.util.XMLWriter; import dr.app.beauti.types.TipDateSamplingType; import dr.app.beauti.generator.BaseComponentGenerator; import dr.app.beauti.options.BeautiOptions; import dr.evolution.util.Taxon; import dr.evolution.util.TaxonList; import dr.evoxml.TaxonParser; import dr.inference.model.ParameterParser; import dr.util.Attribute; import dr.xml.XMLParser; /** * @author Andrew Rambaut * @version $Id$ */ public class TipDateSamplingComponentGenerator extends BaseComponentGenerator { public TipDateSamplingComponentGenerator(final BeautiOptions options) { super(options); } public boolean usesInsertionPoint(final InsertionPoint point) { TipDateSamplingComponentOptions comp = (TipDateSamplingComponentOptions)options.getComponentOptions(TipDateSamplingComponentOptions.class); if (comp.tipDateSamplingType == TipDateSamplingType.NO_SAMPLING) { return false; } switch (point) { case IN_TREE_MODEL: case IN_FILE_LOG_PARAMETERS: return true; case AFTER_TREE_MODEL: case IN_MCMC_PRIOR: return comp.tipDateSamplingType == TipDateSamplingType.SAMPLE_JOINT; default: return false; } } protected void generate(final InsertionPoint point, final Object item, final XMLWriter writer) { TipDateSamplingComponentOptions comp = (TipDateSamplingComponentOptions)options.getComponentOptions(TipDateSamplingComponentOptions.class); TaxonList taxa = comp.getTaxonSet(); if (taxa == null) { // all taxa... taxa = options.taxonList; } switch (point) { case IN_TREE_MODEL: { for (int i = 0; i < taxa.getTaxonCount(); i++) { Taxon taxon = taxa.getTaxon(i); writer.writeOpenTag("leafHeight", new Attribute[]{ new Attribute.Default<String>(TaxonParser.TAXON, taxon.getId()), } ); writer.writeTag(ParameterParser.PARAMETER, new Attribute.Default<String>(XMLParser.ID, "age(" + taxon.getId() + ")"), true); writer.writeCloseTag("leafHeight"); } } break; case AFTER_TREE_MODEL: if (comp.tipDateSamplingType == TipDateSamplingType.SAMPLE_JOINT) { writer.writeOpenTag("compoundParameter", new Attribute[]{ new Attribute.Default<String>(XMLParser.ID, "treeModel.tipDates"), } ); for (int i = 0; i < taxa.getTaxonCount(); i++) { Taxon taxon = taxa.getTaxon(i); writer.writeIDref(ParameterParser.PARAMETER, "age(" + taxon.getId() + ")"); } writer.writeCloseTag("compoundParameter"); } break; case IN_MCMC_PRIOR: if (comp.tipDateSamplingType == TipDateSamplingType.SAMPLE_INDIVIDUALLY) { // nothing to do - individual parameter priors are written automatically } else if (comp.tipDateSamplingType == TipDateSamplingType.SAMPLE_JOINT) { } break; case IN_FILE_LOG_PARAMETERS: if (comp.tipDateSamplingType == TipDateSamplingType.SAMPLE_INDIVIDUALLY) { for (int i = 0; i < taxa.getTaxonCount(); i++) { Taxon taxon = taxa.getTaxon(i); writer.writeIDref(ParameterParser.PARAMETER, "age(" + taxon.getId() + ")"); } } else if (comp.tipDateSamplingType == TipDateSamplingType.SAMPLE_JOINT) { writer.writeIDref(ParameterParser.PARAMETER, "treeModel.tipDates"); } break; default: throw new IllegalArgumentException("This insertion point is not implemented for " + this.getClass().getName()); } } protected String getCommentLabel() { return "Tip date sampling"; } }