/*
* ThreadedCompoundLikelihoodParser.java
*
* Copyright (c) 2002-2015 Alexei Drummond, Andrew Rambaut and Marc Suchard
*
* This file is part of BEAST.
* See the NOTICE file distributed with this work for additional
* information regarding copyright ownership and licensing.
*
* BEAST is free software; you can redistribute it and/or modify
* it under the terms of the GNU Lesser General Public License as
* published by the Free Software Foundation; either version 2
* of the License, or (at your option) any later version.
*
* BEAST is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public
* License along with BEAST; if not, write to the
* Free Software Foundation, Inc., 51 Franklin St, Fifth Floor,
* Boston, MA 02110-1301 USA
*/
package dr.inferencexml.model;
import dr.inference.model.Likelihood;
import dr.inference.model.ThreadedCompoundLikelihood;
import dr.xml.*;
/**
*
*/
public class ThreadedCompoundLikelihoodParser extends AbstractXMLObjectParser {
public static final String THREADED_COMPOUND_LIKELIHOOD = "threadedCompoundLikelihood";
public static final String WEIGHT = "robustWeight";
public String getParserName() {
return THREADED_COMPOUND_LIKELIHOOD;
}
public String[] getParserNames() {
return new String[]{getParserName()};
}
public Object parseXMLObject(XMLObject xo) throws XMLParseException {
ThreadedCompoundLikelihood compoundLikelihood = new ThreadedCompoundLikelihood();
for (int i = 0; i < xo.getChildCount(); i++) {
if (xo.getChild(i) instanceof Likelihood) {
compoundLikelihood.addLikelihood((Likelihood) xo.getChild(i));
} else {
Object rogueElement = xo.getChild(i);
throw new XMLParseException("An element (" + rogueElement + ") which is not a likelihood has been added to a " + THREADED_COMPOUND_LIKELIHOOD + " element");
}
}
double weight = xo.getAttribute(WEIGHT, 0.0);
if (weight < 0)
throw new XMLParseException("Robust weight must be non-negative.");
compoundLikelihood.setWeightFactor(Math.exp(-weight));
return compoundLikelihood;
}
//************************************************************************
// AbstractXMLObjectParser implementation
//************************************************************************
public String getParserDescription() {
return "A likelihood function which is simply the product of its component likelihood functions.";
}
public XMLSyntaxRule[] getSyntaxRules() {
return rules;
}
private final XMLSyntaxRule[] rules = {
new ElementRule(Likelihood.class, 1, Integer.MAX_VALUE),
AttributeRule.newDoubleRule(WEIGHT, true),
};
public Class getReturnType() {
return ThreadedCompoundLikelihood.class;
}
}