Bedtools fisher result
0
0
Entering edit mode
5.3 years ago
sxf710 • 0

I would like to run a fishers exact test using bedtools fisher. I have two files to test the overlap, and a genome file. However, the result does not give me an odds ratio, and the number of possible intervals is the same number as the number of query intervals. I have ran the same code with different files and it works fine, so I am not sure what is going on here. See the result below:

# Number of query intervals: 962714
# Number of db intervals: 18661
# Number of overlaps: 516254
# Number of possible intervals (estimated): 962714
# phyper(516254 - 1, 962714, 962714 - 962714, 18661, lower.tail=F)
# Contingency Table Of Counts
#_________________________________________
#                  |  in -b           | not in -b    |
#     in -a        | 516254           | 446460       |
# not in -a        | 0                | 0            |
#_________________________________________
# p-values for fisher's exact test
left    right   two-tail    ratio
1   1   1   -nan

How would one interpret this result and could there be an issue with the files?

genome bedtools fisher • 1.8k views
ADD COMMENT
0
Entering edit mode

The strange thing is that the counts for "not in -a" (second row of the 2 by 2 contingency table) are zero. I would suggest to check this issue

ADD REPLY

Login before adding your answer.

Traffic: 2550 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6