Mercurial > repos > rliterman > csp2
comparison CSP2/CSP2_env/env-d9b9114564458d9d-741b3de822f2aaca6c6caa4325c4afce/lib/python3.8/site-packages/pybedtools/stats.py @ 68:5028fdace37b
planemo upload commit 2e9511a184a1ca667c7be0c6321a36dc4e3d116d
author | jpayne |
---|---|
date | Tue, 18 Mar 2025 16:23:26 -0400 |
parents | |
children |
comparison
equal
deleted
inserted
replaced
67:0e9998148a16 | 68:5028fdace37b |
---|---|
1 import multiprocessing | |
2 from . import helpers | |
3 import pybedtools | |
4 | |
5 | |
6 def random_jaccard(x, y, genome_fn, shuffle_kwargs, jaccard_kwargs): | |
7 z = x.shuffle(g=genome_fn, **shuffle_kwargs).sort() | |
8 result = z.jaccard(y, **jaccard_kwargs) | |
9 helpers.close_or_delete(z) | |
10 return result | |
11 | |
12 | |
13 def random_intersection(x, y, genome_fn, shuffle_kwargs, intersect_kwargs): | |
14 z = x.shuffle(g=genome_fn, **shuffle_kwargs) | |
15 zz = z.intersect(y, stream=True, **intersect_kwargs) | |
16 result = len(zz) | |
17 helpers.close_or_delete(z, zz) | |
18 return result | |
19 | |
20 | |
21 def random_intersection_bp(x, y, genome_fn, shuffle_kwargs, intersect_kwargs): | |
22 z = x.shuffle(g=genome_fn, **shuffle_kwargs) | |
23 zz = z.intersect(y, stream=True, **intersect_kwargs) | |
24 result = sum(len(i) for i in zz) | |
25 helpers.close_or_delete(z, zz) | |
26 return result |