package dr.evomodelxml.coalescent; import dr.evomodel.coalescent.BMPriorLikelihood; import dr.evomodel.coalescent.VariableDemographicModel; import dr.inference.distribution.ParametricDistributionModel; import dr.inference.model.Parameter; import dr.xml.*; /** */ public class BMPriorLikelihoodParser extends AbstractXMLObjectParser { //TODO Parser is still in a bad state public static final String BM = "BrownianMotion"; public static final String LOG_SPACE = "logUnits"; //static final String NORMALIZE = "normalize"; //static final String DATA = "data"; // static final String TIMES = "times"; // public static final String MEAN = "mean"; public static final String SIGMA = "sigma"; // public static final String LAMBDA = "lambda"; public String getParserDescription() { return ""; } public Class getReturnType() { return BMPriorLikelihood.class; } public String getParserName() { return BM; } private Parameter getParam(XMLObject xo, String name) throws XMLParseException { final XMLObject object = xo.getChild(name); // optional if (object == null) { return null; } final Object child = object.getChild(0); if (child instanceof Parameter) { return (Parameter) child; } double x = object.getDoubleChild(0); return new Parameter.Default(x); } public Object parseXMLObject(XMLObject xo) throws XMLParseException { // Parameter mean = getParam(xo, MEAN); Parameter sigma = getParam(xo, SIGMA); // Parameter lambda = getParam(xo, LAMBDA); final boolean logSpace = xo.getAttribute(LOG_SPACE, false); // final boolean normalize = xo.getAttribute(NORMALIZE, false); VariableDemographicModel m = (VariableDemographicModel) xo.getChild(VariableDemographicModel.class); ParametricDistributionModel popMeanPrior = (ParametricDistributionModel) xo.getChild(ParametricDistributionModel.class); return new BMPriorLikelihood(sigma, m, logSpace, popMeanPrior); } public XMLSyntaxRule[] getSyntaxRules() { return new XMLSyntaxRule[]{ AttributeRule.newBooleanRule(LOG_SPACE, true), new XORRule( new ElementRule(SIGMA, Double.class), new ElementRule(SIGMA, Parameter.class) ), // you can't have a XOR (b AND c), yikes // make all optional and check in parser new ElementRule(VariableDemographicModel.class), // new ElementRule(DATA, new XMLSyntaxRule[]{new ElementRule(Statistic.class)}, true), // new ElementRule(TIMES, new XMLSyntaxRule[]{new ElementRule(Statistic.class)}, true) }; } }