search

JetBrains-Research / span

SPAN Semi-supervised Peak Analyzer
https://doi.org/10.1093/bioinformatics/btab376
MIT License
9 stars 1 forks source link

Stoffer-Liptak test produces Nan #32

Closed olegs closed 2 years ago

olegs commented 2 years ago 
Caused by: java.lang.IllegalStateException: Nan during combining pvalues 4.536832820715596E-24,1.8393263841715094E-27,5.994971392624784E-21,2.06835737356540$
E-42,1.8393263841715094E-27,6.085984087566133E-11,7.803572981068073E-13,3.0593896259456266E-34,5.994971392624784E-21,4.2782295607092716E-33,2.06835737356540$
E-42,1.3611161918897815E-65,4.2782295607092716E-33,2.068357373565406E-42,1.504096577367313E-36,6.066738108307561E-77,3.2379502029748828E-78,3.62203538305709$
E-24,5.534060003552777E-46,7.022847588883002E-12,8.022790630864596E-31,6.003874037766238E-21,0.007552237026621174,0.007552237026621174,8.034673975762777E-31$
5.994971392624784E-21,2.2119764512649164E-7,2.9106729037243954E-26,1.0831485550464778E-15,3.2206438254678727E-4,2.942794634147554E-25,8.117383475641987E-18,$
.686895628392727E-15,6.085984087566133E-11,1.8393263841715094E-27,4.034142241199322E-9,4.034142241199322E-9,2.1589050666765635E-35,3.622035383057091E-24,8.6$
4047594152862E-45,3.2206438254678727E-4,4.731380466089632E-40,5.176489595662274E-22,0.0016398706491104705,6.794465920204226E-20,8.370554457775231E-14,3.0649$
8754606286E-8,5.994971392624784E-21,4.251326388864969E-68,1.2167142708556864E-93,1.2948467168860364E-60,7.755118251507132E-124,2.2698495981315425E-105,3.622$
35383057091E-24,4.3732795626831537E-23,9.090850871918053E-81,6.536893123315908E-57,3.0593896259456266E-34,2.0979958890372006E-74,3.8375216365123104E-58,6.06$
738108307561E-77,5.412172008363167E-108,7.755118251507132E-124,5.105010789765497E-52,3.4976834298945663E-47,3.4976834298945663E-47,3.622035383057091E-24,5.9$
0734275309314E-32,6.536893123315908E-57,4.2524239202042556E-63,6.792263176744292E-20,8.539658036579827E-17,3.622035383057091E-24,8.117383475641987E-18,9.250$
62722468583E-90,8.022790630864596E-31,3.089724755968893E-53,4.722908695981398E-91,1.7324958657514746E-130,2.2368774805061296E-59,8.539658036579827E-17,1.773$
316710600536E-28
        at org.jetbrains.bio.statistics.hypothesis.StofferLiptakTest.combine(StofferLiptakTest.kt:47)
        at org.jetbrains.bio.span.IslandsKt.getChromosomeIslands(Islands.kt:100)
        at org.jetbrains.bio.span.IslandsKt$getIslands$map$1.invoke(Islands.kt:49)
        at org.jetbrains.bio.span.IslandsKt$getIslands$map$1.invoke(Islands.kt:45)
        at org.jetbrains.bio.genome.containers.GenomeMap.lambda-1$lambda-0(GenomeMaps.kt:54)
        at java.util.concurrent.ForkJoinTask$AdaptedCallable.exec(ForkJoinTask.java:1424)
        at java.util.concurrent.ForkJoinTask.doExec(ForkJoinTask.java:289)
        at java.util.concurrent.ForkJoinPool$WorkQueue.runTask(ForkJoinPool.java:1056)
        at java.util.concurrent.ForkJoinPool.runWorker(ForkJoinPool.java:1692)
        at java.util.concurrent.ForkJoinWorkerThread.run(ForkJoinWorkerThread.java:157)
Arguments: [--islands, -t, bams/GSM646371_HepG2_Input_rep2.bam, --chrom.sizes, hg19.chrom.sizes, --peaks, span-islands/GSM646371_HepG2_Input_rep2_200_0.05_10
.islands, --model, span-islands/fit/GSM646371_HepG2_Input_rep2_200.span, --workdir, span-islands, --threads, 4, --bin, 200, --fdr, 0.05, --gap, 10]