In [1]:
# oct4
intersect_lines = sum(1 for line in open('intersection_nexus_chipseq/oct4_intersected.bed'))
intersect_lines
Out[1]:
13950
In [2]:
import gzip
chipseq_lines = sum(1 for line in gzip.open('new_chipseq/oct4/oct4_peaks.bed.gz', 'rb'))
chipseq_lines
Out[2]:
19351
In [3]:
(intersect_lines/chipseq_lines)
Out[3]:
0.7208929771071263
In [4]:
nexus_lines = sum(1 for line in gzip.open('chip-nexus/oct4_peaks.bed.gz', 'rb'))
nexus_lines
Out[4]:
21841
In [5]:
(intersect_lines/nexus_lines)
Out[5]:
0.6387070189093906
In [6]:
# sox2
intersect_lines = sum(1 for line in open('intersection_nexus_chipseq/sox2_intersected.bed'))
intersect_lines
Out[6]:
6436
In [7]:
chipseq_lines = sum(1 for line in gzip.open('new_chipseq/sox2/sox2_peaks.bed.gz', 'rb'))
chipseq_lines
Out[7]:
10200
In [8]:
(intersect_lines/chipseq_lines)
Out[8]:
0.6309803921568627
In [9]:
nexus_lines = sum(1 for line in gzip.open('chip-nexus/sox2_peaks.bed.gz', 'rb'))
nexus_lines
Out[9]:
9396
In [10]:
(intersect_lines/nexus_lines)
Out[10]:
0.6849723286504895
In [11]:
# nanog
intersect_lines = sum(1 for line in open('intersection_nexus_chipseq/nanog_intersected.bed'))
intersect_lines
Out[11]:
17574
In [12]:
chipseq_lines = sum(1 for line in gzip.open('new_chipseq/nanog/nanog_peaks.bed.gz', 'rb'))
chipseq_lines
Out[12]:
40691
In [13]:
(intersect_lines/chipseq_lines)
Out[13]:
0.4318891155292325
In [14]:
nexus_lines = sum(1 for line in gzip.open('chip-nexus/nanog_peaks.bed.gz', 'rb'))
nexus_lines
Out[14]:
18017
In [15]:
(intersect_lines/nexus_lines)
Out[15]:
0.9754121107842593