@@ -0,0 +1,1212 @@
+package test.beast.evolution.speciation;%0A%0Aimport junit.framework.TestCase;%0Aimport beast.core.parameter.RealParameter;%0Aimport beast.evolution.speciation.YuleModel;%0Aimport beast.evolution.tree.Tree;%0Aimport beast.util.TreeParser;%0Aimport org.junit.Test;%0A%0A/**%0A * Created by Denise on 16.11.16.%0A */%0Apublic class YuleModelTest extends TestCase %7B%0A%0A    @Test%0A    public void testYule() throws Exception %7B%0A%0A        YuleModel bdssm =  new YuleModel();%0A%0A        Tree tree1 = new TreeParser(%22((A:1.0,B:1.0):1.0,C:2.0);%22,false);%0A        bdssm.setInputValue(%22tree%22, tree1);%0A%0A        bdssm.setInputValue(%22birthDiffRate%22, new RealParameter(%2210.%22));%0A        bdssm.setInputValue(%22originHeight%22, new RealParameter(%2210.%22));%0A%0A        bdssm.initAndValidate();%0A%0A        double logP1 = bdssm.calculateTreeLogLikelihood(tree1);%0A%0A        Tree tree = new TreeParser(%22((A:1.0,B:1.0):2.0,C:3.0);%22,false);%0A        bdssm.setInputValue(%22tree%22, tree);%0A%0A        bdssm.setInputValue(%22birthDiffRate%22, new RealParameter(%2210.%22));%0A        bdssm.setInputValue(%22originHeight%22, new RealParameter(%2210.%22));%0A%0A%0A        bdssm.initAndValidate();%0A%0A        double logP2 = bdssm.calculateTreeLogLikelihood(tree);%0A%0A        assertEquals(logP1-logP2,10.0);%0A%0A    %7D%0A%0A%7D%0A

add simple yule test for comparison with bdsky