zzh
/
gatk-3.8
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Added another assert to testBinomialProbability() and added a test method for testMultinomialProbability(). 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@544 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
kiran 2009-04-27 14:59:11 +00:00
parent 60fe2e06d1
commit eeb0b78cce
1 changed files with 13 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
java/test/org/broadinstitute/sting/utils/MathUtilsTest.java
View File

@ -23,7 +23,7 @@ public class MathUtilsTest extends BaseTest {
public void testBinomialProbability() {
logger.warn("Executing testBinomialProbability");
Assert.assertTrue(MathUtils.compareDoubles(MathUtils.binomialProbability(2, 3, 0.5), 0.375, 0.001) == 0);
Assert.assertTrue(MathUtils.compareDoubles(MathUtils.binomialProbability(2, 3, 0.5), 0.375, 0.0001) == 0);
Assert.assertTrue(MathUtils.compareDoubles(MathUtils.binomialProbability(10, 100, 0.5), 1.365543e-17, 1e-18) == 0);
Assert.assertTrue(MathUtils.compareDoubles(MathUtils.binomialProbability(73, 217, 0.02), 4.521904e-67, 1e-68) == 0);
Assert.assertTrue(MathUtils.compareDoubles(MathUtils.binomialProbability(100, 300, 0.02), 9.27097e-91, 1e-92) == 0);
@ -31,4 +31,16 @@ public class MathUtilsTest extends BaseTest {
Assert.assertTrue(MathUtils.compareDoubles(MathUtils.binomialProbability(120, 300, 0.98), 3.090054e-221, 1e-222) == 0);
Assert.assertTrue(MathUtils.compareDoubles(MathUtils.binomialProbability(112, 300, 0.98), 2.34763e-236, 1e-237) == 0);
}
/**
* Tests that we get the right values from the multinomial distribution
*/
@Test
public void testMultinomialProbability() {
logger.warn("Executing testMultinomialProbability");
int[] counts = { 10, 20, 30 };
double[] probs = { 0.25, 0.25, 0.50 };
Assert.assertTrue(MathUtils.compareDoubles(MathUtils.multinomialProbability(counts, probs), 0.002870301, 1e-9) == 0);
}
}