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;
}
}