Search in sources :

Example 1 with ApSpInfo

use of dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings.ApSpInfo 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("");
}
Also used : AlloppSpeciesNetworkModel(dr.evomodel.alloppnet.speciation.AlloppSpeciesNetworkModel) AlloppSpeciesBindings(dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings) Individual(dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings.Individual) ApSpInfo(dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings.ApSpInfo)

Example 2 with ApSpInfo

use of dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings.ApSpInfo 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);
}
Also used : AlloppSpeciesBindings(dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings) Individual(dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings.Individual) ApSpInfo(dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings.ApSpInfo) AlloppMulLabTree(dr.evomodel.alloppnet.speciation.AlloppMulLabTree)

Aggregations

AlloppSpeciesBindings (dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings)2 ApSpInfo (dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings.ApSpInfo)2 Individual (dr.evomodel.alloppnet.speciation.AlloppSpeciesBindings.Individual)2 AlloppMulLabTree (dr.evomodel.alloppnet.speciation.AlloppMulLabTree)1 AlloppSpeciesNetworkModel (dr.evomodel.alloppnet.speciation.AlloppSpeciesNetworkModel)1