Examples with AlloppSpeciesBindings dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings used on opensource projects

Example 1 with AlloppSpeciesBindings

use of dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings in project beast-mcmc by beast-dev.

the class AlloppMoveLegsParser method parseXMLObject.

@Override
public Object parseXMLObject(XMLObject xo) throws XMLParseException {
    AlloppSpeciesBindings apsp = (AlloppSpeciesBindings) xo.getChild(AlloppSpeciesBindings.class);
    AlloppSpeciesNetworkModel apspnet = (AlloppSpeciesNetworkModel) xo.getChild(AlloppSpeciesNetworkModel.class);
    final double weight = xo.getDoubleAttribute(MCMCOperator.WEIGHT);
    return new AlloppMoveLegs(apspnet, apsp, weight);
}

Example 2 with AlloppSpeciesBindings

use of dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings in project beast-mcmc by beast-dev.

the class AlloppSequenceReassignmentParser method parseXMLObject.

@Override
public Object parseXMLObject(XMLObject xo) throws XMLParseException {
    AlloppSpeciesBindings apsp = (AlloppSpeciesBindings) xo.getChild(AlloppSpeciesBindings.class);
    AlloppSpeciesNetworkModel apspnet = (AlloppSpeciesNetworkModel) xo.getChild(AlloppSpeciesNetworkModel.class);
    final double weight = xo.getDoubleAttribute(MCMCOperator.WEIGHT);
    return new AlloppSequenceReassignment(apspnet, apsp, weight);
}

Example 3 with AlloppSpeciesBindings

use of dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings in project beast-mcmc by beast-dev.

the class AlloppSpeciesNetworkModelTEST method testLhoodMulLabTree.

public void testLhoodMulLabTree() {
    ApSpInfo[] apspecies = new ApSpInfo[3];
    apspecies[0] = new ApSpInfo("a", 2, new Individual[0]);
    apspecies[1] = new ApSpInfo("b", 2, new Individual[0]);
    apspecies[2] = new ApSpInfo("z", 4, new Individual[0]);
    AlloppSpeciesBindings testASB = new AlloppSpeciesBindings(apspecies);
    AlloppMulLabTree testamlt = new AlloppMulLabTree(testASB);
    double llhood = testamlt.testGeneTreeInMULTreeLogLikelihood();
    assertEquals(llhood, -13.562218135041552713, 1e-10);
    System.out.println(llhood);
}

Example 4 with AlloppSpeciesBindings

use of dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings in project beast-mcmc by beast-dev.

the class AlloppSpeciesNetworkModelTEST method testNetworkToMulLabTree.

public void testNetworkToMulLabTree() {
    ApSpInfo[] apspecies = new ApSpInfo[5];
    apspecies[0] = new ApSpInfo("a", 2, new Individual[0]);
    apspecies[1] = new ApSpInfo("b", 4, new Individual[0]);
    apspecies[2] = new ApSpInfo("c", 4, new Individual[0]);
    apspecies[3] = new ApSpInfo("d", 4, new Individual[0]);
    apspecies[4] = new ApSpInfo("e", 2, new Individual[0]);
    AlloppSpeciesBindings testASB = new AlloppSpeciesBindings(apspecies);
    AlloppSpeciesNetworkModel testASNM = new AlloppSpeciesNetworkModel(testASB);
    System.out.println("Tests of Network To MulLabTree conversion with dips a,e and tets b,c,d");
    String newick;
    newick = testASNM.testExampleNetworkToMulLabTree(1);
    System.out.println(newick + "\n\n");
    assertEquals(0, newick.compareTo("((((b0,c0),d0),a),(((b1,c1),d1),e))"));
    newick = testASNM.testExampleNetworkToMulLabTree(2);
    System.out.println(newick + "\n\n");
    assertEquals(0, newick.compareTo("(((((b0,c0),d0),a),((b1,c1),d1)),e)"));
    newick = testASNM.testExampleNetworkToMulLabTree(3);
    System.out.println(newick + "\n\n");
    assertEquals(0, newick.compareTo("(((((b0,c0),d0),((b1,c1),d1)),a),e)"));
    newick = testASNM.testExampleNetworkToMulLabTree(4);
    System.out.println(newick + "\n\n");
    assertEquals(0, newick.compareTo("((((b0,c0),a),(d0,d1)),((b1,c1),e))"));
    newick = testASNM.testExampleNetworkToMulLabTree(5);
    System.out.println(newick + "\n\n");
    assertEquals(0, newick.compareTo("(((((a,b0),c0),c1),(d0,d1)),(b1,e))"));
    System.out.println("");
    System.out.println("");
}

Example 5 with AlloppSpeciesBindings

use of dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings in project beast-mcmc by beast-dev.

the class AlloppSpeciesNetworkModelParser method parseXMLObject.

@Override
public Object parseXMLObject(XMLObject xo) throws XMLParseException {
    AlloppSpeciesBindings apspb = (AlloppSpeciesBindings) xo.getChild(AlloppSpeciesBindings.class);
    boolean onehyb = xo.getBooleanAttribute(ONEHYBRIDIZATION);
    boolean diprootisroot = xo.getBooleanAttribute(DIPLOIDROOT_ISROOT);
    final XMLObject tippopxo = xo.getChild(TIP_POPULATIONS);
    final double tippopvalue = tippopxo.getAttribute(Attributable.VALUE, 1.0);
    final XMLObject rootpopxo = xo.getChild(ROOT_POPULATIONS);
    final double rootpopvalue = rootpopxo.getAttribute(Attributable.VALUE, 1.0);
    final XMLObject hybpopxo = xo.getChild(HYBRID_POPULATIONS);
    final double hybpopvalue = hybpopxo.getAttribute(Attributable.VALUE, 1.0);
    AlloppSpeciesNetworkModel asnm = new AlloppSpeciesNetworkModel(apspb, tippopvalue, rootpopvalue, hybpopvalue, onehyb, diprootisroot);
    // don't know dimensionality until network created, so replace parameters
    ParameterParser.replaceParameter(tippopxo, asnm.tippopvalues);
    final Parameter.DefaultBounds tippopbounds = new Parameter.DefaultBounds(Double.MAX_VALUE, 0, asnm.tippopvalues.getDimension());
    asnm.tippopvalues.addBounds(tippopbounds);
    ParameterParser.replaceParameter(rootpopxo, asnm.rootpopvalues);
    final Parameter.DefaultBounds rootpopbounds = new Parameter.DefaultBounds(Double.MAX_VALUE, 0, asnm.rootpopvalues.getDimension());
    asnm.rootpopvalues.addBounds(rootpopbounds);
    ParameterParser.replaceParameter(hybpopxo, asnm.logginghybpopvalues);
    final Parameter.DefaultBounds hybpopbounds = new Parameter.DefaultBounds(Double.MAX_VALUE, 0, asnm.logginghybpopvalues.getDimension());
    asnm.logginghybpopvalues.addBounds(hybpopbounds);
    // note hybpopvalues are different and only work for logging.
    return asnm;
}